Application of the DILIsym® Quantitative Systems Toxicology drug-induced liver injury model to evaluate the carcinogenic hazard potential of acetaminophen.

In 2019, the California Office of Environmental Health Hazard Assessment (OEHHA) initiated a review of the carcinogenic hazard potential of acetaminophen. The objective of the analysis herein was to inform this review by assessing whether variability in patient baseline characteristics (e.g. baseline glutathione (GSH) levels, pharmacokinetics, and capacity of hepatic antioxidants) leads to potential differences in carcinogenic hazard potential at different dosing schemes: maximum labeled doses of 4 g/day, repeated doses above the maximum labeled dose (>4-12 g/day), and acute overdoses of acetaminophen (>15 g). This was achieved by performing simulations of acetaminophen exposure in thousands of diverse virtual patients scenarios using the DILIsym® Quantitative Systems Toxicology (QST) model. Simulations included assessments of the dose and exposure response for toxicity and mode of cell death based on evaluations of the kinetics of changes of: GSH, N-acetyl-p-benzoquinone-imine (NAPQI), protein adducts, mitochondrial dysfunction, and hepatic cell death. Results support that, at therapeutic doses, cellular GSH binds to NAPQI providing sufficient buffering capacity to limit protein adduct formation and subsequent oxidative stress. Simulations evaluating repeated high-level supratherapeutic exposures or acute overdoses indicate that cell death precedes DNA damage that could result in carcinogenicity and thus acetaminophen does not present a carcinogenicity hazard to humans at any dose.

A critical review of the acetaminophen preclinical carcinogenicity and tumor promotion data and their implications for its carcinogenic hazard potential.

In 2019 the California Office of Environmental Health Hazard Assessment (OEHHA) initiated a review of the carcinogenic hazard potential of acetaminophen, including an assessment of the long-term rodent carcinogenicity and tumor initiation/promotion studies. The objective of the analysis herein was to inform this review process with a weight-of-evidence assessment of these studies and an assessment of the relevance of these models to humans. In most of the 14 studies, there were no increases in the incidences of tumors in any organ system. In the few studies in which an increase in tumor incidence was observed, there were factors such as absence of a dose response and a rodent-specific tumor supporting that these findings are not relevant to human hazard identification. In addition, we performed qualitative analysis and quantitative simulations of the exposures to acetaminophen and its metabolites and its toxicity profile; the data support that the rodent models are toxicologically relevant to humans. The preclinical carcinogenicity results are consistent with the broader weight of evidence assessment and evaluations of multiple international health authorities supporting that acetaminophen is not a carcinogenic hazard.

Benchmarking renin suppression and blood pressure reduction of direct renin inhibitor imarikiren through quantitative systems pharmacology modeling.

Multiple classes of antihypertensive drugs inhibit components of the renin-angiotensin-aldosterone system (RAAS). The primary physiological effector of the RAAS is angiotensin II (AngII) bound to the AT1 receptor (AT1-bound AngII). There is a strong non-linear feedback from AT1-bound AngII on renin secretion. Since AT1-bound AngII is not readily measured experimentally, plasma renin concentration (PRC) and/or activity (PRA) are typically measured to indicate RAAS suppression. We investigated the RAAS suppression of imarikiren hydrochloride (TAK-272; SCO-272), a direct renin inhibitor currently under clinical development. We employed a previously developed quantitative system pharmacology (QSP) model to benchmark renin suppression and blood pressure regulation with imarikiren compared to other RAAS therapies. A pharmacokinetic (PK) model of imarikiren was linked with the existing QSP model, which consists of a mechanistic representation of the RAAS pathway coupled with a model of blood pressure regulation and volume homeostasis. The PK and pharmacodynamic effects of imarikiren were calibrated by fitting drug concentration, PRA, and PRC data, and trough AT1-bound AngII suppression was simulated. We also prospectively simulated expected mean arterial pressure reduction in a cohort of hypertensive virtual patients. These predictions were benchmarked against predictions for several other (previously calibrated) RAAS monotherapies and dual-RAAS therapies. Our analysis indicates that low doses (5-10 mg) of imarikiren are comparable to current RAAS therapies, and at higher doses (25-200 mg), RAAS suppression may be equivalent to existing dual-RAAS combinations (at registered doses). This study illustrates application of QSP modeling to predict phase II endpoints from phase I data.

