Next Article in Journal
Adsorption Characteristics of Modified Bamboo Charcoal on Cu(II) and Cd(II) in Water
Next Article in Special Issue
Preliminary Toxicological Analysis in a Safe-by-Design and Adverse Outcome Pathway-Driven Approach on Different Silver Nanoparticles: Assessment of Acute Responses in A549 Cells
Previous Article in Journal
Deteriorative Effects of Radiation Injury Combined with Skin Wounding in a Mouse Model
Previous Article in Special Issue
On the In Vitro and In Vivo Hazard Assessment of a Novel Nanomaterial to Reduce the Use of Zinc Oxide in the Rubber Vulcanization Process
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxics
Volume 10
Issue 12
10.3390/toxics10120786
toxics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ambient Benzo[a]pyrene’s Effect on Kinetic Modulation of Amyloid Beta Peptide Aggregation: A Tentative Association between Ultrafine Particulate Matter and Alzheimer’s Disease

by
Samal Kaumbekova
1,
Mehdi Amouei Torkmahalleh
2 and
Dhawal Shah
1,*
1
Department of Chemical and Materials Engineering, School of Engineering and Digital Sciences, Nazarbayev University, Kabanbay Batyr 53, Astana 010000, Kazakhstan
2
Division of Environmental and Occupational Health Sciences, School of Public Health, University of Illinois at Chicago, Chicago, IL 60612, USA
*
Author to whom correspondence should be addressed.
Toxics 2022, 10(12), 786; https://doi.org/10.3390/toxics10120786
Submission received: 9 November 2022 / Revised: 9 December 2022 / Accepted: 13 December 2022 / Published: 14 December 2022
(This article belongs to the Special Issue Nano and Ultrafine Particle Toxicology and Exposure Assessment)
Browse Figures
Versions Notes

Abstract

:
Long-time exposure to ambient ultrafine particles is associated with an increased risk of neurodegenerative diseases such as Alzheimer’s disease (AD), which is triggered by the aggregation of Aβ peptide monomers into toxic oligomers. Among different ultrafine air pollutants, polycyclic aromatic hydrocarbons (PAHs) are known to have a negative neural impact; however, the impact mechanism remains obscure. We herein examined the effect of Benzo[a]Pyrene (B[a]P), one of the typical PAHs on Aβ42 oligomerization using all-atom molecular dynamics simulations. In particular, the simulations were performed using four molecules of Aβ42 in the presence of 5.00 mM, 12.5 mM, and 50.0 mM of B[a]P. The results revealed strong hydrophobic interactions between Aβ42 peptides and B[a]P, which in turn resulted in increased interpeptide electrostatic interactions. Furthermore, 5.00 mM of B[a]P accelerated the kinetics of the formation of peptide tetramer by 30%, and stabilized C-terminus in Aβ42 peptides, suggesting consequent progression of AD in the presence of 5.00 mM B[a]P. In contrast, 12.5 mM and 50.0 mM of B[a]P decreased interpeptide interactions and H-bonding due to the aggregation of numerous B[a]P clusters with the peptides, suppressing oligomerization kinetics of Aβ42 peptides by 13% and 167%, respectively. While the study elucidates the effect of small environmental hydrophobic molecules on the formation of Aβ oligomers, the impact of ambient ultrafine particles on AD in the complex composition of the environmental realm requires further systematic delving into the field.

1. Introduction

Environmental air pollutants are known to have a negative neural impact [1], with the effect ranging from an increased risk of autism spectrum disorders, brain volumetric changes, brain inflammation, cerebrovascular and neuropsychiatric disorders, and the hallmarks of Alzheimer’s disease (AD) [2,3,4,5,6]. The concern about the development of AD-like neurodegenerative diseases increases with the growth of the elderly population globally. AD is associated with the development of amyloid plaques due to the aggregation of amyloid beta (Aβ) peptide monomers into soluble oligomers, and the formation of neurofibrillary tangles in the human brain [7]. The pathogenesis of AD also implicates genetic and environmental factors [8].
Among different air pollutants, Polycyclic Aromatic Hydrocarbons (PAHs) are associated with an increased risk for neurodegeneration [9]. Benzo[a]pyrene (B[a]P, C20H12) is a typical PAH, with five aromatic rings in its structure. B[a]P is usually generated from incomplete combustion of organic material, motor-vehicle exhaust, and from cooking and smoking. According to source apportionment measurements, the average concentration of B[a]P in ambient air might reach up to 1.94 ng/m3 in highly industrialized regions of Italy [10]. Long-term monitoring analysis of ambient air in the Czech Republic revealed high levels of B[a]P (up to 7.7 ng/m3) in industrialized regions, in addition to the elevated concentrations of B[a]P (up to 13.6 ng/m3) observed during the winter season due to the household and local heating [11]. The mean values of average annual concentrations of B[a]P and a total of 24 PAHs in the ultrafine particles in different regions of China measured from October 2012 to September 2013 were 5.02 ng/m3 and 53.5 ng/m3, respectively [12]. Hydrophobic ultrafine molecules, such as PAHs, can diffuse through cell membranes [13] and cross the blood–brain barrier [14]. The concentration of PAHs in the blood of children in Nigeria varied from 535 to 708 µg/L [15]. The effect of B[a]P on the neurobehavioral functions of coke oven workers was recently studied in Taiyuan, China [16]. The concentrations of B[a]P observed in the coke oven’s bottom, side, and top regions were 19.5 ± 13.2, 185.9 ± 38.6, and 1623.5 ± 435.8 ng/m3, respectively. In addition, the authors evaluated the neurobehavioral function of coke oven workers and observed a statistical decline in the exposed group.
The formation of soluble aggregates of Aβ peptides in the human cerebrospinal fluid is related to the progression of AD via synaptic and neuronal loss [17]. The aggregation of Aβ peptides can be observed in various brain regions at distinct stages of AD [18]. According to the amyloid cascade hypothesis, the aggregation of Aβ peptides into β—sheets via fibrillization leads to the formation of amyloid plaques in different areas of the brain [19]. Soluble Aβ peptide oligomers are considered the most toxic form of amyloid aggregates associated with neuronal damage [20]. According to clinical studies, the concentration of Aβ oligomers may reach up to 2647 fg/mL in the human cerebral fluid (CSF) of AD patients with moderate dementia [21]. According to cell binding and toxicity models in cortical neuronal cultures, tetramers possess higher neurotoxicity among different Aβ oligomers with lower molecular weight [22]. Amyloid fibrils and aggregates are formed more rapidly from the Aβ42 peptide variant, with 42 aminoacids in its structure among different Aβ peptide isoforms with 39 to 43 aminoacids [23]. While the first 16 aminoacid residues of the N-terminus of Aβ42 are hydrophilic, the remaining region of the peptide is hydrophobic, including the central region (Aβ12–23) and C-terminus (Aβ30–42) [24,25].
A recent molecular dynamics (MD) study showed that B[a]P molecule decreased the amount of helices, promoting the formation of β-sheets and β-bridges in the secondary structure of Aβ42 peptide monomer [26]. Furthermore, in vitro experiments on the effect of cigarette smoke components on the structure and aggregation of Aβ peptides showed that PAHs increased the aggregation kinetics of Aβ peptides [27]. In particular, the presence of B[a]P decreased the aggregation halftime from 6.5 ± 0.9 h (for 10–20 µM of Aβ40 dissolved in 20 mM of buffer) to 5.1 ± 0.7 h (Aβ dissolved in 100–200 µM of B[a]P) [27]. According to Gao et al. [28,29], chronic exposure to B[a]P increases Aβ42 concentrations in the brain of the zebrafishes. In addition, the study showed that cognitive, memory, and locomotor activities of the fish decreased after the B[a]P exposure. Moreover, Liu et al. [30] observed elevated concentrations of Aβ monomers and oligomers, as well as the formation of Aβ plaques and Aβ fibrils in the brain of mice due to the exposure to B[a]P.
While in vivo and in vitro experiments showed enhanced aggregation of Aβ peptides and the formation of amyloid plaques in the presence of B[a]P, the molecular interactions between Aβ peptides and the pollutant are obscure. In this study, we performed a molecular dynamics study aimed to investigate the effect of B[a]P on the structure of Aβ42 peptide and oligomerization kinetics of four Aβ42 peptide monomers. Furthermore, to investigate the impact of varying concentrations of B[a]P on the oligomerization of Aβ42 peptides, the molecular dynamics simulations, each of 500 ns, were performed in the presence of different concentrations of B[a]P molecules. Moreover, to get statistically relevant results within reasonable computational time in the MD study, the concentrations of the peptides and B[a]P molecules inserted in the simulation box were higher than what would be found in in vitro and in vivo studies [27,28]. In particular, to investigate the impact of B[a]P, we performed simulations in the presence of 4 B[a]P molecules (5 mM), keeping the 1:1 ratio between the peptide and air pollutant, as was used in our previous study, where we simulated systems with 1 monomer and 1 B[a]P molecule [26]. The highest concentration of B[a]P used in the simulations was 50 mM, with the 10:1 ratio between B[a]P and peptides, as was previously used in the experimental study performed by Wallin et al. [27]. The choice of 12.5 mM was taken as an intermediate concentration between the two values. Considering that the simulation box is a simple model of a physical system, the aggregation of peptides occurs faster in the MD simulations, while the results of the MD studies represent the general trend that would be expected to be found in an actual system.

2. Methodology

Atomistic molecular dynamics simulations were performed via Gromacs 2019.6 software [31]. Gromos54a7 forcefield parameters were chosen based on the literature validations [32] and our previous study with the validation of the density of B[a]P at 20 °C [26]. The molecular dynamics simulations were performed in a simulation box with the dimensions of 11 × 11 × 11 nm3. Four Aβ42 monomers (PDB ID: 1Z0Q, with a total charge of −2 on each peptide) were inserted in a simulation box, keeping the concentration of peptides at 5.00 mM. The system was solvated via SPC water model. 0.15 M of NaCl salt was added as a buffer environment. The topology parameters and optimized geometry of B[a]P molecule (C20H12) were taken from the Automated Topology Builder (ATB, Version 3.0) [33]. The concentration of B[a]P varied from 0 mM and 5.00 mM, 12.5 mM, and 50.0 mM, corresponding to 4, 10, and 40 molecules of B[a]P. Due to the limitations of the simulation time and box size, the concentrations of Aβ42 peptide and B[a]P molecules in the simulated systems were comparatively higher than the molecular concentrations normally found in the human blood [15]. However, the relative amount of B[a]P to Aβ, which is 10:1, is similar to those used in the in vitro experiments [27]. The number of molecules used in the simulated systems is shown in Table 1.
Energy minimization was performed with the steepest descent algorithm, applying periodic boundary conditions in all directions, and setting the maximum force for convergence of 100 kJ mol−1 nm−1. After that, NVT—equilibration step with H-bonds constraints was performed for 25 ps with a time step of 0.5 fs. Next, the NPT—equilibration step with all-bonds constraints was performed for 100 ps with a time step of 2 fs, at reference pressure of 1 bar. LINCS (LINear Constraint Solver) constraint algorithm [34] with LINCS iterations of 1 was applied at reference temperature of 298 K. Short-range cutoff distance of 1 nm was applied for both electrostatics and van der Waals interactions with Verlet cut-off scheme algorithm. PME (Particle-mesh Ewald) was set for long-range electrostatic interactions. Molecular dynamics simulations were performed for 500 ns with a time step of 2 fs, considering reaching an equilibrium within a simulated time in all systems under the study, as discussed below. The output coordinates, velocities, and energies were saved every 4 ns.
The oligomerization (and kinetics) of four Aβ42 peptides were investigated by using cluster and intermolecular distance analyses. The intermolecular distances were calculated based on the center of masses (COM) of the residues and taking the average values between the distances between COM of peptides 1 and 2, peptides 1 and 3, peptides 1 and 4, etc. The formation of the clusters was studied for the defined groups of the residues, “four peptides” and “peptides and B[a]P”, to analyze the formation of the interpeptide and peptide-B[a]P clusters, respectively. The cluster analyses were performed setting 0.35 nm as a maximum distance parameter between the center of masses of the residues within the defined group to be defined as a cluster. Time-evolution of solvent accessible surface area (SASA), radius of gyration (RoG), and secondary structure analyses were performed to investigate the deviations in the structure of the peptides within 500 ns of the simulations. The (radial distribution function) rdf analysis was performed to investigate the probability of finding “peptides” and “B[a]P” residues from “peptides” residues using the COM of the residues. H-bonds analyses, root-mean square fluctuations (RMSF), and (rdf) analyses were performed for the last 30 ns of the simulations to investigate the structure of peptides and intermolecular interactions in the end of the simulations. In addition, visual molecular dynamics (VMD) software [35] was used to visualize the systems under the study. Finally, the non-bonded intermolecular interactions were studied for the last 10 ns of the MD simulations. In particular, the short-range (SR) interactions with the cut-off distance of 1 nm, and long-range (LR) interaction energies were calculated separately. In addition, the Lennard-Jones (LJ) potential and Coulombic (Coul) potential energies were quantified, indicating the repulsion-dispersion and electrostatics interactions, respectively.
Finally, an additional short replica of each system was simulated independently for 250 ns, starting from the randomly inserted molecules in the simulation box. The results of the simulations were validated via interpeptide distance and SASA analyses.

