Skip to main content

In silico prediction of novel residues involved in amyloid primary nucleation of human I56T and D67H lysozyme



Amyloidogenic proteins are most often associated with neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, and Huntington’s disease, but there are more than two dozen human proteins known to form amyloid fibrils associated with disease. Lysozyme is an antimicrobial protein that is used as a general model to study amyloid fibril formation. Studies aimed at elucidating the process of amyloid formation of lysozyme tend to focus on partial unfolding of the native state due to the relative instability of mutant amyloidogenic variants. While this is well supported, the data presented here suggest the native structure of the variants may also play a role in primary nucleation.


Three-dimensional structural analysis identified lysozyme residues 21, 62, 104, and 122 as displaced in both amyloidogenic variants compared to wild type lysozyme. Residue interaction network (RIN) analysis found greater clustering of residues 112–117 in amyloidogenic variants of lysozyme compared to wild type. An analysis of the most energetically favored predicted dimers and trimers provided further evidence for a role for residues 21, 62, 104, 122, and 112–117 in amyloid formation.


This study used lysozyme as a model to demonstrate the utility of combining 3D structural analysis with RIN analysis for studying the general process of amyloidogenesis. Results indicated that binding of two or more amyloidogenic lysozyme mutants may be involved in amyloid nucleation by placing key residues (21, 62, 104, 122, and 112–117) in proximity before partial unfolding occurs. Identifying residues in the native state that may be involved in amyloid formation could provide novel drug targets to prevent a range of amyloidoses.


Amyloidoses are a group of diseases defined by the formation of protein aggregates characterized by stacks of cross-beta sheets [1]. There are dozens of different amyloid diseases caused by a variety of both wild type (WT) and mutant forms of proteins [2]. Some of the most well-known amyloidoses are neurodegenerative diseases such as Alzheimer’s disease (involving amyloid-beta peptide) and Parkinson’s disease (involving alpha-synuclein protein). However, not all amyloid diseases affect the brain. Lysozyme amyloidosis is a rare disease characterized by the deposition of amyloid fibrils of the enzyme lysozyme. Lysozyme was discovered by Alexander Fleming in 1922 [3] and is an antimicrobial enzyme synthesized by hepatocytes, cells of the gastrointestinal system, and macrophages [4]. Lysozyme amyloidosis has no known effective treatment and leads to lysozyme amyloid deposits typically concentrated in the liver [5], spleen, gastrointestinal tract [6], and kidneys [7]. Lysozyme amyloidosis is thought to be largely caused by subtle structural changes of the protein caused by genetic mutations that lead to pockets of local instability and a greater likelihood of partial unfolding [8]. The Online Mendelian Inheritance in Man (OMIM) database [9] entry for lysozyme (OMIM ID 153450) reports four lysozyme variants that are associated with the disease: I56T [10], D67H [10], W64R [11], and F57I [12].

Lysozyme has long been used as a model for studying protein structure and folding. Since lysozyme is structurally and functionally well-characterized, the protein provides a useful model for understanding the complex process of amyloid fibril formation [13]. Several studies have investigated the role of amyloidogenic mutations on lysozyme amyloid formation with a focus on the first identified mutations, I56T and D67H. Studies that examined the crystal structure of the WT and I56T variant suggest very little difference in the native structure of these enzymes [14]. The D67H variant, however, destroys the hydrogen bonds that stabilize the beta-domain, leading to the displacement of a long loop of residues [15]. Because the obvious loop displacement between the WT and D67H mutant is not present in the I56T mutant, it is thought that this change is not responsible for amyloidogenesis. Instability of the I56T variant may be caused by subtle changes in bonding between alpha and beta domains of lysozyme; similar bonding changes are also evident in the D67H variant [15].

Since the structure of proteins in amyloid plaques are different from the native structure, amyloidogenic proteins must at least partially unfold during amyloidogenesis. Most studies focus on the unfolding process of lysozyme instead of differences in the native structure. The amyloidogenic proteins likely spend more time partially unfolded, providing more opportunities for unfolded segments to interact and aggregate in the form of amyloid plaques [8]. The importance of the partially unfolded state for lysozyme amyloidosis has been demonstrated in vitro with the use of antibodies that stabilize the protein [16, 17]. Studies have also shown that both I56T and D67H are less stable than WT lysozyme when heated [15], and I56T is also less stable than WT at low pH [14, 15, 18, 19], further supporting a role for instability. However, other factors besides regions of protein instability may be involved in amyloidogenesis. Further studies that examine primary nucleation from different perspectives could provide more insight into this important process that is associated with a variety of diseases.

