Research article Special Issues

The switch point algorithm applied to a harvesting problem

  • In this paper, we investigate an optimal harvesting problem of a spatially explicit fishery model that was previously analyzed. On the surface, this problem looks innocent, but if parameters are set to where a singular arc occurs, two complex questions arise. The first question pertains to Fuller's phenomenon (or chattering), a phenomenon in which the optimal control possesses a singular arc that cannot be concatenated with the bang-bang arcs without prompting infinite oscillations over a finite region. 1) How do we numerically assess whether or not a problem chatters in cases when we cannot analytically prove such a phenomenon? The second question focuses on implementation of an optimal control. 2) When an optimal control has regions that are difficult to implement, how can we find alternative strategies that are both suboptimal and realistic to use? Although the former question does not apply to all optimal harvesting problems, most fishery managers should be concerned about the latter. Interestingly, for this specific problem, our techniques for answering the first question results in an answer to the the second. Our methods involve using an extended version of the switch point algorithm (SPA), which handles control problems having initial and terminal conditions on the states. In our numerical experiments, we obtain strong empirical evidence that the harvesting problem chatters, and we find three alternative harvesting strategies with fewer switches that are realistic to implement and near optimal.

    Citation: Summer Atkins, William W. Hager, Maia Martcheva. The switch point algorithm applied to a harvesting problem[J]. Mathematical Biosciences and Engineering, 2024, 21(5): 6123-6149. doi: 10.3934/mbe.2024269

    Related Papers:

    [1] Ivana Sirangelo, Margherita Borriello, Gaetano Irace, Clara Iannuzzi . Intrinsic blue-green fluorescence in amyloyd fibrils. AIMS Biophysics, 2018, 5(2): 155-165. doi: 10.3934/biophy.2018.2.155
    [2] Alessandro Didonna, Federico Benetti . Post-translational modifications in neurodegeneration. AIMS Biophysics, 2016, 3(1): 27-49. doi: 10.3934/biophy.2016.1.27
    [3] Nick K. Olrichs, J. Bernd Helms . Novel insights into the function of the conserved domain of the CAP superfamily of proteins. AIMS Biophysics, 2016, 3(2): 232-246. doi: 10.3934/biophy.2016.2.232
    [4] Jany Dandurand, Angela Ostuni, Maria Francesca Armentano, Maria Antonietta Crudele, Vincenza Dolce, Federica Marra, Valérie Samouillan, Faustino Bisaccia . Calorimetry and FTIR reveal the ability of URG7 protein to modify the aggregation state of both cell lysate and amylogenic α-synuclein. AIMS Biophysics, 2020, 7(3): 189-203. doi: 10.3934/biophy.2020015
    [5] Nicola Gaetano Gatta, Gaetano Cammarota, Vittorio Gentile . Possible roles of transglutaminases in molecular mechanisms responsible for human neurodegenerative diseases. AIMS Biophysics, 2016, 3(4): 529-545. doi: 10.3934/biophy.2016.4.529
    [6] J. Dandurand, E. Dantras, C. Lacabanne, A. Pepe, B. Bochicchio, V. Samouillan . Thermal and dielectric fingerprints of self-assembling elastin peptides derived from exon30. AIMS Biophysics, 2021, 8(3): 236-247. doi: 10.3934/biophy.2021018
    [7] Daniela Meleleo, Cesare Sblano . Influence of cholesterol on human calcitonin channel formation. Possible role of sterol as molecular chaperone. AIMS Biophysics, 2019, 6(1): 23-38. doi: 10.3934/biophy.2019.1.23
    [8] Marek K. Korzeniowski, Barbara Baird, David Holowka . STIM1 activation is regulated by a 14 amino acid sequence adjacent to the CRAC activation domain. AIMS Biophysics, 2016, 3(1): 99-118. doi: 10.3934/biophy.2016.1.99
    [9] Nicole Shamitko-Klingensmith, Jonathan W. Boyd, Justin Legleiter . Microtubule modification influences cellular response to amyloid-β exposure. AIMS Biophysics, 2016, 3(2): 261-285. doi: 10.3934/biophy.2016.2.261
    [10] Yu-Ting Huang, Hui-Fen Liao, Shun-Li Wang, Shan-Yang Lin . Glycation and secondary conformational changes of human serum albumin: study of the FTIR spectroscopic curve-fitting technique. AIMS Biophysics, 2016, 3(2): 247-260. doi: 10.3934/biophy.2016.2.247
  • In this paper, we investigate an optimal harvesting problem of a spatially explicit fishery model that was previously analyzed. On the surface, this problem looks innocent, but if parameters are set to where a singular arc occurs, two complex questions arise. The first question pertains to Fuller's phenomenon (or chattering), a phenomenon in which the optimal control possesses a singular arc that cannot be concatenated with the bang-bang arcs without prompting infinite oscillations over a finite region. 1) How do we numerically assess whether or not a problem chatters in cases when we cannot analytically prove such a phenomenon? The second question focuses on implementation of an optimal control. 2) When an optimal control has regions that are difficult to implement, how can we find alternative strategies that are both suboptimal and realistic to use? Although the former question does not apply to all optimal harvesting problems, most fishery managers should be concerned about the latter. Interestingly, for this specific problem, our techniques for answering the first question results in an answer to the the second. Our methods involve using an extended version of the switch point algorithm (SPA), which handles control problems having initial and terminal conditions on the states. In our numerical experiments, we obtain strong empirical evidence that the harvesting problem chatters, and we find three alternative harvesting strategies with fewer switches that are realistic to implement and near optimal.



    Alzheimer's disease (AD) is a fatal neurodegenerative disorder that is the most prevalent form of dementia. The neuropathological and neurochemical hallmarks of AD include: Synaptic loss and selective neuronal cell death; decreases in markers for certain neurotransmitters; and abnormalities in neurons and their processes as well as in the extracellular space. Two of the main features associated with AD are neurofibrillary tangles comprised of the protein tau and cerebrovascular, diffuse, and neuritic plaques composed predominantly of the amyloidogenic peptide amyloid-β (Aβ). These proteinaceous deposits of tau and Aβ consist of stable amyloid fibrils, which are β-sheet rich fibrous protein agammaegates. Similar deposition of amyloid is associated with numerous other diseases [1].

    In the early 90's, the amyloid cascade hypothesis was introduced, which postulated that Aβ agammaegation and deposition directly lead to neuronal death, resulting in AD [2]. Yet, the correlation between plaques and cognitive dysfunction in AD has been questionable for years [3][7], and with Aβ-directed therapeutic strategies failing in numerous clinical trials, the role of Aβ in AD progression is being re-evaluated [8][11]. This has led to an enhanced research focus on diffuse, soluble agammaegates of Aβ. Various small, globular agammaegates of Aβ, referred to broadly as Aβ oligomers (AβO), were detected in AD patients a few decades ago [12][14], and these AβO were originally classified as being intermediates toward the formation of amyloid fibrils. Over the years, increasing evidence points toward AβOs playing a central role in AD, as AβOs correlate more strongly with AD progression in patients and animal models [15][19]. For example, AβO formation and buildup occurs early compared with plaque buildup (much earlier than plaque deposition) in the AD brain [20][22]and CSF [23], which has led to extensive efforts to develop assays to detect AβOs for potential application as biomarkers [24][31].

    An extensive body of evidence has linked high levels of AβOs in the brain to a variety of pathogenic consequences associated with AD (summarized in Table 1). As such, there has been significant effort made to characterize AβO formation, structure, and biochemical/biophysical characteristics (such as interactions with other proteins and lipids) in the hope that underlying modes of AβO-related toxicity could be revealed [32][37]. As AβOs are metastable, highly heterogeneous in nature, and can form via a variety of different pathways, this remains a challenging task, yet progress has been made. Here, we review AβOs with a focus on biophysical characterization of their formation and structure.

    Table 1.  Toxic mechanisms associated with AβOs.
    Toxic effects of AβOs Model system(s) References
    Reduction in neural plasticity mice, rat [37][40]
    Stimulation of tau phosphorylation cortical neurons (rat), hippocampal neurons, neuroblastoma, primary neurons, Tg-Mice [41][45]
    Choline acetyltransferase Inhibition cholinergic cell lines [46],[47]
    Oxidative stress cortical neurons, hippocampal neurons, in vitro, neuroblastoma [48][51]
    Endoplasmic reticulum stress cortical astrocytes, cortical neurons (rat), fibroblasts, Tg-mice [43],[52],[53]
    Receptor disturbance cortical neurons, hippocampal neurons [54][57]
    Insulin resistance cortical neurons, hippocampal neurons, rat [40],[42],[57],[58]
    Synapse deterioration hippocampal neurons, pyramidal neurons (rat), Tg-mice [21],[35],[54],[59]
    Axonal transport cortical neurons, hippocampal, in vitro, Tg-mice [60][62]
    Astrocytes/microglia effects cortical astrocytes, Tg-mice [44],[52],[63],[64]
    Cell cycle disruption cortical neurons, Tg-mice [65],[66]
    Selective neuron death Mice [38],[67]
    Inhibition of long-term potentiation Mice [37],[40]
    Calcium dysregulation cortical neurons (rat), hippocampal neurons, primary neurons [43],[45],[48],[52]
    Modulation of metal toxicity neuroblastoma [49]
    Cytoskeleton disruption primary neurons [50]
    Modulation of Receptor/Channel Activity hippocampal neurons [48]

     | Show Table
    DownLoad: CSV

    Aβ is an approximately 4 kDA peptide (typically 40–42 amino acids long) that is derived from the transmembrane portion of the amyloid precursor protein (APP; Figure 1). The production of Aβ is achieved by the sequential cleavage of APP by two membrane-bound endoprotease activities, β- and γ-secretase. The two predominant Aβ peptides produced are 40 and 42 amino acids in length, and these peptides are referred to as Aβ40 and Aβ42 respectively. Aβ is amphipathic in nature (having a predominately hydrophilic N-terminus and a predominately hydrophobic C-terminus), which is thought to drive its agammaegation. As the C-terminal end of Aβ coincides with the transmembrane portion of APP, Aβ42 has a larger hydrophobic domain, making it more fibrillogenic compared to Aβ40 and deposits to a much greater extent in the brain [68][70]. Only about 10% of APP is processed via this Aβ producing pathway. Most APP is cleaved by the α-secretase, generating a series of much more benign peptide fragments.

    The hydrophilic N-terminal region of Aβ can adopt both an α-helical or β-sheet structure dependent on solution conditions, for example pH [71],[72]. The hydrophobic C-terminal end of Aβ has a propensity to adopt β-sheet structure upon agammaegation independent of solvent conditions [71],[72]. Beyond its amphipathic nature, several other domains have been identified in Aβ. The different polyomorphic fibril structures of Aβ are comprised of bundled β-sheets with backbones orthogonal to the fiber axis creating a cross-β structure [73], and two β-strand forming domains (residues 11–21 and 29–39 respectively) that are separated by a turn/bend region (around residues 23–26) identified through various experimental and computational studies [74][77]. The central region of Aβ (residues 16–21), contained within one of the β-strands, has enhanced amyloidogenic properties and represents a hydrophobic core [78].

    Figure 1.  APP processing and subsequent Aβ amino acid sequence with specific domains of interest specified. Proteolytic cleavage of APP is initiated at residue 671 by β-secretase followed by either non-amyloidogenic processing, cleavage at residue 687 by α-secretase, or amyloidogenic processing, cleavage at residue 711, 713, 714, or 726 by γ- secretase. The amino acid sequence below highlights the region of APP from which Aβ is produced with the numbering referring the residues in Aβ. The orange highlighted region represents the intact Aβ42 sequence. Hydropathy indexes of individual residues are color coded according to hydrophobic (blue), slightly hydrophobic (light green), and hydrophilic (red). Regions of interest and cites of secretase activity within the APP and Aβ are indicated.

    The agammaegation of Aβ (and other amyloid-forming proteins) is typically characterized in terms of fibril formation (Figure 2). Aβ fibril formation occurs via a complex agammaegation pathway. Fundamentally, AβOs can be subdivided into species that are intermediates in fibrils formation (referred to as being on pathway) or species that do not directly lead to fibrils (referred to as being off pathway). This is a contributing factor to the immense heterogeneity observed in AβO populations as will be discussed in more detail later. In general, amyloid formation proceeds via a nucleation dependent polymerization mechanism [79][81]. With this mechanism, agammaegation initially occurs via a slow nucleation phase (often called the lag phase) that involves the formation of a thermodynamically unfavorable critical nucleus that is associated with a transition from a native to non-native protein conformation. For Aβ, the critical nucleus is likely a multimeric species [82][84]. Once the critical nucleus has formed, an elongation or growth phase (characterized by a relatively rapid extension of fibril agammaegates) occurs. While numerical models can extract important parameters (lag phase times, elongation rates, critical nucleus size) from experimental data [85], the actual agammaegation pathway toward fibrils can be complicated. For example, other on pathway intermediates, like protofibrils, are also observed in Aβ agammaegation. Protofibrils are amyloid-like, elongated agammaegates with filament-like morphologies and are late-stage intermediate precursors on the agammaegation pathway to fibrils. A key aspect that facilitates on and off pathway agammaegation routes and complicates investigations of AβOs is that they possess structural plasticity and are metastable and transient in nature.

    Figure 2.  Production and aggregation Aβ. Aβ is a cleavage product of APP, a transmembrane protein. Monomeric Aβ transitions between ordered and disordered states. Once dimerization occurs, subsequent aggregation occurs either on pathway or off pathway with respect to fibril formation. The aggregates associated with these different pathways increase in molecular weight from the left to the right of the schematic. Fibrillization can proceed via several potential pathways that can populate various intermediate aggregate states, including oligomers and protofibrils. Off-pathway oligomers of various size may also form. Annular aggregates of Aβ can also form and are thought to potentially be associated with forming pore-like structures.

    Further complicating the issue is the observation that Aβ can agammaegate into a variety of morphologically distinct fibril structures, referred to as polymorphs [86][90]. This phenomenon is predicated on subtle changes of environmental conditions associated with agammaegation, and as a result, preparatory protocols employed in experiments determines the resulting Aβ fibril morphology [87]. While polymorphic agammaegates are readily observed with in vitro studies using synthetic Aβ, polymorphic structures have been observed in amyloids derived from tissue, and it is thought that variations in Aβ agammaegate morphologies may play an important role in AD [91],[92]. For example, polymorphic agammaegates and fibrils may result in distinct biological activities and levels of toxicity that could underlie variations in AD [76], and distinct fibril structures can be directly associated with individual AD patients and clinical phenotype [93],[94].

    A complicating factor in evaluating and comparing studies aimed at elucidating AβO formation, structure, and physiological impact is divergent experimental conditions, such as Aβ preparation protocols, heavily influences experimental outcomes. In terms of oligomers, the emergence of distinct fibril structure strongly suggests that there would also be distinct oligomeric precursors associated with their formation. Indeed, distinct AβO species can be observed within in vitro agammaegation assays under conditions that result in fibril polymorphs [95] (Figure 3). For studies conducted with synthetic peptide, there are a variety of protocols used to solubilize Aβ (Table 2). Typically, these protocols consist of a disagammaegation step and a reconstitution step. The disagammaegation steps usually involves the use of hydrogen bond disrupting solvents, i.e., hexafluoroisopropanol (HFIP) or trifluoroacetic acid (TFA), that break down pre-existing agammaegates within lyophilized stocks of Aβ. These solvents are often removed under vacuum, leaving a peptide film. The reconstitution step involves dissolving these peptide films into a solvent that facilitates dilution into an appropriate buffer. Dimethyl sulfoxide (DMSO) is often used, creating a concentrated stock that is diluted into the desired buffer. Disagammaegation and reconstitution can also be facilitated by changes in pH. Sometimes reconstitution is performed directly into the final buffer. Importantly, variations in preparatory protocols indeed result in different populations of oligomers (Table 2 and Figure 3), and these variations can complicate direct comparisons between different reports within the literature. Further complicating the issue, there are often distinctions observed between studies performed with synthetic Aβ and naturally derived Aβ [96]. In many studies aimed at elucidating activity of AβOs, specific preparations are used to obtain a particular population of oligomer species. These are then directly applied to different model systems, ranging from cell culture to animal models. However, there is often a lack of effort to verify that once added to the model system (which can often be a pronounced change in chemical environment) that these AβO species do not dissociate or agammaegate into a different AβO or Aβ agammaegate. To truly relate specific AβOs to a neurotoxic activity, effort should be invoked to attempt additional controls of this type.

    Table 2.  Representative disaggregation, reconstitution, and miscellaneous protocols for the preparation of Aβ and observed AβOs.
    Classification Disaggregation Reconstitution Miscellaneous Result References
    ADDLs None F12 Media 4 °C Centrifuged 14,000 ×g for 10 mins 5–6 nm by AFM (height) [38]
    ADDLs HFIP DMSO at 5 mM Aβ Sonicate 5 mins, dilute with DMEM/F12 Media A11+ [97]
    Globulomers HFIP DMSO at 5 mM Aβ Sonicate 10 mins, dilute PBS + 0.05% SDS 16–56 kDa by SDS-PAGE [98]
    Globulomers HFIP DMSO at 5 mM Aβ Dilute PBS + 0.05% SDS, Dialyze 38–48 kDa by SDS-PAGE [99]
    Aβ*56 HFIP DMSO at 5 mM Aβ Sonicate 20 min, PBS + 0.2% SDS incubate 6 h, dilute and incubated 18 h; centrifuge 3000 ×g, dialysis 56 kDa by Native-PAGE [100]
    AβOs HFIP DMSO F12 Media incubated at 4 °C for 24 h 1–4 nm by AFM [101],[102]
    AβOs HFIP/NH4OH 10 mM Tris-HCl Addition of Zn2+ 10–12 nm by AFM [97]
    AβOs TFA/HFIP (2Xs) 2 mM NaOH PBS Centrifuge 386,000 ×g 2.5 nm by AFM [87],[95]
    AβOs TFA salt/2 mM NaOH, pH ∼10.5, 1 min sonication HH PBS - 3–12 nm by AFM [103]
    AβOs 10% NH4OH (w/v) sonicated 5 mins, lyophilized 60 mM NaOH - 1–10 nm by DLS [104]
    AβOs LMW - DMSO Sonicate 1 min, centrifuge 16,000 ×g, SEC, PBS, PICUP crosslinking 4–26 kDa by SDS-PAGE [105]

     | Show Table
    DownLoad: CSV
    Figure 3.  Heterogeneity of AβOs. A series of atomic force microscopy images and size analysis, i.e. height histograms, of AβOs formed from synthetic Aβ that had been prepared by some of the protocols described in Table 2. The AβOs were prepared in the following way: (A) protocol provided in the AggreSure β-amyloid kits available from AnaSpec which consists of a reconstitution step directly into Tris buffer with bath sonication; (B) 10% NH4OH disaggregation buffer, followed by additional treatment with HFIP, and reconstitution in 2 mM NaOH (pH > 11) with subsequent dilution into Tris buffer; (C) 10% NH4OH disaggregation buffer, and reconstitution in 60 mM NaOH (pH > 11) with subsequent dilution into HEPES buffer; (D) HFIP disaggregation, reconstitution in DMSO, followed by dilution into PBS; (E) TFA disaggregation with sonication, an additional HFIP disaggregation step, reconstitution with NaOH, and dilution into PBS; (F) No disaggregation step and direct reconstitution in to phosphate buffer.

    Enormous efforts have been extended in identifying the most toxic/disease-relevant AβO species and the relevant underlying structure [9],[106],[107]. This is a daunting endeavor due to the transient nature and extensive heterogeneity associated with AβOs [108][110]. It is possible that different AβO species may activate different deleterious changes associated with AD [9],[106],[107]; however, there may also be relatively benign AβO species or even experimental artifacts [96],[108],[111]. Unraveling the physiological activity of AβOs may require separate analysis of these different species, which may be altered by preparatory protocols. The ability of smaller AβO (like dimers and trimers) to further agammaegate into large, more stable synaptotoxic assemblies [107] must be accounted for in assessing the toxic effects of specific AβO species. That is, upon the addition of preparations of AβO to cell culture or other models, AβOs may further assemble into other higher order species that may influence the associated toxic effects. Despite the inherent complexity of this endeavor, progress has been made. Aβ trimers have been linked to playing a role in inducing pathological conformational changes in tau [112]. However, crosslinked Aβ dimers were shown not to be toxic themselves, but rather contributed to toxicity by further assembling into larger assemblies [113]. A 56 kDa SDS-stable Aβ0 (referred to as Aβ*56) has been identified as a prominent specie in the AD brain [16], CSF [22], and in transgenic mouse models of AD [32],[114]. In terms of an actual biological activity, Aβ*56 interacts with N-methyl-D-aspartate receptors (NMDARs), increasing NMDAR-dependent Ca2+ influx and activation of Ca2+/calmodulin-dependent kinase IIα (CAMKIIα) [115]. Activation of CAMKIIα correlates with enhanced site-specific phosphorylation and mis-sorting of tau [115]. Smaller AβOs, namely dimers and trimers, do not appear to elicit these specific effects.

    There appears to be some common themes emerging with respect to toxicity. Toxic AβOs appear to react with oligomer specific antibodies like A11 (generic for amyloid oligomers in general, [116]) and NU4 (specific for AβOs) [117]; whereas, nontoxic AβOs demonstrate reactivity with anti-fibril antibodies like OC [116]. Importantly, toxic AβOs appear to be unrelated to plaques [116],[118]. AβOs related to amyloid plaques temporally, spatially, and structurally are nontoxic [118]. A number of toxic AβO species are larger than 50 kDa [21],[54],[118], like the previously mentioned Aβ*56 [32]. Smaller AβOs appear to be less toxic or even benign [21],[54],[118],[119], except for their ability to further agammaegated into larger assemblies. This has led to distinguishing between high molecular weight (HMW) and low molecular weight (LMW) oligomers [120]. Agammaegation mechanisms differentiating between the eventual formation of HMW and LMW AβOs already appear to deviate at the dimer stage [121]. Furthermore, LMW and HMW Aβ oligomers differentially impact synapses and memory [122],[123]; although, LMW AβOs are not always associated with memory dysfunction [116],[117]. HMW AβOs are the predominant Aβ species in the native soluble protein fraction of AD brains [124]. These HMW species in the AD brain sometimes appear to be constructed from smaller ∼7 kDa Aβ species [125]. Neurohistopathological and biochemical analyses of AβOs in the temporal cortex of AD brains implicated an Aβ dodecamer (∼55 kDa) [126]. HMW AβOs bind cultured synapses [21],[54],[118], induce reactive oxygen species (ROS) production [123], and disrupt memory function [116],[117]. With respect to the previously defined agammaegation pathways, LMW AβOs are typically on pathway to fibril formation; HMW AβOs are off-pathway [127],[128]. This is consistent with HMW AβOs being potent, as off pathway AβOs appear more toxic [129]. Collectively, these observations point to the complex interplay between different AβO species and their specific activity with respect to neurotoxicity.

    Table 3.  Methods used for analyzing AβOs.
    Technique Features of AβOs Obtained
    Atomic Force Microscopy (AFM) Morphology, population distributions
    Electron Microscopy (EM) Morphology, population distributions
    Ion-Mobility Mass Spectroscopy (IM-MS) Secondary structure, multimeric configurations
    Nuclear Magnetic Resonance 2D (NMR) Secondary structure, multimeric configurations, atomic resolution structure
    Electron Paramagnetic Resonance (EPR) Secondary structure, multimeric configurations, atomic resolution structure
    Powder X-ray Diffraction (PXRD) Atomic resolution structure
    Small Angle X-ray Scattering (SAXS) Atomic resolution structure
    Single Crystal X-ray Diffraction (SCXD) Atomic resolution structure
    SDS-Page Size distribution of multimers

     | Show Table
    DownLoad: CSV

    To facilitate structure/activity analysis, significant efforts have been made to obtain structural details of AβOs using a variety of methods (Table 3). This is required to fully elucidate the modes of interaction of AβOs with other biomolecules and related toxic mechanisms [32],[33],[35][37],[130]. Structural characterization of specific AβOs is challenging due to their transient nature and heterogeneity. AβOs can exhibit conformational plasticity that can be heavily influenced by environmental factors, further complicating such analysis. Despite these inherent challenges, efforts have been made to separate, isolate, and characterize distinct AβO species obtained from synthetic Aβ or from AD brain tissue and cell cultures [18],[32],[35][37],[114],[129],[131][136]. With regard to HMW and LMW AβOs, they are separable in vitro by size exclusion chromatography [137] or ultrafiltration with a 50 kDa molecular weight cutoff [21],[54],[118]. Quick characterization of AβOs can be accomplished by chromatographic techniques, SDS PAGE, and by the use of a variety of oligomer specific antibodies (Table 4).

    Table 4.  Anti-AβO specific antibodies.
    Antibody Epitope References
    A11 Soluble amyloid oligomers [17],[138]
    NAB61 Dimeric, oligomeric, higher order aggregates [139]
    NU-1 ADDLs specific [140]
    NU-2 ADDLs weak binding; no Aβ monomer staining [140]
    NU-4 ADDL trimer, tetramer, and 12–24 mer specific [140]

     | Show Table
    DownLoad: CSV

    One strategy that has been used to overcome obstacles associated with characterizing AβOs has been to stabilize oligomers via cross-linking. Such methods can lead to not only structural but also functional characterization of AβOs. Photo-induced crosslinking (PICUP) was initially used to stabilize and characterize LMW AβOs of Aβ40 and Aβ42 [33],[141][143].While Aβ40 primarily formed an equilibrium of dimers, trimers, and tetramers, Aβ42 agammaegated into pentamers/hexamers that further assembled into protofibrils [143]. Using mutated Aβ42 (F10, Y42), AβOs ranging from dimer to dodecamers were stabilized using PICUP [141], opening up the ability to perform structure/activity analysis on AβOs up to 50 kDa in size. AβOs can also be stabilized using dityrosine crosslinking, which occurs under elevated copper concentrations and oxidative stress [144]. Both of these conditions have been linked to AD, suggesting that the crosslinking associated with this method may be more physiologically relevant [145],[146]. In fact, dityrosine crosslinked proteins are observed in amyloid plaques and CSF from AD patients [145]. Dityrosine linkages are associated with copper-mediated stabilization of AβOs [147]. Copper stabilization is effective enough to allow for 3D structural characterization of AβOs by small-angle x-ray scattering [148], and the copper to Aβ ratio could push agammaegation toward ellipsoidal oligomers of 38 peptides (excess copper) or fibrils (excess Aβ) [119].

    Another method to overcome the metastability of Aβ species and control the agammaegation process is to complex/fuse Aβ sequences within other protein. Aβ sequences are often based on previously identified regions of the peptide that have been identified as playing a role in Aβ agammaegation (Figure 1). Such a strategy has been successful in studying monomeric structure of Aβ sequences [149]. With careful design, this method has been applied to stabilizing oligomers derived from Aβ fragments for structural characterization. Fusion of Aβ18–41 with the CDR3 loop region of a shark Ig new antigen receptor single variable domain antibody resulted in the formation and stabilization of tightly associated Aβ dimers, which could pair to form tetramers [150]. This dimer had a compact structure rather than a β-turn/β-sheet structure. Several engineered peptide macrocycles that incorporate Aβ sequences (Aβ15–23 and Aβ17–36) have been designed to contain agammaegation to the oligomeric state, and these systems have been structurally characterized by X-ray crystallography and NMR spectroscopy [151][153]. Collectively, this macrocycle strategy has revealed numerous dimers, trimers, tetramers, and higher order oligomer species that display a variety of β-sheet based structural heterogeneity.

    Due to the complex nature of Aβ agammaegation, techniques that allow for distinguishing and characterizing distinct morphological features of AβOs within heterogeneous agammaegation reactions are of enormous benefit. Both atomic force microscopy (AFM) and electron microscopy (EM) provide this capability. In particular, AFM has emerged as a particularly useful technique in studying AβO formation and morphology [95],[101],[134],[135],[154][160]. As AFM can be operated in solution, it has the ability to observe and track the behavior of individual AβOs on surfaces under physiological buffer conditions [158],[160]. The surfaces used in AFM experiments have become progressively more biologically relevant and include lipid membranes [135],[161][164]. With regard to Aβ agammaegation on lipid membranes, in solution AFM studies have demonstrated the formation of distinct oligomeric agammaegates associated with point mutation in Aβ [165], the formation of pore-like AβO morphologies [154][156], that preparation history influences AβO formation on bilayers [95], and that mechanical changes occur in bilayers associated with the presence of AβOs [166]. AFM based force spectroscopy has even been used to understand the energetics of AβO formation and stability [167],[168].

    With the recent development of high-speed AFM in solution, insights into the dynamics and fate of individual AβOs has been achieved [128],[169]. Using high-speed AFM to track the dynamics of PICUP-stabilized LMW AβOs demonstrated that these AβOs were highly dynamic in structure, fluctuating between single and multi-globular assemblies [169]. Direct visualization of Aβ agammaegation with high-speed AFM imaging demonstrates that LMW AβOs much more quickly transition to form fibrils with distinct morphologies compared with HMW AβOs [128]. Despite being classified as being off-pathway, HMW AβOs can still contribute to fibril formation by serving as a reservoir of Aβ. That is, HMW AβOs may dissociate into smaller LMW AβO that seed fibrillization [128]. This again points to the necessity to track the fate of AβO species when evaluating their physiological activities. For example, the LMW AβOs that are dissociation products of HMW AβOs may actually be more toxic [170].

    Another technique that is capable of characterizing heterogeneous populations of AβOs is ion mobility mass spectrometry (IM-MS). These IM-MS studies ascertained qualitative differences in AβO structure associated with Aβ40 and Aβ42 [131],[132]. Aβ40 tetramers displayed an enclosed ring-shaped configuration that would inhibit additional contacts required to assemble into larger AβO species [131],[132]. Aβ42 tetramers preferentially had a bent structure that would provide oligomer ends capable of additional contacts and enabling further agammaegation. Indeed, Aβ42 was capable of forming larger donut-shaped dodecamers.

    Another strategy to perform structural analysis of AβOs is utilizing specific conditions to stabilize them via a kinetic trap. Such an approach has been successfully used to enable NMR spectroscopy of AβOs [171],[172]. By incubation of Aβ at 4 °C and freeze-trapping with liquid nitrogen, heterogeneous, spherical AβOs were analyzed with 2D NMR and shown to possess in-register parallel β-strand structure similar to fibrils [172]. Based on NMR analysis, Aβ42 pentamers stabilized at a low 4 °C and 10 mM salt concentration were disordered [171]. These Aβ42 pentamers displayed enhanced toxicity compared with protofibrils or fibrils [171]. An atomic model of Aβ42 oligomers consisting of approximately 15–24 peptides has been proposed from a combination of biophysical techniques [173]. These oligomers were prepared by disagammaegation in HFIP followed by resuspension in dilute ammonium hydroxide, preventing fibril formation. These were not end-stage AβOs, as subsequent dilution in PBS resulted in fibril formation. Powder X-ray diffraction patterns of the AβOs were consistent with helical β-sheet pairs wrapped together into a super-helix. This wrapping results in a hole along the super-helix axis, which is consistent with proposed toxic mechanisms in which Aβ forms pathogenic pores. Another proposed structure, based on site-directed spin labeling and electron paramagnetic resonance, shares similarities with this structure [174]. This study used Aβ42 fused to GroES-ubiquitin that formed stable oligomers that were A11 positive. Based on the EPR data, an AβO model was proposed that consisted of a β-sheet with three antiparallel strands with these strands being arranged head to tail. These sheets are further packed face to back as a group of four.

    Even just tracking AβO formation has been challenging. A number of straight forward spectroscopic assays are well-established to track formation and kinetic parameters of fibrils (e.g., ThT), but methods to easily track AβO formation have been lacking. Recently, the use of 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene (BODIPY) dyes have been utilized to fluorescently detect the AβOs in vitro [175],[176]. These BODIPY fluorescent probes have high quantum yields with the capability of selectively binding to AβOs [177]. Importantly, increases and decreases in the “BD-Oligo” dye fluorescence also correlated with increasing and decreasing intensities of A11 staining, directly relating this signal to an established methods to detect AβOs [175]. In addition, the BD-Oligo dye can be used in parallel with ThT assays that detect fibril formation, allowing for the direct investigation of the correlation between oligomer and fibril formation [175]. A BODIPY-based probe (BAP-1) has also facilitated the direct visualization of Aβ plaques in transgenic mice [177], and rational modifications of BAP-1 also allowed for near-infrared selective detection of tau neurofibrillary tangles [178]. The ability to rationally modify BODIPY dyes for fluorescent detection of specific agammaegate species, the tunability of their spectroscopic properties, and their insensitivity to solvent and pH changes have made BODIPY dyes a recently promising avenue for high throughput studies of AβO formation and stability.

    While the agammaegation of Aβ has been extensively studied, there is still much to understand at the molecular level about AβO formation, structure, and activity. Due to the transient nature of AβOs, their morphological heterogeneity, and the continuing debate concerning specific toxic agammaegate species associated with AD, structural details and physiological activities of the variety of AβOs still need to be fully elucidated. The exact mechanisms associated of how AβOs lead to cellular dysfunction and death have not fully been explained. Understanding these phenomenon may prove crucial in the effectiveness of therapeutic strategies based on manipulating Aβ production, clearance, and agammaegation. Here, we highlighted some specific features of AβOs and techniques that have provided insight into their structure and formation. While this review is far from exhaustive, we hope that collectively they provide a compelling argument toward the importance of understanding the nature of AβOs, highlight some of the intrinsic obstacles associated with studying AβOs, and provide some insight into methods that will play a role in pushing our knowledge of AβOs further.



    [1] M. Demir, S. Lenhart, A spatial food chain model for the Black Sea anchovy, and its optimal fishery, Discrete Contin. Dyn. Syst., 26 (2021), 155–171. https://doi.org/10.3934/dcdsb.2020373 doi: 10.3934/dcdsb.2020373
    [2] R. Hilborn, Overfishing: What Everyone Needs to Know, Oxford University Press, Oxford, 2012.
    [3] M. R. Kelly, Y. Xing, S. Lenhart, Optimal fish harvesting for a population modeled by a nonlinear parabolic partial differential equation, Nat. Resour. Model., 29 (2016), 36–70. https://doi.org/10.1111/nrm.12073 doi: 10.1111/nrm.12073
    [4] G. N. Tuck, H. P. Possingham, Marine protected areas for spatially structured exploited stocks, Marine Ecology Progress Series, 192 (2000), 81–101. https://doi.org/10.3354/meps192089 doi: 10.3354/meps192089
    [5] J. N. Sanchirico, J. E. Wilen, Bioeconomics of spatial exploitation in a patchy environment, J. Environ. Econ. Manage., 37 (1999), 129–150. https://doi.org/10.1006/jeem.1998.1060 doi: 10.1006/jeem.1998.1060
    [6] J. N. Sanchirico, J. E. Wilen, A bioeconomic model of marine reserve creation, J. Environ. Econ. Manage., 42 (2001), 257–276. https://doi.org/10.1006/jeem.2000.1162 doi: 10.1006/jeem.2000.1162
    [7] G. Brown, J. Roughgarden, A metapopulation model with private property and a common pool, Ecol. Econ., 30 (1997), 65–71. https://doi.org/10.1016/S0921-8009(97)00564-8 doi: 10.1016/S0921-8009(97)00564-8
    [8] C. W. Clark, The Worldwide Crises in Fisheries, Cambridge University Press, New York, 2006.
    [9] C. J. Walter, S. Martell, Fisheries Ecology Management, Princeton University Press, Princeton, New Jersey, 2004.
    [10] T. Quinn Ⅱ, Ruminations on the development and future of population dynamic models in fisheries, Nat. Res. Model, 16 (2003), 341–392. https://doi.org/10.1111/j.1939-7445.2003.tb00119.x doi: 10.1111/j.1939-7445.2003.tb00119.x
    [11] J. E. Wilen, Spatial management of fisheries, Marine Resour. Econ., 19 (2004), 7–19. https://doi.org/10.1086/mre.19.1.42629416 doi: 10.1086/mre.19.1.42629416
    [12] M. Neubert, Marine reserves and optimal harvesting, Ecol. Lett., 6 (2003), 843–849. https://doi.org/10.1046/j.1461-0248.2003.00493.x doi: 10.1046/j.1461-0248.2003.00493.x
    [13] M. Neubert, G. Herrera, Triple benefits from spatial resource management, Theor. Ecol., 1 (2008), 5–12. https://doi.org/10.1007/s12080-007-0009-6 doi: 10.1007/s12080-007-0009-6
    [14] W. Ding, S. Lenhart, Optimal harvesting of a spatially explicit fishery model, Nat. Resour. Model., 22 (2009), 173–211. https://doi.org/10.1111/j.1939-7445.2008.00033.x doi: 10.1111/j.1939-7445.2008.00033.x
    [15] H. Joshi, G. Herrera, S. Lenhart, M. Neubert, Optimal dynamic harvest of a mobile renewable resource, Nat. Resour. Model., 22 (2009), 322–343. https:/doi.org/10.1111/j.1939-7445.2008.00038.x doi: https:/doi.org/10.1111/j.1939-7445.2008.00038.x
    [16] P. De Leenheer, Optimal placement of marine protected areas: A trade-off between fisheries goals and conservation efforts, IEEE Trans. Autom. Control, 59 (2014), 1583–1587. https:/doi.org/10.1109/TAC.2013.2292742 doi: https:/doi.org/10.1109/TAC.2013.2292742
    [17] M. R. Kelly, M. G. Neubert, S. Lenhart, Marine reserves and optimal dynamic harvesting when fishing damages habitat, Theor. Ecol., 12 (2019), 131–144. https://doi.org/10.1007/s12080-018-0399-7 doi: 10.1007/s12080-018-0399-7
    [18] M. Demir, S. Lenhart, Optimal sustainable fishery management of the Black Sea anchovy with food chain modeling framework, Nat. Resour. Model., 33 (2020), 1–29. https://doi.org/10.1111/nrm.12253 doi: 10.1111/nrm.12253
    [19] H. Schättler, U. Ledzewicz, Geometric Optimal Control, Springer, 2012. https://doi.org/10.1007/978-1-4614-3834-2
    [20] M. I. Zelikin, V. Borisov, Theory of Chattering Control: With Applications to Astronautics, Robotics, Economics, and Engineering, Birkhäuser, 1994. https://doi.org/10.1007/978-1-4612-2702-1
    [21] H. J. Kelley, R. Kopp, H. G. Moyer, 3 Singular extremals, Math. Sci. Eng., 31 (1967), 63–101. https://doi.org/10.1016/S0076-5392(09)60039-4 doi: 10.1016/S0076-5392(09)60039-4
    [22] R. M. Lewis, Definitions of order and junction conditions in singular optimal control problems, SIAM J. Control Optim., 18 (1980), 21–32. https://doi.org/10.1137/0318002 doi: 10.1137/0318002
    [23] W. W. Hager, Extension of switch point algorithm to boundary-value problems, Comput. Optim. Appl., 86 (2023), 1229–1246. https://doi.org/10.1007/s10589-023-00530-y doi: 10.1007/s10589-023-00530-y
    [24] M. Aghaee, W. W. Hager, The switch point algorithm, SIAM J. Control Optim., 59 (2021), 2570–2593. https://doi.org/10.1137/21M1393315 doi: 10.1137/21M1393315
    [25] M. Schaefer, Some considerations of population dynamics and economics in relation to the management of commercial marine fisheries, J. Fish. Res. Board Can., 14 (1957), 669–681. https://doi.org/10.1139/f57-025 doi: 10.1139/f57-025
    [26] C. W. Clark, Mathematical Bioeconomics: The Optimal Management of Renewable Resources, Wiley, New York, 1990.
    [27] R. A. Fisher, The wave of advance of advantageous genes, Eugenics, 7 (1937), 355–369. https://doi.org/10.1111/j.1469-1809.1937.tb02153.x doi: 10.1111/j.1469-1809.1937.tb02153.x
    [28] A. Cañada, J. L. Gámez, J. A. Montero, Study of an optimal control problem for diffusive nonlinear elliptic equations of logistic type, SIAM J. Control Optim., 36 (1998), 1171–1189. https://doi.org/10.1137/S0363012995293323 doi: 10.1137/S0363012995293323
    [29] J. A. Montero, A uniqueness result for an optimal control problem on a diffusive elliptic Volterra-Lotka type equation, J. Math. Anal. Appl., 243 (2000), 13–31. https://doi.org/10.1006/jmaa.1999.6638 doi: 10.1006/jmaa.1999.6638
    [30] H. Berestycki, P. L. Lions, Some applications of the method of super and subsolutions, in Bifurcation and Nonlinear Eigenvalue Problems. Lecture Notes in Mathematics (eds. C. Bardos, J. M. Lasry and M. Schatzman), Springer Berlin Heidelberg, (1980), 16–41. https://doi.org/10.1007/BFb0090426
    [31] N. V. Krylov, Nonlinear Elliptic and Parabolic Equations of the Second Order, D. Reidel Publishing Company, Dordrecht, Holland, (1987). https://doi.org/10.1007/978-94-010-9557-0
    [32] U. Ledzewicz, H. Schättler, Singular controls and chattering arcs in optimal control problems arising in biomedicine, Control Cybern., 38 (2009), 1501–1523.
    [33] J. P. McDanell, W. F. Powers, Necessary conditions for joining optimal singular and nonsingular subarcs, SIAM J. Control, 9 (1971), 161–173. https://doi.org/10.1137/0309014 doi: 10.1137/0309014
    [34] S. Atkins, Regularization of Singular Control Problems that Arise in Mathematical Biology, PhD thesis, University of Florida, 2021.
    [35] W. W. Hager, H. Zhang, Algorithm 1035: A gradient-based implementation of the polyhedral active set algorithm, ACM Trans. Math. Software, 49 (2023), 1–13. https://doi.org/10.1145/3583559 doi: 10.1145/3583559
    [36] W. W. Hager, H. Zhang, An active set algorithm for nonlinear optimization with polyhedral constraints, Sci. Chin. Math., 59 (2016), 1525–1542. https://doi.org/10.1007/s11425-016-0300-6 doi: 10.1007/s11425-016-0300-6
    [37] S. Atkins, M. Aghaee, M. Martcheva, W. W. Hager, Solving singular control problems in mathematical biology using PASA, in Computational and Mathematical Population Dynamics (eds. N. Tuncer, M. Martcheva, O. Prosper and L. Childs), World Scientific, (2023), 319–419. https://doi.org/10.1142/9789811263033_0009
    [38] M. Caponigro, R. Ghezzi, B. Piccoli, E. Trelat, Regularization of chattering phenomena via bounded variation controls, IEEE Trans. Autom. Control, 63 (2018), 2046–2060. https://doi.org/10.1109/TAC.2018.2810540 doi: 10.1109/TAC.2018.2810540
  • This article has been cited by:

    1. Albert W. Pilkington, Gregory C. Donohoe, Novruz G. Akhmedov, Timothy Ferrebee, Stephen J. Valentine, Justin Legleiter, Hydrogen Peroxide Modifies Aβ–Membrane Interactions with Implications for Aβ40 Aggregation, 2019, 58, 0006-2960, 2893, 10.1021/acs.biochem.9b00233
    2. Albert W. Pilkington, Jane Schupp, Morgan Nyman, Stephen J. Valentine, David M. Smith, Justin Legleiter, Acetylation of Aβ40 Alters Aggregation in the Presence and Absence of Lipid Membranes, 2020, 11, 1948-7193, 146, 10.1021/acschemneuro.9b00483
    3. Kristine L. Teppang, Qilin Zhao, Jerry Yang, Development of fluorophores for the detection of oligomeric aggregates of amyloidogenic proteins found in neurodegenerative diseases, 2023, 11, 2296-2646, 10.3389/fchem.2023.1343118
    4. Akane Iida-Adachi, Hideki Nabika, Changes in Adsorption, Aggregation, and Diffusion Nature of Amyloid β on a Lipid Membrane in an Open System, 2025, 0743-7463, 10.1021/acs.langmuir.4c03663
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1139) PDF downloads(93) Cited by(0)

Figures and Tables

Figures(3)  /  Tables(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog