ACS Publications. Most Trusted. Most Cited. Most Read
Myelin Basic Protein Phospholipid Complexation Likely Competes with Deimination in Experimental Autoimmune Encephalomyelitis Mouse Model
My Activity
  • Open Access
Article

Myelin Basic Protein Phospholipid Complexation Likely Competes with Deimination in Experimental Autoimmune Encephalomyelitis Mouse Model
Click to copy article linkArticle link copied!

  • Anddre Osmar Valdivia
    Anddre Osmar Valdivia
    Department of Ophthalmology, Bascom Palmer Eye Institute, University of Miami, Miami, Florida 33136, United States
    Neuroscience Graduate Program, University of Miami, Miami, Florida 33136, United States
  • Pratul K. Agarwal
    Pratul K. Agarwal
    Department of Biochemistry & Cell and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, United States
    Department of Physiological Sciences  and  High Performance Computing Center and, Oklahoma State University, Stillwater, 106 Math Sciences, Stillwater, Oklahoma 74078-1010, United States
  • Sanjoy K. Bhattacharya*
    Sanjoy K. Bhattacharya
    Department of Ophthalmology, Bascom Palmer Eye Institute, University of Miami, Miami, Florida 33136, United States
    Neuroscience Graduate Program, University of Miami, Miami, Florida 33136, United States
    *Email: [email protected]. Tel: (305) 482-4103. Fax: (305) 326-6547. Department of Ophthalmology, Bascom Palmer Eye Institute, University of Miami, 1638 NW 10th Avenue, Room 707A, Miami, Florida 33136, United States.
Open PDFSupporting Information (1)

ACS Omega

Cite this: ACS Omega 2020, 5, 25, 15454–15467
Click to copy citationCitation copied!
https://doi.org/10.1021/acsomega.0c01590
Published June 16, 2020

Copyright © 2020 American Chemical Society. This publication is licensed under these Terms of Use.

Abstract

Click to copy section linkSection link copied!

Multiple sclerosis has complex pathogenesis encompassing a variety of components (immunologic, genetic, and environmental). The autoimmunogenicity against the host’s myelin basic protein is a major contributor. An increase in myelin basic protein deimination (a post-translational modification) and a change in phospholipid composition have been associated with multiple sclerosis. The interaction of myelin basic protein with phospholipids in the myelin membrane is an important contributor to the stability and maintenance of proper myelin sheath function. The study of this aspect of multiple sclerosis is an area that has yet to be fully explored and that the present study seeks to understand. Several biochemical methods, a capillary electrophoresis coupled system and mass spectrometry, were used in this study. These methods identified four specific phospholipids complexing with myelin basic protein. We show that lysophosphatidylcholine 18:1 provides a robust competitive effect against hyper-deimination. Our data suggest that lysophosphatidylcholine 18:1 has a different biochemical behavior when compared to other phospholipids and lysophosphatidylcholines 14:0, 16:0, and 18:0.

Copyright © 2020 American Chemical Society

1. Introduction

Click to copy section linkSection link copied!

Multiple sclerosis (MS) is a demyelinating disease that manifests with neurological disabilities that include impaired muscle control, tremors, and visual impairments. (1) This is mainly due to the interruption of signal propagation in different neuronal circuits through the disruption of the myelin sheath. (1) MS is a complex disease with various components contributing to pathogenesis; however, immune reactivity against myelin components seems to be a major contributor. (2−5) Early efforts sought to identify myelin components associated with MS. It was found that interferences in the interaction of myelin basic protein (MBP) with the rest of the myelin components produce immune-reactive epitopes. (2,5−7)
The myelin sheath allows for rapid transduction of nerve impulses through the electrical insulating properties of the plasma membrane. (8) Myelin sheath is composed of two major constituents: proteins and lipids. Among proteins, MBP is the second-most abundant component, whose primary role is to maintain a compact myelin sheath, while lipids provide electrical insulation. There are 14 distinct isoforms of MBP that have been reported in the murine central nervous system (CNS), all of which are the results of different exon compositions of a single MBP gene. (9) The classic family of MBP is the one predominantly expressed in the major dense line of the myelin sheath, and isoform 5 (18.5 kDa) is the isoform primarily associated with MS. The most prominent immune-reactive epitopes found in the cerebrospinal fluid of MS patients are against MBP peptides. (10,11) These peptides have been documented to be selectively produced by proteolytic cleavage by a number of metalloproteases. (12−14) The rate of digestion at which metalloproteases produce these immune-reactive peptides has been directly linked to the levels of deimination (a post-translational modification) in MBP.
Deimination is a post-translational modification mediated by the family of peptidyl arginine deiminases (PADs) that converts protein-bound arginine to protein-bound citrulline. (15) This modification effectively neutralizes the positively charged guanidinium group found in arginine by replacing it with a neutral ureido group. (15−17) Deimination is considered to be an irreversible post-translational modification given that the enzyme responsible for the reverse reaction has not been found. (18) In addition, due to the longevity of deimination, this poses the potential for long-lasting effects on protein regulation.
Studies focused on the structure of deiminated MBP have documented that deimination renders MBP to an open structure, which allows proteases access to internal residues, increasing the rate of digestion and production of immune-reactive peptides. (10,19−22) The elevated levels of MBP deimination have been documented to positively correlate with the severity of MS in human subjects, as well as in the experimental autoimmune encephalomyelitis (EAE) mouse. (6,11,23−29) Subsequent studies demonstrated that native MBP takes a C-shape folded structure upon lipid binding, which resulted in an increase in compaction and stability. (30−32) This highlights the importance that MBP–phospholipid complexation has in maintaining myelin stability.
Lipids in the myelin sheath are arranged in an asymmetric pattern. Sphingolipids are the major class of lipids that face the extracellular surface, and phospholipids are the major class facing the cytoplasmic space. (33) Changes in lipid composition have been documented in human MS brains. Members of phospholipid classes [phosphatidylethanolamines (PEs), phosphatidylcholines (PCs), phosphatidylinositols (PIs), and phosphatidylserines (PSs)] have been reported to increase in relative amounts (relative to the total) for specific lipid species and decrease for others. (34−42) MBP is a cytoplasmic side peripheral membrane protein likely to complex with phospholipids (6,43) and help maintain an ordered myelin sheath.
Although studies have demonstrated that MS pathogenicity is associated with changes in lipid profiles and disruption of MBP structural roles, a study exploring the specific phospholipid interactions that are native to MBP remains to be performed. The effects of MBP–phospholipid complexation on protein deimination also remain to be studied. The aim of the study is to address this gap in the scientific literature, that is, to understand the individual phospholipid that complexes with MBP and plays a potential role toward influencing protein deimination.

2. Results and Discussion

Click to copy section linkSection link copied!

Demyelinating diseases such as multiple sclerosis converge diverse components resulting in pathogenicity. The immunological events, despite being most well studied, still warrant further intense investigation. To contribute to our understanding of disease pathogenicity, it is crucial to understand the changes in the basic biochemical conditions associated with the myelin sheath. Autoimmunogenicity against host’s MBP is a hallmark of multiple sclerosis. Understanding the native biochemistry can shed light on how the state of antigenicity becomes perturbed in disease conditions to result in autoimmunogenicity.
The animal models available for the study of MS include the toxin-induced models (cuprizone and lysolecithin, useful for understanding demyelination and remyelination processes), viral-induced model (Theiler’s murine encephalomyelitis, useful for understanding axonal damage and inflammation), and chronic EAE model. (44) The EAE mouse model is one of the most characterized models in the literature. The benefits of this model include chronic demyelination and its effects on the central nervous system, allowing for the study of immunological and neuroinflammatory processes.
We present findings that (1) PAD-mediated hyper-deimination is associated with EAE (Figure 1B), (2) there are broad changes of phospholipid composition in EAE (30 days post injection; Figure 2C), (3) there is a differential complexation of phospholipids with MBP in EAE (Figure 3B,C), and (4) lysophosphatidylcholine (LPC) 18:1 likely competes for hyper-deimination of MBP (Figures 4A,B, 5A, and 7A). Our findings in the EAE model are parallel to the reports of phospholipid deficiency in human multiple sclerosis subjects only in the PC class. (34−42) Although there were overall changes in the general classes of phospholipids, we focused on the specific phospholipids that were deficient in the EAE group, thus enabling the selection of a finite number of species that can be pursued with ease. Specific phospholipid species have been reported to be deficient in human MS samples despite increases in the overall class of phospholipids. (34−42)

Figure 1

Figure 1. Myelin basic protein in EAE is deiminated. (A) Representative sucrose density ultracentrifugation fractions of mouse brains probing for myelin basic protein (Abcam-ab7349). *Sucrose fraction 32 (SF 32) from each group was used for further analysis (coomassie, Western blot, mass spectrometry, and capillary electrophoresis (CE)). EAE, n = 15; sham-injection control, n = 10; and no-injection control, n = 10. Each immunoblot was repeated three times. (B) Coomassie blue stain and Western blot for sucrose fraction 32. Bands were excised from coomassie gel and analyzed by mass spectrometry. Western blot probed for levels of deimination (citrulline, Millipore-MABS54887). (C) Mass spectrometry analysis of excised bands from coomassie blue stain gel. PSM = peptide-spectrum match.

Figure 2

Figure 2. Myelin basic protein isolated from EAE has different phospholipid compositions than controls. (A) Representative capillary electrophoresis (CE) chromatogram of sucrose fraction 32 (SF 32). Arrow points toward the peak collected for protein and phospholipid analysis. (B) Dot immunoblot for MBP in CE fractions, 10 min before the peak, peak (as noted by the arrow), and 10 min after the peak. MBP was probed with an anti-MBP antibody (Abcam, ab7349). (C) Phospholipid analysis of CE peak fractions, three replicates, and bars correspond to standard error mean (SEM) (*p < 0.05).

Figure 3

Figure 3. Confirmation of MBP–phospholipid complexation. (A) Structure of candidate lipids: PC 16:0/22:6, 1-((4Z,7Z,10Z,13Z,16Z,19Z)-docosahexaenoyl)-2-hexadecanoyl-sn-glycero-3-phosphocholine, LMGP01011116; PC 16:1/16:1, 1,2-di-((9Z)-hexadecanoyl)-sn-glycero-3-phosphocholine, LMGP01010684; PI 18:0/20:4, 1-octadecanoyl-2-((5Z,8Z,11Z,14Z)-eicosatetraenoyl)-sn-glycero-3-phospho-(1′-myo-inositol), LMGP06010010; PS 18:0/18:1, 1-octadecanoyl-2-((9Z)-octadecanoyl)-sn-glycero-3-phosphoserine, LMGP03010025; LPC 18:1, 2-((9Z)-octadecanoyl)-sn-glycero-3-phosphocholine, LMGP01050082. (B) Liposome flotation assay (LFA) of candidate phospholipids demonstrating MBP complexation with phospholipids. B, bottom fraction; M, middle fraction; and T, top fraction (MBP, Abcam, ab7349). (C) Protein–lipid overlay assay (PLOA) of candidate phospholipids demonstrating MBP complexation with phospholipids (MBP, Abcam, ab7349).

Figure 4

Figure 4. Complexation of LPC 18:1 with MBP is competitive for MBP hyper-deimination. (A) Dot blot analysis of deimination reaction for MBP (MBP, Abcam-ab7349) and deiminated MBP (citrulline, Millipore-MABS54887) in the presence of candidate phospholipids. PAD, peptidyl arginine deiminase. (B) Deimination fold change analysis of Western blot, three replicates, and densitometry using ImageJ was normalized to MBP (***p < 0.0005, comparison to MBP–PAD). (C) Representative circular dichroism (CD) spectra of deimination reaction for MBP in the presence of candidate phospholipids, three replicates per spectra. (D) CD spectra analysis of secondary structure probability, utilizing CD analysis and plotting tool (CAPITO) software (mean with SEM).

Figure 5

Figure 5. LPC 18:1 competes with hyper-deimination of MBP arginine residues and has differential complexation with MBP. (A) Mass spectrometry analysis of MBP arginine residues that are being deiminated in in vitro deimination assay (Millipore, 13-104, Burlington, MA). The deiminated Myelin basic protein (UniProt accession: P02687) peptides were identified by mass spectrometry (see Table 3). Gray highlight, deiminated arginine; red bold square, arginine not deiminated in LPC 18:1 reaction (residues 63, 112, and 129). (B) Protein–lipid overlay assay of LPC 18:1, LPC 14:0, LPC 16:0, and LPC 18:0 MBP was probed with an anti-MBP antibody (Abcam, ab7349).

2.1. Confirmation of the Experimental Autoimmune Encephalomyelitis (EAE) Mouse Model

To confirm that the multiple sclerosis mouse model utilized for this study (EAE) exhibited its characteristic traits, animals were clinically scored to monitor the progression of EAE (Figure S1A). The clinical score graph demonstrated that animals maintained an elevated clinical score and corresponded to the chronic phase of disease progression (day 30). Additionally, visual impairment was monitored by measuring the inner retinal function (retinal ganglion cell activity) using pattern electroretinogram (PERG). The PERG measurements demonstrated that 30 days post immunization, EAE mice exhibited a decrease in amplitude when compared to the no-injection control and sham-injection control groups (Figure S1B). Flash electroretinogram (FERG) was utilized to determine if the photoreceptor activity was contributing to visual impairment. The FERG demonstrated no change in any of the groups, suggesting that visual impairment was mediated through the inner retina (data not shown). To eliminate the possibility that an aberrant retinal morphology was contributing to visual decline, optical coherence tomography (OCT) was utilized and demonstrated no change in retinal thickness (Figure S1C). Immunohistochemical analysis of optic nerves demonstrated that EAE mice were suffering from optic nerve demyelination as it was observed by a decrease in the MBP signal and topographical distribution (Figure S1D). The blood antibody analysis specific for myelin oligodendrocyte glycoprotein (MOG) further corroborated that experimental animals exhibited characteristic traits of the EAE mouse model (Figure S1E). Thus, our EAE animals demonstrated all clinical neurological deficits, as well as visual function impairment, in parallel to that in human MS.

2.2. Isolation of the MBP–Phospholipid Complex and Deimination in EAE

Several studies have highlighted the importance of MBP–phospholipid complexation; however, the specific lipid species that complex with MBP have not been explored in detail. (10,19−22,30−32) The isolation of intact complex using inert biochemical methods is the first critical step toward the identification of native MBP–phospholipid complex constituents. The inert methods preserve the intermolecular interactions within the complex. From confirmed EAE mice, we collected neuronal tissue (brain) and subjected the sample to homogenization and fractionation through sucrose gradient ultracentrifugation. Sucrose fractions were then evaluated using immunoblotting for the identification of the fractions containing MBP (Figure 1A). The distribution of MBP was rather uniform throughout many different fractions (fractions 14–36; Figure 1A). Sucrose fraction 32, taken as a representative, was analyzed for verification of the presence of MBP and levels of deimination. Coomassie staining demonstrated the same pattern of MBP within the gel as observed in the Western blot analysis (Figure 1B, upper panel); the latter also showed an increase in levels of deimination in the EAE sucrose fraction (Figure 1B, lower panel). Bands from coomassie gel were excised and analyzed, verifying the presence of MBP in the fractions (Figure 1C). Sucrose fraction 32 (SF 32) from all groups was further fractionated through capillary electrophoresis (CE). The CE electropherogram displayed large peaks, as indicated by the arrows. The peaks did not demonstrate a constant retention time; however, they typically eluted between 55 and 65 min (Figure 2A). Peaks were collected as well as the eluents 10 min before and after the peak (each for 10 min duration). They were subjected to dot immunoblot analysis using the antibody against MBP. Results demonstrated the presence of MBP in the peak but not in the other fractions, confirming the presence of MBP and its constituents in the peaks (Figure 2B). We therefore collected MBP complexes at the end of CE fractionation.

2.3. MBP–Phospholipid Complexation

Mass spectrometric lipidomic analysis of the CE peak fractions demonstrated a decrease in phosphatidylcholines (PCs), no change in phosphatidylethanolamines (PEs), and an increase in phosphatidylinositols (PIs) and phosphatidylserines (PSs) in the EAE group when compared to that of control groups (Figure 2C). The following criteria were used to identify candidate phospholipids that potentially complex with MBP: (1) phospholipids had to be present in both the no-injection control and sham-injection control groups but (2) absent in the EAE group. We found five candidate lipids that met these criteria (Table 1 and Figure 3A).
Table 1. Candidate Phospholipids after CEa
lipid moleculecalculated mass
PC 16:0/22:6805.5622
PC 16:1/16:1729.5309
PI 18:0/20:4858.5258
PS 18:0/18:1789.552
LPC 18:1521.3481
a

Candidate phospholipids that are present in the noninject and sham groups but absent in the EAE group.

To corroborate their direct complexation with MBP, the candidate phospholipids were analyzed using the liposome flotation assay (LFA) and protein–lipid overlay assay (PLOA) (Figure 3B,C). For this purpose, liposomes/micelles were prepared using extrusion and were confirmed for being within the appropriate size range, using nanoparticle tracking analysis. (45) The assay revealed that there are various degrees of complexation that exist between MBP and candidate phospholipids, as revealed by the movement of MBP to either the middle or top layer (Figure 3B).
The PLOA analysis demonstrated and confirmed the same pattern observed in the LFA (Figure 3C). However, PC 16:0/22:6 consistently presented difficulties adhering to the membrane during dot immunoblot in the PLOA and was not pursued for further analysis (Table 2). The PLOA demonstrated that LPC 14:0, LPC 16:0, and LPC 18:0 had equal complexation to MBP between the control groups and EAE group. This contrasts to the complexation of LPC 18:1, which showed a difference between MBP derived from controls and EAE (Figure 5B).
Table 2. Curated Phospholipids after LFA and PLOAa
lipid moleculecalculated mass
PC 16:1/16:1729.5309
PI 18:0/20:4858.5258
PS 18:0/18:1789.552
LPC 18:1521.3481
a

Lipids demonstrating direct complexation with MBP. Liposome flotation assay (LFA) and protein–lipid overlay assay (PLOA).

Comparison of phospholipids deficient in the EAE with the control groups helped identify and validate four phospholipids that complex with MBP, of which PC 16:1/16:1 and LPC 18:1 have been reported to be deficient in humans (Figure 3 and Table 2). (40,42)
Lipids in an aqueous environment can attain different geometries due to their hydrophobic tails and hydrophilic head groups. The packing parameters (P) dictate that the geometry lipids will undertake in an aqueous environment. (46,47) Phospholipids (P ∼ 1) and lysophospholipids (P < 1) exhibit a cylindrical and an inverted cone geometry, respectively. In both of these species, the hydrophilic head group will face the outer shell of the liposome (phospholipids) or micelle (lysophospholipids). The myelin sheath mimics the same geometry, that is, the head group is facing the cytoplasmic space and interacting with MBP. Only lipids with a P > 1 (for example, dioleoylphosphocholine) will exhibit an inverted micellar geometry (our study did not identify any such lipid species).
The variable degree of travel observed in the LFA for specific MBP–phospholipid complexes can be accounted for by either the strength of the MBP–phospholipid complexation or the presence of liposomes/micelles of smaller sizes. The smaller-sized liposomes/micelles may fail to move MBP to the top layer due to a lack of attaining sufficient buoyancy. PC 16:0/20:4 was used as a control phospholipid to PC 16:0/22:6. This control was used to determine if reducing the length and saturation of the hydrocarbon chain affected complexation with MBP (Figure S2A,B). LFA and PLOA for PC 16:0/20:4 demonstrated the same pattern observed with PC 16:0/22:6, suggesting that small variations in the length and saturation of the hydrocarbon chain do not affect MBP complexation. Similarly, PC 18:1/18:1 and LPCs (14:0, 16:0, and 18:0) were used as controls for LPC 18:1. These controls were used to determine if the absence of a second hydrocarbon chain and variations in the hydrocarbon chain length affected MBP complexation (Figures S2C,D and 5B). LFA and PLOA demonstrated that PC 18:1/18:1 did not complex with MBP, suggesting that the number of hydrocarbon chains has an effect in complexation with MBP. PLOA of other LPCs demonstrated no preference for MBP complexation derived from controls or EAE group.

2.4. LPC 18:1 Likely Competes for Hyper-deimination of MBP

EAE sucrose fraction demonstrated an increase in deimination (Figure 1B, lower panel); hence, the effects of phospholipid complexation with MBP were analyzed for their competitive effects against hyper-deimination. We present evidence that LPC 18:1 competes for PAD-catalyzed hyper-deimination of arginine residues in MBP when compared to the other phospholipids (Figures 4A,B and 5A). In the present study, LPC 18:1 was identified to be present in healthy, but absent in chronic EAE mice. The absence in EAE suggests that endogenous LPC 18:1 can potentially have important biological functions in healthy animals. The complexation of LPC 18:1 is likely to be beneficial for MBP stability, preventing MBP from further deimination when challenged with PAD. It is likely that LPC 18:1 complexation occurs prior to deimination and thus prevents hyper-deimination. In normal humans, MBP is known to be deiminated in 6 of a total of 19 arginine residues (UniProt accession: P02686-5, in MBP isoform 5; Figure S5). However, deimination in MBP isoform 5 goes beyond these six sites in MS and/or in demyelinating diseases. (48) In other words, hyper-deimination is prevented in normal conditions but occurs only in demyelinating diseases/MS. (6,10,24,28,49) The more severe the demyelinating disease, the more intense the hyper-deimination. (23) In mouse and bovine, four and eight residues are found deiminated in normal or wild-type animals, respectively. The number of sites (beyond four and eight, respectively) increases in mouse and bovine in demyelinating disease conditions consistent with hyper-deimination in humans (Figure S5).
The complexation of LPC 18:1 with MBP is partial in vitro as demonstrated by LFA of MBP in the middle layer and the partial signal in the PLOA, but more complete complexation in vivo cannot be ruled out. When LPC 18:1–MBP complex is challenged by PAD, only nine arginine residues are deiminated (one more residue compared to eight endogenous residues in MBP). In contrast, all other LPCs (14:0, 16:0, and even 18:0) and other phospholipids demonstrated deimination of 13 and 11–13 arginine residues, respectively (Figure 5A and Table 3). This supports the hypothesis that MBP–LPC 18:1 complexation prevents hyper-deimination of arginine residues by PAD (Figure 7B,C). LPC 18:1 complexation with MBP derived from the EAE group (Figure 3B,C) is consistent with the availability of strong binding sites (vacant sites) for LPC 18:1 in hyper-deiminated MBP. In vivo, the EAE animals (or human MS CNS samples) are deficient in LPC 18:1, leaving such binding sites vacant in MBP derived from these sources. This is consistent with our CE experiments, where we found a lack of bound LPC 18:1 with the EAE-derived MBP (Table 1). LPC 18:1 confers protection against hyper-deimination, contrasting other LPCs (14:0, 16:0, and 18:0). These LPCs are the ones frequently used for demyelination (at or beyond 100 μM, they induce toxicity). (50−53) A comparative study of the degree of demyelination caused by various LPC species remains to be performed.
Table 3. Arginine Residues Deiminated in MBP–Lipid Challenge Assaya
lipid moleculetotal deiminated citesprotein coverage (%)arginine residue deiminatedb
MBP + PAD + PC 16:1/16:11385.8041, 47, 52, 63, 78, 96, 106, 112, 129, 158, 161, 168, 169
MBP + PAD + PI 18:0/20:41285.8047, 52, 63, 78, 96, 106, 112, 129, 158, 161, 168, 169
MBP + PAD + PS 18:0/18:11185.8052, 63, 78, 96, 106, 112, 129, 158, 161, 168, 169
MBP + PAD + LPC 18:1994.0847, 52, 78, 96, 106, 158, 161, 168, 169
MBP + PAD + LPC 18:01392.9023, 47, 52, 63, 78, 96, 106, 112, 129, 158, 161, 168, 169
MBP + PAD + LPC 16:01392.9023, 47, 52, 63, 78, 96, 106, 112, 129, 158, 161, 168, 169
MBP + PAD + LPC 14:01392.9023, 47, 52, 63, 78, 96, 106, 112, 129, 158, 161, 168, 169
MBP + PAD + no lipid1294.0847, 52, 63, 78, 96, 106, 112, 129, 158, 161, 168, 169
MBP + No PAD + no lipid894.0852, 78, 96, 106, 158, 161, 168, 169
a

Analysis of arginine deimination of myelin basic protein (MBP). For chymotrypsin-digested MBP (UniProt accession number: P02687; Millipore, 13-104, Burlington, MA), the total number of spectrum matches were 1756 and 20 unique peptides, which were identified by mass spectrometry. The deimination cites had a 99–100% confidence. LPC 18:1 incubation protected residues 63, 112, and 129 from deimination.

b

The sequence positions are in reference to canonical MBP isoform 1 (UniProt accession number: P02687). See Table S1 for more details.

2.5. Structural Analysis of LPC 18:1–MBP Complex

Sequence analysis-based computer methods and full-atomic-level molecular dynamics (MD) simulations are consistent with MBP being an intrinsically disordered protein (IDP). MBP exhibiting a disordered structure (Movie S1 and Figure S3) is supported by previous literature reports. (43,54,55) This was further corroborated by circular dichroism (CD) studies, which demonstrated that despite MBP being an intrinsically disordered protein (IDP), it retains a partial secondary structure. The CD showed altered spectra when MBP was denatured at 105 °C (Figure 4C, top left panel) compared to that at 20 °C. The LPC 18:1 binding does not cause any change in the secondary structure (Figure 4C, bottom left panel). In hyper-deiminated MBP, the loss of secondary structure (Figure 4C, top right panel) is quite pronounced (Figure 4C, top right panel), as expected. However, the loss of secondary structure in hyper-deiminated MBP is somewhat protected (not statistically significant) when complexed with LPC 18:1 (Figure 4C, top right panel and bottom left panel, respectively).
Computer simulations indicated that MBP is a highly flexible protein that shows a variable conformational landscape, a behavior characteristic of IDPs (see Movies S1–S3 in the Supporting Information). MD simulations indicated that the presence of lipids could change the level of deimination of arginine residues due to accessibility by PAD. This change could be due to a combination of two possible factors. The first factor is the presence of the arginine residues on the MBP surface rendering accessibility by PAD. The presence of lipid molecules around the arginine residues could make them inaccessible to PAD. The second factor is the overall conformation of the protein, which is altered by the presence of the lipid. In the extended conformation of MBP, more arginine residues are accessible or present on an accessible surface (Figure S4). The water-only simulations indicated that the arginine residues are solvent accessible (located on the surface) and could be easily accessed by PAD (Figure 6). The alternate case would be when these residues are not accessible to the solvent due to salt linkages or being buried inside the protein. In the LPC 18:1 simulations, the lipid molecules surround the arginine residues, possibly making it difficult for deimination by PAD to occur.

Figure 6

Figure 6. Change in MBP conformations and solvent accessibility of arginine residues. Representative structures are shown from simulations only in water, water with LPC, and water with PC. Helices are shown in red, β-strands in yellow, three N-terminal residues are colored green, and three C-terminal residues are colored cyan. Arginine residues are shown as green/blue sticks, while the lipid is shown as magenta/blue/orange sticks. Water was included in the simulations but omitted here for visualization purposes. In the case of water-only, the arginine residues are exposed on the surface; the presence of LPC protects a number of these residues, and in the presence of PC, MBP adopts an extended conformation, exposing several arginine residues to the solvent.

Figure 7

Figure 7. Hypothesis model of MBP–phospholipid complexation. (A) Working hypothesis for membrane disruption, incorporating key events associated with unstable myelin. The order of events has yet to be determined. (B) Model of phospholipid (PC 16:1/16:1, PI 18:0/20:4, PS 18:0/18:1, LPC 18:0, LPC 16:0, and LPC 14:0) competitive effects against PAD-mediated hyper-deimination and potential cascade of events. (C) Model of LPC 18:1 competitive effects against PAD-mediated hyper-deimination. Δ, “change in”; MBP, myelin basic protein; and PAD, peptidyl arginine deiminase. Open and C-shaped conformations are consistent with the previously published reports.

The behavior in the PC simulations is significantly different. First, in the presence of PC, there are several more arginine residues that appear to be solvent accessible. Second, there is a change in the structure of MBP. The averaged radius of gyration of MBP in the presence of PC is 19.1 Å, compared to 18.2 Å in the case of water-only simulations and 18.7 Å in the presence of LPC 18:1. The more extended conformations in the presence of PC allow more arginine residues to be accessible by other biomolecules, including PAD. This further supports that steric hindrance provided by LPC 18:1 competes against PAD-mediated hyper-deimination.

3. Conclusions

Click to copy section linkSection link copied!

Together, the data presented in this study demonstrate that LPC 18:1 undergoes complexation with MBP. This complexation prevents hyper-deimination when challenged with PAD. The severity of MS has shown strong and positive correlation with hyper-deimination of MBP. (6,10,23,24) Hyper-deiminated MBP is more easily unfolded (10) and becomes vulnerable to autolysis (56) and proteolysis, (22) resulting in the release of MBP peptides, which in turn triggers autoimmunogenicity. (10,28,57) Our findings reported here expands our understanding of the biochemistry associated with the previously recorded molecular events in MS pathogenicity (Figure 7).

4. Experimental Section

Click to copy section linkSection link copied!

4.1. Generation of Experimental Autoimmune Encephalomyelitis (EAE) Model

All animals were housed at the McKnight vivarium, University of Miami, FL, and animals utilized in this study were under the University of Miami Institutional Animal Care and Use Committee-approved protocols (IACUC protocol number: 16-235). Method was adapted from other studies. (58) Immunizing neuroantigen was prepared by the emulsion of 2 mg of myelin oligodendrocyte glycoprotein (MOG35-55) (Biosynthesis Inc., 12668, Lewisville, TX), incomplete Freund’s adjuvant (VWR, 90003-748, Radnor, PA), and 6 mg of Mycobacterium tuberculosis (H37RA, Fisher Scientific, DF3114-33-8, Hampton, NH). This stock solution was used to immunize six animals, in which each individual animal received 0.333 mg of myelin oligodendrocyte glycoprotein (MOG35-55) and 1 mg of M. tuberculosis. Pertussis toxin (Biological Laboratories, 181, Campbell, CA) was prepared at a concentration of 50 μg/mL in sterile phosphate-buffered saline (PBS), pH 7.4, with an intraperitoneal (i.p.) administration of 200 ng per injection. Female C57Bl6/J mice were injected i.p. with pertussis toxin the day before injection. They were then immunized subcutaneously (s.c.) at 2 months of age with 200 μL of neuroantigen (MOG35-55) followed by an i.p. injection of pertussis toxin 2 days post immunization. Three groups were created for this study, an immunized group with MOG35-55 (EAE, 15 animals), an immunized group without MOG35-55 (sham-injection control, 10 animals), and an injection control that received no injections (no-injection control, 10 animals). Animals were monitored for decline in body weight as well as manifestation of clinical signs indicative of encephalomyelitis following the same criteria as previous studies. (58)

4.2. Electroretinogram Recordings

Pattern electroretinogram (PERG) and flash electroretinogram (FERG) were recorded utilizing the Jörvec PERG system (Jörvec PERG Visual Stimulation Box, M014760L, Miami, FL). The experimental setup for FERG was as described in previously published protocols; however, there was a slight modification for the experimental setup for PERG. (59) The modification consisted of the subcutaneous insertion of the recording electrode between the two eyes, the reference electrode in the scalp and the ground electrode in the lower back adjacent to the tail. For both PERG and FERG, animals were anesthetized by the i.p. injection of 100 μL of ketamine and xylazine cocktail (1.5 mg/0.3 mg per 100 μL) per 20 g of body weight. A drop of the balanced saline solution was delivered to the eye to prevent dryness for the duration of the procedure. Animals were situated in a heating pad (Physitemp TCAT-2LV controller) with an anal reference thermometer set at 37 °C to stabilize the internal body temperature. PERG was recorded by three consecutive responses of 600 contrast reversals, and recording settings consisted of 10.0 K gain, 1.0 Hz high pass, 100.0 Hz low pass, and 360.0 μV rejection. FERG was recorded by three consecutive measurements of a flash, with a strength of 20.0 CDS/SM and a frequency of 1.0 Hz. Amplitude for PERG and FERG was measured as the difference in microvolts between the highest peak and the consecutive lowest trough. Latency for PERG and FERG was measured as the time (milliseconds) it took to reach the highest peak (from time zero to the time when the highest peak was recorded).

4.3. Optical Coherence Tomography (OCT)

Animals were anesthetized as described above, followed by the application of one drop of tropicamide ophthalmic solution (Akorn, 17478-102-12) to dilate the pupil and a drop of the balanced saline solution to prevent eye dryness. A 3.5 mm diameter mini contact lens (Ocuscience, ERGACC3) was placed for the visualization of retina. OCT images were taken with the Engineering Spectralis HRA + OCT (Franklin, MA) and analyzed with the Spectralis HRA + OCT software version 6.6.

4.4. Enzyme-Linked Immunosorbent Assay (ELISA)

Animals were anesthetized as described above, followed by exposure of the thoracic cavity. Using a 26G1/2 needle, about 1 mL of blood was aspirated from the left atrium. Animals were immediately decapitated to ensure euthanasia and prevent recovery from anesthesia. The blood was treated with 0.25 mL of anticoagulant citrate dextrose solution per 1 mL of blood. Blood was layered over 2 mL of Ficoll-Paque PLUS, a density of 1.077 g/mL (GE Healthcare, 17-1440-02), and centrifuged at 1800g for 35 min. Plasma was collected and store at −80 °C. Utilizing the AnaSpec Anti-MOG35-55 Quantitative ELISA kit (VWR, AS-54465) and following the manufacturer’s protocol, antibodies against MOG35-55 were detected and normalized to the total protein concentration in isolated blood plasma.

4.5. Immunohistochemistry

Dissected optic nerves were immersed fixed in 4% paraformaldehyde in PBS and incubated overnight at 4 °C. Tissue was washed using PBS, followed by cryoprotection using a gradient of 10% sucrose (overnight at 4 °C), 20% sucrose (overnight at 4 °C), and 30% sucrose (overnight at 4 °C). Tissue was then embedded in an optimum cutting temperature compound (VWR, 25608-930) and stored at −80 °C. Sections were collected at 10 μm thickness and stored at −80 °C. Immunohistochemistry standard protocol was used to detect MBP (Abcam, ab7349) and CD90 (CYMBUS Biotech, CBL1354), and the tissue was mounted using Vectashield with 4′,6-diamidino-2-phenylindole (DAPI) (Vector Laboratories, H-1200).

4.6. Tissue Preparation and Density Ultracentrifugation

Method was adapted from other studies. (60) Brains from 30 days postinjection mice (8–10 brains, ∼3.5 g total weight, when they manifested clinical signs indicative of encephalomyelitis) were rapidly dissected, flash-frozen by liquid nitrogen immersion followed by equilibration on dry ice, and stored at −80 °C. Tissue was homogenized at 4 °C utilizing 0.32 M sucrose solution with a protease inhibitor (5% w/v homogenate). Gradient I was established by layering the homogenate (18 mL) over 18 mL of 0.85 M sucrose at 4 °C and ultracentrifuged at 140 000g at 4 °C for 1 h (Beckman Rotor SW28). The low-density band was discarded, and the main band along with the dispersion band was collected. They were homogenized and diluted to 200 mL, utilizing a 10 mM ethylene glycol tetraacetic acid (EGTA) solution followed by centrifugation at 141 000g at 4 °C for 40 min (Beckman Rotor SW28). The supernatant was discarded, and the pellet was resuspended in 120 mL of 10 mM EGTA solution, homogenized, and ultracentrifuged at 35 000g at 4 °C for 15 min (Beckman Rotor SW28). This step was repeated once. The supernatant was discarded, and the pellet was resuspended in 18 mL with 0.85 M sucrose at 4 °C. Gradient II was established by layering from bottom to top over 2 mL of 1 M sucrose, 18 mL of 0.85 M sucrose-containing pellet, and 18 mL of 0.32 M sucrose at 4 °C followed by ultracentrifugation at 140 000g at 4 °C for 1.5 h (Beckman Rotor SW28). The low-density band was discarded, and the main band along with the dispersion band was collected followed by processing as described above; however, the pellet was finally resuspended in 13 mL with 0.85 M sucrose solution. Gradient III was established by layering from bottom to top over 12 mL of 0.85 M sucrose-containing pellet, 12 mL of 0.75 M sucrose, and 13 mL of 0.32 M sucrose at 4 °C followed by ultracentrifugation at 14 0000g at 4 °C for 16 h (Beckman Rotor SW28). Fractions were collected in separate tubes (1 mL each, 37 fractions total) along with the pellet. Fractions were washed with 2 mL of 2 mM EGTA solution followed by centrifugation at 18 000g 4 °C for 25 min (Beckman Rotor GH-3.8). The supernatant was discarded, and the pellets were resuspended in 1 mL of 50 mM Tris solution with a protease inhibitor and stored at −80 °C.

4.7. Gel Electrophoresis and Immunoblotting

Fractions from density centrifugation were heated to 85 °C for 3 min and diluted using Laemmli loading buffer (Amresco, M337). Equal volume of samples was loaded in 4–20% Tris–glycine gel (Bio-Rad, 5671095) to maximize the protein content in lower fractions and ran using Tris–glycine–sodium dodecyl sulfate (SDS) running buffer at 200 V for 33 min. Poly(vinylidene difluoride) (PVDF) membrane was activated with methanol and washed with ultrapure water and with transfer buffer (TGS-20% methanol) before transfer. Semidry transfer was carried out using the Trans-Blot Turbo Transfer System (Bio-Rad, 1704150). Membrane was blocked with a blocking buffer (Bio-Rad, 170-6404) and probed for either MBP using antibodies (Abcam, ab7349) at a 1:500 dilution or deiminated MBP (Cit-MBP, Millipore, MABS54887) at a 1:500 dilution. Membrane was washed using 25 mM Tris, 2.6 mM KCl, 0.14 M NaCl, 0.2% Tween-20, at pH 8, and chemiluminescence was detected using Western blotting substrate (Thermo, 32106). Dot immunoblotting was performed as described, with the exception that samples were not separated by gel electrophoresis. A modification to this protocol was made when detecting deiminated MBP by modification of the citrulline residues through a reaction with 2,3-butanedione monoxime and antipyrine in a strong acid solution (Millipore, 17-347B), which was performed immediately after transfer, and modification was incubated at 37 °C for 3 h.

4.8. Capillary Electrophoresis (CE)

After confirming the presence of MBP, fraction 32 was selected for further analysis. Fractions were separated using the Agilent Capillary Electrophoresis 7100 coupled with a bare fused silica capillary tube (Agilent, G1600-63311). (61) CE buffer consisted of 100 mM formic acid, 5 mM ammonium acetate, and 100 mM Tris at a pH of 7.3. The sample was prepared at a 1:10 ratio of density centrifugation fraction to CE buffer. Parameters for the CE system included injection at 50 mbar for 10 s, high voltage at 30 kV, 25 μA and 6 W, a cassette temperature of 25 °C, coupled with a high-pressure system of 10 mbar for the duration of the separation (70 min). Signal was detected using UV absorbance at 230 nm (bandwidth 4 nm) and a reference at 350 nm (bandwidth 50 nm). Fractions were collected in separate tubes in intervals of 10 min.

4.9. Lipid and Protein Extraction

Capillary electrophoresis fractions were centrifuged at 18 000g at 4 °C for 1 h (Beckman Microfuge 18). The supernatant was discarded, and the pellet (very small) was resuspended in methanol with 0.516 mg/mL butylated hydroxytoluene (BHT). The resuspended pellet was transferred to a glass tube to minimize contact with plastic, and methyl tert-butyl ether (MTBE) was added, maintaining a 10:3 ratio per volume of MTBE/methanol. Extraction was left shaking overnight (∼16 h) at 4 °C. Extractions were treated with 0.15 M ammonium acetate for phase separation while maintaining a 20:6:5 ratio per volume of MTBE/methanol/ammonium acetate. Extractions were centrifuged at 2000g at 4 °C for 10 min (Thermo Fisher Megafuge 8R), and the upper organic phase was collected in a separate glass tube. Extraction containers were washed with a 20:6:5 ratio per volume of MTBE/methanol/ammonium acetate, centrifuged, and the organic phase was added to the first collection. Organic phase samples were completely dried, utilizing speed vacuum at a temperature of 37 °C; the dried samples were stored at −80 °C and resuspended in 1:1 chloroform/methanol for mass spectrometry analysis. The aqueous phase was used for liposome flotation assay and protein–lipid overlay assay.

4.10. Sample Preparation for Protein Mass Spectrometry

Fifteen micrograms of total protein from the aqueous phase sample was added to four times the volume of acetone at 4 °C and incubated at room temperature for 15 min. Samples were centrifuged at 21 000g at 4 °C for 30 min (Thermo Fisher Megafuge 8R). The supernatant was discarded, and the pellet (very small) was air-dried for 10 min. The pellet was resuspended and reduced with 10 mM dithiothreitol in 100 mM ammonium bicarbonate for 45 min at room temperature. Following reduction, the resuspended pellet was alkylated with 55 mM iodoacetamide in 100 mM ammonium bicarbonate for 30 min while maintained in darkness. Two times the volume of acetone at room temperature was added and centrifuged at 21 000g at room temperature for 15 min (Thermo Fisher Megafuge 8R). The supernatant was discarded, and the pellet (very small) was air-dried for 10 min. The pellet was resuspended in 0.1 μg/μL of chymotrypsin (Promega, V106A) in 15 mM 4-ethylmorpholine and incubated overnight (∼16 h) at 37 °C. Samples were centrifuged at 21 000g at room temperature for 5 min (Thermo Fisher Megafuge 8R) and dried, utilizing speed vacuum at a temperature of 40 °C, leaving about 50 μL of sample. Samples were stored at −80 °C.

4.11. High-Performance Liquid Chromatography (HPLC)–Mass Spectrometry

Mass spectrometry lipidomic chromatography and identification were carried out as described in detail in our previous publication. (62) In brief, samples were run through an Acclaim C30 column (particle size 3.0 μm, 150 × 2.1 mm2 ID; Thermo Fisher Scientific, Waltham, MA). An HPLC Accela instrument (equipped with an autosampler and a 600 pump) was coupled to a Q-Exactive mass spectrometer (Thermo Fisher Scientific). Heated electrospray ionization (HESI) was used as the method of ionization by coupling a HESI probe to the Q-exactive instrument. The instrument was set at full scan with a collision energy of 30 and 19 eV in positive mode. Lipids were identified using LipidSearch software version 4.1, developed by Ryo Taguchi and Mitsui Knowledge Industry Co. (Tokyo, Japan). Mass spectrometry proteomics was performed on a Q-Exactive instrument after fractionation on a coupled Easy nLC 1000 nano-liquid chromatography system (Thermo Fisher Scientific), as described in our other published reports. (63) Peaks were generated using a Thermo Scientific Xcalibur (version 4.1.31.9, released 2017); proteins were identified using Proteome Discoverer 2.2 (version 2.2.0.388, released 2017). UniProt sequence database was used for the identification of proteins (downloaded July 2019). Proteome Discoverer search parameters for chymotrypsin-digested enzymes (max missed cleavage sites: 2, min peptide length: 6, max peptide length: 144), post-translational modification: deimination (+0.984 Da (R), max modification per peptide: 3), precursor mass tolerance: 10 ppm, fragment mass tolerance: 0.02 Da, signal/noise threshold for spectra: 1.5, false discovery rate (FDR): strict for PSMs 0.01, and strict for peptides 0.01. In brief, false discovery rates are calculated as follows: first, the software ascertains whether there are q-values and PEPs available for PSMs. If so, the software uses them and assigns the PSM confidences based on target FDRs for PSMs. Next, the software calculates q-values and PEPs for peptides engaging the quality algorithm. Peptide confidences are then assigned based on target FDRs for peptides. If there are no q-values and PEPs available for PSMs, the PSM confidences are set based on our defined target FDRs for PSM employing the respective search engine scores.

4.12. Liposome Extrusion and Liposome Flotation Assay

Candidate lipids were purchased at Avanti Polar Lipids Inc. (PC 16:0/22:6 cat no. 850461, PC 16:1/16:1 cat no. 850358, PI 18:0/20:4 cat no. 850144, PS 18:0/18:1 cat no. 840039, LPC 18:1 cat no. 855773, LPC 18:0 cat no. 855774, LPC 16:0 cat no. 855675, LPC 14:0 cat no. 855575, PC 16:0/20:4 cat no. 850459, PC 18:1/18:1 cat no. 999989). Lipids were resuspended in chloroform, and 30 μg of lipid was aliquoted in glass vials, desiccated in speed vacuum, and resuspended in 100 μL of PBS using sonication. Lipids were extruded using the NanoSizer MINI extruder kit (T&T Scientific, TT-030-0001) through a 100 nm NanoSizer (T&T Scientific, TT-002-0010). Extruded liposomes (100 μL) were incubated with 2.5 μg (1 μL) of isolated ex vivo mouse protein from the lipid–protein extraction of CE fractions (no-injection control, sham-injection control, and EAE groups; see Section 4.9) for 30 min at room temperature. Liposomes/micelles and protein mixture were then resuspended in 50% sucrose. Sucrose gradient was prepared by layering from bottom to top over 1 mL of liposome/micelle–protein mixture in 50% sucrose, 2 mL of 25% sucrose in PBS, and 1 mL of PBS. The gradient was ultracentrifuged at 114 000g at 4 °C for 3.5 h. Fractions were collected using a mechanical pulp to aspirate 1 mL of the bottom fraction followed by 2 mL of the middle fraction and 1 mL of the top fraction. The collected fractions were washed using PBS followed by centrifugation at 18 000g 4 °C for 25 min (Beckman Rotor GH-3.8). The supernatant was discarded, and the pellet was processed for dot immunoblot analysis, as described above.

4.13. Protein–Lipid Overlay Assay

Candidate lipids listed above were dotted on a PVDF membrane at 30 μg per dot and allowed to completely dry. The membrane was incubated for 1 h at 4 °C with 2.5 μg of isolated ex vivo mouse protein from the lipid–protein extraction of CE fractions (no-injection control, sham-injection control, and EAE groups; see Section 4.9) in PBS, followed by cross-linking at 1200 × 100 μJ for 50 s (Stratagene, La Jolla, CA, Stratalinker UV crosslinker, model 1800). The membrane was blocked, washed, and probed using primary and secondary antibodies, as described in Section 4.7.

4.14. MBP–Lipid Deimination Assay

Candidate lipids (30 μg per lipid) listed above were incubated in the presence of 10 μg of a vendor-purified MBP protein (Millipore, 13-104, Burlington, MA) in PBS and allowed to interact for 30 min at room temperature. Sequence alignment of bovine, mouse, and human MBP demonstrated high conservation between species (Figure S5). After incubation, each reaction was adjusted to a final concentration of 1.66 mM CaCl2 and incubated overnight at 37 °C in the presence of a purified peptidyl arginine deiminase (PAD, Sigma-Aldrich, P1584). Proteins were then precipitated using standard acetone precipitation protocols and analyzed as described in Section 4.7. Densitometry analysis was done using ImageJ program (National Institute of Health, Bethesda, MD). Deiminated MBP (MBP–PAD) densitometry signal was normalized to the MBP signal, and fold change was adjusted to MBP.

4.15. Circular Dichroism

Candidate lipids (PI and LPC, 180 μg per lipid) and 60 μg of a vendor-purified MBP (Millipore, 13-104, Burlington, MA) were separately resuspended in CD buffer (10 mM potassium phosphate, 50 mM ammonium sulfate). Circular dichroism (Jasco, Jasco J-815, Easton, MD) settings included measured range 250–190 nm, data pitch 1 nm, bandwidth 3 nm, scanning speed 20 nm/min, data integration time (D.I.T.) 2 s, and variable temperatures 20 and 105 °C. Baseline measurements were recorded for the CD buffer, and candidate lipids were resuspended in CD buffer before binding assays were performed. Following these measurements, a vendor-purified MBP was added to the corresponding candidate lipids or CD buffer sample and allowed to interact for 30 min at room temperature. CD spectra were recorded for MBP–lipid samples at room temperature and MBP without the presence of any lipids at two temperatures (20 and 105 °C). Baseline CD buffer and candidate lipid spectra were subtracted from their corresponding sample. Samples were then incubated overnight at 37 °C in the presence of purified PAD and 1 mM calcium phosphate. CD spectra were recoded, and baseline spectra were subtracted for each sample. CD analysis was done on the CD analysis and plotting tool (CAPITO) software (mean with SEM). (64)

4.16. Computational Modeling

There is no atomic structure available for the full-length MBP. Previous NMR studies had indicated that a part of the sequence adopts a helical structure, while additional parts are unstructured. (65) Previous studies have suggested that MBP may be an intrinsically disordered protein (IDP). (43,54,55) A check with the PrDOS disordered structure prediction server (66) indicated that MBP sequence matches the criteria for disordered proteins. Robetta server (http://new.robetta.org/) was used to obtain structure prediction (based on the Rosetta method) for MBP. (67) We used five best predictions provided by the server as the starting point and performed MD simulations, allowing MBP to explore the landscape in the presence of water, as well as in the presence of PC and LPC.

4.17. Model Preparation

AMBER’s LEaP module was used for model preparation. Starting with the five structural predictions obtained from the Robetta server (secondary structure depicted in the Supporting Information; Figure S4), five separate systems were prepared. Each structure was immersed in a rectangular box of explicit solvent (extended simple point charge (SPC/E)) (68) such that the protein coordinates are at least 10 Å from the edge of the box, and counterions were added to make the system charge neutral. AMBER’s ff14SB (69) was used for all simulations.
For PC and LPC systems, the coordinates for the protein were surrounded with 25 molecules of LPC or PC, using a model development procedure used and validated in our lab for other protein systems in nonaqueous solvents. (70) These calculations were performed by PackMol (71) such that there was one molecule of MBP for every 25 molecules of lipid, in a cube of 80 Å × 80 Å × 80 Å. The size of this box was selected after trial and error, such that the lipid molecules could easily fit but the cube is smaller in size, possibly to keep the number of atoms in the system at a minimum. Then, the system was immersed in a box of water similar to the case of water-only simulations. The charges for PC and LPC lipids were obtained using the protocol described in AMBER, and ff14SB parameters were used.

4.18. MD Simulations

Each of the 15 systems (five in water, five in PC + water, and five in LPC + water) was equilibrated separately using a protocol developed in our laboratory. (72) Equilibration removes vacuum bubbles and equilibrates the density closer to the density of water and was found to be in the range of 1.05–1.08 g/mL for the 15 model systems. The final temperature was 300 K. Following the equilibration, each of these systems was slowly heated to 305, 310, 315, and 320 K to generate conformational diversity in the MD simulation of 1 ns. These four points were separately cooled to 300 K. Therefore, for each of the 15 systems, there were five starting conformations (the equilibrated conformations at 300 K and four additional ones obtained by heating and cooling) to represent the possible conformational diversity.
For each of these 15 systems with five different starting conformations, production MD simulations were performed for a duration 100 ns under constant energy conditions (NVE ensemble). A total of 2.5 μs aggregated MD sampling was obtained for each system (with water, with water + LPC, and water + PC). NVE was used as this ensemble offers better computational stability and performance for longer MD simulations. (73) The production simulations were performed at a temperature of 300 K and a time step of 2 fs (with SHAKE applied to bonds and angles involving hydrogens). As NVE ensemble was used for production runs, these values correspond to the initial temperature at the start of simulations. A temperature adjusting thermostat was not used in simulations; over the course of 200 ns simulations, the temperature fluctuated around 300 K, with the fluctuations between 2 and 4 K, which is typical for well-equilibrated systems. A total of 100 conformational snapshots collected from each MD run was used for analysis.

4.19. Computational Modeling Analysis

The conformational snapshots from five runs with different starting conformations were collected and analyzed. Therefore, a total of 2500 conformations were analyzed for MBP in water, MBP in PC and water, and MBP in LPC and water. The goal here was to capture as much conformational diversity that MBP samples in the disordered state. The radius of gyration was calculated for the three systems separately, for each of the five structures. Secondary structure was analyzed using AMBER’s PTRAJ program. (74) Previously, we have successfully used the radius of gyration as a measure of conformational flexibility of proteins. (75)

4.20. Statistics

Data are expressed as mean ± standard error (S.E.) of three separate replicates. Statistical significance was assessed using paired (for PERG data) and unpaired (remaining data) t-tests with Prism software version 8 (GraphPad). p < 0.05 was considered statistically significant. Lipid abundance ratios were normalized to the protein concentration found in each samples used for lipidomic analysis. Deimination levels were the measurement of raw densitometry for Cit-MBP normalized to MBP densitometry after membrane stripping. Deimination fold changes reported were all compared and normalized to the endogenous levels of deimination in the MBP sample.

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsomega.0c01590.

  • Confirmation of the EAE model, variation of hydrocarbon chain, computational modeling of MBP, and five model structures of MBP and MBP sequence alignment (Figures S1–S5); mass spectrometry results for arginine residues and confirmation of MBP presence in gel-excised band (Table S1); and computer simulations of MBP (Movies S1–S3) (ZIP)

Accession Codes

MBP (UniProt accession: P02686-5 and P02687). Mass spectrometry data (PRIDE archive accession: PXD015494).

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Author
    • Sanjoy K. Bhattacharya - Department of Ophthalmology, Bascom Palmer Eye Institute, University of Miami, Miami, Florida 33136, United StatesNeuroscience Graduate Program, University of Miami, Miami, Florida 33136, United StatesOrcidhttp://orcid.org/0000-0003-3759-647X Email: [email protected]
  • Authors
    • Anddre Osmar Valdivia - Department of Ophthalmology, Bascom Palmer Eye Institute, University of Miami, Miami, Florida 33136, United StatesNeuroscience Graduate Program, University of Miami, Miami, Florida 33136, United StatesPresent Address: Department of Ophthalmology, Bascom Palmer Eye Institute, University of Miami, 1638 NW 10th Avenue, Room 706A, Miami Florida 33136, United StatesOrcidhttp://orcid.org/0000-0002-8425-9102
    • Pratul K. Agarwal - Department of Biochemistry & Cell and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, United StatesDepartment of Physiological Sciences  and  High Performance Computing Center and, Oklahoma State University, Stillwater, 106 Math Sciences, Stillwater, Oklahoma 74078-1010, United StatesOrcidhttp://orcid.org/0000-0002-3848-9492
  • Author Contributions

    A.O.V. and S.K.B., conceptualization; A.O.V and P.K.A., data acquisition and curation; A.O.V., formal analysis; A.O.V. and S.K.B., supervision; S.K.B., funding acquisition; A.O.V., validation; S.K.B., principal investigator; A.O.V. and S.K.B., methodology; A.O.V., writing-original draft; A.O.V. and S.K.B., project administration; and A.O.V. and S.K.B., writing: review and editing.

  • Funding

    This work was supported by an unrestricted grant to the University of Miami from Research to Prevent Blindness (RPB), Department of Defense Grant WHX81-16-0715 and NIH Grant EY027257, EY14801, and NIH Grant GM105978 (to P.K.A.). This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562. XSEDE computing allocation were awarded to P.K.A. (MCB-180199 and MCB-190044).

  • Notes
    The authors declare no competing financial interest.

Acknowledgments

Click to copy section linkSection link copied!

The authors would like to thank Dr. Jae Lee for his valuable insight into the preparation of the manuscript, Drs. Roberta Brambilla and Haritha Desu for their guidance in the generation of the EAE mouse model, and Dr. David Broyles for his guidance during circular dichroism experiments, EY14801and NIH Grant GM105978 (to PKA). This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562. XSEDE computing allocation were awarded to P.K.A. (MCB-180199 and MCB-190044).

Abbreviations

Click to copy section linkSection link copied!

MS

multiple sclerosis

MBP

myelin basic protein

PE

phosphatidylethanolamines

PC

phosphatidylcholines

PI

phosphatidylinositols

PS

phosphatidylserines

EAE

experimental autoimmune encephalomyelitis

MOG

myelin oligodendrocyte glycoprotein

PERG

pattern electroretinogram

FERG

flash electroretinogram

BHT

butylated hydroxytoluene

MTBE

methyl tert-butyl ether

OCT

optical coherence tomography

CE

capillary electrophoresis

PAD

peptidyl arginine deiminase

LFA

liposome floatation assay

PLOA

protein–lipid overlay assay

CD

circular dichroism

MD

molecular dynamic

IDP

intrinsic-disordered protein

CNS

central nervous system

References

Click to copy section linkSection link copied!

This article references 75 other publications.

  1. 1
    Chen, L.; Gordon, L. K. Ocular manifestations of multiple sclerosis. Curr. Opin. Ophthalmol. 2005, 16, 315320,  DOI: 10.1097/01.icu.0000179804.49842.e2
  2. 2
    Dendrou, C. A.; Fugger, L.; Friese, M. A. Immunopathology of multiple sclerosis. Nat. Rev. Immunol. 2015, 15, 545558,  DOI: 10.1038/nri3871
  3. 3
    Nylander, A.; Hafler, D. A. Multiple sclerosis. J. Clin. Invest. 2012, 122, 1180,  DOI: 10.1172/JCI58649
  4. 4
    Lassmann, H. What drives disease in multiple sclerosis: Inflammation or neurodegeneration?. Clin. Exp. Neuroimmunol. 2010, 1, 211,  DOI: 10.1111/j.1759-1961.2009.00003.x
  5. 5
    Bielekova, B.; Sung, M.-H.; Kadom, N.; Simon, R.; McFarland, H.; Martin, R. Expansion and Functional Relevance of High-Avidity Myelin-Specific CD4+ T Cells in Multiple Sclerosis. J. Immunol. 2004, 172, 38933904,  DOI: 10.4049/jimmunol.172.6.3893
  6. 6
    Moscarello, M. A.; Wood, D. D.; Boulias, C.; Ackerley, C. Myelin in multiple sclerosis is developmentally immature. J. Clin. Invest. 1994, 94, 146154,  DOI: 10.1172/JCI117300
  7. 7
    Kursula, P. Structural properties of proteins specific to the myelin sheath. Amino Acids 2008, 34, 175185,  DOI: 10.1007/s00726-006-0479-7
  8. 8
    Hartline, D. K.; Colman, D. R. Rapid conduction and the evolution of giant axons and myelinated fibers. Curr. Biol. 2007, 17, R29R35,  DOI: 10.1016/j.cub.2006.11.042
  9. 9
    Valdivia, A. O.; Farr, V.; Bhattacharya, S. K. A novel myelin basic protein transcript variant in the murine central nervous system. Mol. Biol. Rep. 2019, 46, 25472553,  DOI: 10.1007/s11033-019-04635-8
  10. 10
    Musse, A. A.; Boggs, J. M.; Harauz, G. Deimination of membrane-bound myelin basic protein in multiple sclerosis exposes an immunodominant epitope. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 44224427,  DOI: 10.1073/pnas.0509158103
  11. 11
    Kidd, B. A.; Ho, P. P.; Sharpe, O.; Zhao, X.; Tomooka, B. H.; Kanter, J. L.; Steinman, L.; Robinson, W. H. Epitope spreading to citrullinated antigens in mouse models of autoimmune arthritis and demyelination. Arthritis Res. Ther. 2008, 10, R119  DOI: 10.1186/ar2523
  12. 12
    Cuzner, M. L.; Davison, A. N. Changes in cerebral lysosomal enzyme activity and lipids in multiple sclerosis. J. Neurol. Sci. 1973, 19, 2936,  DOI: 10.1016/0022-510X(73)90053-1
  13. 13
    Einstein, E. R.; Csejtey, J.; Dalal, K. B.; Adams, C. W.; Bayliss, O. B.; Hallpike, J. F. Proteolytic activity and basic protein loss in and around multiple sclerosis plaques: combined biochemical and histochemical observations. J. Neurochem. 1972, 19, 653662,  DOI: 10.1111/j.1471-4159.1972.tb01382.x
  14. 14
    Richards, P. T.; Cuzner, M. L. Proteolytic activity in CSF. Adv. Exp. Med. Biol. 1978, 100, 521527,  DOI: 10.1007/978-1-4684-2514-7_38
  15. 15
    Dreyton, C. J.; Knuckley, B.; Jones, J. E.; Lewallen, D. M.; Thompson, P. R. Mechanistic studies of protein arginine deiminase 2: evidence for a substrate-assisted mechanism. Biochemistry 2014, 53, 4426,  DOI: 10.1021/bi500554b
  16. 16
    Vossenaar, E. R.; Zendman, A. J.; van Venrooij, W. J.; Pruijn, G. J. PAD, a growing family of citrullinating enzymes: genes, features and involvement in disease. BioEssays 2003, 25, 11061118,  DOI: 10.1002/bies.10357
  17. 17
    Tarcsa, E.; Marekov, L. N.; Mei, G.; Melino, G.; Lee, S. C.; Steinert, P. M. Protein unfolding by peptidylarginine deiminase. Substrate specificity and structural relationships of the natural substrates trichohyalin and filaggrin. J. Biol. Chem. 1996, 271, 30709,  DOI: 10.1074/jbc.271.48.30709
  18. 18
    Takahara, H. A History of Deimination Research in Japan: The Founding Fathers. In Protein Deimination in Human Health and Disease; Nicholas, A. P.; Bhattacharya, S. K.; Thompson, P. R., Eds.; Springer International Publishing: Cham, 2017; pp 110.
  19. 19
    Cao, L.; Goodin, R.; Wood, D.; Moscarello, M. A.; Whitaker, J. N. Rapid release and unusual stability of immunodominant peptide 45-89 from citrullinated myelin basic protein. Biochemistry 1999, 38, 61576163,  DOI: 10.1021/bi982960s
  20. 20
    Whitaker, J. N. Myelin encephalitogenic protein fragments in cerebrospinal fluid of persons with multiple sclerosis. Neurology 1977, 27, 911920,  DOI: 10.1212/WNL.27.10.911
  21. 21
    Whitaker, J. R.; Granum, P. E. An absolute method for protein determination based on difference in absorbance at 235 and 280 nm. Anal. Biochem. 1980, 109, 156159,  DOI: 10.1016/0003-2697(80)90024-X
  22. 22
    Pritzker, L. B.; Joshi, S.; Gowan, J. J.; Moscarello, M. A.; Harauz, G. Deimination of myelin basic protein. 1. Effect of deimination of arginyl residues of myelin basic protein on its structure and susceptibility to digestion by cathepsin D. Biochemistry 2000, 39, 53745381,  DOI: 10.1021/bi9925569
  23. 23
    Wood, D. D.; Bilbao, J. M.; O’Connors, P.; Moscarello, M. A. Acute multiple sclerosis (Marburg type) is associated with developmentally immature myelin basic protein. Ann. Neurol. 1996, 40, 1824,  DOI: 10.1002/ana.410400106
  24. 24
    Mastronardi, F. G.; Wood, D. D.; Mei, J.; Raijmakers, R.; Tseveleki, V.; Dosch, H. M.; Probert, L.; Casaccia-Bonnefil, P.; Moscarello, M. A. Increased citrullination of histone H3 in multiple sclerosis brain and animal models of demyelination: a role for tumor necrosis factor-induced peptidylarginine deiminase 4 translocation. J. Neurosci. 2006, 26, 1138711396,  DOI: 10.1523/JNEUROSCI.3349-06.2006
  25. 25
    Nicholas, A. P.; Whitaker, J. N. Preparation of a monoclonal antibody to citrullinated epitopes: its characterization and some applications to immunohistochemistry in human brain. Glia 2002, 37, 328336,  DOI: 10.1002/glia.10039
  26. 26
    Raijmakers, R.; Vogelzangs, J.; Croxford, J. L.; Wesseling, P.; van Venrooij, W. J.; Pruijn, G. J. Citrullination of central nervous system proteins during the development of experimental autoimmune encephalomyelitis. J. Comp. Neurol. 2005, 486, 243253,  DOI: 10.1002/cne.20529
  27. 27
    Cao, L.; Sun, D.; Whitaker, J. N. Citrullinated myelin basic protein induces experimental autoimmune encephalomyelitis in Lewis rats through a diverse T cell repertoire. J. Neuroimmunol. 1998, 88, 2129,  DOI: 10.1016/S0165-5728(98)00063-0
  28. 28
    Carrillo-Vico, A.; Leech, M. D.; Anderton, S. M. Contribution of myelin autoantigen citrullination to T cell autoaggression in the central nervous system. J. Immunol. 2010, 184, 28392846,  DOI: 10.4049/jimmunol.0903639
  29. 29
    Quinn, T. A.; Dutt, M.; Shindler, K. S. Optic neuritis and retinal ganglion cell loss in a chronic murine model of multiple sclerosis. Front. Neurol. 2011, 2, 50  DOI: 10.3389/fneur.2011.00050
  30. 30
    Beniac, D. R.; Luckevich, M. D.; Czarnota, G. J.; Tompkins, T. A.; Ridsdale, R. A.; Ottensmeyer, F. P.; Moscarello, M. A.; Harauz, G. Three-dimensional structure of myelin basic protein. I. Reconstruction via angular reconstitution of randomly oriented single particles. J. Biol. Chem. 1997, 272, 42614268,  DOI: 10.1074/jbc.272.7.4261
  31. 31
    Ridsdale, R. A.; Beniac, D. R.; Tompkins, T. A.; Moscarello, M. A.; Harauz, G. Three-dimensional structure of myelin basic protein. II. Molecular modeling and considerations of predicted structures in multiple sclerosis. J. Biol. Chem. 1997, 272, 42694275,  DOI: 10.1074/jbc.272.7.4269
  32. 32
    Haas, H.; Oliveira, C. L.; Torriani, I. L.; Polverini, E.; Fasano, A.; Carlone, G.; Cavatorta, P.; Riccio, P. Small angle x-ray scattering from lipid-bound myelin basic protein in solution. Biophys. J. 2004, 86, 455460,  DOI: 10.1016/S0006-3495(04)74122-3
  33. 33
    Jahn, O.; Tenzer, S.; Werner, H. Myelin Proteomics: Molecular Anatomy of an Insulating Sheath. Mol. Neurobiol. 2009, 40, 5572,  DOI: 10.1007/s12035-009-8071-2
  34. 34
    Fewster, M. E.; Hirono, H.; Mead, J. F. Lipid composition of myelin in multiple sclerosis. J. Neurol. 1976, 213, 119131,  DOI: 10.1007/BF00313273
  35. 35
    Wilson, R.; Tocher, D. R. Lipid and fatty acid composition is altered in plaque tissue from multiple sclerosis brain compared with normal brain white matter. Lipids 1991, 26, 915,  DOI: 10.1007/BF02544017
  36. 36
    Alling, C.; Vanier, M. T.; Svennerholm, L. Lipid alterations in apparently normal white matter in multiple sclerosis. Brain Res. 1971, 35, 325336,  DOI: 10.1016/0006-8993(71)90478-1
  37. 37
    Woelk, H.; Borri, P. Lipid and fatty acid composition of myelin purified from normal and MS brains. Eur. Neurol. 1973, 10, 250260,  DOI: 10.1159/000114281
  38. 38
    Gerstl, B.; Eng, L. F.; Tavaststjerna, M.; Smith, J. K.; Kruse, S. L. Lipids and proteins in multiple sclerosis white matter. J. Neurochem. 1970, 17, 677689,  DOI: 10.1111/j.1471-4159.1970.tb00547.x
  39. 39
    Arnetoli, G.; Pazzagli, A.; Amaducci, L. Fatty acid and aldehyde changes in choline- and ethanolamine-containing phospholipids in the white matter of multiple sclerosis brains. J. Neurochem. 1969, 16, 461463,  DOI: 10.1111/j.1471-4159.1969.tb10387.x
  40. 40
    Göpfert, E.; Pytlik, S.; Debuch, H. 2′,3′-Cyclic nucleotide 3′-phosphohydrolase and lipids of myelin from multiple sclerosis and normal brains. J. Neurochem. 1980, 34, 732739,  DOI: 10.1111/j.1471-4159.1980.tb11205.x
  41. 41
    Neu, I.; Woelk, H. Investigations of the lipid metabolism of the white matter in multiple sclerosis Changes in glycero-phosphatides and lipid-splitting enzymes. Neurochem. Res. 1982, 7, 727735,  DOI: 10.1007/BF00965525
  42. 42
    Del Boccio, P.; Pieragostino, D.; Di Ioia, M.; Petrucci, F.; Lugaresi, A.; De Luca, G.; Gambi, D.; Onofrj, M.; Di Ilio, C.; Sacchetta, P.; Urbani, A. Lipidomic investigations for the characterization of circulating serum lipids in multiple sclerosis. J. Proteomics 2011, 74, 28262836,  DOI: 10.1016/j.jprot.2011.06.023
  43. 43
    Boggs, J. M. Myelin Basic Protein; Nova Science Publishers, Incorporated: New York, 2008.
  44. 44
    Procaccini, C.; De Rosa, V.; Pucino, V.; Formisano, L.; Matarese, G. Animal models of Multiple Sclerosis. Eur. J. Pharmacol. 2015, 759, 182191,  DOI: 10.1016/j.ejphar.2015.03.042
  45. 45
    Niroomand, H.; Pamu, R.; Mukherjee, D.; Khomami, B. Microenvironment alterations enhance photocurrents from photosystem I confined in supported lipid bilayers. J. Mater. Chem. A 2018, 6, 1228112290,  DOI: 10.1039/C8TA00898A
  46. 46
    Israelachvili, J. N.; Mitchell, D. J.; Ninham, B. W. Theory of self-assembly of lipid bilayers and vesicles. Biochim. Biophys. Acta, Biomembr. 1977, 470, 185201,  DOI: 10.1016/0005-2736(77)90099-2
  47. 47
    Israelachvili, J. N. Intermolecular and Surface Forces; Academic Press, 2015.
  48. 48
    Wood, D. D.; Moscarello, M. A. The isolation, characterization, and lipid-aggregating properties of a citrulline containing myelin basic protein. J. Biol. Chem. 1989, 264, 51215127
  49. 49
    Harauz, G.; Musse, A. A Tale of Two Citrullines—Structural and Functional Aspects of Myelin Basic Protein Deimination in Health and Disease. Neurochem. Res. 2007, 32, 137158,  DOI: 10.1007/s11064-006-9108-9
  50. 50
    Matsumoto, T.; Kobayashi, T.; Kamata, K. Role of lysophosphatidylcholine (LPC) in atherosclerosis. Curr. Med. Chem. 2007, 14, 32093220,  DOI: 10.2174/092986707782793899
  51. 51
    Muramatsu, R.; Kuroda, M.; Matoba, K.; Lin, H.; Takahashi, C.; Koyama, Y.; Yamashita, T. Prostacyclin prevents pericyte loss and demyelination induced by lysophosphatidylcholine in the central nervous system. J. Biol. Chem. 2015, 290, 1151511525,  DOI: 10.1074/jbc.M114.587253
  52. 52
    Plemel, J. R.; Michaels, N. J.; Weishaupt, N.; Caprariello, A. V.; Keough, M. B.; Rogers, J. A.; Yukseloglu, A.; Lim, J.; Patel, V. V.; Rawji, K. S.; Jensen, S. K.; Teo, W.; Heyne, B.; Whitehead, S. N.; Stys, P. K.; Yong, V. W. Mechanisms of lysophosphatidylcholine-induced demyelination: A primary lipid disrupting myelinopathy. Glia 2018, 66, 327347,  DOI: 10.1002/glia.23245
  53. 53
    Birgbauer, E.; Rao, T. S.; Webb, M. Lysolecithin induces demyelination in vitro in a cerebellar slice culture system. J. Neurosci. Res. 2004, 78, 157166,  DOI: 10.1002/jnr.20248
  54. 54
    Boggs, J. M. Myelin basic protein: a multifunctional protein. Cell. Mol. Life Sci. 2006, 63, 19451961,  DOI: 10.1007/s00018-006-6094-7
  55. 55
    Harauz, G.; Ishiyama, N.; Hill, C. M.; Bates, I. R.; Libich, D. S.; Fares, C. Myelin basic protein-diverse conformational states of an intrinsically unstructured protein and its roles in myelin assembly and multiple sclerosis. Micron 2004, 35, 503542,  DOI: 10.1016/j.micron.2004.04.005
  56. 56
    D’Souza, C. A.; Wood, D. D.; She, Y.-M.; Moscarello, M. A. Autocatalytic cleavage of myelin basic protein: An alternative to molecular mimicry. Biochemistry 2005, 44, 1290512913,  DOI: 10.1021/bi051152f
  57. 57
    Zhou, S. R.; Moscarello, M. A.; Whitaker, J. N. The effects of citrullination or variable amino-terminus acylation on the encephalitogenicity of human myelin basic protein in the PL/J mouse. J. Neuroimmunol. 1995, 62, 147152,  DOI: 10.1016/0165-5728(95)00112-3
  58. 58
    Madsen, P. M.; Motti, D.; Karmally, S.; Szymkowski, D. E.; Lambertsen, K. L.; Bethea, J. R.; Brambilla, R. Oligodendroglial TNFR2 Mediates Membrane TNF-Dependent Repair in Experimental Autoimmune Encephalomyelitis by Promoting Oligodendrocyte Differentiation and Remyelination. J. Neurosci. 2016, 36, 51285143,  DOI: 10.1523/JNEUROSCI.0211-16.2016
  59. 59
    Porciatti, V.; Saleh, M.; Nagaraju, M. The pattern electroretinogram as a tool to monitor progressive retinal ganglion cell dysfunction in the DBA/2J mouse model of glaucoma. Invest. Ophthalmol. Visual Sci. 2007, 48, 745751,  DOI: 10.1167/iovs.06-0733
  60. 60
    Menon, K.; Rasband, M. N.; Taylor, C. M.; Brophy, P.; Bansal, R.; Pfeiffer, S. E. The myelin–axolemmal complex: biochemical dissection and the role of galactosphingolipids. J. Neurochem. 2003, 87, 9951009,  DOI: 10.1046/j.1471-4159.2003.02075.x
  61. 61
    Valdivia, A. O.; Myer, C.; Suarez, M. F.; Bhattacharya, S. K. Protein–Lipid Complex Separation Utilizing a Capillary Electrophoresis System. In Metabolomics: Methods and Protocols; Bhattacharya, S. K., Ed.; Springer: New York, 1996; pp 95100.
  62. 62
    Chauhan, M. Z.; Valencia, A. K.; Piqueras, M. C.; Enriquez-Algeciras, M.; Bhattacharya, S. K. Optic Nerve Lipidomics Reveal Impaired Glucosylsphingosine Lipids Pathway in Glaucoma. Invest. Ophthalmol. Visual Sci. 2019, 60, 17891798,  DOI: 10.1167/iovs.18-25802
  63. 63
    Travers, T. S.; Harlow, L.; Rosas, I. O.; Gochuico, B. R.; Mikuls, T. R.; Bhattacharya, S. K.; Camacho, C. J.; Ascherman, D. P. Extensive Citrullination Promotes Immunogenicity of HSP90 through Protein Unfolding and Exposure of Cryptic Epitopes. J. Immunol. 2016, 197, 19261936,  DOI: 10.4049/jimmunol.1600162
  64. 64
    Wiedemann, C.; Bellstedt, P.; Gorlach, M. CAPITO--a web server-based analysis and plotting tool for circular dichroism data. Bioinformatics 2013, 29, 17501757,  DOI: 10.1093/bioinformatics/btt278
  65. 65
    Ahmed, M. A.; De Avila, M.; Polverini, E.; Bessonov, K.; Bamm, V. V.; Harauz, G. Solution nuclear magnetic resonance structure and molecular dynamics simulations of a murine 18.5 kDa myelin basic protein segment (S72-S107) in association with dodecylphosphocholine micelles. Biochemistry 2012, 51, 74757487,  DOI: 10.1021/bi300998x
  66. 66
    Ishida, T.; Kinoshita, K. PrDOS: prediction of disordered protein regions from amino acid sequence. Nucleic Acids Res. 2007, 35, W460W464,  DOI: 10.1093/nar/gkm363
  67. 67
    Song, Y.; DiMaio, F.; Wang, R. Y.; Kim, D.; Miles, C.; Brunette, T.; Thompson, J.; Baker, D. High-resolution comparative modeling with RosettaCM. Structure 2013, 21, 17351742,  DOI: 10.1016/j.str.2013.08.005
  68. 68
    Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The missing term in effective pair potentials. J. Phys. Chem. A 1987, 91, 62696271,  DOI: 10.1021/j100308a038
  69. 69
    Maier, J. A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K. E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 36963713,  DOI: 10.1021/acs.jctc.5b00255
  70. 70
    Duff, M. R., Jr.; Borreguero, J. M.; Cuneo, M. J.; Ramanathan, A.; He, J.; Kamath, G.; Chennubhotla, S. C.; Meilleur, F.; Howell, E. E.; Herwig, K. W.; Myles, D. A. A.; Agarwal, P. K. Modulating Enzyme Activity by Altering Protein Dynamics with Solvent. Biochemistry 2018, 57, 42634275,  DOI: 10.1021/acs.biochem.8b00424
  71. 71
    Martínez, L.; Andrade, R.; Birgin, E. G.; Martinez, J. M. PACKMOL: a package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 21572164,  DOI: 10.1002/jcc.21224
  72. 72
    Agarwal, P. K. Cis/trans isomerization in HIV-1 capsid protein catalyzed by cyclophilin A: insights from computational and theoretical studies. Proteins 2004, 56, 449463,  DOI: 10.1002/prot.20135
  73. 73
    Beck, D. A.; Daggett, V. Methods for molecular dynamics simulations of protein folding/unfolding in solution. Methods 2004, 34, 112120,  DOI: 10.1016/j.ymeth.2004.03.008
  74. 74
    Roe, D. R.; Cheatham, T. E., 3rd PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 30843095,  DOI: 10.1021/ct400341p
  75. 75
    Shukla, S.; Bafna, K.; Gullett, C.; Myles, D. A. A.; Agarwal, P. K.; Cuneo, M. J. Differential Substrate Recognition by Maltose Binding Proteins Influenced by Structure and Dynamics. Biochemistry 2018, 57, 58645876,  DOI: 10.1021/acs.biochem.8b00783

Cited By

Click to copy section linkSection link copied!

This article is cited by 9 publications.

  1. Anwaar S. Chaudary, Yanglin Guo, Yuri N. Utkin, Maryam Barancheshmeh, Ruben K. Dagda, Edward S. Gasanoff. Sphingomyelin Inhibits Hydrolytic Activity of Heterodimeric PLA2 in Model Myelin Membranes: Pharmacological Relevance. The Journal of Membrane Biology 2024, 133 https://doi.org/10.1007/s00232-024-00327-y
  2. Majid Dejbakht, Morteza Akhzari, Sajad Jalili, Fouziyeh Faraji, Mahdi Barazesh. Multiple Sclerosis: New Insights into Molecular Pathogenesis and Novel Platforms for Disease Treatment. Current Drug Research Reviews 2024, 16 (2) , 175-197. https://doi.org/10.2174/2589977516666230915103730
  3. Faith Christine Harvey, Anddre Osmar Valdivia, Colin Hayter, Sanjoy K. Bhattacharya. Isobaric Incorporation of C13-Histidine for the Assessment of Remyelination. 2023, 169-175. https://doi.org/10.1007/978-1-0716-2699-3_17
  4. Anddre Osmar Valdivia, Carolyn Cray, Sanjoy K. Bhattacharya. Capillary Electrophoresis Plasma Fractionation for Lipids Analysis. 2023, 175-181. https://doi.org/10.1007/978-1-0716-2966-6_16
  5. Irina A. Ishina, Maria Y. Zakharova, Inna N. Kurbatskaia, Azad E. Mamedov, Alexey A. Belogurov, Alexander G. Gabibov. MHC Class II Presentation in Autoimmunity. Cells 2023, 12 (2) , 314. https://doi.org/10.3390/cells12020314
  6. Mariam F. Farid, Yara S. Abouelela, Noha A. E. Yasin, Mohamed R. Mousa, Marwa A. Ibrahim, Abdelbary Prince, Hamdy Rizk. A novel cell-free intrathecal approach with PRP for the treatment of spinal cord multiple sclerosis in cats. Inflammation and Regeneration 2022, 42 (1) https://doi.org/10.1186/s41232-022-00230-w
  7. Juan Sebastian Cruz-Méndez, María Paula Herrera-Sánchez, Ángel Enrique Céspedes-Rubio, Iang Schroniltgen Rondón-Barragán. Molecular characterization of myelin basic protein a (mbpa) gene from red-bellied pacu (Piaractus brachypomus). Journal of Genetic Engineering and Biotechnology 2022, 20 (1) , 8. https://doi.org/10.1186/s43141-022-00296-6
  8. Magdalena Kołodziej, Karolina Chrabąszcz, Ewa Pięta, Natalia Piergies, Julia Rudnicka-Czerwiec, Halina Bartosik-Psujek, Czesława Paluszkiewicz, Marian Cholewa, Wojciech M. Kwiatek. Spectral signature of multiple sclerosis. Preliminary studies of blood fraction by ATR FTIR technique. Biochemical and Biophysical Research Communications 2022, 593 , 40-45. https://doi.org/10.1016/j.bbrc.2022.01.046
  9. Anddre Osmar Valdivia, Sanjoy K. Bhattacharya. Lyso-Lipid-Induced Oligodendrocyte Maturation Underlies Restoration of Optic Nerve Function. eneuro 2022, 9 (1) , ENEURO.0429-21.2022. https://doi.org/10.1523/ENEURO.0429-21.2022

ACS Omega

Cite this: ACS Omega 2020, 5, 25, 15454–15467
Click to copy citationCitation copied!
https://doi.org/10.1021/acsomega.0c01590
Published June 16, 2020

Copyright © 2020 American Chemical Society. This publication is licensed under these Terms of Use.

Article Views

2592

Altmetric

-

Citations

Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Myelin basic protein in EAE is deiminated. (A) Representative sucrose density ultracentrifugation fractions of mouse brains probing for myelin basic protein (Abcam-ab7349). *Sucrose fraction 32 (SF 32) from each group was used for further analysis (coomassie, Western blot, mass spectrometry, and capillary electrophoresis (CE)). EAE, n = 15; sham-injection control, n = 10; and no-injection control, n = 10. Each immunoblot was repeated three times. (B) Coomassie blue stain and Western blot for sucrose fraction 32. Bands were excised from coomassie gel and analyzed by mass spectrometry. Western blot probed for levels of deimination (citrulline, Millipore-MABS54887). (C) Mass spectrometry analysis of excised bands from coomassie blue stain gel. PSM = peptide-spectrum match.

    Figure 2

    Figure 2. Myelin basic protein isolated from EAE has different phospholipid compositions than controls. (A) Representative capillary electrophoresis (CE) chromatogram of sucrose fraction 32 (SF 32). Arrow points toward the peak collected for protein and phospholipid analysis. (B) Dot immunoblot for MBP in CE fractions, 10 min before the peak, peak (as noted by the arrow), and 10 min after the peak. MBP was probed with an anti-MBP antibody (Abcam, ab7349). (C) Phospholipid analysis of CE peak fractions, three replicates, and bars correspond to standard error mean (SEM) (*p < 0.05).

    Figure 3

    Figure 3. Confirmation of MBP–phospholipid complexation. (A) Structure of candidate lipids: PC 16:0/22:6, 1-((4Z,7Z,10Z,13Z,16Z,19Z)-docosahexaenoyl)-2-hexadecanoyl-sn-glycero-3-phosphocholine, LMGP01011116; PC 16:1/16:1, 1,2-di-((9Z)-hexadecanoyl)-sn-glycero-3-phosphocholine, LMGP01010684; PI 18:0/20:4, 1-octadecanoyl-2-((5Z,8Z,11Z,14Z)-eicosatetraenoyl)-sn-glycero-3-phospho-(1′-myo-inositol), LMGP06010010; PS 18:0/18:1, 1-octadecanoyl-2-((9Z)-octadecanoyl)-sn-glycero-3-phosphoserine, LMGP03010025; LPC 18:1, 2-((9Z)-octadecanoyl)-sn-glycero-3-phosphocholine, LMGP01050082. (B) Liposome flotation assay (LFA) of candidate phospholipids demonstrating MBP complexation with phospholipids. B, bottom fraction; M, middle fraction; and T, top fraction (MBP, Abcam, ab7349). (C) Protein–lipid overlay assay (PLOA) of candidate phospholipids demonstrating MBP complexation with phospholipids (MBP, Abcam, ab7349).

    Figure 4

    Figure 4. Complexation of LPC 18:1 with MBP is competitive for MBP hyper-deimination. (A) Dot blot analysis of deimination reaction for MBP (MBP, Abcam-ab7349) and deiminated MBP (citrulline, Millipore-MABS54887) in the presence of candidate phospholipids. PAD, peptidyl arginine deiminase. (B) Deimination fold change analysis of Western blot, three replicates, and densitometry using ImageJ was normalized to MBP (***p < 0.0005, comparison to MBP–PAD). (C) Representative circular dichroism (CD) spectra of deimination reaction for MBP in the presence of candidate phospholipids, three replicates per spectra. (D) CD spectra analysis of secondary structure probability, utilizing CD analysis and plotting tool (CAPITO) software (mean with SEM).

    Figure 5

    Figure 5. LPC 18:1 competes with hyper-deimination of MBP arginine residues and has differential complexation with MBP. (A) Mass spectrometry analysis of MBP arginine residues that are being deiminated in in vitro deimination assay (Millipore, 13-104, Burlington, MA). The deiminated Myelin basic protein (UniProt accession: P02687) peptides were identified by mass spectrometry (see Table 3). Gray highlight, deiminated arginine; red bold square, arginine not deiminated in LPC 18:1 reaction (residues 63, 112, and 129). (B) Protein–lipid overlay assay of LPC 18:1, LPC 14:0, LPC 16:0, and LPC 18:0 MBP was probed with an anti-MBP antibody (Abcam, ab7349).

    Figure 6

    Figure 6. Change in MBP conformations and solvent accessibility of arginine residues. Representative structures are shown from simulations only in water, water with LPC, and water with PC. Helices are shown in red, β-strands in yellow, three N-terminal residues are colored green, and three C-terminal residues are colored cyan. Arginine residues are shown as green/blue sticks, while the lipid is shown as magenta/blue/orange sticks. Water was included in the simulations but omitted here for visualization purposes. In the case of water-only, the arginine residues are exposed on the surface; the presence of LPC protects a number of these residues, and in the presence of PC, MBP adopts an extended conformation, exposing several arginine residues to the solvent.

    Figure 7

    Figure 7. Hypothesis model of MBP–phospholipid complexation. (A) Working hypothesis for membrane disruption, incorporating key events associated with unstable myelin. The order of events has yet to be determined. (B) Model of phospholipid (PC 16:1/16:1, PI 18:0/20:4, PS 18:0/18:1, LPC 18:0, LPC 16:0, and LPC 14:0) competitive effects against PAD-mediated hyper-deimination and potential cascade of events. (C) Model of LPC 18:1 competitive effects against PAD-mediated hyper-deimination. Δ, “change in”; MBP, myelin basic protein; and PAD, peptidyl arginine deiminase. Open and C-shaped conformations are consistent with the previously published reports.

  • References


    This article references 75 other publications.

    1. 1
      Chen, L.; Gordon, L. K. Ocular manifestations of multiple sclerosis. Curr. Opin. Ophthalmol. 2005, 16, 315320,  DOI: 10.1097/01.icu.0000179804.49842.e2
    2. 2
      Dendrou, C. A.; Fugger, L.; Friese, M. A. Immunopathology of multiple sclerosis. Nat. Rev. Immunol. 2015, 15, 545558,  DOI: 10.1038/nri3871
    3. 3
      Nylander, A.; Hafler, D. A. Multiple sclerosis. J. Clin. Invest. 2012, 122, 1180,  DOI: 10.1172/JCI58649
    4. 4
      Lassmann, H. What drives disease in multiple sclerosis: Inflammation or neurodegeneration?. Clin. Exp. Neuroimmunol. 2010, 1, 211,  DOI: 10.1111/j.1759-1961.2009.00003.x
    5. 5
      Bielekova, B.; Sung, M.-H.; Kadom, N.; Simon, R.; McFarland, H.; Martin, R. Expansion and Functional Relevance of High-Avidity Myelin-Specific CD4+ T Cells in Multiple Sclerosis. J. Immunol. 2004, 172, 38933904,  DOI: 10.4049/jimmunol.172.6.3893
    6. 6
      Moscarello, M. A.; Wood, D. D.; Boulias, C.; Ackerley, C. Myelin in multiple sclerosis is developmentally immature. J. Clin. Invest. 1994, 94, 146154,  DOI: 10.1172/JCI117300
    7. 7
      Kursula, P. Structural properties of proteins specific to the myelin sheath. Amino Acids 2008, 34, 175185,  DOI: 10.1007/s00726-006-0479-7
    8. 8
      Hartline, D. K.; Colman, D. R. Rapid conduction and the evolution of giant axons and myelinated fibers. Curr. Biol. 2007, 17, R29R35,  DOI: 10.1016/j.cub.2006.11.042
    9. 9
      Valdivia, A. O.; Farr, V.; Bhattacharya, S. K. A novel myelin basic protein transcript variant in the murine central nervous system. Mol. Biol. Rep. 2019, 46, 25472553,  DOI: 10.1007/s11033-019-04635-8
    10. 10
      Musse, A. A.; Boggs, J. M.; Harauz, G. Deimination of membrane-bound myelin basic protein in multiple sclerosis exposes an immunodominant epitope. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 44224427,  DOI: 10.1073/pnas.0509158103
    11. 11
      Kidd, B. A.; Ho, P. P.; Sharpe, O.; Zhao, X.; Tomooka, B. H.; Kanter, J. L.; Steinman, L.; Robinson, W. H. Epitope spreading to citrullinated antigens in mouse models of autoimmune arthritis and demyelination. Arthritis Res. Ther. 2008, 10, R119  DOI: 10.1186/ar2523
    12. 12
      Cuzner, M. L.; Davison, A. N. Changes in cerebral lysosomal enzyme activity and lipids in multiple sclerosis. J. Neurol. Sci. 1973, 19, 2936,  DOI: 10.1016/0022-510X(73)90053-1
    13. 13
      Einstein, E. R.; Csejtey, J.; Dalal, K. B.; Adams, C. W.; Bayliss, O. B.; Hallpike, J. F. Proteolytic activity and basic protein loss in and around multiple sclerosis plaques: combined biochemical and histochemical observations. J. Neurochem. 1972, 19, 653662,  DOI: 10.1111/j.1471-4159.1972.tb01382.x
    14. 14
      Richards, P. T.; Cuzner, M. L. Proteolytic activity in CSF. Adv. Exp. Med. Biol. 1978, 100, 521527,  DOI: 10.1007/978-1-4684-2514-7_38
    15. 15
      Dreyton, C. J.; Knuckley, B.; Jones, J. E.; Lewallen, D. M.; Thompson, P. R. Mechanistic studies of protein arginine deiminase 2: evidence for a substrate-assisted mechanism. Biochemistry 2014, 53, 4426,  DOI: 10.1021/bi500554b
    16. 16
      Vossenaar, E. R.; Zendman, A. J.; van Venrooij, W. J.; Pruijn, G. J. PAD, a growing family of citrullinating enzymes: genes, features and involvement in disease. BioEssays 2003, 25, 11061118,  DOI: 10.1002/bies.10357
    17. 17
      Tarcsa, E.; Marekov, L. N.; Mei, G.; Melino, G.; Lee, S. C.; Steinert, P. M. Protein unfolding by peptidylarginine deiminase. Substrate specificity and structural relationships of the natural substrates trichohyalin and filaggrin. J. Biol. Chem. 1996, 271, 30709,  DOI: 10.1074/jbc.271.48.30709
    18. 18
      Takahara, H. A History of Deimination Research in Japan: The Founding Fathers. In Protein Deimination in Human Health and Disease; Nicholas, A. P.; Bhattacharya, S. K.; Thompson, P. R., Eds.; Springer International Publishing: Cham, 2017; pp 110.
    19. 19
      Cao, L.; Goodin, R.; Wood, D.; Moscarello, M. A.; Whitaker, J. N. Rapid release and unusual stability of immunodominant peptide 45-89 from citrullinated myelin basic protein. Biochemistry 1999, 38, 61576163,  DOI: 10.1021/bi982960s
    20. 20
      Whitaker, J. N. Myelin encephalitogenic protein fragments in cerebrospinal fluid of persons with multiple sclerosis. Neurology 1977, 27, 911920,  DOI: 10.1212/WNL.27.10.911
    21. 21
      Whitaker, J. R.; Granum, P. E. An absolute method for protein determination based on difference in absorbance at 235 and 280 nm. Anal. Biochem. 1980, 109, 156159,  DOI: 10.1016/0003-2697(80)90024-X
    22. 22
      Pritzker, L. B.; Joshi, S.; Gowan, J. J.; Moscarello, M. A.; Harauz, G. Deimination of myelin basic protein. 1. Effect of deimination of arginyl residues of myelin basic protein on its structure and susceptibility to digestion by cathepsin D. Biochemistry 2000, 39, 53745381,  DOI: 10.1021/bi9925569
    23. 23
      Wood, D. D.; Bilbao, J. M.; O’Connors, P.; Moscarello, M. A. Acute multiple sclerosis (Marburg type) is associated with developmentally immature myelin basic protein. Ann. Neurol. 1996, 40, 1824,  DOI: 10.1002/ana.410400106
    24. 24
      Mastronardi, F. G.; Wood, D. D.; Mei, J.; Raijmakers, R.; Tseveleki, V.; Dosch, H. M.; Probert, L.; Casaccia-Bonnefil, P.; Moscarello, M. A. Increased citrullination of histone H3 in multiple sclerosis brain and animal models of demyelination: a role for tumor necrosis factor-induced peptidylarginine deiminase 4 translocation. J. Neurosci. 2006, 26, 1138711396,  DOI: 10.1523/JNEUROSCI.3349-06.2006
    25. 25
      Nicholas, A. P.; Whitaker, J. N. Preparation of a monoclonal antibody to citrullinated epitopes: its characterization and some applications to immunohistochemistry in human brain. Glia 2002, 37, 328336,  DOI: 10.1002/glia.10039
    26. 26
      Raijmakers, R.; Vogelzangs, J.; Croxford, J. L.; Wesseling, P.; van Venrooij, W. J.; Pruijn, G. J. Citrullination of central nervous system proteins during the development of experimental autoimmune encephalomyelitis. J. Comp. Neurol. 2005, 486, 243253,  DOI: 10.1002/cne.20529
    27. 27
      Cao, L.; Sun, D.; Whitaker, J. N. Citrullinated myelin basic protein induces experimental autoimmune encephalomyelitis in Lewis rats through a diverse T cell repertoire. J. Neuroimmunol. 1998, 88, 2129,  DOI: 10.1016/S0165-5728(98)00063-0
    28. 28
      Carrillo-Vico, A.; Leech, M. D.; Anderton, S. M. Contribution of myelin autoantigen citrullination to T cell autoaggression in the central nervous system. J. Immunol. 2010, 184, 28392846,  DOI: 10.4049/jimmunol.0903639
    29. 29
      Quinn, T. A.; Dutt, M.; Shindler, K. S. Optic neuritis and retinal ganglion cell loss in a chronic murine model of multiple sclerosis. Front. Neurol. 2011, 2, 50  DOI: 10.3389/fneur.2011.00050
    30. 30
      Beniac, D. R.; Luckevich, M. D.; Czarnota, G. J.; Tompkins, T. A.; Ridsdale, R. A.; Ottensmeyer, F. P.; Moscarello, M. A.; Harauz, G. Three-dimensional structure of myelin basic protein. I. Reconstruction via angular reconstitution of randomly oriented single particles. J. Biol. Chem. 1997, 272, 42614268,  DOI: 10.1074/jbc.272.7.4261
    31. 31
      Ridsdale, R. A.; Beniac, D. R.; Tompkins, T. A.; Moscarello, M. A.; Harauz, G. Three-dimensional structure of myelin basic protein. II. Molecular modeling and considerations of predicted structures in multiple sclerosis. J. Biol. Chem. 1997, 272, 42694275,  DOI: 10.1074/jbc.272.7.4269
    32. 32
      Haas, H.; Oliveira, C. L.; Torriani, I. L.; Polverini, E.; Fasano, A.; Carlone, G.; Cavatorta, P.; Riccio, P. Small angle x-ray scattering from lipid-bound myelin basic protein in solution. Biophys. J. 2004, 86, 455460,  DOI: 10.1016/S0006-3495(04)74122-3
    33. 33
      Jahn, O.; Tenzer, S.; Werner, H. Myelin Proteomics: Molecular Anatomy of an Insulating Sheath. Mol. Neurobiol. 2009, 40, 5572,  DOI: 10.1007/s12035-009-8071-2
    34. 34
      Fewster, M. E.; Hirono, H.; Mead, J. F. Lipid composition of myelin in multiple sclerosis. J. Neurol. 1976, 213, 119131,  DOI: 10.1007/BF00313273
    35. 35
      Wilson, R.; Tocher, D. R. Lipid and fatty acid composition is altered in plaque tissue from multiple sclerosis brain compared with normal brain white matter. Lipids 1991, 26, 915,  DOI: 10.1007/BF02544017
    36. 36
      Alling, C.; Vanier, M. T.; Svennerholm, L. Lipid alterations in apparently normal white matter in multiple sclerosis. Brain Res. 1971, 35, 325336,  DOI: 10.1016/0006-8993(71)90478-1
    37. 37
      Woelk, H.; Borri, P. Lipid and fatty acid composition of myelin purified from normal and MS brains. Eur. Neurol. 1973, 10, 250260,  DOI: 10.1159/000114281
    38. 38
      Gerstl, B.; Eng, L. F.; Tavaststjerna, M.; Smith, J. K.; Kruse, S. L. Lipids and proteins in multiple sclerosis white matter. J. Neurochem. 1970, 17, 677689,  DOI: 10.1111/j.1471-4159.1970.tb00547.x
    39. 39
      Arnetoli, G.; Pazzagli, A.; Amaducci, L. Fatty acid and aldehyde changes in choline- and ethanolamine-containing phospholipids in the white matter of multiple sclerosis brains. J. Neurochem. 1969, 16, 461463,  DOI: 10.1111/j.1471-4159.1969.tb10387.x
    40. 40
      Göpfert, E.; Pytlik, S.; Debuch, H. 2′,3′-Cyclic nucleotide 3′-phosphohydrolase and lipids of myelin from multiple sclerosis and normal brains. J. Neurochem. 1980, 34, 732739,  DOI: 10.1111/j.1471-4159.1980.tb11205.x
    41. 41
      Neu, I.; Woelk, H. Investigations of the lipid metabolism of the white matter in multiple sclerosis Changes in glycero-phosphatides and lipid-splitting enzymes. Neurochem. Res. 1982, 7, 727735,  DOI: 10.1007/BF00965525
    42. 42
      Del Boccio, P.; Pieragostino, D.; Di Ioia, M.; Petrucci, F.; Lugaresi, A.; De Luca, G.; Gambi, D.; Onofrj, M.; Di Ilio, C.; Sacchetta, P.; Urbani, A. Lipidomic investigations for the characterization of circulating serum lipids in multiple sclerosis. J. Proteomics 2011, 74, 28262836,  DOI: 10.1016/j.jprot.2011.06.023
    43. 43
      Boggs, J. M. Myelin Basic Protein; Nova Science Publishers, Incorporated: New York, 2008.
    44. 44
      Procaccini, C.; De Rosa, V.; Pucino, V.; Formisano, L.; Matarese, G. Animal models of Multiple Sclerosis. Eur. J. Pharmacol. 2015, 759, 182191,  DOI: 10.1016/j.ejphar.2015.03.042
    45. 45
      Niroomand, H.; Pamu, R.; Mukherjee, D.; Khomami, B. Microenvironment alterations enhance photocurrents from photosystem I confined in supported lipid bilayers. J. Mater. Chem. A 2018, 6, 1228112290,  DOI: 10.1039/C8TA00898A
    46. 46
      Israelachvili, J. N.; Mitchell, D. J.; Ninham, B. W. Theory of self-assembly of lipid bilayers and vesicles. Biochim. Biophys. Acta, Biomembr. 1977, 470, 185201,  DOI: 10.1016/0005-2736(77)90099-2
    47. 47
      Israelachvili, J. N. Intermolecular and Surface Forces; Academic Press, 2015.
    48. 48
      Wood, D. D.; Moscarello, M. A. The isolation, characterization, and lipid-aggregating properties of a citrulline containing myelin basic protein. J. Biol. Chem. 1989, 264, 51215127
    49. 49
      Harauz, G.; Musse, A. A Tale of Two Citrullines—Structural and Functional Aspects of Myelin Basic Protein Deimination in Health and Disease. Neurochem. Res. 2007, 32, 137158,  DOI: 10.1007/s11064-006-9108-9
    50. 50
      Matsumoto, T.; Kobayashi, T.; Kamata, K. Role of lysophosphatidylcholine (LPC) in atherosclerosis. Curr. Med. Chem. 2007, 14, 32093220,  DOI: 10.2174/092986707782793899
    51. 51
      Muramatsu, R.; Kuroda, M.; Matoba, K.; Lin, H.; Takahashi, C.; Koyama, Y.; Yamashita, T. Prostacyclin prevents pericyte loss and demyelination induced by lysophosphatidylcholine in the central nervous system. J. Biol. Chem. 2015, 290, 1151511525,  DOI: 10.1074/jbc.M114.587253
    52. 52
      Plemel, J. R.; Michaels, N. J.; Weishaupt, N.; Caprariello, A. V.; Keough, M. B.; Rogers, J. A.; Yukseloglu, A.; Lim, J.; Patel, V. V.; Rawji, K. S.; Jensen, S. K.; Teo, W.; Heyne, B.; Whitehead, S. N.; Stys, P. K.; Yong, V. W. Mechanisms of lysophosphatidylcholine-induced demyelination: A primary lipid disrupting myelinopathy. Glia 2018, 66, 327347,  DOI: 10.1002/glia.23245
    53. 53
      Birgbauer, E.; Rao, T. S.; Webb, M. Lysolecithin induces demyelination in vitro in a cerebellar slice culture system. J. Neurosci. Res. 2004, 78, 157166,  DOI: 10.1002/jnr.20248
    54. 54
      Boggs, J. M. Myelin basic protein: a multifunctional protein. Cell. Mol. Life Sci. 2006, 63, 19451961,  DOI: 10.1007/s00018-006-6094-7
    55. 55
      Harauz, G.; Ishiyama, N.; Hill, C. M.; Bates, I. R.; Libich, D. S.; Fares, C. Myelin basic protein-diverse conformational states of an intrinsically unstructured protein and its roles in myelin assembly and multiple sclerosis. Micron 2004, 35, 503542,  DOI: 10.1016/j.micron.2004.04.005
    56. 56
      D’Souza, C. A.; Wood, D. D.; She, Y.-M.; Moscarello, M. A. Autocatalytic cleavage of myelin basic protein: An alternative to molecular mimicry. Biochemistry 2005, 44, 1290512913,  DOI: 10.1021/bi051152f
    57. 57
      Zhou, S. R.; Moscarello, M. A.; Whitaker, J. N. The effects of citrullination or variable amino-terminus acylation on the encephalitogenicity of human myelin basic protein in the PL/J mouse. J. Neuroimmunol. 1995, 62, 147152,  DOI: 10.1016/0165-5728(95)00112-3
    58. 58
      Madsen, P. M.; Motti, D.; Karmally, S.; Szymkowski, D. E.; Lambertsen, K. L.; Bethea, J. R.; Brambilla, R. Oligodendroglial TNFR2 Mediates Membrane TNF-Dependent Repair in Experimental Autoimmune Encephalomyelitis by Promoting Oligodendrocyte Differentiation and Remyelination. J. Neurosci. 2016, 36, 51285143,  DOI: 10.1523/JNEUROSCI.0211-16.2016
    59. 59
      Porciatti, V.; Saleh, M.; Nagaraju, M. The pattern electroretinogram as a tool to monitor progressive retinal ganglion cell dysfunction in the DBA/2J mouse model of glaucoma. Invest. Ophthalmol. Visual Sci. 2007, 48, 745751,  DOI: 10.1167/iovs.06-0733
    60. 60
      Menon, K.; Rasband, M. N.; Taylor, C. M.; Brophy, P.; Bansal, R.; Pfeiffer, S. E. The myelin–axolemmal complex: biochemical dissection and the role of galactosphingolipids. J. Neurochem. 2003, 87, 9951009,  DOI: 10.1046/j.1471-4159.2003.02075.x
    61. 61
      Valdivia, A. O.; Myer, C.; Suarez, M. F.; Bhattacharya, S. K. Protein–Lipid Complex Separation Utilizing a Capillary Electrophoresis System. In Metabolomics: Methods and Protocols; Bhattacharya, S. K., Ed.; Springer: New York, 1996; pp 95100.
    62. 62
      Chauhan, M. Z.; Valencia, A. K.; Piqueras, M. C.; Enriquez-Algeciras, M.; Bhattacharya, S. K. Optic Nerve Lipidomics Reveal Impaired Glucosylsphingosine Lipids Pathway in Glaucoma. Invest. Ophthalmol. Visual Sci. 2019, 60, 17891798,  DOI: 10.1167/iovs.18-25802
    63. 63
      Travers, T. S.; Harlow, L.; Rosas, I. O.; Gochuico, B. R.; Mikuls, T. R.; Bhattacharya, S. K.; Camacho, C. J.; Ascherman, D. P. Extensive Citrullination Promotes Immunogenicity of HSP90 through Protein Unfolding and Exposure of Cryptic Epitopes. J. Immunol. 2016, 197, 19261936,  DOI: 10.4049/jimmunol.1600162
    64. 64
      Wiedemann, C.; Bellstedt, P.; Gorlach, M. CAPITO--a web server-based analysis and plotting tool for circular dichroism data. Bioinformatics 2013, 29, 17501757,  DOI: 10.1093/bioinformatics/btt278
    65. 65
      Ahmed, M. A.; De Avila, M.; Polverini, E.; Bessonov, K.; Bamm, V. V.; Harauz, G. Solution nuclear magnetic resonance structure and molecular dynamics simulations of a murine 18.5 kDa myelin basic protein segment (S72-S107) in association with dodecylphosphocholine micelles. Biochemistry 2012, 51, 74757487,  DOI: 10.1021/bi300998x
    66. 66
      Ishida, T.; Kinoshita, K. PrDOS: prediction of disordered protein regions from amino acid sequence. Nucleic Acids Res. 2007, 35, W460W464,  DOI: 10.1093/nar/gkm363
    67. 67
      Song, Y.; DiMaio, F.; Wang, R. Y.; Kim, D.; Miles, C.; Brunette, T.; Thompson, J.; Baker, D. High-resolution comparative modeling with RosettaCM. Structure 2013, 21, 17351742,  DOI: 10.1016/j.str.2013.08.005
    68. 68
      Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The missing term in effective pair potentials. J. Phys. Chem. A 1987, 91, 62696271,  DOI: 10.1021/j100308a038
    69. 69
      Maier, J. A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K. E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 36963713,  DOI: 10.1021/acs.jctc.5b00255
    70. 70
      Duff, M. R., Jr.; Borreguero, J. M.; Cuneo, M. J.; Ramanathan, A.; He, J.; Kamath, G.; Chennubhotla, S. C.; Meilleur, F.; Howell, E. E.; Herwig, K. W.; Myles, D. A. A.; Agarwal, P. K. Modulating Enzyme Activity by Altering Protein Dynamics with Solvent. Biochemistry 2018, 57, 42634275,  DOI: 10.1021/acs.biochem.8b00424
    71. 71
      Martínez, L.; Andrade, R.; Birgin, E. G.; Martinez, J. M. PACKMOL: a package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 21572164,  DOI: 10.1002/jcc.21224
    72. 72
      Agarwal, P. K. Cis/trans isomerization in HIV-1 capsid protein catalyzed by cyclophilin A: insights from computational and theoretical studies. Proteins 2004, 56, 449463,  DOI: 10.1002/prot.20135
    73. 73
      Beck, D. A.; Daggett, V. Methods for molecular dynamics simulations of protein folding/unfolding in solution. Methods 2004, 34, 112120,  DOI: 10.1016/j.ymeth.2004.03.008
    74. 74
      Roe, D. R.; Cheatham, T. E., 3rd PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 30843095,  DOI: 10.1021/ct400341p
    75. 75
      Shukla, S.; Bafna, K.; Gullett, C.; Myles, D. A. A.; Agarwal, P. K.; Cuneo, M. J. Differential Substrate Recognition by Maltose Binding Proteins Influenced by Structure and Dynamics. Biochemistry 2018, 57, 58645876,  DOI: 10.1021/acs.biochem.8b00783
  • Supporting Information

    Supporting Information


    The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsomega.0c01590.

    • Confirmation of the EAE model, variation of hydrocarbon chain, computational modeling of MBP, and five model structures of MBP and MBP sequence alignment (Figures S1–S5); mass spectrometry results for arginine residues and confirmation of MBP presence in gel-excised band (Table S1); and computer simulations of MBP (Movies S1–S3) (ZIP)

    Accession Codes

    MBP (UniProt accession: P02686-5 and P02687). Mass spectrometry data (PRIDE archive accession: PXD015494).


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.