Residue interaction networks (RIN) abstract protein structure into a network of likely side-chain interactions with residues represented as nodes and interactions represented as edges, the connections between the nodes [20]. Several metrics are available for studying networks and identifying subnetworks of interest [21,22,23]. Many network features have been associated with and applied to protein structural and functional characteristics [22, 24,25,26], demonstrating the relevance of RINs to structural biology. Clusters are particularly interesting in RIN analysis because they identify areas with many chemical interactions, suggesting structural rigidity or functional importance [27]. RIN analysis is most useful when combined with 3D structural analysis [28]. This study uses the two-pronged approach of combining 3D structural analysis with RIN analysis to identify residues in the native structure that are likely involved in amyloid formation.


Three-dimensional structure visualization and structure comparison

The Online Mendelian Inheritance in Man (OMIM) database [9] was searched for mutations in lysozyme that have been associated with amyloidosis. The Protein Data Bank (PDB) [29] was then searched for human lysozyme structures with these mutations, resulting in a dataset of WT human lysozyme (PDB ID: 1REX, [30]) and two amyloidogenic variants, I56T (PDB ID: 1LOZ, [15]) and D67H (PDB ID: 1LYY, [15]). The 3D protein structure coordinates were downloaded from the PDB and visualized using UCSF Chimera v1.11.2 software [31]. The 3D structures were overlapped using the MatchMaker application [32] in UCSF Chimera with default settings, and the root mean square deviations (RMSDs) for the full residues from the wild type (1REX) structure of lysozyme were calculated in the Multialign Viewer [32]. Side chains that had different locations when compared to WT lysozyme in both amyloidogenic variants were selected for further study. Because the resolution of the PDB files used for the comparison was less than or equal to 1.8 Å, only residues with a RMSD from WT greater than or equal to 1.9 Å were considered.

Generating residue interaction networks and calculating clusters and metrics

To detect network clusters in the proteins, the PDB files were converted to GML format using the Protein Graph Converter software from the Protein Graph Repository (PGR) [33]. Each alpha carbon was considered a node, and an edge was drawn between every alpha carbon within seven angstroms of another. The GML files were then analyzed for clusters using the MCODE application in the network analysis software Cytoscape v3.4.0 [34]. Only clusters with MCODE scores greater than or equal to 5.00 were selected for further analysis. The residues involved in the clusters in the amyloidogenic mutants were compared to those identified in WT lysozyme. As with the 3D structure comparison, cluster changes that are in common between the amyloidogenic mutants and different from WT lysozyme were selected for further analysis. UCSF Chimera was used to calculate the average residue B-factor for each of the clusters. The Pearson correlation of MCODE scores for the selected clusters of each of the lysozyme PDB files with the average B-factor for the clusters was calculated using GraphPad Prism 7, and a p-value of < 0.05 was considered significant.

Generating predicted dimer and trimer structures and calculating interprotein bond number and energies

ClusPro v2.0 software [35] was used to generate predicted structures of homodimers, homotrimers, and heterodimers for WT lysozyme and the amyloidogenic variants. After multimer generation, the resulting PDB files for the top predicted dimer and trimer structures were edited so that each lysozyme protein was given a unique name. Next, the edited PDB files were uploaded to the Residue Interaction Network Generator (RING) v2.0 software [36], and residue interaction networks were created using a strict distance threshold between the closest atoms of residues separated by at least two other residues. Multiple edges per residue pair were allowed but only one edge per interaction type. The resulting graph files were analyzed for the number of interactions and overall bonding energy occurring between lysozyme proteins in dimers or trimers. The number and bond energies of the interprotein interactions were analyzed for the whole complexes and the residues of interest to provide information about the relative importance of the residues of interest to the formation of dimers and trimers.


Three-dimensional structural comparison of lysozyme

Three-dimensional structure overlaps from the MatchMaker software revealed residues of the amyloidogenic variants that differed from WT lysozyme. The D67H variant diverges more from WT lysozyme than the I56T variant. However, there are only four residues with a RMSD greater than or equal to 1.9 Å that were shared by both the I56T and the D67H amyloidogenic lysozyme variants: residues 21, 62, 104 and 122. In all four cases, the structural changes are in proximity to each other in the 3D structure (Fig. 1).

Fig. 1

Three-dimensional overlap of lysozyme structures 1REX (WT, tan), 1LOZ (I56T, orange), and 1LYY (D67H, blue). The side chains of residues 21, 62, 104, 112–117, and 122 are shown and outlined in green

