METABOLOMICS

The Metabolome and Metabolomics

The metabolome is the complete set of metabolites—the small molecule substrates,
intermediates, and products of metabolism—within a cell, tissue, or biological sample at a
specific time point. It is inherently dynamic, with small molecules constantly being
absorbed, synthesized, and degraded, interacting within and between biological systems and
with the environment. The metabolome is considered the most complex and informative of
all the "omics" techniques, as it is the only true reflection of the phenotype at any given
moment. Any changes in gene expression, protein expression, or the environment can be
directly observed through the metabolome.
Objectives of Metabolomics: Metabolomics aims to identify and quantify these small
molecules to understand biological function. By studying metabolites, researchers can
understand cellular processes like growth, defense, and stimulation, and how they are
impacted by factors such as age, disease, and drugs.
Importance: Metabolomics is a non-invasive tool that is closely linked to the phenotype,
making it ideal for the pharmaceuticals, preventive healthcare, and agricultural industries. It
has enabled informed decision-making in areas like biomarker discovery and drug safety
screens. In the future, it has the potential to enable personalized medicine and more effective,
individualized treatment strategies.
Metabolic Reactions
Metabolic reactions are the chemical reactions required to sustain life. They are organized
into metabolic pathways, which are a series of enzyme-catalyzed reactions where the product
of one reaction becomes the substrate for the next.
Types: There are two fundamental types of metabolic reactions:
 Catabolic reactions break down complex compounds to release energy. An example is
the process of converting carbohydrates from food into glucose, which is then used to
generate ATP (adenosine triphosphate), the cell's energy currency.
 Anabolic reactions build up (synthesize) compounds. For instance, amino acids
combine to form proteins.
Examples of Metabolic Reactions:
The breakdown of drugs by the body for easier excretion.
Bacteria in the large intestine breaking down fiber into short-chain fatty acids that serve as
an energy source for intestinal cells.
Metabolomics Profiling and Techniques
Metabolomics involves profiling the metabolome using various analytical techniques to
capture and analyze data.
Profiling Approaches:
Untargeted metabolomics is a global approach used to identify metabolites in a sample set to
generate hypotheses.
Targeted metabolomics is used to quantify a specific set of known metabolites to test
hypotheses and better understand the findings of untargeted studies.
Technologies: The two main technologies for obtaining metabolite data are Nuclear
Magnetic Resonance (NMR) and Mass Spectrometry (MS).
 Nuclear Magnetic Resonance (NMR): This technique identifies compounds by exposing
samples to a magnetic field and radio frequency pulse. Its key features are that it is non-
invasive, meaning it doesn't destroy the sample, which can be reused. It is also highly
reproducible, quantitative, and a powerful method for identifying unknown metabolites,
making it great for novel biomarker discovery.
 Liquid Chromatography-Mass Spectrometry (LC-MS): This method measures small
molecules by ionizing them and detecting them as they move through a spectrometer. It
uses liquid chromatography to physically separate molecules prior to ionization, allowing
researchers to work with very small sample sizes. LC-MS is more sensitive than NMR
and is better suited for targeted metabolomics, as it can measure thousands of metabolites
at a time and achieve a wider coverage of the metabolome.
Comparison: While NMR is non-invasive and simpler in sample preparation, LC-MS is more
sensitive and can analyze a wider range of metabolites.
Implementation: A metabolomics project can be executed either in-house or outsourced to a
service provider. In-house projects offer direct control and immediate processing but require
significant financial investment in instruments, software, and trained personnel. Outsourcing
provides access to experts, high-quality technology, and validated methods but may not be
feasible for urgent analyses due to shipping times and budget constraints.
Databases and Multi-omics Approaches
To fully understand the human metabolome, comprehensive databases are needed to guide
hypothesis generation.
 The Human Metabolome Database (HMDB): This database was first introduced in 2007
and contains records for thousands of endogenous metabolites, drugs, and food
components. It is a curated collection of human metabolite and metabolism data compiled
from literature and experimental analyses.
 The Virtual Metabolic Human (VMH) Database: Released in 2019, this database captures
information on human and gut microbial metabolism and links this data to diseases and
nutritional information. It is an open-source project that reveals previously hidden
connections between human and gut microbial metabolic processes.
 Multi-omics: This biological analysis approach integrates data from multiple "omics" such
as the genome, transcriptome, proteome, and metabolome to better understand the
complexity of biological systems. By combining metabolomics datasets with other
"omics," researchers can complete the puzzle of biological processes. Multi-omics studies
are used in clinical research to identify patient subgroups, predict drug resistance, and
discover predictive biomarkers. It also helps overcome complexities caused by genetic and
phenotypic heterogeneity, which can hinder the understanding of diseases like cancer.

Proteomics
Analytical Proteomics Approach
Analytical proteomics is a field that studies the complete set of proteins in a biological
sample. Most analytical proteomics problems begin with a complex mixture of proteins with
varying molecular weights, modifications, and solubilities. To analyze this mixture, the
proteins must first be separated into simpler fractions and then cleaved into peptides. This is
necessary because the mass spectrometers used in proteomics cannot directly measure or
obtain sequences from intact proteins. The document outlines the essential elements of this
approach.
Protein Extraction: The process starts with a biological sample (e.g., tissue, cultured cells,
bacteria, or a leaf). This sample is disrupted to create a soup from which proteins are then
extracted. The goal is to recover as much protein as possible while minimizing
contamination from other biomaterials like lipids, nucleic acids, and carbohydrates.
Aids for Extraction: Various agents are used to facilitate protein extraction:
 Detergents (e.g., SDS, CHAPS) help solubilize membrane proteins and separate them
from lipids.
 Reductants (e.g., DTT, mercaptoethanol) reduce disulfide bonds and prevent protein
oxidation.
 Denaturing agents (e.g., urea, acids) disrupt protein-protein interactions and
secondary/tertiary structures by altering pH and ionic strength.
 Enzymes (e.g., DNase, RNase) digest contaminating nucleic acids, carbohydrates, and
lipids.
Separation: Proteins can be separated first, then cleaved, or cleaved first, then separated.
Separating the complex mixture into simpler fractions allows a greater number of
components to be analyzed by the mass spectrometer.
Mass Spectrometry (MS) Analysis: The separated peptides are then analyzed by one of two
types of mass spectrometers. The resulting data is matched with databases using specialized
software to identify peptides and their sequences, which ultimately establishes the identity of
the proteins in the original sample. This matching process is often done without direct
interpretation of the peptide sequences from the data.
Analytical Protein and Peptide Separations
Separating proteins and peptides is a crucial step in analytical proteomics. These separations
can be done using gel-based or gel-free methods.
Gel-Based Methods
These methods are used to separate intact proteins before they are digested into peptides.
One-Dimensional SDS-PAGE (1D-SDS-PAGE): This is the most widely used analytical
separation in protein chemistry. The method is based on the principle that SDS (a detergent)
binds to proteins, imparting a net negative charge proportional to their molecular weight.
When a high voltage is applied, the protein-SDS complexes migrate through a
polyacrylamide gel, and the proteins are resolved into bands in order of their molecular
weight.
Two-Dimensional SDS-PAGE (2D-SDS-PAGE): This method has become synonymous
with proteomics and is considered the best method for resolving highly complex protein
mixtures. It combines two different types of separations:
First Dimension: Proteins are resolved based on their isoelectric point (pI) using isoelectric
focusing (IEF).
Second Dimension: The focused proteins are further resolved based on their molecular
weight by electrophoresis on a polyacrylamide gel.
Preparative Isoelectric Focusing (IEF): This technique is analogous to the first step of 2D-
SDS-PAGE. The separation is performed on a gel strip or in a solution using soluble
ampholytes, which create a stable pH gradient when voltage is applied. The proteins then
separate based on their isoelectric point.
Gel-Free Methods
These methods utilize mass spectrometers for the analysis of proteins and peptides directly,
without the use of gels. The two main types of instruments are MALDI-TOF and ESI-tandem
MS.
Mass Spectrometers (MS)
Mass spectrometers are instruments used to analyze components of a mixture based on their
mass-to-charge (m/z) ratio.
 Essential Components: Every mass spectrometer has three essential parts:
 Source: Produces ions from the sample.
 Mass Analyzer: Resolves ions based on their m/z ratio.
 Detector: Detects the ions resolved by the mass analyzer.
 Data Quality: For proteomics, high-quality MS data needs to have three characteristics:
 Sensitivity: The instruments must be able to obtain data from very small quantities of
peptides, ideally in the femtomole (10 −15mole) range or less.
 Resolution: The ability to distinguish ions with very similar m/z values.
 Mass Accuracy: The instrument's ability to measure the mass of an ion precisely.
MALDI-TOF Mass Spectrometers
MALDI-TOF stands for Matrix-Assisted Laser Desorption Ionization-Time of Flight. This
type of instrument is primarily used to measure the masses of peptides.
How the MALDI Source Works:
A sample containing peptides is mixed with a chemical matrix (e.g., 2,5-dihydroxybenzoic
acid) and spotted onto a plate.
The solvent evaporates, forming a crystal lattice with the peptides integrated into it.
A laser beam fires at the target. The matrix absorbs the photons, gets electronically excited,
and transfers this energy to the peptides.
The peptides are then ejected into the gas phase, where they acquire a single proton to form
singly charged positive ions, with an m/z value that is roughly their mass plus one.
How the TOF Mass Analyzer Works:
The TOF analyzer measures the time it takes for ions to travel from one end of a flight tube
to a detector.
The speed of the ions is proportional to their m/z value; ions with a greater m/z fly faster.
Resolution can be improved using a technique called delayed extraction, which introduces a
slight delay between the laser pulse and the ion extraction. This ensures that all ions of the
same m/z value start their flight at the same time and hit the detector simultaneously.
Pros and Cons of MALDI-TOF:
Pros:
User-Friendly: MALDI-TOF instruments are among the easiest MS instruments to operate
and are compatible with "walk-up" or "open access" formats, making them suitable for
shared facilities. A single instrument can handle hundreds of analyses per day.
High-Throughput Compatibility: They are compatible with robotic sample preparation
devices, which increases the speed and reproducibility of analyses.
Accuracy: The accuracy and resolution of modern TOF analyzers are well-suited for the
demanding task of peptide mass mapping.
High Sensitivity: They can routinely deliver quality data on low femtomole quantities of
peptides, with some instruments reaching attomole (10 −18 mole) sensitivity.
Cons: The document doesn't explicitly mention cons for MALDI-TOF but implies that no
single MS instrument is perfect for all analytical proteomics work.

Electrospray Ionization (ESI)-tandem MS Instruments

ESI-tandem MS instruments are primarily used to obtain sequence data for peptides. The
document highlights their role in obtaining data that describes peptide fragmentation, which
is necessary for sequencing.

PHYLOGENETIC ANALYSIS
Principles of Phylogenetic Analysis
Phylogenetic analysis is a method used to infer the evolutionary relationships among
organisms or genes. The core principles guiding this analysis are based on the assumption
that the sequences being compared are homologous.
Homology: This refers to a similarity between sequences due to a common ancestry.
Homologous sequences evolved and diverged from a common ancestral sequence over time.
Independent Evolution: A key assumption is that each position within a sequence evolves
independently.
Steps in Phylogenetic Tree Construction
Constructing a phylogenetic tree is a multi-step process:
[Link] of Molecular Markers: The first step is to choose appropriate molecular markers
—the biological information (e.g., nucleic acid or protein sequences) used to infer
evolutionary relationships.
[Link] Sequence Alignment (MSA): This is a critical step that aligns multiple sequences
to identify blocks of conserved residues. A good alignment has fewer or shorter gaps, which
indicate sequences that were gained or lost during evolution.
[Link] of an Evolutionary Model: An evolutionary model of sequence data is a model of
nucleotide or amino acid substitution and the consequent divergence of sequences. This is
important for analyzing molecular sequence data, as it filters the complexity of the biological
mutation process into simpler, predictable patterns.
[Link] of the Phylogenetic Tree: Trees are constructed using different methods,
including distance-based and character-based methods.
[Link] of Tree Reliability: The final step is to assess the reliability of the constructed
tree.
Sequence Homology, Identity, and Similarity
These terms are fundamental to understanding phylogenetic relationships.
Sequence Homology: Sequences are homologous if they have a common evolutionary origin
from a common ancestral sequence. The document emphasizes that sequences are either
homologous or not—there is no quantification of homology. Problems arise in concluding
homology when the similarity is low. However, two proteins can be considered homologous
despite low similarity if the similarity is statistically significant over a long stretch of
sequence, the same pattern of residues is seen in multiple sequences, or the sequence pattern
reflects the similarity between their experimentally determined structures.
Sequence Identity: This is a direct measure of the same residues (amino acids or bases) being
present at corresponding positions in two sequences. It is calculated as the number of
identical amino acids divided by the total number of amino acids, multiplied by 100.
Sequence Similarity: This measures the presence of similar residues at corresponding
positions. For nucleic acids, sequence similarity is the same as sequence identity. For
proteins, it involves amino acids with similar physicochemical and functional properties. It is
calculated as the sum of identical and similar substitutions, divided by the total number of
amino acids, and multiplied by 100.
Molecular Markers
A molecular marker is the biological information (e.g., DNA, protein) used to infer
evolutionary relationships.
Nucleic Acid vs. Protein: The document recommends using protein sequences for coding
sequences because amino acid substitution matrices are more sophisticated than nucleotide
matrices. Also, the existence of codon bias can artificially inflate nucleotide sequence
variation.
When to use Nucleotide Sequences: Nucleotide sequences are preferred when comparing
genes with highly conserved regions, studying the evolution of genes within geographically
separated populations of a species, or when using slowly evolving genes for distant species
or rapidly evolving genes for recently evolved ones.
Key Markers:
Small Subunit Ribosomal RNA (rRNA): This is an ancient molecule that is vertically
transferred with a low mutation rate. The 16S rDNA is used for prokaryotes, and the 18S
rDNA is used for eukaryotes. The SSU sequence has multiple variable regions flanked by
conserved blocks.
Rare Genomic Sequences: These markers include indels (insertions and deletions), sequence
signatures, and changes in amino acid composition.
Multiple Sequence Alignment (MSA)
MSA is the core of phylogenetic analysis, as it identifies blocks of conserved residues among
multiple sequences.
Principles: The basic premise is that every residue in each column of the alignment is
homologous, meaning it evolved from the same position in a common ancestral sequence
without insertions or deletions.
Progressive Sequence Alignment: This is a method for constructing MSA. It involves two
main steps:
[Link] the Guide Tree: This is the first step, which determines the order in which
sequences are added to the growing alignment.
[Link] Addition: Sequences are added sequentially based on the order from the guide
tree.
The "Cardinal Rule": The cardinal rule of progressive sequence alignment is "once a gap,
always a gap," meaning gaps can only be added or enlarged, never moved or removed.
Phylogenetic Tree Construction Methods
The methods for calculating phylogenetic trees fall into two general categories.
[Link]-Matrix Methods: These are also known as clustering or algorithmic methods,
such as UPGMA and Neighbour-Joining.
Process: A single statistic—the distance, or percent sequence difference—is calculated for all
pairwise combinations of sequences, and these distances are then assembled into a tree.
Features: Distance methods are much faster than discrete data methods but yield little
information other than the tree itself.
[Link] are also known as tree-searching methods, such as Parsimony and Maximum
Likelihood.
Process: These methods examine each column of the alignment separately and look for the
tree that best accommodates all of the information.
Features: Discrete data analyses are information-rich because they generate a hypothesis for
every column in the alignment, allowing for the tracing of evolution at specific sites.
Evolutionary Models
Evolutionary models of sequence data are crucial for phylogenetic analysis because they
filter the complexity of mutations into simpler patterns that can be described and predicted.
Problem: A simple count of the differences between two sequences can underestimate the
amount of evolution that has occurred, especially when multiple mutations have happened at
the same site.
Features: Substitution models predict the rate of substitution for nucleotides and amino acids
at a given site and the distribution of substitutions across the entire sequence. The differential
rate of substitution is called rate heterogeneity.
Types:
Jukes-Cantor Model: This is a simple model that assumes all nucleotide substitutions are
equally likely.
Kimura Model: This model is more sophisticated because it accounts for transitions (purine-
to-purine or pyrimidine-to-pyrimidine changes) and transversions (purine-to-pyrimidine
changes) separately. Transitions are generally more common than transversions during
evolution.
Assessing Tree Reliability
After a tree is constructed, its reliability must be assessed.
Bootstrapping: This is the simplest test of phylogenetic accuracy. It determines whether the
entire dataset supports the tree or if it is just a marginal winner among many alternatives.
Process: Bootstrapping involves taking random subsamples of the original dataset, building a
tree from each subsample, and then calculating the frequency with which the parts of the tree
are reproduced across all subsamples. The results are displayed as a bootstrap value (as a
percentage), with values of 50% or higher typically reported.
Homologous Genes
Homologous genes share a common evolutionary origin. There are two types of homologous
genes:
Orthologues: These genes only duplicate when their host divides. They are strictly vertically
transmitted (from parent to offspring), so their phylogeny directly traces the lineage of their
host.
Paralogues: These genes are members of multigene families and arise from gene duplication
events within a single species.
Tree Types and Visualization
Types of Trees:
1. Monophyletic (or Holophyletic): A group that consists of a common ancestor and all
of its descendants.
2. Paraphyletic: A group that consists of a common ancestor and some, but not all, of its
descendants.
3. Polyphyletic: A group that consists of two or more distantly related lineages that do
not share an immediate common ancestor.
Visualization: Trees can be represented in various layouts, including rectangular, circular,
radial, and slanted. Standard tree file formats include New Hampshire/Newick, NHX, and
Nexus. Tree visualization tools like Dendroscope and MEGA are used to display and
annotate trees. Branch lengths on a tree are almost always drawn to scale and provide a
general impression of the relative rates of change.

Systems Biology -1
Systems Biology: Principles and Core Concepts
Systems biology is an interdisciplinary field that uses a holistic approach to study the
complex interactions within biological systems. Instead of studying isolated components, it
aims to understand how components like genes, proteins, and metabolites interact
dynamically to produce cellular functions and emergent properties. It integrates biology with
mathematics and computational science.
Guiding Principles:
The "whole is more than the sum of its parts" is a central theme, as it focuses on emergent
properties that arise from the collective behavior of a system's components.
It studies biological systems as networks of interacting components.
It uses quantitative data and computational models to simulate and predict biological
functions, which are then iteratively refined through experimental validation.
It investigates biology at multiple levels, from molecules to cells to organisms.
Key Components and Disciplines
Systems biology relies on the integration of various disciplines and technologies to build
predictive models and answer biological questions.
1. Data Generation and Integration: This involves collecting large-scale experimental data
from high-throughput technologies known as omics. These include:
 Genomics: The study of the complete DNA sequence and genetic variation.
 Transcriptomics: The examination of RNA transcripts to understand gene expression.
 Proteomics: The analysis of the entire set of proteins and their functions.
 Metabolomics: The study of small molecules and metabolites that reflect biochemical
activity.
 Epigenomics: The assessment of chemical modifications on DNA or histones that
regulate gene activity.
 Data is standardized to allow for coherent integration and sharing.
2. Mathematical and Computational Modeling: Mathematical models are constructed to
represent interactions among biological components. Computational tools and
 simulations are used to predict the behavior and dynamics of these systems, with models
being iteratively refined through experimental validation.
3. Network and Systems Analysis: This involves analyzing molecular interaction networks
and identifying functional modules within them. This analysis helps explore emergent
properties that are not predictable from studying isolated parts.
4. Experimental Validation: This component involves designing system-wide experiments
based on modeling predictions. Omics data and other experimental results are used to test
hypotheses generated from computational analyses.
Approaches and Methodologies
Systems biology uses different approaches and methodologies to study biological systems,
integrating both experimental and computational frameworks
 Bottom-Up Approach: This approach builds a system by combining small functional
units. It starts from the level of genes, proteins, and metabolites and works up to
cellular and tissue functions.
 Top-Down Approach: This approach starts with a global view of the system and then
studies smaller subsystems to understand its behavior. It begins at the organism or
organ system level and moves down to the molecular level.
Modeling Workflow: The process of creating a systems biology model is iterative and
involves several steps:
1. Problem Formulation: Clearly state the questions, background, and hypotheses to be
addressed by the model.
2. Verification of Available Information: Review existing data about the system's
structure.
3. Selection of Model Structure: Determine the general type of model, such as
microscopic or macroscopic, deterministic or stochastic, or discrete or continuous.
4. Establishing a Simple Model: Create an initial model in words, a scheme, or a
mathematical formulation.
5. Sensitivity Analysis: Verify the choice of parameters, as mathematical simulations can
be highly sensitive to parameter changes.
6. Experimental Tests: Conduct experiments to test the model's predictions.
7. Stating Agreements and Divergences: Compare model results with experimental
behavior and check for false assumptions if they don't agree.
8. Iterative Refinement: Continuously refine the model until it accurately represents the
biological system.
Core Concepts of Systems Biology
Systems biology is built upon several core concepts that distinguish it from the traditional
reductionist approach.
  Nonlinearity: The relationship between a system's inputs and outputs is often not
linear. A small change in a component can have a large, disproportionate effect on the
system's behavior.
 Feedback Loops: Biological systems are governed by feedback loops, where the
output of a process can either inhibit (negative feedback) or stimulate (positive
feedback) the process itself. For example, a metabolic pathway's end product can bind
to and inactivate an enzyme at an earlier step, regulating the pathway.
 Open Systems: Living systems are open systems that dissipate energy and exchange
matter with their environment.
 Memory: Biological systems can have a form of memory, where their response is
dependent on past history. For example, a new protein may remain in a cell after an
initial response, shifting the reaction rate the next time the cell is exposed to a
chemical.
 Modularity: Biological networks are not random; they are organized into functional
modules or smaller sets of interacting components that perform common functions.
These modules are frequently reused in different networks. For example, a
chemosensory module can detect extracellular molecules, a transduction module can
process the signal, and an actuator module like a flagellum can execute a response .
 Robustness: This refers to a system's insensitivity to variations in its parameters or
components.
Types of Models in Systems Biology
Systems biology employs different types of models to capture various aspects of a biological
system.
1. Biochemical Models: These are mathematical models that describe chemical reactions
and molecular interactions, focusing on the transformation of substrates to products via
enzymes. They often use ordinary differential equations (ODEs) to represent time-
dependent changes in molecule concentrations and can be used for flux balance analysis
(FBA) to predict metabolite flow and identify bottlenecks in metabolic engineering. An
example is the reconstruction of the metabolic network in Escherichia coli.
2. Genetic Models: These models represent how genes and their interactions control
biological functions, development, and inheritance. They can be statistical (describing
inheritance) or mechanistic (capturing regulatory relationships). A common tool is
Boolean networks, which represent genes as on/off states and their interactions as logical
rules. An example is the segment polarity gene network in Drosophila melanogaster.
3. Systems Models: These are integrative models that combine genetic, biochemical, and
other data to simulate entire biological systems or subsystems, revealing emergent
properties and predicting behavior under various stimuli. They can be deterministic,
stochastic, or hybrid. An example is a digital twin of a human immune system, which
integrates molecular, cellular, and organ-level interactions to simulate personalized
disease progression and treatment effects.
Evolution of Systems Biology
Systems biology emerged from the convergence of several key scientific and technological
advances over the past century.
 Early Foundations: The field's roots can be traced back to early concepts like non-
equilibrium thermodynamics (1931) and the discovery of feedback regulation in
metabolism (1952). Mathematical models of biological processes, such as the Hodgkin-
Huxley model of a neuron, also laid the groundwork.
 Molecular Revolution: The unraveling of DNA structure (1953) was a pivotal moment,
leading to the development of recombinant DNA technology (1973) and the Polymerase
Chain Reaction (PCR) (1983).
 The Data Explosion: The first bacterial genome was sequenced in 1995, followed by the
completion of the Human Genome Project in 2003. This massive increase in available
data,coupled with advances in high-throughput omics technologies and high-performance
computing, made a systems-level analysis not just possible, but critical.

Systems Biology -2
Equilibrium Binding and Cooperativity
Binding is the process where two or more molecules recognize and attach to each other
through intermolecular interactions. These interactions are typically non-covalent and are
characterized by molecular complementarity. A ligand is a small molecule that binds to a
specific receptor, which is a larger, relatively stationary molecule. The binding energy
causes a structural change in the receptor, which is associated with a molecular
communication event.

Studying Binding: The two general methods to study ligand-receptor interactions are
equilibrium thermodynamics and association/dissociation kinetics.

Association Constant (Kf): This constant describes the equilibrium in favor of the
complex, where a large value indicates a high affinity of the receptor for the ligand.

[RL]
Kf =
[ R][ L]

Dissociation Constant (Kd): This is the inverse of the association constant and is the
typical way to describe receptor affinity. A small K d value indicates a high affinity for the
ligand, as the equilibrium favors the complex.
1 k −1
Kd= =
K f k +1
Michaelis-Menten Equation
The Michaelis-Menten (MM) model is a classic biochemical model used to
mathematically describe enzyme-catalyzed reactions. In this model, an enzyme (E) binds
to a substrate (S) to form an enzyme-substrate complex (ES), which then either dissociates
or forms a product (P). This model is fundamental for understanding enzyme efficiency
and for analyzing metabolic pathways.
Equation: The MM equation describes the reaction velocity (v) as a function of substrate
concentration [S]:

Assumptions: The model assumes that the substrate concentration is much larger than the
enzyme concentration, that the system is in a steady state where the ES complex
concentration is constant, and that the initial velocity is measured before product
accumulation.
Cooperativity
Cooperativity is a phenomenon in which the binding of one substrate molecule to an
enzyme's subunit changes the shape of neighboring subunits, affecting their binding
affinity for subsequent substrate molecules.
 Positive Cooperativity: The binding of the first substrate molecule makes it easier for
the second to bind.
 Negative Cooperativity: The binding of the first substrate makes it more difficult for
the second to bind.
 Hill Equation: This is a classic biochemical model that quantifies cooperative binding.
The Hill coefficient (n) indicates the type of cooperativity:
n=1: No cooperativity.
n>1: Positive cooperativity.
n<1: Negative cooperativity.

Biochemical Models and Central Dogma

Biochemical models are mathematical models that describe chemical reactions and
molecular interactions in biological systems. They focus on the transformation of
substrates to products via enzymes and often represent these as networks of biochemical
reactions.
 Example: A classic biochemical model is the reconstruction of the metabolic network in E.
coli, which integrates all known metabolic reactions to simulate cellular growth and
metabolic fluxes.
 Metabolic Flux Analysis (MFA): This technique caused a shift in the paradigm of the
central dogma, which traditionally describes the flow of genetic information from DNA
to RNA to protein. MFA demonstrates that by modifying enzymes (proteins), the
metabolic flux and the final product can be changed, a concept that is distinct from the
one-way flow of information in the central dogma.
 Metabolic Models of Microorganisms: These models represent the growth of a
microorganism as a reaction, using the set of metabolic reactions occurring within it.
They can be used to model growth behaviors and predict growth rates based on factors
like glucose availability. A biomass objective function (BOF), which is an equation that
represents the cellular biomass composition, is a key component of these models.

Systems Biology -3
Cis-Regulatory Elements (CREs)
Cis-regulatory elements (CREs) are segments of DNA that control gene expression,
orchestrating processes like tissue identity, developmental timing, and responses to stimuli.
They are considered "DNA switches" that can precisely turn genes on and off.
Types of CREs:
 Promoters: These elements initiate the transcription process by providing a binding site
for transcription factors and RNA polymerase.
 Enhancers: These elements boost gene expression by facilitating transcription factor
interactions at the promoter.
 Silencers: These elements reduce gene expression by inhibiting transcription factor
binding or their interaction with the promoter.
Engineered CREs: A platform has been developed using deep neural networks to design and
validate synthetic CREs capable of driving gene expression with programmed cell-type
specificity. These synthetic sequences have been shown to be more effective at driving cell-
type-specific expression in three cell lines compared to natural sequences from the human
genome.
Biological Oscillators and Rhythms
Biological oscillators are the internal, often molecular, mechanisms that generate and
regulate biological rhythms. Biological rhythms are the observable patterns, while oscillators
are the underlying systems that drive them. Rhythms originate from feedback loops within
regulatory networks and control every aspect of cell physiology.
Biological Rhythms: These include a wide range of cycles, such as circadian rhythms
(following a ~24-hour cycle), diurnal, ultradian, and infradian rhythms, as well as heartbeat,
respiration, and the cell cycle.
Types of Oscillators:
 Cytoplasmic Oscillators: These are formed by protein enzymes and small molecule
metabolites interacting within the cytoplasm. Glycolytic oscillations in yeast cells are a
classic example, where metabolite concentrations in the glycolysis pathway fluctuate in
a periodic manner, driven by the cyclic regulation of phosphofructokinase.
 Membrane Oscillators: These are formed by regulatory ion channels restricted to the
cell membrane, such as action potentials in neural and cardiac pacemaker cells.
 Genetic Oscillators: These are composed of genes that regulate each other inside cells,
causing rhythmic changes in gene expression and protein concentrations due to
transcriptional/translational feedback loops.
Genetic Oscillators
A genetic oscillator involves rhythmic changes in gene expression driven by regulatory
feedback loops. The core mechanism is a negative feedback loop where a gene product
inhibits its own production.
Requirements for Oscillation:
1. Negative Feedback: A negative feedback loop is necessary to carry the system back to
its "starting point" of oscillation.
2. Time Delay: The feedback signal must be sufficiently delayed so that the reactions do
not settle into a stable steady state. This delay can be caused by transcription, translation,
protein modification, or the time needed for molecules to be transported between cellular
compartments.
3. Nonlinearity: The kinetic rate laws of the reaction mechanism must be sufficiently
nonlinear to destabilize the steady state. This is often modeled by Hill-type repression.
Examples:
 PER/CRY loop: The core of the circadian clock is a negative feedback loop. CLOCK and
BMAL1 proteins activate PER and CRY gene transcription. The resulting PER/CRY
proteins then form a complex that inhibits CLOCK/BMAL1, suppressing their own
transcription and creating a ~24-hour cycle.
 The Repressilator: The first synthetic genetic oscillator, built in E. coli in 2000. It is a
ring of three genes, where each gene's repressor protein inhibits the next gene in the loop,
creating a cyclic negative feedback loop that leads to oscillating protein levels.
 Goodwin Oscillator: A theoretical model of a single gene with negative autoregulation,
where the gene's product represses its own transcription after a delay, allowing for
periodic expression.
Feedback Loops and Related Concepts
Feedback loops are a core concept in systems biology that enable regulatory networks to
function.
 Negative Feedback Loops: The output of a system inhibits its own production, which
stabilizes the system and can generate oscillations.
 Positive Feedback Loops: The output enhances its own production, which amplifies
signals, reinforces states, and can lead to bistability (on/off states). In the cell cycle, a
double-negative feedback loop (where Cdk1 inhibits its inhibitor, Wee1) acts as a
functional positive feedback loop, leading to a switch-like, irreversible activation of
Cdk1.
 Coupled Feedback Loops: Many biological systems, like circadian clocks, use multiple
feedback loops that interact with each other to enhance stability and precision. The
presence of positive feedback loops can increase the amplitude and robustness of
oscillations.
 Bistability and Hysteresis: Bistability is a state where a reaction network has two stable
steady states. A bistable switch is a key consequence of a positive feedback loop that
can delay a negative feedback signal, preventing the system from damping into a single
stable steady state. Hysteresis is a property of bistable systems where the state of the
system depends on its recent history.

Synchronization and Cell Cycle Coupling

The synchronization of spatially distributed oscillations is vital for the proper functioning of
many biological systems.
 Synchronization: The process where two interacting oscillating systems with different
periods assume the same period. This is seen in pacemaker cells of the heart, β-cells in the
pancreas, and in the synchronized flashing of male fireflies. The proper functioning of
many systems relies on this collective activity, and its disruption, as in a heart arrhythmia,
can be life-threatening.
 Cell Cycle and Circadian Clock Coupling: The circadian clock can regulate the
progression of the cell cycle by influencing regulators of cyclin/Cdk activities, such as
Wee1 and p21. This coupling can result in circadian gating, where cell cycle progression is
allowed only at certain circadian phases. Alternatively, the circadian clock can entrain
(synchronize) the cell cycle into oscillations with a specific period (e.g., 24 or 48 hours),
depending on the strength of the regulation.

Systems Biology -4
 Gene Regulation: Prokaryotes vs. Eukaryotes
Gene regulation is the process by which a cell controls which genes are expressed and at
what level. This process differs significantly between prokaryotic and eukaryotic
organisms.

Prokaryotic Gene Regulation:

 Structure: Prokaryotes lack a nucleus, so their DNA is found in the cytoplasm.
 Process: Transcription and translation occur almost simultaneously in the cytoplasm.
 Regulation: Gene expression is regulated primarily at the transcriptional level. The
operon model, such as the lac and trp operons, is a classic example of this type of
regulation.
Eukaryotic Gene Regulation:
 Structure: Eukaryotes have a nucleus that contains their DNA.
 Process: Transcription occurs in the nucleus, followed by translation in the cytoplasm,
physically separating the two processes.
 Regulation: Gene expression is regulated at many different levels, including epigenetic,
transcriptional, post-transcriptional, translational, and post-translational control.

Operon Model of Prokaryotic Gene Regulation

Operons are clusters of genes that are transcribed together and regulated by a single
promoter. The trp and lac operons in E. coli are key examples.
The trp Operon: This operon contains the five genes needed to synthesize the amino acid
tryptophan.
 Absence of Tryptophan: The repressor protein is inactive and does not bind to the
operator, allowing RNA polymerase to bind to the promoter and begin transcription.
RNA synthesis proceeds, and the genes for tryptophan synthesis are expressed.
 Presence of Tryptophan: Tryptophan acts as a corepressor, binding to the repressor
protein and activating it. The activated repressor then binds to the operator, blocking
RNA polymerase from transcribing the genes and stopping RNA synthesis.
The lac Operon: This operon contains genes needed to metabolize lactose.
 Absence of Lactose: The lac repressor binds to the operator, blocking RNA polymerase
and preventing transcription.
 Presence of Lactose: Lactose binds to the repressor, causing it to be released from the
operator. Transcription can then proceed, but at a low rate.
 Catabolite Repression: Transcription of the lac operon is only highly efficient when
glucose is absent. In the absence of glucose, a cAMP-CAP complex forms and binds to
the promoter, stimulating RNA polymerase activity and significantly increasing RNA
synthesis. However, even with the cAMP-CAP complex, transcription is still blocked if
the repressor is bound to the operator.

Eukaryotic Gene Regulation Mechanisms

 Eukaryotes regulate gene expression at multiple levels, allowing for more complex control
than in prokaryotes.
1. Epigenetic Control: This involves controlling how tightly DNA is packed, influencing
the accessibility of genes for transcription.
2. Transcriptional Control: This is the primary level of regulation and determines whether
or not a gene is transcribed into pre-mRNA. This is achieved through the binding of
transcription factors to promoter and enhancer regions. Enhancers, which can be located
far upstream from the gene, bind to activator proteins, causing the DNA to bend and
facilitate the binding of RNA polymerase to the promoter to initiate transcription.
3. Post-Transcriptional Control: This occurs after transcription but before translation. It
includes processes like alternative splicing, where a single pre-mRNA molecule can be
spliced in different ways to create multiple different protein products.
4. Translational Control: This involves regulating the translation of mRNA into protein. The
stability of the mRNA molecule is influenced by RNA-binding proteins that attach to the
5' or 3' untranslated regions (UTRs).
5. Post-Translational Control: This occurs after a protein has been synthesized, regulating
its activity, stability, or localization through modifications like phosphorylation.

Transcription Networks
A transcription network describes the regulatory interactions within a cell, where nodes
represent genes and directed edges represent the transcriptional regulation of one gene by
the protein product of another. These networks are a core part of the cell's information
processing function, determining the production rate of each protein in response to
environmental signals.
Transcription Factors (TFs): These are special proteins that act as mediators between
environmental signals and gene expression. They rapidly switch between active and
inactive states, modulated by a specific environmental signal. When active, a TF can bind
to DNA to regulate the transcription rate of its target genes.
Activators vs. Repressors:
 An activator is a transcription factor that increases the transcription rate when it binds
to the promoter.
 A repressor is a transcription factor that decreases the transcription rate when it binds
to the promoter. Some genes still show a weak (basal) transcription even when a
repressor is bound. In E. coli and yeast, activation (positive control) is more common
than repression (negative control).
Regulatory Circuits: These are representations of the interactions and logic within a
system. They are dynamic, often involving nonlinear functions and feedback loops.
Examples include negative autoregulation, feed-forward loops, multi-component loops,
and regulator chains.
 Modeling Gene Regulation
Mathematical modeling is an essential tool for studying gene regulatory circuits.
Production and Degradation: A simple model of gene expression can be described by a
single ordinary differential equation (ODE), where the protein concentration at steady
state ( X st ) is the ratio of its production rate (β) to its degradation rate (γ). The cell can
control protein levels by modulating transcription, translation, or the degradation rates of
mRNA and protein.

Input Functions: The relationship between the concentration of an active transcription

factor ¿ ¿ ) and the rate of gene production (Y) is described by an input function.
 Hill Equation: The Hill input function is a nonlinear model that describes activation
and repression.
 Hill Coefficient (n): This parameter governs the steepness of the input function.
Higher values of n create a more switch-like (sigmoidal) response, while n=1 results in
a hyperbolic (Michaelis-Menten-like) response. The Hill coefficient provides a
measure of cooperativity of substrate binding to the protein.

Negative vs. Positive Control: These concepts are defined by the operon's response when
no regulator protein is present.
 Negative Control: Genes are expressed by default and are switched off by a repressor
protein. This acts as a "fail-safe" mechanism.
 Positive Control: Genes are expressed only when an active regulator protein (an
activator) is present. This mechanism requires an extrinsic or intrinsic event to trigger
it.

Systems Biology -5
Network Building & Analysis
Network building is a fundamental concept in systems biology that involves representing
cellular systems as graphs to understand their organization and function.
Graph Theory: Graphs are mathematical structures made of vertices (nodes) and
connections (edges). In cell biology, nodes represent cellular components (like genes or
proteins), and edges represent the interactions, reactions, or relationships between them. This
approach is useful for understanding the topology of biological systems.
Modeling Requirements: To create effective biological models, several requirements must
be met:
i. Network Visualization: Models should be intuitive and easily interpreted by people.
ii. Topological Analysis: A significant portion of the work is based on analyzing the
topology of biological networks, such as degree distribution and network motifs
iii. Modularity and Hierarchy: The cell is organized as a set of connected modules with
specific functions, and this structure must be captured.
iv. Standardization: Biological models should be represented in a common format for
exchange between different tools. The Systems Biology Markup Language (SBML) has
become the de facto standard for this community.

Types of Biological Graph Networks

Biological networks are classified according to their structure and the information they
convey.
 Undirected Graphs: These are the most basic network representations, where nodes are
linked by lines without a specified direction. They indicate that a relationship exists
between nodes but provide no information about the nature of the interaction (e.g.,
activation or repression). An example is a protein-protein interaction between scaffold
proteins.
 Directed Graphs: In these graphs, nodes are connected by directed edges (arrows) that
restrict signal flow and show a cause-and-effect relationship. Directed edges can
represent reversible reactions where the flow is dominant in one direction.
 Sign-Specified Directed Graphs: These are the most useful type of directed graph for
developing dynamical models. Nodes are connected by black arrows for stimulatory
(activating) relationships and red plungers for inhibitory (repressing) relationships. In E.
coli and yeast, 60-80% of interactions are positive (activating).
 Weighted Graphs: In these graphs, the link between two nodes has a relative value or
weight. The weight can represent the strength of a relationship, such as a weak or strong
interaction between proteins, and can be determined by methods like NMR or X-ray
crystallography.

Types of Networks Relevant to Cellular Systems Biology

The document classifies biological networks according to their biological function.
 Signaling Networks: These networks involve signal transduction, a communication
process where the cell receives and responds to external stimuli from other cells and the
 environment. This process affects most of the cell's basic control mechanisms, such as
differentiation and apoptosis.
 Gene Regulatory Networks (GRN): GRNs control the expression of genes and,
consequently, all cellular functions. They are essential for the survival of the cell,
ensuring that only selected functions are active at a given moment. In a GRN, nodes are
genes, and edges represent transcriptional regulation of one gene by the protein product
of another.
 Metabolic Networks: These networks are composed of a set of biochemical reactions that
a cell uses to sustain its growth and energy requirements. They include catabolic and
anabolic pathways of enzyme-catalyzed reactions that transform substrates into energy
and building blocks, forming complex interconnected networks.
 Protein-Protein Interaction Networks: These networks describe the direct interactions
between proteins within a cell. The yeast-2-hybrid (Y2H) method is a tool for detecting
these direct interactions.

Topological Properties of Networks

The topological properties of a network describe its overall structure and organization. Three
main types of networks are observed in complex systems.

 Random Networks: In a random network, the degree distribution of nodes follows a

Poisson distribution, meaning most nodes have an intermediate number of connections.
This model is typical in network engineering but rarely explains the architecture of
biological networks.
 Scale-Free Networks: In a scale-free network, the degree distribution follows a power-
law, implying the existence of nodes with a significantly higher number of connections.
These are called hubs. Most biological networks, including metabolic, genetic, and
protein-protein networks, are scale-free.
 Hierarchical Networks: These networks are also scale-free but contain a modular or
hierarchical structure. In E. coli, for example, some transcription factors interact with
many genes and are called global regulators.

Centrality Measurements
Centrality measurements are used to analyze a network's graph properties and identify
distinctive nodes with specific connectivity properties. The most common measures are
degree, betweenness, closeness, and clustering.
Degree Centrality: This measures the number of links connected to a node.
 1. In-degree: The number of links pointing to a node, highlighting nodes that receive
more incoming connections.
2. Out-degree: The number of links starting from a node, highlighting nodes that emit
more connections.
Betweenness Centrality: This measures how critical a node is as a bridge of information
across the network. It estimates the proportion of nodes over which a given node has
influence and indicates which nodes spread or control information flux better. Blocking or
removing nodes with high betweenness centrality can effectively disrupt communications in
the network.
Closeness Centrality: This measures the average number of link steps from a given node to
all other nodes in the network.
Clustering Coefficient: This measures the ratio of the number of triangles incident to a node
relative to the maximum possible number of triangles in the network.

Biological Databases
Major classes of biological databases are used to store and organize biological data, which is
essential for network analysis and modeling.
Genes/Genomics:
GenBank: A database for DNA sequences and sequence analysis.
GEO (Gene Expression Omnibus): A genomics database supported by the NCBI that
contains microarray and sequence-based data.
OMIM (Online Mendelian Inheritance in Man): A catalog of human genes and genetic
disorders.
ENCODE (Encyclopedia of DNA Elements): Contains information on all functional
elements in the human DNA sequence.
Proteins:
PDB (Protein Data Bank): A database for protein structures.
UniProtKB: A database for protein characteristics.
Pathways and Systems:
KEGG (Kyoto Encyclopedia of Genes and Genomes): A database of biological functions and
systems, including pathways
Pathway Commons: A database of biological pathways from multiple organisms.