Multiscale Mathematical Model of Drug-Induced Proximal Tubule Injury: Linking Urinary Biomarkers to Epithelial Cell Injury and Renal Dysfunction.

Drug-induced nephrotoxicity is a major cause of acute kidney injury, and thus detecting the potential for nephrotoxicity early in the drug development process is critical. Various urinary biomarkers exhibit different patterns following drug-induced injury, which may provide greater information than traditional biomarkers like serum creatinine. In this study, we developed a multiscale quantitative systems pharmacology model relating drug exposure to proximal tubule (PT) epithelial cell injury and subsequently to expression of multiple urinary biomarkers and organ-level functional changes. We utilized urinary kidney injury molecule-1 (Kim-1), alpha glutathione S-transferase, albumin (αGST), glucose, and urine volume time profiles as well as serum creatinine and histopathology data obtained from rats treated with the nephrotoxicant cisplatin to develop the model. Although the model was developed using single-dose response to cisplatin, the model predicted the serum creatinine response to multidose cisplatin regimens. Further, using only the urinary Kim-1 response to gentamicin (a nephrotoxicant with a distinctly different injury time course than cisplatin), the model detected and predicted mild to moderate PT injury, as confirmed with histopathology, even when serum creatinine was unchanged. Thus, the model is generalizable, and can be used to deconvolute the underlying degree and time course of drug-induced PT injury and renal dysfunction from a small number of urinary biomarkers, and may provide a tool to determine optimal dosing regimens that minimize renal injury.

Primary proximal tubule hyperreabsorption and impaired tubular transport counterregulation determine glomerular hyperfiltration in diabetes: a modeling analysis.

Glomerular hypertension and hyperfiltration in early diabetes are associated with development and progression of diabetic kidney disease. The tubular hypothesis of diabetic hyperfiltration proposes that it is initiated by a primary increase in sodium (Na) reabsorption in the proximal tubule (PT) and the resulting tubuloglomerular feedback (TGF) response and lowering of Bowman space pressure (PBow). Here we utilized a mathematical model of the human kidney to investigate over acute and chronic timescales the mechanisms responsible for the magnitude of the hyperfiltration response. The model implicates that the primary hyperreabsorption of Na in the PT produces a Na imbalance that is only partially restored by the hyperfiltration induced by TGF and changes in PBow Thus secondary adaptations are needed to restore Na balance. This may include neurohumoral transport regulation and/or pressure-natriuresis (i.e., the decrease in Na reabsorption in response to increased renal perfusion pressure). We explored the role of each tubular segment in contributing to this compensation and the consequences of impairment in tubular compensation. The simulations indicate that impaired secondary downregulation of transport potentiated the rise in glomerular hypertension and hyperfiltration needed to restore Na balance at a given level of primary PT hyperreabsorption. Therefore, we propose for the first time that both the extent of primary PT hyperreabsorption and the degree of impairment of the distal tubular responsiveness to regulatory signals determine the level of glomerular hypertension and hyperfiltration in the diabetic kidney, thereby extending the tubule-centric concept of diabetic hyperfiltration and potential therapeutic approaches beyond the proximal tubule.

Conformational properties of interacting neurofilaments: Monte Carlo simulations of cylindrically grafted apposing neurofilament brushes.

Neurofilaments are essential cytoskeletal filaments that impart mechanical stability to axons. They are mostly assembled from three neurofilament proteins that form the core of the filament and its sidearms. Adjacent neurofilaments interact with each other through their apposing sidearms and attain unique conformations depending on the ionic condition, phosphorylation state, and interfilament separations. To understand the conformational properties of apposing sidearms under various conditions and gain insight into interfilament interactions, we performed Monte Carlo simulations of neurofilament pairs. We employed a sequence-based coarse-grained model of apposing NF sidearms that are end-tethered to cylindrical geometries according to the stoichiometry of the three neurofilament subunits. Monte Carlo simulations were conducted under different conditions such as phosphorylation state, ionic condition, and interfilament separations. Under salt-free conditions, apposing sidearms are found to adopt mutually excluding stretched but bent away conformations that are reminiscent of a repulsive type of interaction. Under physiological conditions, apposing sidearms are found to be in a coiled conformation, suggesting a short-range steric repulsive type of interaction. Increased sidearm mutual interpenetration and a simultaneous decrease in the individual brush heights were observed as the interfilament separation was reduced from 60 to 40 nm. The observed conformations suggest entropic interaction as a likely mechanism for sidearm-mediated interfilament interactions under physiological conditions.

Effects of molecular model, ionic strength, divalent ions, and hydrophobic interaction on human neurofilament conformation.

The present study examines the effects of the model dependence, ionic strength, divalent ions, and hydrophobic interaction on the structural organization of the human neurofilament (NF) brush, using canonical ensemble Monte Carlo (MC) simulations of a coarse-grained model with the amino-acid resolution. The model simplifies the interactions between the NF core and the sidearm or between the sidearms by the sum of excluded volume, electrostatic, and hydrophobic interactions, where both monovalent salt ions and solvents are implicitly incorporated into the electrostatic interaction potential. Several important observations are made from the MC simulations of the coarse-grained model NF systems. First, the mean-field type description of monovalent salt ions works reasonably well in the NF system. Second, the manner by which the NF sidearms are arranged on the surface of the NF backbone core has little influence on the lateral extension of NF sidearms. Third, the lateral extension of the NF sidearms is highly affected by the ionic strength of the system: at low ionic strength, NF-M is most extended but at high ionic strength, NF-H is more stretched out because of the effective screening of the electrostatic interaction. Fourth, the presence of Ca(2+) ions induces the attraction between negatively charged residues, which leads to the contraction of the overall NF extension. Finally, the introduction of hydrophobic interaction does not change the general structural organization of the NF sidearms except that the overall extension is contracted.

Computational investigation of the key factors affecting the second stage activation mechanisms of domain II m-calpain.

The unique conformation of the active site in calpains along with the implication of their role in several diseases has prompted widespread research interest in the scientific community. Structural studies devoted to m- and µ-calpains have proposed a two-stage calcium-dependent activation mechanism for calpains. In this work, we performed conventional and targeted molecular dynamics simulations to investigate global and local changes in the structure of the protease core of m-calpain upon calcium binding. Simulations were performed on the protease core of calcium free (pdbid: 1kfu) and calcium bound (pdbid: 3df0) m-calpain with and without the presence of calcium ions. Our results indicate that the inactive, open conformation of the protease core does not transform into the active, closed conformation simply upon removal of constraints from the neighbor domains. The role of other factors, including calcium binding and the subsequent formation of an Arg94-Glu305 inter-domain salt bridge and the change in the orientation of Trp288 side chain, in the activation of the protease core is elicited.

Expansion of neurofilament medium C terminus increases axonal diameter independent of increases in conduction velocity or myelin thickness.

Maturation of the peripheral nervous system requires specification of axonal diameter, which, in turn, has a significant influence on nerve conduction velocity. Radial axonal growth initiates with myelination, and is dependent upon the C terminus of neurofilament medium (NF-M). Molecular phylogenetic analysis in mammals suggested that expanded NF-M C termini correlated with larger-diameter axons. We used gene targeting and computational modeling to test this new hypothesis. Increasing the length of NF-M C terminus in mice increased diameter of motor axons without altering neurofilament subunit stoichiometry. Computational modeling predicted that an expanded NF-M C terminus extended farther from the neurofilament core independent of lysine-serine-proline (KSP) phosphorylation. However, expansion of NF-M C terminus did not affect the distance between adjacent neurofilaments. Increased axonal diameter did not increase conduction velocity, possibly due to a failure to increase myelin thickness by the same proportion. Failure of myelin to compensate for larger axonal diameters suggested a lack of plasticity during the processes of myelination and radial axonal growth.

Neurofilament stoichiometry simulations during neurodegeneration suggest a remarkable self-sufficient and stable in vivo protein structure.

BACKGROUND: Neurofilaments (Nfs) are protein biomarkers of neurodegeneration in human disease. There is in vivo evidence of changes of the Nf stoichiometry in the cerebrospinal fluid (CSF) of patients. The protein-structural implications of these findings are not known but may be assessed indirectly using simulations studies. METHODS: Monte Carlo simulations were performed using a coarse-grained model of a Nf brush. Based on the published in vivo CSF data the tested Nf stoichiometries (NfL:NfM:NfH) were 16:11:4 for multiple system atrophy (MSA), 24:5:2 for relapsing remitting multiple sclerosis (RRMS), and 30:0:1 for clinically isolated syndromes (CIS). Simulations were performed in a wide range of ionic strength (1 mM-100 mM) for dephosphorylated and phosphorylated NF isoforms. RESULTS: At lower ionic strengths (1 mM, 10 mM), NfM is the main determinant for the radius of gyration (R(g)) ranging from ≈15 nm in the dephosphorylated state at 10 mM ionic strength to ≈27 nm at 1mM ionic strength if fully phosphorylated. At high ionic strength (100mM) NfH becomes the main determinant with R(g) of 14.8±0.2 nm if dephosphorylated and 15±0.2 nm if phosphorylated. There was no significant difference in the structures of the three Nf sidearms for MSA, RRMS or CIS. CONCLUSION: Large changes of the in vivo Nf stoichiometry have only little effect on the simulated structure of Nf sidearms independent of phosphorylation and ionic strength. This suggests that the axonal cytoskeleton is remarkably stable, possibly relying on NfL which forms a dense brush around the Nf backbone and virtually excludes NfM and NfH from the core region, such that the dropout of NfM and NfH can be dealt with structurally.