Residue interaction network clustering analysis of lysozyme structures

The 3D structure of lysozyme and its resulting PGR residue interaction network (RIN) representation are shown in Fig. 2. The MCODE application in Cytoscape revealed four clusters in each of the lysozyme structures that had an MCODE score of greater than 5.00. Some of the clusters in each of the lysozyme variants involved similar or identical sets of residues. As before, we focused on the differences from the WT clusters that were present in both amyloidogenic proteins. The most robust and consistent difference was the cluster around residues 112–117. The MCODE score for the cluster containing these residues in WT lysozyme was 5.11, but it increased in both amyloidogenic variants to 6.00. Residues 104 and 106–108 are clustered with 112–117 in WT lysozyme but not in I56T and D67H. The average B-factors of the clusters containing residues 112–117 also decreased in both amyloidogenic variants compared to WT lysozyme. These results are shown in Table 1. Residues 112–117 are shown in Fig. 1. There was a statistically significant negative correlation (r2 = 0.44, p = 0.0184) between MCODE scores for the top four clusters in each of the three lysozyme network structures and the average B-factor for each cluster. The residues are in proximity to the other residues of interest as shown in Fig. 1.

Fig. 2

a Three-dimensional WT lysozyme, PDB 1REX. b Network representation of WT lysozyme with amino acids represented as nodes and edges drawn between alpha carbons within seven angstroms of each other. Residue numbering only for this image starts at zero

Table 1 RIN clusters in WT and amyloidogenic variants of lysozyme

Interprotein bonds involving residues 21, 62, 104, 112–117, and 122 in predicted lysozyme dimers and trimers

The top-rated ClusPro models of dimer and trimer structures for each of the variants of lysozyme are shown in Fig. 3 (homodimers), Fig. 4 (homotrimers), and Fig. 5 (heterodimers). The number of interprotein residue interactions and the strength of the bonding energy for these interactions for all residues and for the residues of interest were quantified using RING 2.0 software, and the results are shown in Table 2 (homodimers), Table 3 (homotrimers), and Table 4 (heterodimers).

Fig. 3

Predicted 3D structures of WT or mutant lysozyme homodimers. a WT: WT, (b) I56T: I56T, and (c) D67H: D67H. The side chains of residues 21, 62, 104, 112–117, and 122 are shown and outlined in green

Fig. 4

Predicted 3D structures WT or mutant lysozyme homotrimers. a WT: WT: WT, (b) I56T: I56T: I56T, and (c) D67H: D67H: D67H. The side chains of residues 21, 62, 104, 112–117, and 122 are shown and outlined in green

Fig. 5

Predicted 3D structures of heterodimers of each of the amyloidogenic lysozyme variants with WT lysozyme. a WT: I56T, (b) WT: D67H. The side chains of residues 21, 62, 104, 112–117, and 122 are shown and outlined in green

Table 2 ClusPro predicted interprotein binding energies for WT and mutant lysozyme homodimers
Table 3 ClusPro predicted interprotein binding energies for WT and mutant lysozyme homotrimers
Table 4 ClusPro predicted interprotein binding energies for WT and mutant lysozyme heterodimers

For homodimers, residues of interest made up 36% of the number of interprotein residue interactions in WT lysozyme, contributing to 39.1% of the predicted interprotein bonding energy. The overall number of interprotein residue interactions increased to 49.3% of the total in I56T, but it decreased to 29.5% of the total number in D67H. However, the total percentage of interprotein bonding energy contributed by the residues of interest increased for both amyloidogenic variants even though the total interprotein bonding energy for dimers was less than WT dimers.

Trimers showed different trends from the homodimers (Table 3). The total number of interprotein residue interactions increased for I56T (133 interactions) and D67H (108 interactions) compared to WT (89 interactions). The contribution of the residues of interest to the total number of interprotein residue interactions increased from 2.2% for WT to 55.6% for I56T and 22.2% for D67H. Interprotein bonding energy showed similar trends. The total interprotein bonding energy increased from 824.8 kJ/mol in the WT to 1392.8 kJ/mol in I56T and 1060.8 kJ/mol in D67H. Residues of interest contributed 1.5% of the total interprotein bonding energy in the WT but 58.7% in I56T and 30.2% in D67H.

