Multi-scale modelling of amyloid formation from unfolded proteins using a set of theory derived rate constants

Biophysical Reviews' "Meet the Editors Series": a profile of Damien Hall

This piece introduces Damien Hall, Chief Editor of the Biophysical Reviews journal since 2019. Currently working as an Assistant Professor at Kanazawa University, the author describes his association with the journal along with some parts of his family history and academic journey.

MIL-CELL: a tool for multi-scale simulation of yeast replication and prion transmission

The single-celled baker's yeast, Saccharomyces cerevisiae, can sustain a number of amyloid-based prions, the three most prominent examples being [URE3], [PSI+], and [PIN+]. In the laboratory, haploid S. cerevisiae cells of a single mating type can acquire an amyloid prion in one of two ways (i) spontaneous nucleation of the prion within the yeast cell, and (ii) receipt via mother-to-daughter transmission during the cell division cycle. Similarly, prions can be lost due to (i) dissolution of the prion amyloid by its breakage into non-amyloid monomeric units, or (ii) preferential donation/retention of prions between the mother and daughter during cell division. Here we present a computational tool (Monitoring Induction and Loss of prions in Cells; MIL-CELL) for modelling these four general processes using a multiscale approach describing both spatial and kinetic aspects of the yeast life cycle and the amyloid-prion behavior. We describe the workings of the model, assumptions upon which it is based and some interesting simulation results pertaining to the wave-like spread of the epigenetic prion elements through the yeast population. MIL-CELL is provided as a stand-alone GUI executable program for free download with the paper. MIL-CELL is equipped with a relational database allowing all simulated properties to be searched, collated and graphed. Its ability to incorporate variation in heritable properties means MIL-CELL is also capable of simulating loss of the isogenic nature of a cell population over time. The capability to monitor both chronological and reproductive age also makes MIL-CELL potentially useful in studies of cell aging.

Amyloid fibril formation, structure and domain swapping of acyl-coenzyme A thioesterase-7

Acyl-coenzyme A thioesterase (Acot) enzymes are involved in a broad range of essential intracellular roles including cell signalling, lipid metabolism, inflammation and the opening of ion channels. Dysregulation in lipid metabolism has been linked to neuroinflammatory and neurological disorders such as Alzheimer's and Parkinson's diseases. Structurally, Acot enzymes adopt a circularised trimeric arrangement with each monomer containing an N- and a C-terminal hotdog domain. Acot7 spontaneously forms amyloid fibrils in vitro under physiological conditions. The resultant amyloid fibrillar structures were characterised by dye-binding fluorescence assays, far-UV circular dichroism spectroscopy, transmission electron microscopy and X-ray fibre diffraction. Acot7 has an unusual mechanism of aggregation with no lag phase. The initial phase (~ 18 h) of aggregation involves conformational rearrangement within the oligomers to form species of enhanced β-sheet character. The subsequent loss of α-helical structure is accompanied by large-scale amyloid fibril formation. The crystal structure of Acot7 revealed an unexpected arrangement of the two domains within the circularised trimeric structure, which is the basis for a proposed mechanism of amyloid fibril formation involving domain swapping during the initial phase of aggregation. Acot7 formed fibrils in the presence of its substrate arachidonoyl-CoA and its inhibitors and maintained its enzyme activity during fibril assembly. It is proposed that the Acot7 fibrillar form acts as functional amyloid.

A fortunate period of overlap with Prof. Haruki Nakamura

The author recounts a period of overlap with Prof. Haruki Nakamura that stretched from 2007 till the present day. Starting as a short-term research fellow in his laboratory, the author has also been a coauthor, academic colleague, and joint journal editorial board member of Prof. Nakamura.

Condition-dependent structural collapse in the intrinsically disordered N-terminal domain of prion protein

