Quantitative bias analysis of the association between subclinical thyroid disease and two perfluoroalkyl substances in a single study.

Exposure to perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) has been associated with the occurrence of thyroid disease in some epidemiologic studies. We hypothesized that in a specific epidemiologic study based on the National Health and Nutrition Examination Survey, the association of subclinical thyroid disease with serum concentration of PFOA and PFOS was due to reverse causality. Thyroid hormone affects glomerular filtration, which in turn affects excretion of PFOA and PFOS. We evaluated this by linking a model of thyroid disease status over the lifetime to physiologically based pharmacokinetic models of PFOA and PFOS. Using Monte Carlo methods, we simulated the target study population and analyzed the data using multivariable logistic regression. The target and simulated populations were similar with respect to age, estimated glomerular filtration rate, serum concentrations of PFOA and PFOS, and prevalence of subclinical thyroid disease. Our findings suggest that in the target study the associations with subclinical hypothyroidism were overstated and the results for subclinical hyperthyroidism were, in general, understated. For example, for subclinical hypothyroidism in men, the reported odds ratio per ln(PFOS) increase was 1.98 (95% CI 1.19-3.28), whereas in the simulated data the bias due to reverse causality gave an odds ratio of 1.19 (1.16-1.23). Our results provide evidence of bias due to reverse causality in a specific cross-sectional study of subclinical thyroid disease with exposure to PFOA and PFOS among adults.

Structurally Linked Dynamics in Lactate Dehydrogenases of Evolutionarily Distinct Species.

We present new findings about how primary and secondary structure affects the role of fast protein motions in the reaction coordinates of enzymatic reactions. Using transition path sampling and committor distribution analysis, we examined the difference in the role of these fast protein motions in the reaction coordinate of lactate dehydrogenases (LDHs) of Apicomplexa organisms Plasmodium falciparum and Cryptosporidium parvum. Having evolved separately from a common malate dehydrogenase ancestor, the two enzymes exhibit several important structural differences, notably a five-amino acid insertion in the active site loop of P. falciparum LDH. We find that these active site differences between the two organisms' LDHs likely cause a decrease in the contribution of the previously determined LDH rate-promoting vibration to the reaction coordinate of P. falciparum LDH compared to that of C. parvum LDH, specifically in the coupling of the rate-promoting vibration and the hydride transfer. This effect, while subtle, directly shows how changes in structure near the active site of LDH alter catalytically important motions. Insights provided by studying these alterations would prove to be useful in identifying LDH inhibitors that specifically target the isozymes of these parasitic organisms.

Systematic review and meta-analysis of epidemiologic data on vaccine response in relation to exposure to five principal perfluoroalkyl substances.

BACKGROUND: Epidemiologic studies of serum per- and polyfluoroalkyl substances (PFAS) and antibody response to vaccines have suggested an adverse association, but the consistency and magnitude of this association remain unclear. OBJECTIVE: The goal of this systematic review was to determine the size of the association between a doubling in perfluoroalkyl substances (PFAS) serum concentration and difference in loge antibody concentration following a vaccine, with a focus on five PFAS: perfluorooctanoic acid (PFOA), perfluorooctane sulfonate (PFOS), perfluorohexane sulfonate (PFHxS), perfluorononanoic acid (PFNA), and perfluorodecanoic acid (PFDA). DATA SOURCE: We conducted online searches of PubMed and Web of Science through May 17, 2022 and identified 14 eligible reports published from 2012 to 2022. STUDY ELIGIBILITY CRITERIA, PARTICIPANTS, AND INTERVENTIONS: We included studies conducted in humans, including mother-child pairs, which examined serum PFAS concentration in relation to serum concentration of antibody to a specific antigen following a vaccine. STUDY APPRAISAL AND SYNTHESIS METHODS: We used the risk of bias assessment for non-randomized studies of exposure and certainty assessment method proposed by Morgan et al. (2019). Using a multilevel meta-regression model, we quantitatively synthesized the data. RESULTS: The 14 reports represented 13 unique groups of subjects; the frequency of studies of a given antibody was Tetanus (n = 7); followed by Diphtheria (6); Measles (4); Rubella (3); Haemophilus influenzae type b and Influenza A H1N1 (2 each); and Hepatitis A, Hepatitis B, Influenza A H2N3, Influenza B, and Mumps (1 each). There were approximately 4,830 unique participants included in the analyses across the 14 reports. The models of coefficients between antibody concentration and the five principal PFAS showed homogeneity of associations across antibody types for each principal PFAS. In the models with all antibodies treated as one type, evidence of effect modification by life stage was present for PFOA and PFOS, and for consistency, all associations were evaluated for all ages and for children. The summary associations (coefficients for difference in loge[antibody concentration] per doubling of serum PFAS) with 95% confidence intervals that excluded zero ("statistical support"), and certainty of evidence ratings were as follows: for PFOA and all antibodies treated as one type in all ages, -0.06 (-0.10, -0.01; moderate) and in children, -0.10 (-0.16, -0.03; moderate); for Diphtheria in children, -0.12 (-0.23, -0.00; high); for Rubella in all ages, -0.09 (-0.17, -0.01; moderate), and for Tetanus in children, -0.12 (-0.24, -0.00; moderate). For PFOS the summary associations were, for all antibodies treated as one type in all ages, -0.06 (-0.11, -0.01; moderate) and in children, -0.10 (-0.18, -0.03; moderate); for Rubella in all ages, -0.09 (-0.15, -0.03; high) and in children, -0.12 (-0.20, -0.04; high). For PFHxS the summary associations were, for all antibodies treated as one type in all ages, -0.03 (-0.06, -0.00; moderate) and in children, -0.05 (-0.09, -0.00; low); and for Rubella in children, -0.07 (-0.11, -0.02; high). Summary associations for PFNA and PFDA did not have statistical support, but all PFAS studied tended to have an inverse association with antibody concentrations. LIMITATIONS AND CONCLUSIONS: Epidemiologic data on immunosuppression and five principal PFAS suggest an association, with support across antibodies against multiple types of antigens. Data on Diphtheria, Rubella, and Tetanus were more supportive of an association than for other antibodies, and support was greater for associations with PFOA, PFOS, and PFHxS, than for PFNA or PFDA. The data on any specific antibody were scarce. Confounding factors that might account for the relation were not identified. Nearly all studies evaluated were judged to have a low or moderate risk of bias.

A Generic Pharmacokinetic Model for Quantifying Mother-to-Offspring Transfer of Lipophilic Persistent Environmental Chemicals.

Lipophilic persistent environmental chemicals (LPECs) can accumulate in a woman's body and transfer to her developing child across the placenta and via breast milk. To assess health risks associated with developmental exposures to LPECs, we developed a pharmacokinetic (PK) model that quantifies mother-to-offspring transfer of LPECs during pregnancy and lactation and facilitates internal dosimetry calculations for offspring. We parameterized the model for mice, rats, and humans using time-varying functions for body mass and milk consumption rates. The only required substance-specific parameter is the elimination half-life of the LPEC in the animal species of interest. We used the model to estimate whole-body concentrations in mothers and offspring following maternal exposures to hexachlorobenzene (HCB) and 2,2',4,4',5,5'-hexachlorobiphenyl (PCB 153) and compared these with measured concentrations from animal studies. We also compared estimated concentrations for humans to those generated using a previously published human LPEC PK model. Finally, we compared human equivalent doses (HEDs) calculated using our model and an allometric scaling method. Estimated and observed whole-body concentrations of HCB and PCB 153 in offspring followed similar trends and differed by less than 60%. Simulations of human exposure yielded concentration estimates comparable to those generated using the previously published model, with concentrations in offspring differing by less than 12%. HEDs calculated using our PK model were about 2 orders of magnitude lower than those generated using allometric scaling. Our PK model can be used to calculate internal dose metrics for offspring and corresponding HEDs and thus informs assessment of developmental toxicity risks associated with LPECs.

Microsolvation of uracil and its conjugate bases: a DFT study of the role of solvation on acidity.

The effect of microsolvation on the deprotonation energies of uracil was examined using DFT. The structures of uracil and its N(1) and N(3) conjugate bases were optimized with zero to six associated water molecules. Multiple configurations (upward of 93) of these hydrated clusters were located at PBE1PBE/6-311+G(d,p). Trends in these geometries are discussed, with the waters generally forming chains with small numbers of waters (one-three), rings (three-five waters), or cages (five-six waters). The difference in energy between the N1 and N3 conjugate bases is 13 kcal mol(-1) in the gas phase, and it decreases with each added water up to four. At this point the energy difference has been halved, but addition of a fifth or sixth water has little effect on the energy difference. This is understood in terms of the water structures and their ability to stabilize the negatively charged atoms in the conjugate bases.

The concentration of several perfluoroalkyl acids in serum appears to be reduced by dietary fiber.

Fiber-rich food intake has been associated with lower serum concentrations of perfluoroalkyl substances (PFAS) in some studies and dietary fiber was related to lower serum PFAS in a recent study. Given the previous epidemiologic data suggesting that fiber might decrease serum PFAS concentrations, we examined the relation of serum PFAS concentrations to intake of dietary fiber in National Health and Nutrition Examination Survey (NHANES) data. We examined the PFAS-fiber association among 6482 adults who participated in the NHANES, 2005-2016. Fiber intake was estimated based on two 24-hour diet recalls. We adjusted the models for determinants of PFAS and potentially confounding factors such as intake of foods reported to increase PFAS exposure. Results were expressed as the percent difference in PFAS concentration per interquartile range (IQR) increase in fiber (and 95 percent confidence interval), and the NHANES sampling parameters were used to make the results generalizable to the U.S. The adjusted percent difference in perfluorooctanoic acid (PFOA) per IQR increase in fiber was -3.64 (-6.15, -1.07); for perfluorooctane sulfonic acid (PFOS) was -6.69 (-9.57, -3.73), and for perfluorononanoic acid (PFNA) was -8.36 (-11.33, -5.29). These results suggest that dietary fiber increases the gastrointestinal excretion of PFOA, PFOS, and PFNA. Because fiber also lowers serum cholesterol, in some studies of the serum cholesterol-PFAS relationship confounding by fiber may be worth evaluating.

Birth weight and perfluorooctane sulfonic acid: a random-effects meta-regression analysis.

Perfluorooctane sulfonic acid (PFOS) is a ubiquitous environmental contaminant. Most people in developed countries have detectable serum concentrations. Lower birth weight has been associated with serum PFOS in studies world-wide, many of which have been published only recently. METHODS: To facilitate a causal assessment of the birth weight and PFOS association, we updated previous meta-analyses of the association and employed a method that facilitated inclusion of all available data in one analysis. Our analysis was based on observations from 29 studies. RESULTS: The random effects summary was -3.22 g/ng/ml (95% confidence interval [CI] = -5.11, -1.33). In a subgroup analysis stratified by when in pregnancy the PFOS concentration was measured, the summary for the early group was -1.35 (95% CI = -2.33, -0.37) and for the later group was -7.17 (95% CI = -10.93, -3.41). In a meta-regression model including a term for timing of blood draw, the intercept was slightly positive but essentially zero (0.59 g/ng/ml, 95% CI = -1.94, 3.11). In other words, the model indicated that when blood was drawn at the very beginning of pregnancy, there was essentially no relation of birth weight to PFOS. The results from the subgroup analyses differed from those from the model because the average gestational age at blood draw in the early group was 14 weeks, when bias would still be expected. A stronger inverse association in Asian studies was not completely explained by their blood draws being from later in pregnancy. CONCLUSIONS: The evidence was weakly or not supportive of a causal association.

Pharmacokinetic bias analysis of an association between clinical thyroid disease and two perfluoroalkyl substances.

Exposure to perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) has been associated with the occurrence of thyroid disease in some epidemiologic studies. We hypothesized that in a specific epidemiologic study based on the National Health and Nutrition Examination Survey, the association of clinical thyroid disease with serum concentration of PFOA and PFOS was due to reverse causality. Thyroid hormone affects glomerular filtration, which in turn affects excretion of PFOA and PFOS. We evaluated this by linking a model of thyroid disease status over the lifetime to a physiologically based pharmacokinetic model of PFOA and PFOS. Using Monte Carlo methods, we simulated the target study population and analyzed the data using multivariable logistic regression. The target and simulated populations were similar with respect to age, estimated glomerular filtration rate, serum concentrations of PFOA and PFOS, and prevalence of clinical thyroid disease. The analysis showed little or no evidence of bias from the hypothesized mechanism. The largest bias was for the fourth quartile of PFOA in females, with an odds ratio of 0.93 (95% CI, 0.90, 0.97). The reported odds ratio of clinical thyroid disease for this group was 1.63 (1.07, 2.47), and if it were corrected for the bias would have been 1.74 (1.14, 2.65). Our results suggest that little of the reported association in the target study was due to reverse causality.

Development of an in vitro approach to point-of-contact inhalation toxicity testing of volatile compounds, using organotypic culture and air-liquid interface exposure.

In vitro chemical risk assessment using human cells is emerging as an alternative to in vivo animal testing with reduced costs, fewer animal welfare concerns, and the possibility of greater human health relevance. In vitro inhalation toxicity testing of volatile compounds poses particular challenges. Here we report our efforts to establish a testing protocol in our own lab using the EpiAirway bronchial epithelium cell culture model and the Vitrocell 12/12 system for air-liquid interface (ALI) exposures. For purposes of method development, we used methyl iodide (MeI) as a test compound. We examined viability, cytotoxicity, and epithelial integrity responses. Dose-dependent, reproducible responses were observed with all assays. EpiAirway and BEAS-2B cytotoxicity responses to acute exposure were roughly similar, but EpiAirway was more resistant than BEAS-2B by the viability measurement, suggesting a proliferative response at low MeI concentrations. If wells were sealed to prevent evaporation, in-solution MeI concentration-response could be used to predict the response to MeI vapor within 2-fold by converting from the media- to the air-concentration at equilibrium using the blood:air partition coefficient for MeI. The long-term stability of EpiAirway cultures enabled repeated exposures over a 5-d period, which produced responses at lower concentrations than did acute exposure.

A model of functional thyroid disease status over the lifetime.

Mathematical models of the natural history of disease can predict incidence rates based on prevalence data and support simulations of populations where thyroid function affects other aspects of physiology. We developed a Markov chain model of functional thyroid disease status over the lifetime. Subjects were in one of seven thyroid disease states at any given point in their lives [normal, subclinical hypothyroidism, overt hypothyroidism, treated thyroid disease (ever), subclinical hyperthyroidism, overt hyperthyroidism, and reverted to normal thyroid status]. We used a Bayesian approach to fitting model parameters. A priori probabilities of changing from each disease state to another per unit time were based on published data and summarized using meta-analysis, when possible. The probabilities of changing state were fitted to observed prevalence data based on the National Health and Nutrition Examination Survey 2007-2012. The fitted model provided a satisfactory fit to the observed prevalence data for each disease state, by sex and decade of age. For example, for males 50-59 years old, the observed prevalence of ever having treated thyroid disease was 4.4% and the predicted value was 4.6%. Comparing the incidence rates of treated disease predicted from our model with published values revealed that 82% were within a 4-fold difference. The differences seemed to be systematic and were consistent with expectation based on national iodine intakes. The model provided new and comprehensive estimates of functional thyroid disease incidence rates for the U.S. Because the model provides a reasonable fit to national prevalence data and predicts thyroid disease status over the lifetime, it is suitable for simulating populations, thereby making possible quantitative bias analyses of selected epidemiologic data reporting an association of thyroid disease with serum concentrations of environmental contaminants.

Targeting a Rate-Promoting Vibration with an Allosteric Mediator in Lactate Dehydrogenase.

We present a new type of allosteric modulation in which a molecule bound outside the active site modifies the chemistry of an enzymatic reaction through rapid protein dynamics. As a test case for this type of allostery, we chose an enzyme with a well-characterized rate-promoting vibration, lactate dehydrogenase; identified a suitable small molecule for binding; and used transition path sampling to obtain ensembles of reactive trajectories. We found that the small molecule significantly affected the reaction by changing the position of the transition state and, through applying committor distribution analysis, showed that it removed the protein component from the reaction coordinate. The ability of a small-molecule to disrupt enzymatic reactions through alteration of subpicosecond protein motion opens the door for new experimental studies on protein motion coupled to enzymatic reactions and possibly the design of drugs to target these enzymes.

Another Look at the Mechanisms of Hydride Transfer Enzymes with Quantum and Classical Transition Path Sampling.

The mechanisms involved in enzymatic hydride transfer have been studied for years, but questions remain due, in part, to the difficulty of probing the effects of protein motion and hydrogen tunneling. In this study, we use transition path sampling (TPS) with normal mode centroid molecular dynamics (CMD) to calculate the barrier to hydride transfer in yeast alcohol dehydrogenase (YADH) and human heart lactate dehydrogenase (LDH). Calculation of the work applied to the hydride allowed for observation of the change in barrier height upon inclusion of quantum dynamics. Similar calculations were performed using deuterium as the transferring particle in order to approximate kinetic isotope effects (KIEs). The change in barrier height in YADH is indicative of a zero-point energy (ZPE) contribution and is evidence that catalysis occurs via a protein compression that mediates a near-barrierless hydride transfer. Calculation of the KIE using the difference in barrier height between the hydride and deuteride agreed well with experimental results.