Heterodimers (WT: I56T and WT: D67H) showed fewer and less energetic interprotein interactions than homodimers (Table 4). The WT: I56T heterodimer had 48 interprotein interactions and 481 kJ/mol interprotein bonding energy, and the WT: D67H heterodimer had 49 interprotein interactions and 482.6 kJ/mol bonding energy. In both cases, there was a high reliance on the residues of interest for the bonding. The residues of interest made up more than 50% of the number of interprotein bonding interactions and nearly 50% of the interprotein bonding energies of both heterodimers.


Three-dimensional structural comparison suggests residues 21, 62, 104, and 122 may be involved in lysozyme amyloidogenesis

Comparing 3D protein structures revealed that both amyloidogenic variants of lysozyme analyzed differ from WT lysozyme in the location of residues 21, 62, 104, and 122. Following the reasoning of Booth and colleagues [15], because these differences are common to both amyloidogenic variants, they suggest these residues may play a role in the formation of amyloid fibrils. The four residues (21, 62, 104, and 122) are in proximity to each other in all three variants (Fig. 1). As discussed above, the prevailing hypothesis for lysozyme amyloidosis is that the mutations disrupt the hydrogen bonds near the residues between alpha and beta domains, leading to partial unfolding followed by fibril formation [8, 19, 37]; interprotein interactions between the native structures are not thought to play a large role. However, the consistency of the native structure changes observed in both amyloidogenic variants hints at a role for the residues in the native structure in amyloidosis. We hypothesized that these residues may facilitate an interprotein interaction between native state amyloidogenic lysozyme proteins, contributing to the first steps of amyloidosis. To gather further support for this hypothesis, we examined the structures for network cluster changes.

Residues 112–117 may also be involved in lysozyme amyloidogenesis

In addition to the 3D structural changes described above, amyloidogenic variants of lysozyme were associated with changes in network clusters. The network cluster consisting of residues 112–117 of I56T and D67H had the greatest and most consistent cluster changes, showing an increased MCODE score in both variants compared to WT (Table 1). Residues 104 and 106–108 were included in the cluster containing residues 112–117 in WT lysozyme only, so the loss of this part of the cluster may also have structural implications. More clustering has been shown to be associated with greater structural stability [27, 38], so we hypothesized that residues 112–117 have greater structural stability in the amyloidogenic variants compared to WT lysozyme. To test this hypothesis, we used UCSF Chimera to calculate the average B-factor for each of the residues in each of the PDB files. The B-factor is a measure of flexibility where a lower B-factor indicates greater stability [39]. Even in our relatively small data set of three lysozyme structures, we found a statistically significant negative correlation (r2 = 0.44, p = 0.0184) between MCODE scores and average B-factors for the top four clusters of each of the PDB files shown in Table 1. Consistent with the hypothesis of greater cluster stability, residues 112–117 had smaller average B-factors in amyloidogenic variants compared to the cluster containing these residues in WT lysozyme. Because most studies focus on the instability caused by amyloidogenic mutations, disrupting this cluster while stabilizing other regions may provide a novel therapeutic approach. The side chains of residues 112–117 are shown in Fig. 1 along with the other residues of interest (21, 62, 104, and 122) from 3D structural comparison. Residues 112–117 are in proximity to the residues identified through 3D structural comparison. Therefore, residues 112–117 may also be involved in facilitating interactions between different lysozyme molecules and possibly contribute to primary nucleation of amyloid fibrils. To test this hypothesis, we simulated intermolecular interactions between native state structures.

Predicted dimer and trimer structures provide further evidence for the involvement of residues 21, 62, 104, 112–117, and 122 in lysozyme amyloidogenesis

To test the hypothesis that residues 21, 62, 104, 112–117, and 122 in the structure of amyloidogenic lysozyme variants are involved in primary nucleation, we used ClusPro docking software to predict the 3D structure of dimers and trimers of lysozyme for the PDB files. While ClusPro generates the structures of many predicted dimers and trimers, we only analyzed the top-ranked structures. Visual inspection of the predicted homodimers (Fig. 3) suggested some of the residues of interest may be involved in dimer interactions. To test this hypothesis, we examined the number and strength of interprotein residue interactions (Table 2). The number of interprotein bonds was less for amyloidogenic variant homodimers, and there was no consistent trend with the percent contribution of the residues of interest to the number of interprotein bonds in the dimers. Both amyloidogenic variant homodimers had overall less bonding energy. However, the residues of interest contributed to a greater degree to the interprotein bonding energy compared to WT lysozyme. Overall, the contribution of the residues of interest to lysozyme homodimer formation in amyloidogenic variants was not as convincing as the evidence for their role in simulated trimer formation.

Visual inspection of the predicted trimer structures of lysozyme suggested a greater role for the residues of interest in interprotein interactions of amyloidogenic variants than WT lysozyme (Fig. 4). When examined quantitatively, both amyloidogenic variants had stronger and more numerous bonds between proteins compared to WT (Table 3). This suggests the residues of interest may be largely facilitating the predicted sharp increase in interprotein interactions, supporting the important role predicted for these residues by 3D structural analysis and network analysis.

A role for mutant native structure in enhancing lysozyme amyloid fibril formation?

Taken together, these data support the hypothesis that residues 21, 62, 104, 112–117, and 122 in the mutant native states are involved in lysozyme amyloid primary nucleation. While the instability of amyloidogenic variants of lysozyme is almost certainly the most important factor for fibril formation [40], it may not be the only factor involved. The greater interprotein interactions predicted to occur between trimeric amyloidogenic mutants compared to WT lysozyme may lead to amyloidogenic variants oligomerizing more readily before unfolding. We hypothesize these aggregates are composed of mostly mutant proteins because heterodimers with WT lysozyme have less interprotein bonding energies than mutant homodimers. This suggests when one of the lysozyme variants partially unfolds, it may already be in proximity to or bound to another mutant molecule, leading to a greater probability of amyloid nucleation. The sequence of events of the unfolding process for WT lysozyme and the I56T variant are consistent with our hypothesis. It has been demonstrated that lysozyme alpha helices A, B, and D [41] are some of the last regions to unfold [42]. All the residues identified in this report except residue 62 are in or near regions of the protein that unfold later in the process, so these positions have a greater chance of maintaining their structure to facilitate intermolecular interactions during this process.

Furthermore, Ahn and colleagues concluded that mutations in the alpha domain of lysozyme are less likely to influence the formation of transient intermediate states compared to mutations in the beta domain [43]. The resistance of changes in the alpha domain to form transient intermediates is largely consistent with our findings which suggest subtle changes in the position and interactions of residues 21, 62, 104, 112–117, and 122 (largely in the alpha domain) play a role in the early steps of primary nucleation before significant unfolding of the alpha domain occurs. Sequence-based predictions of “hot spots” of aggregation using Aggrescan [44] and TANGO [45] suggest sequences in both the alpha and beta domains may be important for amyloid formation. While we also report a role for residues in the alpha domain, the residues identified here are not all included in the results of the sequence-based predictions. This may be because we predict residues 21, 62, 104, 112–117 and 122 to play a role in amyloidogenesis prior to the formation of unfolded intermediates by positioning lysozyme mutants in or near the native state close to each other. Because this is predicted to occur before significant unfolding and fibrillization of lysozyme, it is understandable that sequence-based algorithms did not identify the same residues. Fibrillization is a multistep process, and it is likely that different residues may be involved in different steps.

Residues 21, 62, 104, and 122 do not appear to be involved in increasing the flexibility of lysozyme in the amyloidogenic mutants because the total predicted bond energies from RING analysis for monomers do not suggest consistent and structurally important differences when compared to WT. Furthermore, the average B-factors for these residues in mutant lysozyme do not consistently differ from WT. It has been suggested that residues not present in the partially unfolded region can be altered without affecting the process of amyloid fibril formation [46], but our data challenge this suggestion. It has been noted that lysozyme amyloid plaques are nearly entirely composed of the amyloidogenic variant free from WT protein [18]. This may be due not only to the greater instability of amyloidogenic proteins, but also due to the predicted favored intermolecular interactions of homodimers and homotrimers compared to heterodimers (Table 4). A better understanding of the process of amyloidogenesis for lysozyme could yield insights into treatments for many different types of amyloidoses.

Study limitations and future studies

Data from various computational approaches used in this study support a role for lysozyme residues 21, 62, 104, 112–117, and 122 in lysozyme amyloidosis. However, this study has several limitations. The most obvious limitation is the small sample size of 3D structures used. Unfortunately, the study is limited by the availability of PDB files of human amyloidogenic lysozyme variants. We draw our conclusions from three PDB files, making this a preliminary study. The findings reported here may not be shared by all amyloidogenic variants of lysozyme. To increase confidence in our conclusions, further studies should explore the structures of WT and amyloidogenic lysozyme under various experimental conditions. In addition, it was found that different software programs predicted different lowest-energy oligomeric structures that influence results. Furthermore, these hypotheses need to be experimentally tested to verify the importance of these residues for amyloidogenesis. Lysozyme mutants with smaller side chains or nonpolar side chains at the residues of interest could be created and the fibril formation kinetics studied. Data from studies of double mutants may also be useful. Future studies could be performed where residues 21, 62, 104, 112–117, or 122 are mutated in combination with I56T or D67H and tested for altered fibril formation kinetics. However, it may be most useful to mutate residue 21 and a residue in 112–117 or 122 because residues 32 through 108 have been shown to form the core of the lysozyme fibril [47]. Therefore, mutating these residues may interfere with amyloid formation after unfolding.