Phosphorylation-mediated conformational changes in the mouse neurofilament architecture: insight from a neurofilament brush model.

Neurofilaments (NFs) are important cytoskeletal filaments that consist of long flexible C-terminal tails that are abundant with charges. The tails attain additional negative charges through serine phosphorylation of Lys-Ser-Pro (KSP) repeat motifs that are particularly found in neurofilament heavy (NF-H) and neurofilament medium (NF-M) proteins. These side-arm protrusions mediate the interaction between neighboring filaments and maintain axonal diameter. However, the precise role of NF proteins and their phosphorylation in regulating interfilament distances and axonal diameter still remains unclear. In this regard, a recent gene replacement study revealed that the phosphorylation of mouse NF-M KSP repeats does not affect axonal cytoarchitecture, challenging the conventional viewpoint on the role of NF phosphorylation. To better understand the effect of phosphorylation, particularly NF-M phosphorylation, we applied a computational method to reveal phosphorylation-mediated conformational changes in mouse NF architecture. We employed a three-dimensional sequence-based coarse-grained NF brush model to perform Monte Carlo simulations of mouse NF by using the sequence and stoichiometry of mouse NF proteins. Our result shows that the phosphorylation of mouse NF-M does not change the radial extension of NF-M side arms under a salt-free condition and in ionic solution, highlighting a structural factor that supports the notion that NF-M KSP phosphorylation has no effect on the axonal diameter of mouse. On the other hand, significant phosphorylation-mediated conformational changes were found in NF-H side arms under the salt-free condition, while the changes in ionic solution are not significant. However, NF-H side arms are found at the periphery of mouse NF architecture, implying a role in linking neighboring filaments.

Structural properties of neurofilament sidearms: sequence-based modeling of neurofilament architecture.

Neurofilaments (NFs) are essential cytoskeletal filaments that impart mechanical integrity to nerve cells. They are assembled from three distinct molecular mass proteins that bind to each other to form a 10-nm-diameter filamentous rod with sidearm extensions. The sidearms are considered to play a critical role in modulating interfilament spacing and axonal caliber. However, the precise mechanism by which NF protrusions regulate axonal diameter remains to be well understood. In particular, the role played by individual NF protrusions in specifying interfilament distances is yet to be established. To gain insight into the role of individual proteins, we investigated the structural organization of NF architecture under different phosphorylation conditions. To this end, a physically motivated sequence-based coarse-grain model of NF brush has been developed based on the three-dimensional architecture of NFs. The model incorporates the charge distribution of sidearms, including charges from the phosphorylation sites corresponding to Lys-Ser-Pro repeat motifs. The model also incorporates the proper grafting of the real NF sidearms based on the stoichiometry of the three subunits. The equilibrium structure of the NF brush is then investigated under different phosphorylation conditions. The phosphorylation of NF modifies the structural organization of sidearms. Upon phosphorylation, a dramatic change involving a transformation from a compact conformation to an extended conformation is found in the heavy NF (NF-H) protein. However, in spite of extensive phosphorylation sites present in the NF-H subunit, the tails of the medium NF subunit are found to be more extended than the NF-H sidearms. This supports the notion that medium NF protrusions are critical in regulating NF spacings and, hence, axonal caliber.

Intrinsic bending and structural rearrangement of tubulin dimer: molecular dynamics simulations and coarse-grained analysis.

Microtubules are long polymers of alphabeta-tubulin heterodimers. They undergo a process known as dynamic instability, in which the ends of a microtubule switch stochastically between phases of slow growth and rapid shrinkage. The molecular mechanisms inducing the depolymerization of microtubules were attributed to the hydrolysis of the guanosine triphosphate (GTP) nucleotide bound to the beta-tubulin. The hydrolysis of GTP is thought to cause microtubule instability by promoting outward curving of the protofilaments constituting the microtubule lattice. The bending of protofilaments is associated with the structural transformation of a tubulin dimer from straight to curved conformations. However, the nature of intrinsic bending of the dimer remains elusive. This study uses molecular dynamics (MD) simulations and coarse-grained analysis to reveal the intrinsic bending, as well as the local structural rearrangements, of the unassembled tubulin dimer as the dimer relaxes from its lattice-constrained, straight conformation of a zinc-induced tubulin sheet. The effect of the nucleotide state on dimer-bending is investigated by the introduction of gamma-phosphate into the beta-tubulin to form GTP-bound tubulin. In agreement with recent experimental studies that proposed nucleotide-independent curved conformations, both guanosine diphosphate (GDP)-bound and GTP-bound tubulin dimers were found to have curved conformations, but with a tendency toward smaller bending in the GTP-tubulin than in the GDP-tubulin. The perturbation induced through the introduction of gamma-phosphate is posited to play a role in straightening the intradimer bending. The local structural rearrangements of GDP-tubulin because of the bending mode of motion of the dimer reveal that one of the three functional domains, the intermediate domain, exhibits significantly lower bending deformation compared with the others, signifying a dynamic connection to the functionally defined domains.

Mesoscopic modeling of bacterial flagellar microhydrodynamics.

A particle-based hybrid method of elastic network model and smooth-particle hydrodynamics has been employed to describe the propulsion of bacterial flagella in a viscous hydrodynamic environment. The method explicitly models the two aspects of bacterial propulsion that involve flagellar flexibility and long-range hydrodynamic interaction of low-Reynolds-number flow. The model further incorporates the molecular organization of the flagellar filament at a coarse-grained level in terms of the 11 protofilaments. Each of these protofilaments is represented by a collection of material points that represent the flagellin proteins. A computational model of a single flexible helical segment representing the filament of a bacterial flagellum is presented. The propulsive dynamics and the flow fields generated by the motion of the model filament are examined. The nature of flagellar deformation and the influence of hydrodynamics in determining the shape of deformations are examined based on the helical filament.

Spatially heterogeneous dynamics and the Adam-Gibbs relation in the Dzugutov liquid.

We perform molecular dynamics simulations of a one-component glass-forming liquid and use the inherent structure formalism to test the predictions of the Adam-Gibbs (AG) theory and to explore the possible connection between these predictions and spatially heterogeneous dynamics. We calculate the temperature dependence of the average potential energy of the equilibrium liquid and show that it obeys the Rosenfeld-Tarazona T(3/5) law for low temperature T, while the average inherent structure energy is found to be inversely proportional to temperature at low T, consistent with a Gaussian distribution of potential energy minima. We investigate the shape of the basins around the local minima in configuration space via the average basin vibrational frequency and show that the basins become slightly broader upon cooling. We evaluate the configurational entropy S(conf), a measure of the multiplicity of potential energy minima sampled by the system, and test the validity of the AG relation between S(conf) and the bulk dynamics. We quantify the dynamically heterogeneous motion by analyzing the motion of particles that are mobile on short and intermediate time scales relative to the characteristic bulk relaxation time. These mobile particles move in one-dimensional "strings", and these strings form clusters with a well-defined average cluster size. The AG approach predicts that the minimum size of cooperatively rearranging regions (CRR) of molecules is inversely proportional to S(conf), and recently (Phys. Rev. Lett. 2003, 90, 085506) it has been shown that the mobile-particle clusters are consistent with the CRR envisaged by Adam and Gibbs. We test the possibility that the mobile-particle strings, rather than clusters, may describe the CRR of the Adam-Gibbs approach. We find that the strings also follow a nearly inverse relation with S(conf). We further show that the logarithm of the time when the strings and clusters are maximum, which occurs in the late-beta-relaxation regime of the intermediate scattering function, follows a linear relationship with 1/TS(conf), in agreement with the AG prediction for the relationship between the configurational entropy and the characteristic time for the liquid to undergo a transition to a new configuration. Since strings are the basic elements of the clusters, we propose that strings are a more appropriate measure of the minimum size of a CRR that leads to configurational transitions. That the cluster size also has an inverse relationship with S(conf) may be a consequence of the fact that the clusters are composed of strings.