3. Results and Discussion

The kinetics of the tetramerization of Aβ42 peptides was initially studied in the absence of B[a]P and in the presence of 4 B[a]P molecules via interpeptide distance analysis (Figure 1A). Similarly, the effect of B[a]P concentration on the oligomerization was investigated by interpeptide distance analysis (Figure 1B), followed by the peptide—B[a]P distance analysis (Figure 1C).
The analysis of the interpeptide distances showed that in the absence of B[a]P molecules, four Aβ42 peptides aggregated within the simulated time; in particular, within 100 ns of the simulation (Figure 1A). A different trend was observed by Brown et al. [36], wherein the tetramerization of Aβ42 was seen within 250 ns of the simulation. The difference in the tetramerization time could be the usage of older forcefield parameters (GROMOS96 53A6), higher temperature of simulations (310 K), and large simulation box (with dimensions of 12.7 nm) by Brown and Bevan [36]. According to Figure 1A, the presence of 4 B[a]P molecules significantly decreased the time required for the oligomerization. The peptides aggregated within 50 ns of the simulation, which was the fastest aggregation among the systems under the study (Figure 1B). The results obtained from the interpeptide distance analysis were validated via an additional short replica for 250 ns of each system. Additional plots for each system under the study are available in the Supplementary Material (Figure S1). While quantitatively the values were different, qualitatively the trends remained the same.
With increasing the concentration of B[a]P molecules, i.e., using 10 and 40 molecules, the time required for tetramerization increased, and became 100 ns and 200 ns for the two systems, respectively. In addition, according to Figure 1C, the distance analyses between peptides and B[a]P molecules showed that B[a]P molecules were also bound to the peptides within 70 ns of the simulations in all systems under the study, indicating the formation of agglomerate encompassing the PAHs and the peptides. Furthermore, the aggregation kinetics of four Aβ42 peptides was analyzed in terms of the growth of the interpeptide clusters (Figure 2A), as well as the formation of the clusters of B[a]P and peptides (Figure 2B).
According to Figure 2A, the rate of tetramerization of four peptides was fastest in the presence of 4 B[a]P molecules (in ~52 ns of the MD run). In comparison, in the absence of B[a]P molecules, a stable single cluster of four peptides was produced in ~75 ns. In the presence of 10 B[a]P molecules, the cluster was formed in ~85 ns of the simulation, while in the presence of 40 B[a]P, formation of a stable cluster of four peptides was completed in ~200 ns of the production run. According to Figure 2B, the fastest formation of a stable peptide—B[a]P cluster occurred in ~40 ns of the simulation in the presence of 4 B[a]P molecules. In comparison, additional time was required for the complete aggregation of peptides and B[a]P molecules in the systems with 10 B[a]P molecules (~85 ns of the simulation) and 40 B[a]P molecules (~155 ns of the simulation).
Overall, according to cluster analyses, the presence of 5.00 mM B[a]P enhanced the oligomerization kinetics, based on the time required to make one stable cluster of four Aβ42 peptide monomers, by 30%, in comparison to the aggregation kinetics of peptides in the absence of B[a]P. In contrast, the presence of 12.5 mM B[a]P and 50.0 mM B[a]P decreased the aggregation kinetics of peptides by 13% and 167%, respectively. In addition, in the system with 10 B[a]P molecules, the formation of interpeptide clusters and formation of peptide—B[a]P clusters occurred simultaneously, while in the systems with 4 and 40 B[a]P molecules, the formation of a stable interpeptide cluster occurred only after all B[a]P molecules were bound to peptides. The results indicated that oligomerization was essentially driven by the hydrophobic interactions between B[a]P and Aβ42 peptides [37], rather than the intrinsic interpeptide hydrophobic interactions [38]. The secondary structure of the interpeptide clusters was further analyzed for the last 30 ns of the simulations (Figure 3), wherein the tetramerized, clustered structure was stable. In addition, the intermolecular clusters observed at the end of the simulations, along with the initial monomer structure, were visualized via VMD (Figure 4).
The results of the secondary structure analysis obtained from our previous study [26] showed that in the presence of 1 B[a]P molecule, high amounts of the coils, bends, and β-content with low amount of helices were observed in the secondary structure of the peptide monomer. The results of this study (Figure 3) showed that the secondary structure of the oligomers would also depend on the concentration of the B[a]P molecules present in the simulation box. In particular, the analysis of the time-averaged secondary structure of Aβ42 peptides showed that in the presence of 4 B[a]P molecules, the coil region was predominant (38%). In comparison, in the presence of 10 and 40 B[a]P molecules, the helix region was prevalent (34% and 36%, respectively). Typically, the progression of Alzheimer’s disease is associated with the formation of β-sheet fibrils [39,40], however, more simulation time will be required for the formation of the stable β-sheets from the random coils in the peptide oligomers [41]. Interestingly, the β-content did not significantly vary with different B[a]P concentrations, with 6–9% of the β-content observed in all systems under the study. Nonetheless, the formation of the compact coil structures in Aβ oligomers is also an important step toward the formation of β-sheets [42]. In particular, in the presence of 4 B[a]P molecules, the high number of coils in the secondary structure of the peptides was correlated with enhanced oligomerization kinetics.
According to Figure 4A, in the beginning of the simulation, Aβ42 monomers consisted of helixes (Aβ10–23, Aβ28–31), turn, bend and coil regions (N- and C- terminuses), with no β—sheets observed in their secondary structures. At the end of 500 ns of the simulation system with no B[a]P had β-sheets in the regions of Aβ31–36 and Aβ40 of the C-terminus (Figure 4B). According to Figure 4C, in the presence of 4 B[a]P, β-sheets were observed in the regions of Aβ28–32 and Aβ35–39 at the end of the simulation. In the system with 10 B[a]P, the formation of the β-sheets in occurred in the regions of Aβ2–5 and Aβ30–41 (Figure 4D), whereas with 40 B[a]P, β-sheets were formed in the regions of Aβ28–32 and Aβ35–39 (Figure 4E). These observations showed that β-sheets were mainly formed in the regions of C-terminus, considered as a possible nucleation site for the self-aggregation, as was highlighted in the earlier literature [43,44]. In addition, in the presence of 4 B[a]P, the tetramerization occurred with the formation of oligomers with increased coil region, that could be expected to form β-sheets further [45]. In comparison, in the systems with 10 and 40 B[a]P molecules, the binding of numerous B[a]P to peptides via aromatic π-π interactions [46] resulted in the formation of B[a]P—peptide agglomerates. Consequently, typical oligomerization was suppressed due to the presence of numerous B[a]P molecules and their clusters, resulted in steric hindrance [47], as shown in Figure 4D,E. To investigate the deviations in the peptide structure with the formation of intermolecular clusters, RMSF, RoG, SASA, and H-bonds analyses were further performed (Figure 5, Table 2).
According to RMSF analyses of the last 30 ns of the simulations, the fluctuations in the positions of aminoacid residues in Aβ42 were in the range of ~0.1–0.22 nm, in the absence of B[a]P molecules (Figure 5A). In system 2 with 4 B[a]P molecules, in general, the RMSF values of the aminoacid residues were low in the Aβ3–42 segment (RMSF values of ~0.06–0.12 nm) with enhanced fluctuations in the N-terminus (RMSF values of ~0.36 nm in the Aβ1–2 region). In addition, the results showed that the presence of 4 B[a]P molecules stabilized the central hydrophobic core (RMSF values of ~0.09 nm in the Aβ15–20 region) and the beta-sheets region of the C-terminus (RMSF values of ~0.09 nm in the Aβ26–28 and Aβ38–41 regions), promoting the formation of tetramers. In the presence of 10 B[a]P molecules, the RMSF values of the aminoacid residues were low in Aβ3–36 segment (RMSF values of ~0.06–0.1 nm), indicating a similar effect of the tetramer stabilization. Moreover, in this system under study, the enhanced fluctuations were observed in the N- and C-terminuses, which were correlated with the suppression of oligomerization (RMSF values up to ~0.23 nm for Aβ1–2, and ~0.14 nm for Aβ38–42). Similarly, enhanced fluctuations in the positions of the amino acids of the N-terminus (Aβ1–2, RMSF values of ~0.28 nm), and C-terminus (Aβ37–41 segment RMSF values of ~0.14–0.21 nm) were observed in the system with 40 B[a]P molecules in the end of the simulation, indicating increased movement of the aminoacids in the peptide terminuses due to the interactions with large amounts of B[a]P molecules and the hydrophobic interactions with the PAHs, as was discussed previously.
This observation was further corroborated with the time-evolution of the radius of gyration (RoG) of the peptide. The RoG values, averaged between four Aβ42 peptides, decreased from the initial value of 1.6 nm to ~1.1 nm within the simulation time in the systems with no B[a]P and in the presence of 4 B[a]P molecules (Figure 5B). In contrast, in the presence of 10 B[a]P and 40 B[a]P molecules, the radius of gyration deviated more significantly, up to ~1.3 nm and ~1.4 nm, respectively, due to the interactions with large amounts of B[a]P molecules.
In agreement with cluster and intermolecular distance analyses, for the systems with no B[a]P molecules, with 4 B[a]P, and 10 B[a]P molecules, the SASA values significantly decreased within first 100 ns of the MD simulations, indicating the interpeptide binding with consequent loss of peptide surface area available for solvent (Figure 5C). In addition, according to Figure 5C and Table 2, the final values of SASA at 500 ns were comparatively high in the systems with high number of B[a]P molecules in the simulation box. In particular, in the end of the simulations, the SASA values of peptides were 94.6 nm2, 102.7 nm2, 110.0 nm2, and 139.2 nm2 in the systems with no B[a]P, with 4 B[a]P, 10 B[a]P, and 40 B[a]P molecules, respectively (Table 2). This observation indicates that binding of large amount of B[a]P to Aβ42 peptides would consequently inhibit the oligomerization. Moreover, it was also noted that in the first 20 ns of the simulations the aggregation of peptides was the fastest in the system with 10 B[a]P molecules (Figure 5D). However, the initial rate of aggregation of peptides, within 20 ns of the simulation, was facilitated in the presence of 10 B[a]P molecules, and the consequent aggregation of B[a]P molecules to peptides inhibited oligomerization of four peptides. The results obtained from the SASA analysis were validated via an additional short replica of 250 ns performed for all systems under the study. Additional plots for each system under the study are available in the Supplementary Material (Figure S2). While quantitatively the values were different, qualitatively the trends remained the same.
The analysis of H-bonds in the last 30 ns of the simulations showed increased amounts of H-bonds in the presence of 4 B[a]P molecules (119 ± 5 bonds), in comparison to the system with 40 B[a]P molecules (109 ± 5 bonds), indicating that the presence of a large amount of B[a]P molecules inhibited the formation of interpeptide H-bonds (Table 2). Overall, RoG, SASA, and H-bonds analyses showed that the presence of 10 and 40 B[a]P molecules could interfere with the oligomerization process by binding to Aβ42 peptides, leading to elevated radius of gyration, SASA of Aβ42 peptides, and suppressed interpeptide H-bonding.
Furthermore, the radial distribution function (rdf) analysis showed that high interpeptide interactions were observed in the systems with no B[a]P molecules (maximum peaks of ~55 at 0.5 nm and ~53 at 0.6 nm on Figure 6A), and in the presence of 4 B[a]P molecules (maximum peak of ~53 at 0.6 nm on Figure 6A). In comparison, decreased interpeptide interactions in the presence of 10 B[a]P molecules (maximum peak value of ~50 at 0.6 nm on Figure 6A), and in the presence of 40 B[a]P molecules (highest peaks of ~45 at 0.4 nm and 0.6 nm on Figure 6A), indicated that interpeptide interactions were decreased with the addition of higher amounts of B[a]P molecules. According to the rdf analysis of peptide—B[a]P interactions (Figure 6B), strong peptide—B[a]P interactions were observed in the system with 4 B[a]P molecules (maximum peaks of ~30 at 0.47 nm and ~60 at 0.72 nm on Figure 6B). In comparison, low peptide—B[a]P interactions were observed in the systems with 10 B[a]P molecules (maximum peaks of ~17 at 0.47 nm and ~30 at 0.72 nm on Figure 6B) and 40 B[a]P molecules (maximum peaks of ~10 at 0.4 nm and ~20 at 0.76 nm on Figure 6B).
The energy analysis was further performed for the last 10 ns of the simulations to investigate the type of molecular interactions. According to Table 3, in the presence of B[a]P molecules, short-range electrostatic interpeptide interactions increased (Coul–SR), while short-range Lennard-Jones interactions (LJ–SR) between peptides decreased. These observations were also correlated to the enhanced hydrophobic interactions between peptides and B[a]P (LJ–SR), related to the aggregation of peptides and B[a]P molecules.

