Next generation risk assessment for occupational chemical safety - A real world example with sodium-2-hydroxyethane sulfonate.

Next Generation Risk Assessment (NGRA) is an exposure-led approach to safety assessment that uses New Approach Methodologies (NAMs). Application of NGRA has been largely restricted to assessments of consumer use of cosmetics and is not currently implemented in occupational safety assessments, e.g. under EU REACH. By contrast, a large proportion of regulatory worker safety assessments are underpinned by toxicological studies using experimental animals. Consequently, occupational safety assessment represents an area that would benefit from increasing application of NGRA to safety decision making. Here, a workflow for conducting NGRA under an occupational safety context was developed, which is illustrated with a case study chemical; sodium 2-hydroxyethane sulphonate (sodium isethionate or SI). Exposures were estimated using a standard occupational exposure model following a comprehensive life cycle assessment of SI and considering factory-specific data. Outputs of this model were then used to estimate internal exposures using a Physiologically Based Kinetic (PBK) model, which was constructed with SI specific Absorption, Distribution, Metabolism and Excretion (ADME) data. PBK modelling indicated a worst-case plasma maximum concentration (Cmax) of 0.8 µM across the SI life cycle. SI bioactivity was assessed in a battery of NAMs relevant to systemic, reproductive, and developmental toxicity; a cell stress panel, high throughput transcriptomics in three cell lines (HepG2, HepaRG and MCF-7 cells), pharmacological profiling and specific assays relating to developmental toxicity (Reprotracker and devTOX quickPredict). Points of Departure (PoDs) for SI ranged from 104 to 5044 µM. Cmax values obtained from PBK modelling of occupational exposures to SI were compared with PoDs from the bioactivity assays to derive Bioactivity Exposure Ratios (BERs) which demonstrated the safety for workers exposed to SI under current levels of factory specific risk management. In summary, the tiered and iterative workflow developed here represents an opportunity for integrating non animal approaches for a large subset of substances for which systemic worker safety assessment is required. Such an approach could be followed to ensure that animal testing is only conducted as a "last resort" e.g. under EU REACH.

Quantitative phosphoproteomics reveal cellular responses from caffeine, coumarin and quercetin in treated HepG2 cells.

Protein phosphorylation is the most common type of post-translational modification where serine, threonine or tyrosine are reversibly bound to the phosphate group of ATP in a reaction catalyzed by protein kinases. Phosphorylation plays an important role in regulation of cell homeostasis, including but not limited to signal perception and transduction, gene expression and function of proteins. Protein phosphorylation happens on a fast time scale and represents an energy-efficient way for the cell to adapt to exposure to chemical stressors. To understand the cascade of cellular signaling induced by exposure to chemicals, we have exposed HepG2 cells to three chemicals with different modes of action, namely, caffeine, coumarin, and quercetin in a concentration and time response manner. Significantly upregulated and downregulated phosphosites were screened to analyze the activation/deactivation of signaling pathways by protein kinases. In total, 69, 44 and 12 signaling pathways were found enriched in caffeine, coumarin and quercetin treated cells, respectively, of which 9 pathways were co-enriched with 11 jointly responded kinases. Among identified co-responded kinases, CDK1, MAPK1 and MAPK3 play important roles in cell cycle and insulin signaling pathways. Quantitative phosphoproteomics can sensitively distinguish the effects of different chemicals on cells, allowing the assessment of chemical safety through changes in substrates and metabolic pathways at the cellular level, which is important for the development of non-animal approaches for chemical safety assessment.

Proteomic analysis of the cellular response to a potent sensitiser unveils the dynamics of haptenation in living cells.

Haptenation of model nucleophiles, representing the key MIE in skin sensitisation, is routinely measured in chemico to provide data for skin allergy risk assessment. Better understanding of the dynamics of haptenation in human skin could provide the metrics required to improve determination of sensitiser potency for risk assessment of chemicals. We have previously demonstrated the applicability and sensitivity of the dual stable isotope labelling approach to detect low level haptenation in complex mixtures of proteins. In the present study, we investigated haptenation in a relevant living cell model over time at a subtoxic concentration. DNCB, an extremely potent sensitiser, caused minimal changes in overall protein differential expression in HaCaT cells and haptenated approximately 0.25 % of all available nucleophiles when applied at a subtoxic concentration (10µM) for 4 h. The data shows that the maximum level of haptenation occurs at 2 h and that DNCB, whilst being a promiscuous hapten, shows a preference for Cys residues, despite the considerably higher concentration of amine-based nucleophiles. Although a proportion of highly abundant proteins were haptenated, numerous haptenated sites were also detected on low abundant proteins. Certain proteins were modified at residues buried deep inside the protein structure which are less accessible to haptenation compared with surface exposed nucleophiles. The microenvironment of the buried residues may be a result of several factors influencing the reactivity of both the target nucleophile and the hapten.

The free energy landscape of the oncogene protein E7 of human papillomavirus type 16 reveals a complex interplay between ordered and disordered regions.

When present, structural disorder makes it very challenging to characterise the conformational properties of proteins. This is particularly the case of proteins, such as the oncogene protein E7 of human papillomavirus type 16, which contain both ordered and disordered domains, and that can populate monomeric and oligomeric states under physiological conditions. Nuclear magnetic resonance (NMR) spectroscopy is emerging as a powerful method to study these complex systems, most notably in combination with molecular dynamics simulations. Here we use NMR chemical shifts and residual dipolar couplings as structural restraints in replica-averaged molecular dynamics simulations to determine the free energy landscape of E7. This landscape reveals a complex interplay between a folded but highly dynamical C-terminal domain and a disordered N-terminal domain that forms transient secondary and tertiary structures, as well as an equilibrium between a high-populated (98%) dimeric state and a low-populated (2%) monomeric state. These results provide compelling evidence of the complex conformational heterogeneity associated with the behaviour and interactions of this disordered protein associated with disease.

The RNF168 paralog RNF169 defines a new class of ubiquitylated histone reader involved in the response to DNA damage.

Site-specific histone ubiquitylation plays a central role in orchestrating the response to DNA double-strand breaks (DSBs). DSBs elicit a cascade of events controlled by the ubiquitin ligase RNF168, which promotes the accumulation of repair factors such as 53BP1 and BRCA1 on the chromatin flanking the break site. RNF168 also promotes its own accumulation, and that of its paralog RNF169, but how they recognize ubiquitylated chromatin is unknown. Using methyl-TROSY solution NMR spectroscopy and molecular dynamics simulations, we present an atomic resolution model of human RNF169 binding to a ubiquitylated nucleosome, and validate it by electron cryomicroscopy. We establish that RNF169 binds to ubiquitylated H2A-Lys13/Lys15 in a manner that involves its canonical ubiquitin-binding helix and a pair of arginine-rich motifs that interact with the nucleosome acidic patch. This three-pronged interaction mechanism is distinct from that by which 53BP1 binds to ubiquitylated H2A-Lys15 highlighting the diversity in site-specific recognition of ubiquitylated nucleosomes.

Molecular Dynamics Simulation of Tau Peptides for the Investigation of Conformational Changes Induced by Specific Phosphorylation Patterns.

The Tau protein plays an important role due to its biomolecular interactions in neurodegenerative diseases. The lack of stable structure and various posttranslational modifications such as phosphorylation at various sites in the Tau protein pose a challenge for many experimental methods that are traditionally used to study protein folding and aggregation. Atomistic molecular dynamics (MD) simulations can help around deciphering relationship between phosphorylation and various intermediate and stable conformations of the Tau protein which occur on longer timescales. This chapter outlines protocols for the preparation, execution, and analysis of all-atom MD simulations of a 21-amino acid-long phosphorylated Tau peptide with the aim of generating biologically relevant structural and dynamic information. The simulations are done in explicit solvent and starting from nearly extended configurations of the peptide. The scaled MD method implemented in AMBER14 was chosen to achieve enhanced conformational sampling in addition to a conventional MD approach, thereby allowing the characterization of folding for such an intrinsically disordered peptide at 293 K. Emphasis is placed on the analysis of the simulation trajectories to establish correlations with NMR data (i.e., chemical shifts and NOEs). Finally, in-depth discussions are provided for commonly encountered problems.

Mapping the Protein Fold Universe Using the CamTube Force Field in Molecular Dynamics Simulations.

It has been recently shown that the coarse-graining of the structures of polypeptide chains as self-avoiding tubes can provide an effective representation of the conformational space of proteins. In order to fully exploit the opportunities offered by such a 'tube model' approach, we present here a strategy to combine it with molecular dynamics simulations. This strategy is based on the incorporation of the 'CamTube' force field into the Gromacs molecular dynamics package. By considering the case of a 60-residue polyvaline chain, we show that CamTube molecular dynamics simulations can comprehensively explore the conformational space of proteins. We obtain this result by a 20 µs metadynamics simulation of the polyvaline chain that recapitulates the currently known protein fold universe. We further show that, if residue-specific interaction potentials are added to the CamTube force field, it is possible to fold a protein into a topology close to that of its native state. These results illustrate how the CamTube force field can be used to explore efficiently the universe of protein folds with good accuracy and very limited computational cost.