Prion protein is composed of a structure-unsolved N-terminal domain and a globular C-terminal domain. Under limited trypsin digestion, mouse recombinant prion protein can be cleaved into two parts at residue Lys105. Here, we termed these two fragments as the N-domain (sequence 23-105) and the C-domain (sequence 106-230). In this study, the structural properties of the N-domain, the C-domain, and the full-length protein were explored using small-angle X-ray scattering, analytical ultracentrifugation, circular dichroism spectroscopy, and the 8-anilino-1-naphthalenesulfonic acid binding assay. The conformation and size of the prion protein were found to change sensitively under the solvent conditions. The positive residues in the sequence 23-99 of the N-domain were found to be responsible for the enhanced flexibility with the salt concentration reduced below 5 mM. The C-domain containing a hydrophobic patch tends to unfold and aggregate during a salt-induced structural collapse. The N-domain collapsed together with the C-domain at pH 5.2, whereas it collapsed independently at pH 4.2. The positively charged cluster (sequence 100-105) in the N-domain contributed to protecting the exposed hydrophobic surface of the C-domain.

A simple method for modeling amyloid kinetics featuring position biased fiber breakage

A mathematical model of amyloid fiber formation is described that is able to simply specify different rates of fiber breakage at internal versus end regions. This Note presents the derivation of the relevant equations and provides results showing the dramatic effects of position biased fiber breakage on the kinetics of amyloid growth.

Structure-based functional analysis of BRCA1 RING domain variants: Concordance of computational mutagenesis, experimental assay, and clinical data

A significant impediment to the improvement of clinical outcomes in treating breast and ovarian cancers rests with the lack of available interpretations for BRCA1 variants of unknown significance. Two research groups recently implemented large-scale functional assays for quantifying effects of single missense mutations on homology-directed DNA repair activity of BRCA1 variants, which is critical for tumor suppression and strongly correlates with cancer risk, and their results are significantly concordant with each other as well as with known pathogenic and benign variant clinical data. In this work, we implemented an established computational mutagenesis procedure to characterize structural impacts of single residue replacements to the BRCA1 RING domain. The computational data showed similarly strong concordance with known clinical data as well as with experimental data from both functional assays. Predictions made by models trained on our computational data offer a complementary and orthogonal approach for classifying all remaining unexplored BRCA1 RING domain variants.

Mathematical Models to Shed Light on Amyloid-Beta and Tau Protein Dependent Pathologies in Alzheimer's Disease

The number of patients suffering from dementia due to Alzheimer's disease (AD) is constantly rising worldwide. This has accordingly resulted in huge burdens on the health systems and involved families. Lack of profound understanding of neural networking in normal brain and their interruption in AD makes the treatment of this neurodegenerative multifaceted disease a challenging issue. In recent years, mathematical and computational methods have paved the way towards a better understanding of the brain functional connectivity. Thus, much attention has been paid to this matter from both basic science researchers and clinicians with an interdisciplinary approach to determine what is not functioning properly in AD patients and how this malfunctioning can be addressed. In this review, a number of AD-related articles and well-studied pathophysiologic topics (e.g., amyloid-beta, neurofibrillary tangles, Ca²⁺ dysregulation, and synaptic plasticity alterations) has been literally surveyed from a computational and systems biology point of view. The neural networks were discussed from biological and mathematical point of views and their alterations in recent findings were further highlighted. Application of the graph theoretical analysis in the brain imaging was reviewed, depicting the relations between brain structure and function, without diving into mathematical details. Moreover, differential rate equations were briefly articulated, emphasizing the potential use of these equations in simplifying complex processes in relevance to pathologies of AD. Comprehensive insights were given into the AD progression from neural networks perspective, which may lead us towards potential strategies for early diagnosis and effective treatment of AD.

On the nature of the optimal form of the holdase-type chaperone stress response

The holdase paradigm of chaperone action involves preferential binding by the chaperone to the unfolded protein state, thereby preventing it from either, associating with other unstable proteins (to form large dysfunctional aggregates), or being degraded by the proteolytic machinery of the cell/organism. In this paper, we examine the necessary physical constraints imposed upon the holdase chaperone response in a cell-like environment and use these limitations to comment on the likely nature of the optimal form of chaperone response in vivo.

Unified theoretical description of the kinetics of protein aggregation