4. Conclusions

To conclude, molecular dynamics simulations revealed that with the addition of B[a]P molecules, electrostatic interactions between Aβ42 peptides increased due to the enhanced B[a]P—peptide hydrophobic interactions and binding of Aβ42 peptides and B[a]P. Although the presence of 4 B[a]P molecules did not significantly affect interpeptide H-bonding and interpeptide interactions, the presence of 10 and 40 B[a]P molecules suppressed H-bonding and decreased interpeptide interactions, observed in the end of the simulations. Overall, according to the results of our study, the presence of 5.00 mM B[a]P, i.e., 4 molecules, accelerated the formation of Aβ42 peptide tetramers by 30% and stabilized C-terminus of peptides, suggesting consequent progression of Alzheimer’s disease. Nonetheless, the effect of ultrafine air pollutants on AD would be more complicated due to the more complex composition of pollutants present in the environment.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxics10120786/s1, Figure S1: Time-evolution of average interpeptide distances in the systems with additional replicas; Figure S2: Time-evolution of Solvent Accessible Surface Area (SASA) of Aβ42 peptides in the systems with additional replicas.

Author Contributions

Conceptualization, S.K.; Investigation, S.K.; Methodology, D.S.; Software, D.S.; Supervision, M.A.T. and D.S.; Writing—original draft, S.K.; Writing—review & editing, M.A.T. and D.S. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the support provided by the Nazarbayev University under the project number 11022021FD2905 entitled “Efficient thermal valorization of municipal sewage sludge in fluidized bed systems: Advanced experiments with process modeling” and the project number 11022021CRP1503 Exposure to Cooking Ultrafine Particles and Neurodegenerative Disease: Clinical Exposure Studies and Computer Modeling.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated and analyzed during the current study are available from the corresponding author on a reasonable request.

