Relative genotoxicity of polycyclic aromatic hydrocarbons inferred from free energy perturbation approaches.

Utilizing molecular dynamics and free energy perturbation, we examine the relative binding affinity of several covalent polycyclic aromatic hydrocarbon - DNA (PAH-DNA) adducts at the central adenine of NRAS codon-61, a mutational hotspot implicated in cancer risk. Several PAHs classified by the International Agency for Research on Cancer as probable, possible, or unclassifiable as to carcinogenicity are found to have greater binding affinity than the known carcinogen, benzo[a]pyrene (B[a]P). van der Waals interactions between the intercalated PAH and neighboring nucleobases, and minimal disruption of the DNA duplex drive increases in binding affinity. PAH-DNA adducts may be repaired by global genomic nucleotide excision repair (GG-NER), hence we also compute relative free energies of complexation of PAH-DNA adducts with RAD4-RAD23 (the yeast ortholog of human XPC-RAD23) which constitutes the recognition step in GG-NER. PAH-DNA adducts exhibiting the greatest DNA binding affinity also exhibit the least RAD4-RAD23 complexation affinity and are thus predicted to resist the GG-NER machinery, contributing to their genotoxic potential. In particular, the fjord region PAHs dibenzo[a,l]pyrene, benzo[g]chrysene, and benzo[c]phenanthrene are found to have greater binding affinity while having weaker RAD4-RAD23 complexation affinity than their respective bay region analogs B[a]P, chrysene, and phenanthrene. We also find that the bay region PAHs dibenzo[a,j]anthracene, dibenzo[a,c]anthracene, and dibenzo[a,h]anthracene exhibit greater binding affinity and weaker RAD4-RAD23 complexation affinity than B[a]P. Thus, the study of PAH genotoxicity likely needs to be substantially broadened, with implications for public policy and the health sciences. This approach can be broadly applied to assess factors contributing to the genotoxicity of other unclassified compounds.

Regularization of least squares problems in CHARMM parameter optimization by truncated singular value decompositions.

We examine the use of the truncated singular value decomposition and Tikhonov regularization in standard form to address ill-posed least squares problems Ax = b that frequently arise in molecular mechanics force field parameter optimization. We illustrate these approaches by applying them to dihedral parameter optimization of genotoxic polycyclic aromatic hydrocarbon-DNA adducts that are of interest in the study of chemical carcinogenesis. Utilizing the discrete Picard condition and/or a well-defined gap in the singular value spectrum when A has a well-determined numerical rank, we are able to systematically determine truncation and in turn regularization parameters that are correspondingly used to produce truncated and regularized solutions to the ill-posed least squares problem at hand. These solutions in turn result in optimized force field dihedral terms that accurately parameterize the torsional energy landscape. As the solutions produced by this approach are unique, it has the advantage of avoiding the multiple iterations and guess and check work often required to optimize molecular mechanics force field parameters.

Occupational Years of Service and Leukocyte Epigenetic Aging: Relationships in United States Firefighters.

OBJECTIVE: The aim of the study is to examine associations between years of firefighting service and eight chronological age-adjusted measures of blood leukocyte epigenetic age acceleration: Horvath, Hannum, SkinBloodClock, Intrinsic, Extrinsic, PhenoAge, GrimAge, and DNAm telomere length. METHODS: The study used a repeated measures analysis of data from 379 incumbent firefighters from eight career departments and 100 recruit firefighters from two of the departments, across the United States. RESULTS: Incumbent firefighters had on average greater epigenetic age acceleration compared with recruit firefighters, potentially due to the cumulative effect of occupational exposures. However, among incumbent firefighters, additional years of service were associated with epigenetic age deceleration, particularly for GrimAge, a strong predictor of mortality. CONCLUSIONS: Long-term studies with more specific occupational exposure classification are needed to better understand the relationship between years of service and aging biomarkers.