Solution conditions chosen for the production of amyloid can also promote formation of significant extents of amorphous protein aggregate. In one interpretation, the amyloid and amorphous aggregation pathways are considered to be in competition with each other. An alternative conceptualization involves considering amorphous aggregation as an obligatory intermediate process of the amyloid formation pathway. Here, we review recently developed macroscopic-level theories of protein aggregation that unify these two competing models into a single paradigm. Key features of the unified model included (1) a description of the amorphous aggregate as a second liquid phase with the degree of liquid-like character determined by the mobility of the monomer within it, and (2) heterogeneous growth pathways based on nucleation, growth, and fragmentation of amyloid occurring within different phases and at their interfacial boundary. Limiting-case behaviors of the protein aggregation reaction, either singly involving amyloid or amorphous aggregate production, and mixed-case behaviors, involving competitive and/or facilitated growth of amorphous and amyloid species, are presented and reviewed in context. This review principally describes an approach developed by Hirota and Hall 2019 (Hirota, N. and Hall, D. 2019. Protein Aggregation Kinetics: A Unified Theoretical Description. Chapter 7 of 'Protein Solubility and Amorphous Aggregation: From Academic Research to Applications in Drug Discovery and Bioindustry' edited by Y. Kuroda and F. Arisaka. CMC Publishers). Sections of that work are translated from the original Japanese and republished here with the full permission of CMC Publishing Corporation.

Ionic liquids and protein folding-old tricks for new solvents

One important aspect of the green chemistry revolution has been the use of ionic liquids as the solvent in liquid-phase enzymatic catalysis. An essential requirement for protein enzyme function is the correct folding of the polypeptide chain into its functional "native" state. Quantitative assessment of protein structure may be carried out either empirically, or by using model-based characterization procedures, in which the parameters are defined in terms of a standard reference state. In this short note, we briefly outline the nature of the parameters associated with different empirical and model-based characterization procedures and point out factors which affect their interpretation when using a base solvent different from water. This review principally describes arguments developed by Wakayama et al., Protein Solubility and Amorphous Aggregation: From Academic Research to Applications in Drug Discovery and Bioindustry, 2019, edited by Y. Kuroda and F. Arisaka; CMC Publishing House. Sections of that work are translated from the original Japanese and republished here with the full permission of CMC Publishing Corporation.

The functional roles of the unstructured N- and C-terminal regions in αB-crystallin and other mammalian small heat-shock proteins

Small heat-shock proteins (sHsps), such as αB-crystallin, are one of the major classes of molecular chaperone proteins. In vivo, under conditions of cellular stress, sHsps are the principal defence proteins that prevent large-scale protein aggregation. Progress in determining the structure of sHsps has been significant recently, particularly in relation to the conserved, central and β-sheet structured α-crystallin domain (ACD). However, an understanding of the structure and functional roles of the N- and C-terminal flanking regions has proved elusive mainly because of their unstructured and dynamic nature. In this paper, we propose functional roles for both flanking regions, based around three properties: (i) they act in a localised crowding manner to regulate interactions with target proteins during chaperone action, (ii) they protect the ACD from deleterious amyloid fibril formation and (iii) the flexibility of these regions, particularly at the extreme C-terminus in mammalian sHsps, provides solubility for sHsps under chaperone and non-chaperone conditions. In the eye lens, these properties are highly relevant as the crystallin proteins, in particular the two sHsps αA- and αB-crystallin, are present at very high concentrations.

Measurement of amyloid formation by turbidity assay-seeing through the cloud

Detection of amyloid growth is commonly carried out by measurement of solution turbidity, a low-cost assay procedure based on the intrinsic light scattering properties of the protein aggregate. Here, we review the biophysical chemistry associated with the turbidimetric assay methodology, exploring the reviewed literature using a series of pedagogical kinetic simulations. In turn, these simulations are used to interrogate the literature concerned with in vitro drug screening and the assessment of amyloid aggregation mechanisms.

Fibril formation from the amyloid-β peptide is governed by a dynamic equilibrium involving association and dissociation of the monomer

Here I review the molecular mechanisms by which water-soluble monomeric amyloid-β (Aβ) peptides are transformed into well-organized supramolecular complexes called amyloid fibrils. The mechanism of amyloid formation is considered theoretically on the basis of experimental results, and the structural and mechanistic similarities of amyloid fibrils to three-dimensional crystals are highlighted. A number of important results from the literature are described. These include the observation that a correct ratio of monomer association and dissociation rate constants is key for formation of well-organized amyloid fibrils. The dynamic nature of the amyloid-β structure is discussed, along with the possibly obligate requirement of the transient formation of a hairpin-like fold prior to its incorporation into amyloid fibrils. Many rounds of monomer association and dissociation events may be present during an apparently silent lag-period. Amongst these association/dissociation events, interaction between the C-terminal regions of the Aβ peptide seems to be more favored. Such association and dissociation events occurring in a "trial-and-error" fashion may be an important requirement for the formation of well-organized amyloid fibrils.

A multi-pathway perspective on protein aggregation: implications for control of the rate and extent of amyloid formation

The nucleation-growth model has been used extensively for characterizing in vitro amyloid fibril formation kinetics and for simulating the relationship between amyloid and disease. In the majority of studies amyloid has been considered as the dominant, or sole, aggregation end product, with the presence of other competing non-amyloid aggregation processes, for example amorphous aggregate formation, being largely ignored. Here, we examine possible regulatory effects that off-pathway processes might exert on the rate and extent of amyloid formation - in particular their potential for providing false positives and negatives in the evaluation of anti-amyloidogenic agents. Furthermore, we investigate how such competing reactions might influence the standard interpretation of amyloid aggregation as a two-state system. We conclude by discussing our findings in terms of the general concepts of supersaturation and system metastability - providing some mechanistic insight as to how these empirical phenomena may manifest themselves in the amyloid arena.

Computational modeling of the relationship between amyloid and disease

Amyloid is a title conferred upon a special type of linear protein aggregate that exhibits a common set of structural features and dye binding capabilities. The formation of amyloid is associated with over twenty-seven distinct human diseases which are collectively referred to as the amyloidoses. Although there is great diversity amongst the amyloidoses with regard to the polypeptide monomeric precursor, targeted tissues and the nature and time course of disease development, the common underlying link of a structurally similar amyloid aggregate has prompted the search for a unified theory of disease progression in which amyloid production is the central element. Computational modeling has allowed the formulation and testing of scientific hypotheses for exploring this relationship. However, the majority of computational studies on amyloid aggregation are pitched at the atomistic level of description, in simple ideal solution environments, with simulation time scales of the order of microseconds and system sizes limited to a hundred monomers (or less). The experimental reality is that disease related amyloid aggregation processes occur in extremely complex reaction environments (i.e. the human body), over time-scales of months to years with monitoring of the reaction achieved using extremely coarse or indirect experimental markers that yield little or no atomistic insight. Clearly a substantial gap exists between computational and experimental communities with a deficit of 'useful' computational methodology that can be directly related to available markers of disease progression. This Review will place its focus on the development of these latter types of computational models and discuss them in relation to disease onset and progression.

A model of amyloid's role in disease based on fibril fracture

Although the correlative evidence relating the presence of amyloid fibrils and certain disease states (e.g. Alzheimer's disease and Type 2 Diabetes) is overwhelming, a direct causative role for amyloid has proved harder to establish. Current thinking links a narrow region of the aggregate protein size distribution, the so called 'early aggregate' domain to cellular toxicity. A troubling feature of this theory however is that the nucleated reaction mechanism by which amyloid formation is believed to occur results in a very low number concentration of early aggregates which are rapidly extended to form amyloid fibrils. This situation means that the concentration of early aggregates is very low at the time when they are supposedly at their most toxic. Here we adopt a novel explicit simulation strategy to examine a kinetic regime involving nucleated growth combined with fibril fragmentation under which this situation can be reversed so as to produce a high number concentration of small on-pathway toxic aggregates. Dependent upon the rate of fragmentation, the time scale for generation of toxic early aggregates may be coupled, uncoupled or disassociated from the time scale for the appearance of amyloid fibrils. Furthermore the model presents itself as a biochemical 'switch' transitioning between modes of amyloid induced cell death dependent upon either specific amyloid toxicity or non-specific solid mass induced tissue damage.