Conflicts of Interest

The authors have no relevant financial or non-financial interest to disclose.

References

  1. Kim, H.; Kim, W.-H.; Kim, Y.-Y.; Park, H.-Y. Air Pollution and Central Nervous System Disease: A Review of the Impact of Fine Particulate Matter on Neurological Disorders. Front. Public Health 2020, 8, 575330. [Google Scholar] [CrossRef] [PubMed]
  2. Hahad, O.; Lelieveld, J.; Birklein, F.; Lieb, K.; Daiber, A.; Münzel, T. Ambient Air Pollution Increases the Risk of Cerebrovascular and Neuropsychiatric Disorders through Induction of Inflammation and Oxidative Stress. Int. J. Mol. Sci. 2020, 21, 4306. [Google Scholar] [CrossRef] [PubMed]
  3. Brockmeyer, S.; D’Angiulli, A. How air pollution alters brain development: The role of neuroinflammation. Transl. Neurosci. 2016, 7, 24–30. [Google Scholar] [CrossRef] [PubMed]
  4. Costa, L.G.; Cole, T.B.; Dao, K.; Chang, Y.-C.; Coburn, J.; Garrick, J.M. Effects of air pollution on the nervous system and its possible role in neurodevelopmental and neurodegenerative disorders. Pharmacol. Ther. 2020, 210, 107523. [Google Scholar] [CrossRef] [PubMed]
  5. Cipriani, G.; Danti, S.; Carlesi, C.; Borin, G. Danger in the Air: Air Pollution and Cognitive Dysfunction. Am. J. Alzheimer’s Dis. Other Dementiasr. 2018, 33, 333–341. [Google Scholar] [CrossRef] [PubMed]
  6. Heusinkveld, H.J.; Wahle, T.; Campbell, A.; Westerink, R.H.S.; Tran, L.; Johnston, H.; Stone, V.; Cassee, F.R.; Schins, R.P.F. Neurodegenerative and neurological disorders by small inhaled particles. Neurotoxicology 2016, 56, 94–106. [Google Scholar] [CrossRef] [PubMed]
  7. Mathis, C.A.; Wang, Y.; Klunk, W.E. Imaging beta-amyloid plaques and neurofibrillary tangles in the aging human brain. Curr. Pharm. Des. 2004, 10, 1469–1492. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, H.; Yang, F.; Zhang, S.; Xin, R.; Sun, Y. Genetic and environmental factors in Alzheimer’s and Parkinson’s diseases and promising therapeutic intervention via fecal microbiota transplantation. NPJ Park. Dis. 2021, 7, 70. [Google Scholar] [CrossRef]
  9. Cho, J.; Sohn, J.; Noh, J.; Jang, H.; Kim, W.; Cho, S.-K.; Seo, H.; Seo, G.; Lee, S.-K.; Noh, Y.; et al. Association between exposure to polycyclic aromatic hydrocarbons and brain cortical thinning: The Environmental Pollution-Induced Neurological EFfects (EPINEF) study. Sci. Total Environ. 2020, 737, 140097. [Google Scholar] [CrossRef]
  10. Gianelle, V.; Colombi, C.; Caserini, S.; Ozgen, S.; Galante, S.; Marongiu, A.; Lanzani, G. Benzo(a)pyrene air concentrations and emission inventory in Lombardy region, Italy. Atmos. Pollut. Res. 2013, 4, 257–266. [Google Scholar] [CrossRef]
  11. Schreiberová, M.; Vlasáková, L.; Vlček, O.; Šmejdířová, J.; Horálek, J.; Bieser, J. Benzo[a]pyrene in the Ambient Air in the Czech Republic: Emission Sources, Current and Long-Term Monitoring Analysis and Human Exposure. Atmosphere 2020, 11, 955. [Google Scholar] [CrossRef]
  12. Yu, Q.; Ding, X.; He, Q.; Yang, W.; Zhu, M.; Li, S.; Zhang, R.; Shen, R.; Zhang, Y.; Bi, X.; et al. Nationwide increase of polycyclic aromatic hydrocarbons in ultrafine particles during winter over China revealed by size-segregated measurements. Atmos. Chem. Phys. 2020, 20, 14581–14595. [Google Scholar] [CrossRef]
  13. Leikauf, G.D.; Kim, S.-H.; Jang, A.-S. Mechanisms of ultrafine particle-induced respiratory health effects. Exp. Mol. Med. 2020, 52, 329–337. [Google Scholar] [CrossRef] [PubMed]
  14. Gao, W.; Liu, Y.; Jing, G.; Li, K.; Zhao, Y.; Sha, B.; Wang, Q.; Wu, D. Rapid and efficient crossing blood-brain barrier: Hydrophobic drug delivery system based on propionylated amylose helix nanoclusters. Biomaterials 2017, 113, 133–144. [Google Scholar] [CrossRef] [PubMed]
  15. Wirnkor, V.A.; Ngozi, V.E.; Ajero, C.M.; Charity, L.K.; Ngozi, O.S.; Ebere, E.C.; Emeka, A.C. Biomonitoring of concentrations of polycyclic aromatic hydrocarbons in blood and urine of children at playgrounds within Owerri, Imo State, Nigeria. Environ. Anal. Health Toxicol. 2019, 34, e2019011. [Google Scholar] [CrossRef] [Green Version]
  16. Niu, Q.; Zhang, H.; Li, X.; Li, M. Benzo[a]pyrene-induced neurobehavioral function and neurotransmitter alterations in coke oven workers. Occup. Environ. Med. 2009, 67, 444–448. [Google Scholar] [CrossRef]
  17. De, S.; Whiten, D.R.; Ruggeri, F.S.; Hughes, C.; Rodrigues, M.; Sideris, D.I.; Taylor, C.G.; Aprile, F.A.; Muyldermans, S.; Knowles, T.P.J.; et al. Soluble aggregates present in cerebrospinal fluid change in size and mechanism of toxicity during Alzheimer’s disease progression. Acta Neuropathol. Commun. 2019, 7, 120. [Google Scholar] [CrossRef] [Green Version]
  18. I Sideris, D.; Danial, J.S.H.; Emin, D.; Ruggeri, F.S.; Xia, Z.; Zhang, Y.P.; Lobanova, E.; Dakin, H.; De, S.; Miller, A.; et al. Soluble amyloid beta-containing aggregates are present throughout the brain at early stages of Alzheimer’s disease. Brain Commun. 2021, 3, fcab147. [Google Scholar] [CrossRef]
  19. Barage, S.H.; Sonawane, K.D. Amyloid cascade hypothesis: Pathogenesis and therapeutic strategies in Alzheimer’s disease. Neuropeptides 2015, 52, 1–18. [Google Scholar] [CrossRef]
  20. Walsh, D.M.; Selkoe, D.J. A beta oligomers—A decade of discovery. J. Neurochem. 2007, 101, 1172–1184. [Google Scholar] [CrossRef]
  21. Hölttä, M.; Hansson, O.; Andreasson, U.; Hertze, J.; Minthon, L.; Nägga, K.; Andreasen, N.; Zetterberg, H.; Blennow, K. Evaluating Amyloid-β Oligomers in Cerebrospinal Fluid as a Biomarker for Alzheimer’s Disease. PLoS ONE 2013, 8, e66381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Jana, M.K.; Cappai, R.; Pham, C.L.L.; Ciccotosto, G.D. Membrane-bound tetramer and trimer Aβ oligomeric species correlate with toxicity towards cultured neurons. J. Neurochem. 2016, 136, 594–608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sun, X.; Chen, W.; Wang, Y. Beta-Amyloid: The key peptide in the pathogenesis of Alzheimer’s disease. Front. Pharmacol. 2015, 6, 221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rauk, A. Why is the amyloid beta peptide of Alzheimer’s disease neurotoxic? Dalton Trans. 2008, 10, 1273–1282. [Google Scholar] [CrossRef]
  25. Kepp, K.P. Bioinorganic Chemistry of Alzheimer’s Disease. Chem. Rev. 2012, 112, 5193–5239. [Google Scholar] [CrossRef] [Green Version]
  26. Kaumbekova, S.; Torkmahalleh, M.A.; Sakaguchi, N.; Umezawa, M.; Shah, D. Effect of ambient polycyclic aromatic hydrocarbons and nicotine on the structure of Aβ42 protein. Front. Environ. Sci. Eng. 2022, 17, 15. [Google Scholar] [CrossRef]
  27. Wallin, C.; Sholts, S.B.; Österlund, N.; Luo, J.; Jarvet, J.; Roos, P.M.; Ilag, L.; Gräslund, A.; Wärmländer, S.K. Alzheimer’s disease and cigarette smoke components: Effects of nicotine, PAHs, and Cd(II), Cr(III), Pb(II), Pb(IV) ions on amyloid-β peptide aggregation. Sci. Rep. 2017, 7, 14423. [Google Scholar] [CrossRef] [Green Version]
  28. Gao, D.; Wu, M.; Wang, C.; Wang, Y.; Zuo, Z. Chronic exposure to low benzo[a]pyrene level causes neurodegenerative disease-like syndromes in zebrafish (Danio rerio). Aquat. Toxicol. 2015, 167, 200–208. [Google Scholar] [CrossRef]
  29. Gao, D.; Wang, C.; Xi, Z.; Zhou, Y.; Wang, Y.; Zuo, Z. Early-Life Benzo[a]Pyrene Exposure Causes Neurodegenerative Syndromes in Adult Zebrafish (Danio rerio) and the Mechanism Involved. Toxicol. Sci. 2017, 157, 74–84. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, D.; Zhao, Y.; Qi, Y.; Gao, Y.; Tu, D.; Wang, Y.; Gao, H.-M.; Zhou, H. Benzo(a)pyrene exposure induced neuronal loss, plaque deposition, and cognitive decline in APP/PS1 mice. J. Neuroinflammation 2020, 17, 258. [Google Scholar] [CrossRef]
  31. Abraham, M.J.; Murtola, T.; Schulz, R.; Pall, S.; Smith, J.C.; Hess, B.; Lindahl, E. Gromacs: High performance molecular simulations through multi-level parallelism from laprtops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  32. Gerben, S.; Lemkul, J.; Brown, A.M.; Bevan, D. Comparing atomistic molecular mechanics force fields for a difficult target: A case study on the Alzheimer’s amyloid β-peptide. J. Biomol. Struct. Dyn. 2013, 32, 1817–1832. [Google Scholar] [CrossRef] [PubMed]
  33. Malde, A.K.; Zuo, L.; Breeze, M.; Stroet, M.; Poger, D.; Nair, P.C.; Oostenbrink, C.; Mark, A.E. An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0. J. Chem. Theory Comput. 2011, 7, 4026–4037. [Google Scholar] [CrossRef]
  34. Hess, B.; Bekker, H.; Berendsen, H.J.; Fraaije, J.G. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  35. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  36. Brown, A.M.; Bevan, D.R. Molecular Dynamics Simulations of Amyloid β -Peptide (1-42): Tetramer Formation and Membrane Interactions. Biophys. J. 2016, 111, 937–949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Östersund, N.; Kulkarni, Y.S.; Misiaszek, A.D.; Wallin, C.; Krüger, D.M.; Liao, Q.; Mashayekhy Rad, F.; Jarvet, J.; Strodel, B.; Wärmländer, S.K.; et al. Amyloid-β Peptide Interactions with Amphiphilic Surfactants: Electrostatic and Hydrophobic Effects. ACS Chem. Neurosci. 2018, 9, 1680–1692. [Google Scholar] [CrossRef]
  38. Berhanu, W.M.; Hansmann, U.H.E. Structure and Dynamics of Amyloid-β Segmental Polymorphisms. PLoS ONE 2012, 7, e41479. [Google Scholar] [CrossRef] [Green Version]
  39. Jokar, S.; Erfani, M.; Bavi, O.; Khazaei, S.; Sharifzadeh, M.; Hajiramezanali, M.; Beiki, D.; Shamloo, A. Design of peptide-based inhibitor agent against amyloid-β aggregation: Molecular docking, synthesis and in vitro evaluation. Bioorganic Chem. 2020, 102, 104050. [Google Scholar] [CrossRef]
  40. Grasso, G.; Lionello, C.; Stojceski, F. Highlighting the effect of amyloid beta assemblies on the mechanical properties and conformational stability of cell membrane. J. Mol. Graph. Model. 2020, 100, 107670. [Google Scholar] [CrossRef]
  41. Paul, A.; Samantray, S.; Anteghini, M.; Khaled, M.; Strodel, B. Thermodynamics and kinetics of the amyloid-β peptide revealed by Markov state models based on MD data in agreement with experiment. Chem. Sci. 2021, 12, 6652–6669. [Google Scholar] [CrossRef] [PubMed]
  42. Chen, G.F.; Xu, T.H.; Yan, Y.; Zhou, Y.R.; Jiang, Y.; Melcher, K.; Xu, H.E. Amyloid beta: Structure, biology and structure-based therapeutic development. Acta Pharmacol. Sin. 2017, 38, 1205–1235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Tung, N.T.; Derreumaux, P.; Vu, V.V.; Nam, P.C.; Ngo, S.T. C-Terminal Plays as the Possible Nucleation of the Self-Aggregation of the S-Shape Aβ. ACS Omega 2019, 4, 11066–11073. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Nguyen, H.L.; Krupa, P.; Hai, N.M.; Linh, H.Q.; Li, M.S. Structure and Physicochemical Properties of the Aβ42 Tetramer: Multiscale Molecular Dynamics Simulations. J. Phys. Chem. B 2019, 123, 7253–7269. [Google Scholar] [CrossRef]
  45. Fatafta, H.; Khaled, M.; Owen, M.C.; Sayyed-Ahmad, A.; Strodel, B. Amyloid-β peptide dimers undergo a random coil to β-sheet transition in the aqueous phase but not at the neuronal membrane. Proc. Natl. Acad. Sci. USA 2021, 118, e2106210118. [Google Scholar] [CrossRef]
  46. Aitken, J.F.; Loomes, K.M.; Konarkowska, B.; Cooper, G.J.S. Suppression by polycyclic compounds of the conversion of human amylin into insoluble amyloid. Biochem. J. 2003, 374 Pt 3, 779–784. [Google Scholar] [CrossRef] [Green Version]
  47. Lu, Q.; Tang, Q.; Xiong, Y.; Qing, G.; Sun, T. Protein/Peptide Aggregation and Amyloidosis on Biointerfaces. Materials 2016, 9, 740. [Google Scholar] [CrossRef]
