ENGINE SCHEMATIC

The complexity is ours.
The clarity is yours.

You give us your protocol and your bloodwork. The engine does the rest.

A standard blood test is one-dimensional — one marker, one range. The TRT Plus engine reads the clinical relationships between markers, weighing ratios and dynamic interactions, not just absolute values.

Inputs
You provide
The Engine
We built
Outputs
You receive
Inputs
You provide
MY
PROTOCOL

Your individual protocol — dose, frequency, all compounds.

Testosterone
+Test Cyp💪🏼140 mg/W (D)
Compounds
+GH🧬1 IU (D)
+Anavar🐂20 mg (3xW)
+BPC-157🐺250 mcg (D)
+Fish Oil🐟6 g (D)
MY
BIOMARKERS

Comprehensive 15+ biomarker panels analysed against optimized reference ranges.

IGF-1
L204ng
Total T
1278ng/dL
E2
41pg/mL
Hct
54%
BP
120/79
ApoB
0.7g/L
Body Fat
8.8%
TSH
1.3mIU/L
CRP
1.4mg/L
SHBG
35nmol/L
DHT
47ng/dL
PSA
0.78ug/L
HbA1c
4.9%
GGT
22U/L
CysC
0.68mg/L
The Engine
We built

Deep Analysis of
Real-World TRT Patients

De-identified results form the foundation of every model and reference range.

Data Foundation

HIPAA-compliant · De-identified · Encrypted · Secure

Clinical Framework

Evidence-based protocols from real clinician input.

AI-Enhanced

Pattern recognition across thousands of cases.

Continuous Learning

Models refine as new data & studies flow in.

Healthspan Advocacy

Promotes real bloodwork follow-through and whole-health vigilance.

Outputs
You receive
BIOMARKER SCORE
& RATIOS
37:1
T:SHBG
Optimal
28:1
A:E
Balanced
PERSONALIZED BIOMARKERS
REPORT
  • Recommendations
  • Hormonal Balance
  • Muscle-Building Potential
  • Dose Response
  • Compound Analysis
  • Healthspan Assessment
  • Cardiovascular Risk Profile
  • Overall Health
The TestimatorTM
OPTIMIZED STEADY-STATE
Read your response curve — dial in your protocol.
Under the Hood

Real engineering.
Real math.

Every score, range, and insight in TRT Plus is produced by a deterministic, programmatically engineered model — not approximated by a generic AI. Clinical-grade analytical depth.

Multi-Route PK

Multi-Route Compound Plotter

Steady-State Concentration  (Multi-Route):Ctotal(t)=i=1nRi(tti)11ekiτiRoute-Specific Release Functions:R(t)=(AEkαroute)ekt(depot, oral, peptide)αroute{αdepot,  αaqueous,  αpeptide,  αoral}R(t)=Dabskakekake(eketekat)(transdermal)Dabs=VmaxDappliedKm+Dapplied(Michaelis-Menten absorption)Exponential Decay:C(t)=C0ekt,k=ln(2)t1/2Parameters:n=total doses,ti=administration time,τ=dosing intervalA=amount,E=effectiveness,k=decay constantαroute=absorption rate modifier (route-dependent)ka=transdermal absorption rate constant (site-specific)ke=elimination rate constantDabs=absorbed dose per administrationVmax,Km=Michaelis-Menten saturation parameters\begin{gathered} \text{Steady-State Concentration}\;\text{(Multi-Route):} \\[0.3em] C_{total}(t) = \sum_{i=1}^{n} R_i(t - t_i) \cdot \frac{1}{1 - e^{-k_i\tau_i}} \\[1.5em] \text{Route-Specific Release Functions:} \\[0.3em] R(t) = (A \cdot E \cdot k \cdot \alpha_{route}) \cdot e^{-kt} \quad \text{(depot, oral, peptide)} \\[0.3em] \alpha_{route} \in \{\alpha_{depot}, \; \alpha_{aqueous}, \; \alpha_{peptide}, \; \alpha_{oral}\} \\[0.8em] R(t) = \frac{D_{abs} \cdot k_a \cdot k_e}{k_a - k_e} \cdot \left(e^{-k_e t} - e^{-k_a t}\right) \quad \text{(transdermal)} \\[0.3em] D_{abs} = \frac{V_{max} \cdot D_{applied}}{K_m + D_{applied}} \quad \text{(Michaelis-Menten absorption)} \\[1.5em] \text{Exponential Decay:} \\[0.3em] C(t) = C_0 \cdot e^{-kt}, \quad k = \frac{\ln(2)}{t_{1/2}} \\[1em] {\scriptstyle \text{Parameters:}} \\[0.2em] {\scriptstyle n = \text{total doses}, \quad t_i = \text{administration time}, \quad \tau = \text{dosing interval}} \\[0.3em] {\scriptstyle A = \text{amount}, \quad E = \text{effectiveness}, \quad k = \text{decay constant}} \\[0.3em] {\scriptstyle \alpha_{route} = \text{absorption rate modifier (route-dependent)}} \\[0.2em] {\scriptstyle k_a = \text{transdermal absorption rate constant (site-specific)}} \\[0.2em] {\scriptstyle k_e = \text{elimination rate constant}} \\[0.2em] {\scriptstyle D_{abs} = \text{absorbed dose per administration}} \\[0.2em] {\scriptstyle V_{max},\, K_m = \text{Michaelis-Menten saturation parameters}} \end{gathered}
Powers PK charting across every compound

A unified pharmacokinetic engine that handles injectable oils, orals, peptides, and testosterone cream — each modelled accurately for its specific delivery route, with route-specific absorption and elimination kinetics.

Androgenic Load

Androgen:Estrogen Asymmetric Additive Formula

Total Androgenic Load (A):A=max ⁣(1,  T+β(DHTDHTref)ϕ)ϕ={1if DHT<DHTref(DHTDHTref)k1if DHTDHTrefA:E Ratio:A:E=AE2Parameters:T=Total Testosterone (ng/dL)E2=Estradiol (pg/mL)DHT=Dihydrotestosterone (ng/dL)DHTref=calibrated AR reference midpointβ=combined AR potency=αaff×αresαaff=AR binding affinity relative to Tαres=AR receptor residence time relative to Tk=paracrine amplification exponent (above DHTref)\begin{gathered} \text{Total Androgenic Load (A):} \\[0.3em] A = \max\!\left(1,\; T + \beta \cdot (DHT - DHT_{ref}) \cdot \phi\right) \\[0.8em] \phi = \begin{cases} 1 & \text{if } DHT < DHT_{ref} \\[0.3em] \left(\dfrac{DHT}{DHT_{ref}}\right)^{k-1} & \text{if } DHT \geq DHT_{ref} \end{cases} \\[1.5em] \text{A:E Ratio:} \\[0.3em] \text{A:E} = \dfrac{A}{E_2} \\[1.5em] {\scriptstyle \text{Parameters:}} \\[0.2em] {\scriptstyle T = \text{Total Testosterone (ng/dL)}} \\[0.2em] {\scriptstyle E_2 = \text{Estradiol (pg/mL)}} \\[0.2em] {\scriptstyle DHT = \text{Dihydrotestosterone (ng/dL)}} \\[0.2em] {\scriptstyle DHT_{ref} = \text{calibrated AR reference midpoint}} \\[0.2em] {\scriptstyle \beta = \text{combined AR potency} = \alpha_{aff} \times \alpha_{res}} \\[0.2em] {\scriptstyle \alpha_{aff} = \text{AR binding affinity relative to T}} \\[0.2em] {\scriptstyle \alpha_{res} = \text{AR receptor residence time relative to T}} \\[0.2em] {\scriptstyle k = \text{paracrine amplification exponent (above } DHT_{ref}\text{)}} \end{gathered}
Powers your A:E ratio analysis

Incorporates DHT's ~10× androgen receptor potency and paracrine amplification — beyond T:E alone. Handles suppressed and elevated DHT non-linearly, producing a truer picture of total androgenic load.

Cardio-Metabolic Risk

The Cardio-Metabolic Cascade Model

Cardio-Metabolic Risk Gate:Ω=mM1 ⁣[χm=Cgreen]M={SBP,  DBP,  ApoB,  CRP,  BF%,  HbA1c}Ψ=1 ⁣[mM:σ(χm)σred]Dynamic Hematocrit Reference Bounds:Hgreen=Hgreenstd+ΩΔHgreenHamber=Hamberstd+ΩΔHamberFull Cascade:M    Ω,Ψ    (Hgreen,Hamber)    f()    χHct    σ-cap    χTM=set of six cardiovascular biomarkersΩ=elite CV gate(1m:χm=green)Ψ=red override flag(1m:χm=red)Hgreen,Hamber=dynamic Hct upper bounds, gated by ΩΔHgreen,ΔHamber=elite threshold incrementsTcaplow,Tcaphigh=testosterone buffer zone boundariesχT(0)=raw testosterone colour prior to cascade\begin{gathered} \text{Cardio-Metabolic Risk Gate:} \\[0.3em] \Omega = \prod_{m \in \mathcal{M}} \mathbf{1}\!\left[\chi_m = \mathcal{C}_{\text{green}}\right] \\[0.8em] \mathcal{M} = \{\text{SBP},\; \text{DBP},\; \text{ApoB},\; \text{CRP},\; \text{BF\%},\; \text{HbA1c}\} \\[0.8em] \Psi = \mathbf{1}\!\left[\exists\, m \in \mathcal{M} : \sigma(\chi_m) \geq \sigma_{\text{red}}\right] \\[1.5em] \text{Dynamic Hematocrit Reference Bounds:} \\[0.3em] H_{\text{green}}^{*} = H_{\text{green}}^{\text{std}} + \Omega \cdot \Delta H_{\text{green}} \\[0.5em] H_{\text{amber}}^{*} = H_{\text{amber}}^{\text{std}} + \Omega \cdot \Delta H_{\text{amber}} \\[1.5em] \text{Full Cascade:} \\[0.3em] \mathcal{M} \;\xrightarrow{\;\Omega,\,\Psi\;}\; \left(H_{\text{green}}^{*},\, H_{\text{amber}}^{*}\right) \;\xrightarrow{\;f(\cdot)\;}\; \chi_{\text{Hct}} \;\xrightarrow{\;\sigma\text{-cap}\;}\; \chi_T \\[1.5em] {\scriptstyle \mathcal{M} = \text{set of six cardiovascular biomarkers}} \\[0.2em] {\scriptstyle \Omega = \text{elite CV gate} \;(\mathbf{1} \iff \forall\, m: \chi_m = \text{green})} \\[0.2em] {\scriptstyle \Psi = \text{red override flag} \;(\mathbf{1} \iff \exists\, m: \chi_m = \text{red})} \\[0.2em] {\scriptstyle H_{\text{green}}^{*},\, H_{\text{amber}}^{*} = \text{dynamic Hct upper bounds, gated by } \Omega} \\[0.2em] {\scriptstyle \Delta H_{\text{green}},\, \Delta H_{\text{amber}} = \text{elite threshold increments}} \\[0.2em] {\scriptstyle T_{\text{cap}}^{\text{low}},\, T_{\text{cap}}^{\text{high}} = \text{testosterone buffer zone boundaries}} \\[0.2em] {\scriptstyle \chi_T^{(0)} = \text{raw testosterone colour prior to cascade}} \end{gathered}
Powers your biomarker score & T range

Six markers — Blood Pressure, ApoB, CRP, Body Fat, HbA1c — cascade into your hematocrit threshold, which sets your personalised testosterone reference range. Lean, metabolically healthy users earn a higher ceiling.

Cream Absorption

Michaelis-Menten Saturation Model

Michaelis-Menten absorption:A=VmaxDKm+DEster effectiveness adjustment:E=Aεˉcyp/enth=1εˉcyp/enthVmaxDKm+Dreff(D)=VmaxKm+D0asDlimD0reff=VmaxKm,ED=Km=Vmax2εˉcyp/enthA=mg testosterone base absorbedD=weekly applied dose (mg/W)E=injectable ester equivalent (mg/W)εˉcyp/enth=mean ester effectivenessVmax=site-specific absorption ceilingKm=site-specific half-saturation dosereff=effective absorption rate\begin{gathered} \text{Michaelis-Menten absorption:} \\[0.8em] A = \dfrac{V_{\max} \cdot D}{K_m + D} \\[1.8em] \text{Ester effectiveness adjustment:} \\[0.8em] E = \dfrac{A}{\bar{\varepsilon}_{\,\text{cyp/enth}}} = \dfrac{1}{\bar{\varepsilon}_{\,\text{cyp/enth}}} \cdot \dfrac{V_{\max} \cdot D}{K_m + D} \\[1.8em] r_{\text{eff}}(D) = \dfrac{V_{\max}}{K_m + D} \longrightarrow 0 \quad \text{as} \quad D \to \infty \\[1.8em] \lim_{D \to 0}\, r_{\text{eff}} = \dfrac{V_{\max}}{K_m} \quad,\qquad E\big|_{D=K_m} = \dfrac{V_{\max}}{2\,\bar{\varepsilon}_{\,\text{cyp/enth}}} \\[1.8em] {\scriptstyle A = \text{mg testosterone base absorbed}} \\[0.2em] {\scriptstyle D = \text{weekly applied dose (mg/W)}} \\[0.2em] {\scriptstyle E = \text{injectable ester equivalent (mg/W)}} \\[0.2em] {\scriptstyle \bar{\varepsilon}_{\,\text{cyp/enth}} = \text{mean ester effectiveness}} \\[0.2em] {\scriptstyle V_{\max} = \text{site-specific absorption ceiling}} \\[0.2em] {\scriptstyle K_m = \text{site-specific half-saturation dose}} \\[0.2em] {\scriptstyle r_{\text{eff}} = \text{effective absorption rate}} \end{gathered}
Powers transdermal testosterone modelling

Cream absorption is non-linear — skin transport proteins and follicular shunts saturate at higher doses. Modelled with site-specific bioavailability for scrotum, abdomen, thighs, and shoulders.

Free Testosterone

Zakharov Nonlinear Allosteric Model

Complete System Equation:f([T]free)=[T]total[T]freeKalb[Alb][T]free[SHBG]α+2αβ1+2α+αβSHBG Allosteric Binding:θSHBG=α+2αβ1+2α+αβα=Ka1[T]free,β=Ka2[T]freeNewton-Raphson Solution:[T]free(n+1)=[T]free(n)+f([T]free(n))f([T]free(n))\begin{gathered} \text{Complete System Equation:} \\[0.3em] f([T]_{free}) = \\[0.5em] [T]_{total} - [T]_{free} - K_{alb}[Alb][T]_{free} - [SHBG]\dfrac{\alpha + 2\alpha\beta}{1 + 2\alpha + \alpha\beta} \\[1.2em] \text{SHBG Allosteric Binding:} \\[0.3em] \theta_{SHBG} = \dfrac{\alpha + 2\alpha\beta}{1 + 2\alpha + \alpha\beta} \\[0.5em] \alpha = K_{a1}[T]_{free}, \quad \beta = K_{a2}[T]_{free} \\[1.2em] \text{Newton-Raphson Solution:} \\[0.3em] [T]_{free}^{(n+1)} = [T]_{free}^{(n)} + \dfrac{f([T]_{free}^{(n)})}{f'([T]_{free}^{(n)})} \end{gathered}
Powers your Free T calculation

Reflects how SHBG dimers actually behave under elevated androgenic load. Significantly more accurate than Vermeulen — the population-level standard — at TRT and supraphysiological levels.

Notation simplified for display. Full derivations and parameters documented in the in-app knowledge base.

>Engine readout complete

The engine lives inside TRT Plus.