Conclusions and testable hypotheses generated

This preliminary study used a combination of 3D structural and residue interaction network analyses to support roles for residues 21, 62, 104, 112–117, and 122 in lysozyme amyloidosis. By comparing two amyloidogenic variants to WT lysozyme, we were able to identify network and 3D structural changes that were shared between the amyloidogenic variants. Modeling dimer and trimer interactions further supported a role for these residues. These residues appear to be especially important for trimer formation. This study generates several hypotheses that can be experimentally tested. 1) Residues 112–117 are less flexible in amyloidogenic variants of lysozyme than in the WT. 2) Residues 21, 62, 104, 112–117, and 122 are involved in the primary nucleation of lysozyme by facilitating intermolecular interactions between mutant lysozyme molecules. 3) Trimers of lysozyme are more stable than dimers. 4) Lysozyme mutant molecules favor self-interactions over interactions with WT molecules. 5) Interprotein interactions in or close to the native state likely play a larger role in amyloid formation in general than previously hypothesized. This study demonstrates the utility of combining 3D structural and network analysis for understanding amyloid formation. Furthermore, it provides insight into lysozyme amyloid formation that may be applicable to the study of many other amyloidoses.



Online mendelian inheritance in man


Protein data bank


Protein graph repository


Residue interaction network


Residue interaction network generator


Root mean square deviation


Wild type


  1. 1.

    Greenwald J, Riek R. Biology of amyloid: structure, function, and regulation. Structure. 2010;18:1244–60.

    Article  PubMed  CAS  Google Scholar 

  2. 2.

    Knowles TP, Vendruscolo M, Dobson CM. The amyloid state and its association with protein misfolding diseases. Nat Rev Mol Cell Biol. 2014;15:384–96.

    Article  PubMed  CAS  Google Scholar 

  3. 3.

    Fleming A. On a remarkable bacteriolytic element found in tissues and secretions. Proc R Soc B Biol Sci The Royal Society. 1922;93:306–17.

    Article  CAS  Google Scholar 

  4. 4.

    Granel B, Valleix S, Serratrice J, Chérin P, Texeira A, Disdier P, et al. Lysozyme amyloidosis: report of 4 cases and a review of the literature. Medicine (Baltimore). 2006;85:66–73.

    Article  Google Scholar 

  5. 5.

    Harrison RF, Hawkins PN, Roche WR, MacMahon RF, Hubscher SG, Buckels JA. “Fragile” liver and massive hepatic haemorrhage due to hereditary amyloidosis. Gut. 1996;38:151–2.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. 6.

    Granel B, Serratrice J, Valleix S, Grateau G, Droz D, Lafon J, et al. A family with gastrointestinal amyloidosis associated with variant lysozyme. Gastroenterology. 2002;123:1346–9.

    Article  PubMed  CAS  Google Scholar 

  7. 7.

    Gillmore JD, Booth DR, Madhoo S, Pepys MB, Hawkins PN. Hereditary renal amyloidosis associated with variant lysozyme in a large English family. Nephrol Dial Transplant. 1999;14:2639–44.

    Article  PubMed  CAS  Google Scholar 

  8. 8.

    Canet D, Sunde M, Last AM, Miranker A, Spencer A, Robinson CV, et al. Mechanistic studies of the folding of human lysozyme and the origin of amyloidogenic behavior in its disease-related variants. Biochemistry. 1999;38:6419–27.

    Article  PubMed  CAS  Google Scholar 

  9. 9.

    Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. Online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015;43:D789–98.

  10. 10.

    Pepys MB, Hawkins PN, Booth DR, Vigushin DM, Tennent GA, Soutar AK, et al. Human lysozyme gene mutations cause hereditary systemic amyloidosis. Nature. 1993;362:553–7.

    Article  PubMed  CAS  Google Scholar 

  11. 11.

    Valleix S, Drunat S, Philit J-B, Adoue D, Piette J-C, Droz D, et al. Hereditary renal amyloidosis caused by a new variant lysozyme W64R in a French family. Kidney Int. 2002;61:907–12.

    Article  PubMed  CAS  Google Scholar 

  12. 12.

    Yazaki M, Farrell SA, Benson MD. A novel lysozyme mutation Phe57Ile associated with hereditary renal amyloidosis. Kidney Int. 2003;63:1652–7.

    Article  PubMed  CAS  Google Scholar 

  13. 13.

    Merlini G, Bellotti V. Lysozyme: a paradigmatic molecule for the investigation of protein structure, function and misfolding. Clin Chim Acta. 2005;357:168–72.

    Article  PubMed  CAS  Google Scholar 

  14. 14.

    Funahashi J, Takano K, Ogasahara K, Yamagata Y, Yutani K. The structure, stability, and folding process of amyloidogenic mutant human lysozyme. J Biochem. 1996;120:1216–23.

    Article  PubMed  CAS  Google Scholar 

  15. 15.

    Booth DR, Sunde M, Bellotti V, Robinson CV, Hutchinson WL, Fraser PE, et al. Instability, unfolding and aggregation of human lysozyme variants underlying amyloid fibrillogenesis. Nature. 1997;385:787–93.

    Article  PubMed  CAS  Google Scholar 

  16. 16.

    Dumoulin M, Last AM, Desmyter A, Decanniere K, Canet D, Larsson G, et al. A camelid antibody fragment inhibits the formation of amyloid fibrils by human lysozyme. Nature. 2003;424:783–8.

    Article  PubMed  CAS  Google Scholar 

  17. 17.

    Chan P-H, Pardon E, Menzer L, De Genst E, Kumita JR, Christodoulou J, et al. Engineering a camelid antibody fragment that binds to the active site of human lysozyme and inhibits its conversion into amyloid fibrils. Biochemistry. 2008;47:11041–54.

    Article  PubMed  CAS  Google Scholar 

  18. 18.

    Morozova-Roche LA, Zurdo J, Spencer A, Noppe W, Receveur V, Archer DB, et al. Amyloid fibril formation and seeding by wild-type human lysozyme and its disease-related mutational variants. J Struct Biol. 2000;130:339–51.

    Article  PubMed  CAS  Google Scholar 

  19. 19.

    Buell AK, Dhulesia A, Mossuto MF, Cremades N, Kumita JR, Dumoulin M, et al. Population of nonnative states of lysozyme variants drives amyloid fibril formation. J Am Chem Soc. 2011;133:7737–43.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. 20.

    Di Paola L, De Ruvo M, Paci P, Santoni D, Giuliani A. Protein contact networks: an emerging paradigm in chemistry. Chem Rev. 2013;113:1598–613.

    Article  PubMed  CAS  Google Scholar 

  21. 21.

    Junker B, Schreiber F, editors. Analysis of Biological Networks. Hoboken: Wiley; 2008.

  22. 22.

    Scardoni G, Laudanna C. Centralities Based Analysis of Complex Networks. In: Zhang Y, editor. New Frontiers in Graph Theory. London: InTech; 2012. p. 323–48.

  23. 23.

    Estrada E. The Structure of Complex Networks: Theory and Applications. Oxford: Oxford University Press; 2012.

  24. 24.

    Liu R, Hu J. Computational prediction of heme-binding residues by exploiting residue interaction network. Uversky VN, editor. PLoS One 2011;6:e25560.

  25. 25.

    Emerson IA, Louis PT. Detection of active site residues in bovine rhodopsin using network analysis. Trends Bioinforma. 2015;8:63–74.

    Article  Google Scholar 

  26. 26.

    Hu G, Yan W, Zhou J, Shen B. Residue interaction network analysis of Dronpa and a DNA clamp. J Theor Biol. 2014;348:55–64.

    Article  PubMed  CAS  Google Scholar 

  27. 27.

    Vishveshwara S, Brinda KV, Kannan N. Protein structure: insights from graph theory. J Theor Comput Chem. 2002;1:0–0.

    Article  CAS  Google Scholar 

  28. 28.

    Amitai G, Shemesh A, Sitbon E, Shklar M, Netanely D, Venger I, et al. Network analysis of protein structures identifies functional residues. J Mol Biol. 2004;344:1135–46.

    Article  PubMed  CAS  Google Scholar 

  29. 29.

    Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res. Oxford University Press. 2000;28:235–42.

  30. 30.

    Muraki M, Harata K, Sugita N, Sato K. Origin of carbohydrate recognition specificity of human lysozyme revealed by affinity labeling. Biochemistry. 1996;35:13562–7.

    Article  PubMed  CAS  Google Scholar 

  31. 31.

    Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem. 2004;25:1605–12.

  32. 32.

    Meng EC, Pettersen EF, Couch GS, Huang CC, Ferrin TE. Tools for integrated sequence-structure analysis with UCSF Chimera. BMC Bioinf. 2006;7:339.

  33. 33.

    Dhifli W, Diallo AB. PGR: a graph repository of protein 3D-structures. In: arXiv; 2016.

    Google Scholar 

  34. 34.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. 35.

    Kozakov D, Hall DR, Xia B, Porter KA, Padhorny D, Yueh C, et al. The ClusPro web server for protein–protein docking. Nat. Protoc. Nat Res. 2017;12:255–78.

    CAS  Google Scholar 

  36. 36.

    Piovesan D, Minervini G, Tosatto SE. The RING 2.0 web server for high quality residue interaction networks. Nucleic Acids Res. 2016;44:W367–74.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. 37.

    Johnson RJK, Christodoulou J, Dumoulin M, Caddy GL, Alcocer MJC, Murtagh GJ, et al. Rationalising lysozyme amyloidosis: insights from the structure and solution dynamics of T70N lysozyme. J Mol Biol. 2005;352:823–36.

    Article  PubMed  CAS  Google Scholar 

  38. 38.

    Kannan N, Vishveshwara S. Aromatic clusters: a determinant of thermal stability of thermophilic proteins. Protein Eng. 2000;13:753–61.

    Article  PubMed  CAS  Google Scholar 

  39. 39.

    Yuan Z, Bailey TL, Teasdale RD. Prediction of protein B-factor profiles. Proteins Struct Funct Bioinf. 2005;58:905–12.

    Article  CAS  Google Scholar 

  40. 40.

    Dumoulin M, Canet D, Last AM, Pardon E, Archer DB, Muyldermans S, et al. Reduced global cooperativity is a common feature underlying the amyloidogenicity of pathogenic lysozyme mutations. J Mol Biol. 2005;346:773–88.

    Article  PubMed  CAS  Google Scholar 

  41. 41.

    Moraitakis G, Goodfellow JM. Simulations of human lysozyme: probing the conformations triggering amyloidosis. Biophys J. 2003;84:2149–58.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. 42.

    Dhulesia A, Cremades N, Kumita JR, Hsu S-TD, Mossuto MF, Dumoulin M, et al. Local cooperativity in an amyloidogenic state of human lysozyme observed at atomic resolution. J Am Chem Soc American Chemical Society. 2010;132:15580–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. 43.

    Ahn M, Hagan CL, Bernardo-Gancedo A, De Genst E, Newby FN, Christodoulou J, et al. The significance of the location of mutations for the native-state dynamics of human lysozyme. Biophys J. 2016;111:2358–67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. 44.

    Conchillo-Solé O, de Groot NS, Avilés FX, Vendrell J, Daura X, Ventura S. AGGRESCAN: a server for the prediction and evaluation of "hot spots" of aggregation in polypeptides. BMC Bioinf. 2007;8:65.

  45. 45.

    Fernandez-Escamilla A-M, Rousseau F, Schymkowitz J, Serrano L. Prediction of sequence-dependent and mutational effects on the aggregation of peptides and proteins. Nat Biotechnol. 2004;22:1302–6.

    Article  PubMed  CAS  Google Scholar 

  46. 46.

    Ahn M, De Genst E, Kaminski-Schierle GS, Erdelyi M, Kaminski CF, Dobson CM, et al. Analysis of the native structure, stability and aggregation of biotinylated human lysozyme. Massiah M, editor PLoS One 2012;7:e50192.

  47. 47.

    Frare E, Mossuto MF, Polverino de Laureto P, Dumoulin M, Dobson CM, Fontana A. Identification of the core structure of lysozyme amyloid fibrils by proteolysis. J Mol Biol. 2006;361:551–61.

    Article  PubMed  CAS  Google Scholar 

Download references


Funding for this study came from East Tennessee State University.

Availability of data and materials

Data analyzed in this study are available from the Protein Data Bank (PDB IDs: 1REX. 1LOZ, and 1LYY).

Author information




JG designed the study, performed research, analyzed data, and wrote the manuscript. PB significantly contributed to data analysis and writing the manuscript. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Jeddidiah W. D. Griffin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Griffin, J.W.D., Bradshaw, P.C. In silico prediction of novel residues involved in amyloid primary nucleation of human I56T and D67H lysozyme. BMC Struct Biol 18, 9 (2018).

Download citation


  • Amyloidosis
  • Lysozyme
  • Residues interaction networks
  • Native structure
  • B-factor