Toxics 10 00786 g001 550
Figure 1. Time-evolution of (A,B) average interpeptide distances during 500 ns of the simulation, (C) average distances between Aβ42 peptides and B[a]P molecules during 500 ns of the simulation.
Figure 1. Time-evolution of (A,B) average interpeptide distances during 500 ns of the simulation, (C) average distances between Aβ42 peptides and B[a]P molecules during 500 ns of the simulation.
Toxics 10 00786 g001
Toxics 10 00786 g002 550
Figure 2. Time-evolution of (A) Formation of interpeptide clusters during 500 ns of the simulation, (B) formation of clusters of Aβ42 peptides and B[a]P molecules during 500 ns of the simulation.
Figure 2. Time-evolution of (A) Formation of interpeptide clusters during 500 ns of the simulation, (B) formation of clusters of Aβ42 peptides and B[a]P molecules during 500 ns of the simulation.
Toxics 10 00786 g002
Toxics 10 00786 g003 550
Figure 3. Composition of the secondary structure of the peptides averaged among the last 30 ns of the simulations.
Figure 3. Composition of the secondary structure of the peptides averaged among the last 30 ns of the simulations.
Toxics 10 00786 g003
Toxics 10 00786 g004 550
Figure 4. Representative snapshots of (A) four Aβ42 monomers before the simulation; (B) interpeptide cluster of four Aβ42 peptide monomers in the end of the simulation with no B[a]P; intermolecular cluster of four Aβ42 peptide monomers and B[a]P in the end of the simulations with (C) 4 B[a]P molecules; (D) 10 B[a]P molecules; (E) 40 B[a]P molecules. Color index: 1. Secondary structure: beta sheet = yellow, bridge − beta = tan, alpha helix = purple, 3_10_Helix = blue, Pi-Helix = red, turn, bend = cyan, coil = white, 2. B[a]P molecule = grey.
Figure 4. Representative snapshots of (A) four Aβ42 monomers before the simulation; (B) interpeptide cluster of four Aβ42 peptide monomers in the end of the simulation with no B[a]P; intermolecular cluster of four Aβ42 peptide monomers and B[a]P in the end of the simulations with (C) 4 B[a]P molecules; (D) 10 B[a]P molecules; (E) 40 B[a]P molecules. Color index: 1. Secondary structure: beta sheet = yellow, bridge − beta = tan, alpha helix = purple, 3_10_Helix = blue, Pi-Helix = red, turn, bend = cyan, coil = white, 2. B[a]P molecule = grey.
Toxics 10 00786 g004
Toxics 10 00786 g005 550
Figure 5. (A) RMSF of Aβ42 peptide residues, averaged among four peptides in the systems under the study, in the last 30 ns of the simulations, (B) Time-evolution of Radius of Gyration (RoG) of Aβ42 peptides, averaged among four peptides in the systems under the study, (C) Time-evolution of Solvent Accessible Surface Area (SASA) of Aβ42 peptides within 500 ns of the simulation, (D) Time-evolution of SASA of Aβ42 peptides within first 20 ns of the simulation in the systems with no B[a]P, with 4 B[a]P and 10 B[a]P molecules.
Figure 5. (A) RMSF of Aβ42 peptide residues, averaged among four peptides in the systems under the study, in the last 30 ns of the simulations, (B) Time-evolution of Radius of Gyration (RoG) of Aβ42 peptides, averaged among four peptides in the systems under the study, (C) Time-evolution of Solvent Accessible Surface Area (SASA) of Aβ42 peptides within 500 ns of the simulation, (D) Time-evolution of SASA of Aβ42 peptides within first 20 ns of the simulation in the systems with no B[a]P, with 4 B[a]P and 10 B[a]P molecules.
Toxics 10 00786 g005
Toxics 10 00786 g006 550
Figure 6. Radial distribution function (rdf) of (A) interpeptide interactions, (B) peptide—B[a]P ineractions in the systems under the study.
Figure 6. Radial distribution function (rdf) of (A) interpeptide interactions, (B) peptide—B[a]P ineractions in the systems under the study.
Toxics 10 00786 g006
Table
Table 1. Number of molecules in the simulated systems.
Table 1. Number of molecules in the simulated systems.
SystemAβ42B[a]PB[a]P ConcentrationH2ONa+Cl−
4 Aβ42400 mM42,296128120
4 Aβ42 + 4 B[a]P445.00 mM42,270128120
4 Aβ42 + 10 B[a]P41012.5 mM42,179128120
4 Aβ42 + 40 B[a]P44050.0 mM41,798128120
Table
Table 2. Solvent Accessible Surface Area (SASA) of peptides in the beginning (SASA0 ns), and in the end (SASA500 ns) of the simulation, minimum SASA (SASAmin) values of peptides, and an average number of interpetide H—bonds observed in the last 30 ns of the simulation.
Table 2. Solvent Accessible Surface Area (SASA) of peptides in the beginning (SASA0 ns), and in the end (SASA500 ns) of the simulation, minimum SASA (SASAmin) values of peptides, and an average number of interpetide H—bonds observed in the last 30 ns of the simulation.
SASA0 ns (nm2)SASA500 ns (nm2)SASAmin (nm2)H-Bonds (Last 30 ns)
4 Aβ42171.694.6087.30115 ± 5
4 Aβ42 + 4 B[a]P169.0102.795.10119 ± 5
4 Aβ42 + 10 B[a]P169.6110.0101.5113 ± 5
4 Aβ42 + 40 B[a]P166.5139.2132.7109 ± 5
Table
Table 3. Short-range (SR) and long-range (LR) Coulombic and Lennard-Jones potential between peptides, and peptide—B[a]P, observed in the last 10 ns of the simulations in the systems under the study.
Table 3. Short-range (SR) and long-range (LR) Coulombic and Lennard-Jones potential between peptides, and peptide—B[a]P, observed in the last 10 ns of the simulations in the systems under the study.
SystemPeptide–Peptide (kJ/mol)Peptide–B[a]P (kJ/mol)
Coul–SRLJ–SRCoul–LRLJ–LRCoul–SRLJ–SR
4 Aβ42−44,048 ± 26−3712 ± 1127,804 ± 8−38.5 ± 3.4--
4 Aβ42 + 4 B[a]P−44,333 ± 10−3577 ± 728,052 ± 12−44.1 ± 1.9−49.2 ± 2.2−423 ± 2
4 Aβ42 + 10 B[a]P−43,963 ± 40−3536 ± 928,042 ± 7−7.9 ± 3−60.6 ± 2.8−762 ± 16
4 Aβ42 + 40 B[a]P−43,945 ± 22−3016 ± 728,043 ± 69.8 ± 3−191.7 ± 3−2020 ± 8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kaumbekova, S.; Torkmahalleh, M.A.; Shah, D. Ambient Benzo[a]pyrene’s Effect on Kinetic Modulation of Amyloid Beta Peptide Aggregation: A Tentative Association between Ultrafine Particulate Matter and Alzheimer’s Disease. Toxics 2022, 10, 786. https://doi.org/10.3390/toxics10120786

AMA Style

Kaumbekova S, Torkmahalleh MA, Shah D. Ambient Benzo[a]pyrene’s Effect on Kinetic Modulation of Amyloid Beta Peptide Aggregation: A Tentative Association between Ultrafine Particulate Matter and Alzheimer’s Disease. Toxics. 2022; 10(12):786. https://doi.org/10.3390/toxics10120786

Chicago/Turabian Style

Kaumbekova, Samal, Mehdi Amouei Torkmahalleh, and Dhawal Shah. 2022. "Ambient Benzo[a]pyrene’s Effect on Kinetic Modulation of Amyloid Beta Peptide Aggregation: A Tentative Association between Ultrafine Particulate Matter and Alzheimer’s Disease" Toxics 10, no. 12: 786. https://doi.org/10.3390/toxics10120786

APA Style

Kaumbekova, S., Torkmahalleh, M. A., & Shah, D. (2022). Ambient Benzo[a]pyrene’s Effect on Kinetic Modulation of Amyloid Beta Peptide Aggregation: A Tentative Association between Ultrafine Particulate Matter and Alzheimer’s Disease. Toxics, 10(12), 786. https://doi.org/10.3390/toxics10120786

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop