Deconvolution of the MBP-Bri2 Interaction by a Yeast Two Hybrid System and Synergy of the AlphaFold2 and High Ambiguity Driven Protein-Protein Docking

: Myelin basic protein (MBP) is one of the key proteins in the development of multiple sclerosis (MS). However, very few intracellular MBP partners have been identiﬁed up to now. In order to ﬁnd proteins interacting with MBP in the brain, an expression library from the human brain was screened using a yeast two-hybrid system. Here we showed that MBP interacts with the C-terminal 24-residue peptide of Integral transmembrane protein II associated with familial British and Danish dementia (ITM2B/Bri2 or Bri2). This peptide (Bri23R) was one residue longer than the known Bri23 peptide, which is cleaved from the C-terminus of Bri2 during its maturation in the Golgi and has physiological activity as a modulator of amyloid precursor protein processing. Since the spatial structures for both MBP and Bri2 were not known, we used computational methods of structural biology including an artiﬁcial intelligence system AlphaFold2 and high ambiguity driven protein-protein docking (HADDOCK 2.1) to gain a mechanistic explanation of the found protein-protein interaction and elucidate a possible structure of the complex of MBP with Bri23R peptide. As expected, MBP was mostly unstructured, although it has well-deﬁned α -helical regions, while Bri23R forms a stable β -hairpin. Simulation of the interaction between MBP and Bri23R in two different environments, as parts of the two-hybrid system fusion proteins and in the form of single polypeptides, showed that MBP twists around Bri23R. The observed interaction results in the adjustment of the size of the internal space between MBP α -helices to the size of the β -hairpin of Bri23R. Since Bri23 is known to inhibit aggregation of amyloid oligomers, and the association of MBP to the inner leaﬂet of the membrane bilayer shares features with amyloid ﬁbril formation, Bri23 may serve as a peptide chaperon for MBP, thus participating in myelin membrane assembly.


Introduction
In 1971, the structural biology community established the single worldwide archive for macromolecular structure data-the Protein Data Bank (PDB). Today, hundreds of data resources and millions of researchers, which explore fundamental biology and biomedicine, use PDB. The tremendous amount of structural information collected over 50 years was used in the creation and testing of AlphaFold2-an artificial intelligence system developed by DeepMind (https://deepmind.com, accessed on 10 October 2021) to predict the spatial GAG (MBP pGBKT7 rev). PCR amplification of MBP-encoding sequence was performed through the following PCR program: 94 • C for 4 min, then 25 cycles each at 94 • C for 30 s/60 • C for 30 s/72 • C for 30 s, followed by incubation at 72 • C for 10 min and cooling to 4 • C. Vector plasmid pGBKT7 was digested using restriction enzymes NdeI and SalI (NewEnglandBiolabs, Ipswich, MA, USA), and MBP PCR fragment was digested using NdeI and XhoI (NewEnglandBiolabs, Ipswich, MA, USA). Both fragments were gel purified and ligated using T4 DNA ligase (ThermoScientific, Waltham, MA, USA). Ligate was transformed into XL-1 cells, which were spread onto kanamycin plates.

Preparation of Competent Yeast Cells
To prepare the yeast competent cells we followed the protocol provided with the Yeastmaker Yeast Transformation System 2 (Takara Bio USA, Inc., San Jose, CA, USA). Briefly, Y187 or Y2H strains from a frozen yeast stock were streaked on YPDA agar plates. The plates were incubated at 30 • C until colonies appeared (approximately 3 days). After that, one fresh 2-3 mm colony was inoculated into 3 mL of YPDA medium containing 50 µg/mL kanamycin in a sterile 15 mL culture tube. The tubes were incubated at 30 • C with shaking at 250 rpm for 8-12 h. Then the 5 µL of the culture were transferred to 50 mL of YPDA in a 250 mL flask, and the culture was incubated with shaking until the OD 600 reached 0.3 (20 h). Then the cells were collected by centrifugation at 700 g for 5 min at room temperature. The pellet was resuspended in 100 mL of fresh YPDA and incubated at 30 • C until the OD 600 reached 0.5 (3 h). The cells from the culture were collected at 700 g for 5 min at room temperature into two 50 mL sterile Falcon conical tubes resuspended in 30 mL sterile, deionized H 2 O. The cells were again centrifuged at 700 g for 5 min at room temperature, the supernatant was discarded, and each pellet was resuspended in 1.5 mL of 1.1× TE/LiAc. Cell suspensions were transferred into 1.5 mL tubes, and the cells were pelleted at high speed for 15 s. Finally, each cell pellet was resuspended in 600 µL of 1.1× TE/LiAc, and that suspension was used for transformation with plasmid DNA.

Transformation of Yeast Cells
The transformation of yeast cells was performed according to the small-scale protocol provided with the Yeastmaker Yeast Transformation System 2 (Takara Bio USA, Inc., San Jose, CA, USA). Briefly, in a pre-chilled 1.5 mL tube, 100 ng of appropriate plasmid DNA was combined and mixed with 50 µg of denatured Yeastmaker Carrier DNA and 50 µL freshly prepared yeast competent cells. Then, 500 µL of PEG/LiAc was added to the plasmid/cells mixture. The resulting volume was gently mixed and incubated at 30 • C for 30 min with inverting the tube every 10 min. The 20 µL of DMSO was added to the volume and mixed by inverting the tube. Then the volume was gently incubated at 42 • C for 15 min with inverting the tube every 5 min. After incubation, the cells were collected by centrifugation at high speed for 15 s and resuspended into 1 mL of YPD Plus Medium, and the suspension was incubated at 30 • C for 30 min with shaking. The cells were collected by centrifugation at high speed for 15 s and resuspended into 1 mL of 0.9% NaCl solution. Then the cells were seeded to the appropriate selection media.

Obtaining of Control Mated Strains
To obtain control strains containing two plasmids, a procedure for mating yeast cells was performed. One colony for each individual strain to mate was selected for use in the mating procedure. Both colonies were placed in one 1.5 mL centrifuge tube containing 500 µL of 2× YPDA and vortexed. The mixture was incubated with shaking at 200 rpm at 30 • C overnight for 20 h. From the mated culture (0.5 mL), 100 µL of 1/10, 1/100 and 1/1000 dilutions were plated on each of the following Petri dishes with agar: SD/-Trp; SD/-Leu; SD/-Leu/-Trp (=DDO); SD/-Leu/-Trp /X-α-Gal/AbA (=DDO/X/A). Cells were grown for 3-5 days at 30 • C.

Two-Hybrid Library Screening Using Yeast Mating
For library screening, the concentrated bait culture was mixed with 1 mL of the normalized Mate & Plate library (Takara Bio USA, Inc., San Jose, CA, USA) and incubated overnight before plating on DDO/X/A selective medium. A concentrated overnight culture of the Y2HGold-MBP bait strain was prepared. For this, one fresh large (2-3 mm) colony of the bait strain was seeded in 50 mL of SD/-Trp liquid medium. The cells were incubated with shaking (250-270 rpm) at 30 • C until OD 600 reached 0.8 (20 h). The cells were precipitated by centrifugation (at 1000× g for 5 min), after which the supernatant was removed. The pellet was resuspended to a cell density of >1 × 10 8 cells/mL in SD/-Trp liquid media (4-5 mL). Then the library strain was combined with the bait strain. For this, a 1 mL aliquot of the library strain was thawed in a water bath at room temperature. 10 µL was left for titration on 100 mm plates with solid SD/-Leu medium. 1 mL of the Mate & Plate Library was mixed with 4-5 mL of the bait strain in a sterile 2 L flask, and 45 mL of 2× YPDA liquid medium (containing 50 µg/mL of kanamycin) was added. Cells from the library flask were washed twice with 1 mL of 2× YPDA, and washes were added to the 2 L flask. The cells were incubated at 30 • C for 20-24 h with minimal shaking (30-50 rpm). After 20 h, a drop of the culture was checked under a phase contrast microscope (40×) for the presence of zygotes. Since they were present in large numbers, which indicated a successful crossing, the cells were harvested by centrifugation at 1000× g for 10 min. The flask was washed twice with 50 mL 0.5× YPDA (with 50 µg/mL kanamycin), the washes were combined and used to resuspend the pelleted cells. The cells were centrifuged at 1000× g for 10 min and the supernatant was removed. The pelleted cells were resuspended in 10 mL of 0.5× YPDA/Kan liquid medium. The total volume of the cell suspension was measured. To count the number of clones resulting from the mated culture, 100 µL of 1/10, 1/100, 1/1000 and 1/10,000 dilutions were plated on each of the following 100 mm agar plates and incubated at 30 • C for 3-5 days: SD/-Trp; SD/-Leu; SD/-Leu/-Trp. The rest of the cell culture was plated by adding 100 µL onto 120 100 mm plates with DDO/X/A agar medium. The library was allowed to grow at 30 • C for 3-5 days. All blue colonies that grew on DDO/X/A medium were streaked onto plates with QDO/X/A agar (SD/-Ade/-His/-Leu/-Trp/X-α-Gal/AbA) with a higher stringency using sterile yellow pipette tips.

Analysis of Positive Interactions by PCR
For primary positive colony screening by PCR, DNA from blue yeast colonies on QDO/X/A was isolated using a DNA-Sorb-AM kit (Federal Budget Institution of Science «Central Research Institute of Epidemiology» of the Federal Service for Surveillance on Consumer Rights Protection and Human Wellbeing, Moscow, Russia) and used as a template in PCR using primers MM_screen_For and MM_screen_Rev, flanking the expressed inserts in pGADT7 plasmid. The primer sequences are 5 -CTATTCGATGATGAAGATACCCCACCA AACC-3 for MM_screen_For and 5 -GTGAACTTGCGGGGTTTTTCAGTATCTACGATT-3 for MM_screen_Rev. The PCR conditions were the following: 300 s at 95 • C, 30 cycles of 420 s at 72 • C and 15 s at 95 • C with a final extension of 600 s at 72 • C. PCR products were resolved in 1% agarose gel containing 0.5 µg/mL of ethidium bromide. Positive bands were excised from the gel, purified by QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany), and their nucleotide sequences were determined by Sanger sequencing with T7 Sequencing Primer.

Confirmation of MBP-Bri23R Positive Interactions in Y2HGold Yeast Strain
Using the small-scale transformation procedure, 100 ng of each of the following pairs of vectors (pGBKT7/MBP + Bri23R in pGADT7; empty pGBKT7 + Bri23R in pGADT7) were cotransformed into Y2HGold competent cells. 100 µL of 1/10 and 1/100 dilutions of the transformation mixes were spread on the DDO/X and QDO/X/A plates. The results of the genuine interaction of MBP with Bri23R were represented by blue colonies grown on DDO/X and QDO/X/A plates for pGBKT7/MBP + Bri23R in pGADT7 pair and white colonies on DDO/X and no colonies on QDO/X/A plates for the empty pGBKT7 + Bri23R in pGADT7 pair (see Supplementary Figure S1). False positives would give blue colonies on any plate for both pairs.

Additional Confirmation of MBP-Bri23R Positive Interactions in Y2HGold Yeast Strain Using Reciprocal Approach
MBP and Bri23R coding sequences were reversely cloned in pGADT7 and pGBKT7 plasmids, respectively. For this, the NdeI/XhoI MBP fragment from p.2.1 of Methods above was cloned in pGADT7 digested with NdeI/XhoI, and Bri23R coding fragment was amplified from Bri23R/pGADT7 using MM_screen_For and MM_screen_Rev primers (300 s at 95 • C, 20 cycles of 60 s at 72 • C and 15 s at 95 • C with a final extension of 420 s at 72 • C), digested with NdeI/XhoI and cloned in pGBKT7 digested with NdeI/SalI.
Then the reciprocal analysis for interaction was performed the same as described in p.2.8 of Methods above for direct mode interaction, except the following pairs of plasmids were used: pGBKT7/Bri23R + MBP/pGADT7; pGBKT7/Bri23R + pGADT7. The results were shown in Supplementary Figure S6, where pGBKT7/Bri23R + MBP/pGADT7 plasmid pair gave blue colonies on both DDO/X and QDO/X/A plates, and pGBKT7/Bri23R + pGADT7 plasmid pair gave white colonies and no colonies on DDO/X and QDO/X/A plates, respectively, confirming true MBP-Bri23R interaction (see Supplementary Figure S1).
The qualities of the models were validated using the AlphaFold lDDT-Cα metric (Supplementary Figure S3), which produces a chart of per-residue estimated confidence on a scale from 0 to 100. Regions with pLDDT > 90 are expected to have very high confidence, regions with pLDDT between 70 and 90 are expected to have good confidence, regions with pLDDT between 50 and 70 are of low confidence, and pLDDT < 50 is a strong predictor of IDP suggesting that such region is either unstructured in physiological conditions or only structured as part of a complex.

Docking Calculations
The 3D models of the complexes were obtained with high ambiguity driven proteinprotein docking (HADDOCK 2.1) [25,26]. The docking algorithm has three stages. At first, rigid body docking was performed, in which interacting molecules were rotated and translated randomly in turn to minimize intermolecular energy. It was followed by the annealing stage, in which annealing of torsion angle space was performed to refine the orientation of the molecules and the side chains and/or backbones of the interface residues. The last stage was solvent refinement in which the structures were further refined in explicit solvent layers. The 3D coordinates of six AlphaFold2-built models were used in the calculation.
HADDOCK presupposes the use of the defined intermolecular ambiguous interaction constraints (AIR) in the calculations. The AIR residues are assumed to have no less than 40% solvent accessible surface area. Surface accessibilities were assessed with the program FreeSASA [27]. Unless otherwise specified, solvent accessible residues of Bri23R and MBP were suggested as AIR. The 1000 structures of the complexes were calculated with the rigid docking protocol. The 200 structures with the lowest AIR violations were energy minimized with the side chains left flexible. Then the 30 best structures were minimized again in an 8 Å shell of explicit TIP3P water [28]. These best structures were grouped into clusters based on the FCC clustering algorithm [29]. The buried surface area in the complexes was analyzed and visualized using PyMOL and PDBePisa [30].

Interaction of MBP with a C-Terminal Peptide of Bri2 Was Discovered Using a Yeast Two Hybrid System
To identify the spectrum of proteins interacting with MBP, Matchmaker Gold Y2HS was used. In a Matchmaker GAL4-based two-hybrid assay, a bait protein is expressed as a fusion to the Gal4 DNA-binding domain (Gal4BD), while libraries of prey proteins are expressed as fusions to the Gal4 activation domain (Gal4AD) [23]. In the Y2HS, when the bait and library (prey) fusion proteins interact, the Gal4BD and Gal4AD are brought into proximity to activate transcription of four independent reporter genes: AUR1-C, ADE2, HIS3, and MEL1 ( Figure 1A).
in explicit solvent layers. The 3D coordinates of six AlphaFold2-built models were used in the calculation.
HADDOCK presupposes the use of the defined intermolecular ambiguous interaction constraints (AIR) in the calculations. The AIR residues are assumed to have no less than 40% solvent accessible surface area. Surface accessibilities were assessed with the program FreeSASA [27]. Unless otherwise specified, solvent accessible residues of Bri23R and MBP were suggested as AIR.
The 1000 structures of the complexes were calculated with the rigid docking protocol. The 200 structures with the lowest AIR violations were energy minimized with the side chains left flexible. Then the 30 best structures were minimized again in an 8 Å shell of explicit TIP3P water [28]. These best structures were grouped into clusters based on the FCC clustering algorithm [29]. The buried surface area in the complexes was analyzed and visualized using PyMOL and PDBePisa [30].

Interaction of MBP with a C-Terminal Peptide of Bri2 Was Discovered Using a Yeast Two Hybrid System
To identify the spectrum of proteins interacting with MBP, Matchmaker Gold Y2HS was used. In a Matchmaker GAL4-based two-hybrid assay, a bait protein is expressed as a fusion to the Gal4 DNA-binding domain (Gal4BD), while libraries of prey proteins are expressed as fusions to the Gal4 activation domain (Gal4AD) [23]. In the Y2HS, when the bait and library (prey) fusion proteins interact, the Gal4BD and Gal4AD are brought into proximity to activate transcription of four independent reporter genes: AUR1-C, ADE2, HIS3, and MEL1 ( Figure 1A).  In Clontech's Saccharomyces cerevisiae Y2H Gold strain, which is used to detect twohybrid interactions, the aforementioned integrated reporter genes are under the control of three distinct Gal4-responsive promoters. AUR1-C is a dominant mutant version of the AUR1 gene that encodes the enzyme inositol phosphoryl ceramide synthase. The expression of AUR1-C confers strong resistance to the otherwise highly toxic drug Aureobasidin A. This drug reporter is preferable to nutritional reporters alone, due to lower background activity. HIS3: Y2HGold is unable to synthesize histidine and is, therefore, unable to grow on media that lacks this essential amino acid. Gal4-responsive HIS3 expression permits cells to biosynthesize histidine and grow on His-minus minimal medium. ADE2: Y2HGold is also unable to grow on minimal media that does not contain adenine. However, when the Ade2 expression is activated, the cells can grow on minimal medium without adenine. MEL1 encodes α-galactosidase, an enzyme occurring naturally in many yeast strains, which is expressed and secreted by yeast, and yeast colonies that express Mel1 turn blue in the presence of the chromogenic substrate X-α-Gal.
To produce bait fusion protein Gal4BD-MBP, containing the Gal4BD at its N-terminal end and MBP at its C-terminus, we inserted human MBP coding sequence in the ORF of GAL4BD into the pGBKT7 and generated pGBKT7-MBP ( Figure 1B).
Before the screening process, we obtained the Y2HGold strain transformed with the pGBKT7-MBP bait plasmid and performed all the necessary check-ups for the GAL4-MBP fusion potential toxicity and autoactivation properties according to the manufacturer's instructions. To generate positive and negative control strains we used plasmid constructs supplied in the kit (pGBKT7-53 Control vector, which encodes the Gal4BD fused with murine p53; pGADT7-T Control vector, encodes the Gal4AD fused with SV40 large Tantigen; and pGBKT7-Lam Control vector, which encodes the Gal4BD fused with lamin). pGBKT7-53 and pGADT7-T plasmids are used to create a positive control strain, as well as pGBKT7-Lam and pGADT7-T plasmids that are used to create negative control strain.
As prey, we used Takara's Mate & Plate Library-Human Brain (Normalized). This yeast two-hybrid library was constructed from mRNA isolated from human brain tissue and transformed into yeast strain Y187. The cDNA was normalized prior to library construction to reduce the copy number of abundant cDNAs derived from highly represented mRNAs, thereby increasing the representation of low copy number transcripts.
For screening, the obtained yeast strain Y2HGold-MBP was mated with Mate & Plate Library in Y187 strain. Primary interactions between respective proteins were identified by growth on double dropout selective media (QDO/X/A) and blue color development. All primary signals were streaked onto plates with quadruple dropout selective media (QDO/X/A) agar providing higher selection stringency. Colonies growing and maintaining the blue color on QDO/X/A agar were subjected to PCR analysis of the inserts in the pGADT7 plasmid. In total, about 200 primary potentially positive clones were selected for analysis, which number was reduced to 56 in QDO/X/A and finally to the 17 unique clones after removing duplicates by PCR and sequence analysis.
Finally, of the identified 17 unique positive signals, only one clone corresponded to the coding region of the protein (Supplementary Table S1, the coding clone highlighted in bold). This protein was identified as ITM2B/Bri2 (Integral transmembrane protein II associated with familial British and Danish dementia (FBD and FDD)), with Genebank ID NM_021999.5 and Uniprot ID Q9Y287. The found clone encoded C-terminal 24-residue peptide of Bri2 (see Figure 2B for illustration).  Figure 2 and Supplementary Figure S2. were fused at the C-termini of Gal4BD and Gal4AD, respectively, or vice versa. Coloring is according to the model confidence (left panels) and domain organization (right panels): Gal4AD-dark gray, Gal4BD-light gray; in MBP, salmon and cyan colors correspond coils and α-helices, respectively; in Bri23R, magenta and green colors correspond β-strands and coils, respectively. See Table  1 and pictograms on the right side for illustration.  Figure 2 and Supplementary Figure S2.  Table 1 and pictograms on the right side for illustration.
Bri2 is a 266 amino acid long transmembrane protein of type 2 consisting of a cytoplasmic domain, helical transmembrane domain and luminal domain [31]. It has a suggested role in neuronal maturation and differentiation [32]. Two different autosomal dominant mutations in the ITM2B gene are associated with the aforementioned forms of early-onset dementia [31]. Disease-associated mutations caused the expression of Cterminally extended 277 amino acid long Bri2 proteins, whose furin-dependent proteolysis generated 34 amino acid long peptides (ABri in FBD and ADan in FDD) that accumulated into neurotoxic amyloid aggregates.
To date, the exact physiological function of Bri2 has not been fully elucidated. However, its interaction with Amyloid precursor protein (APP) and involvement in the regulation of APP processing and inhibition of beta-amyloid production was shown [33]. This function as a modulator of APP processing is associated with Bri2 maturation, which started with furindependent cleavage of the C-terminal 23 amino acid long Bri23 peptide from immature Bri2 during its transport in the Golgi [34]. This Bri23 peptide itself prevents the aggregation of APP amyloid-beta protein 42 into toxic oligomers [35]. Truncated Bri2 undergoes further proteolytic processing by ADAM10, which cuts out the extracellular 100 amino acid long evolutionary conserved Brichos domain, and by Signal Peptidase-Like 2B, which cleaves the remaining membrane-bound polypeptide [36].
The 24-residue peptide found by means of the two-hybrid system represents Bri23 peptide with R residue attached at the N-terminus (Bri23R). Interaction between Bri23R and MBP was confirmed in yeasts on selective media using the yeast cotransformation procedure in direct and reciprocal modes where protein-coding inserts were switched between the plasmids coding for Gal4BD and -AD domains (procedure details are described in pp.2.7 and 2.8 of Methods, and results are presented in Figure 1C and Supplementary Figure S1). It should be noted, that in 2018, Rebelo's group conducted a study of the Bri2 interactome in the brain [37]. MBP was identified among the proteins co-precipitated with Bri2 from brain tissue. However, further studies on the specificity of their interaction have not been carried out. These data raised the question what form or domain of Bri2 is actually involved in the interaction with MBP.

Molecular Modeling and Characterization of the Putative Complexes between MBP and the C-Terminal β-Hairpin of Bri2
The 3D models of human MBP and Bri23R as well as those of four fusion proteins, described in Table 1, were prepared by means of AlphaFold2. The sequences of the fusion proteins were identical to those of the interacting protein components in the yeast twohybrid system (Supplementary Figure S2).
All obtained 3D models are shown in Figure 2 with the confidence (pLDDT) level colormapped in the 3D structures. As stated earlier, MBP is an intrinsic disordered protein (IDP), although three α-helices (α1, α2 and α3), predicted with the confidence level from good to high, were clearly visible in the structures (Figure 2A,C and Supplementary Figure S3). The 24 C-terminal amino acids of Bri2 (Bri23R-the Bri23 peptide [36] with one additional R residue at the N-terminus) always form a stable β-hairpin predicted with good confidence ( Figure 2B,C and Supplementary Figure S3). In the case of the luminal domain of native Bri2, this β-hairpin participates in the formation of a curved β-sheet ( Figure 2B), in which it is partially buried; but, as we explained above, the Bri23 peptide is cleaved from Bri2 during its maturation in the Golgi [34,35]. Both Gal4 domains contain α-helical and disordered areas. The α-helical core of Gal4BD is predicted with high accuracy, while Gal4AD is mostly disordered with a small α-helical region predicted with moderate accuracy ( Figure 2C and Supplementary Figure S3). Superimposition at the Cα-atoms of MBP and Gal4 domains from the fusion proteins with those in their native prototypes showed that the secondary structure composition and general structural topologies were similar (data not shown).
Using a data-driven docking program HADDOCK, three types of complexes were designed: Gal4BD-MBP/Gal4AD-Bri23R, Gal4AD-MBP/Gal4BD-Bri23R and MBP/Bri23R (a complex of MBP with the 24 C-terminal amino acid of Bri2). The first two complexes correspond to two variants of the fusion proteins, whose interactions were detected by the Y2HS (see Supplementary Figure S1), while in the third complex neither MBP nor Bri23R was fused to the Gal4 domains and thus, represented native proteins actually existing in cells of higher eukaryotes.
The characteristics of the best clusters generated from the 30 final complexes in all three cases are summarized in Table 2. The lowest energy structures representing the obtained clusters are shown in Figure 3 for the Gal4BD-MBP/Gal4AD-Bri23R and Gal4AD-MBP/Gal4BD-Bri23R complexes and in Figure 4 for the MBP/Bri23R complex.   Table 2 for description) obtained for the best 30 HADDOCK-generated complexes are shown as ribbon-mesh pictograms generated by PDBePISA. MBP-containing fusion proteins are colored gold, while Bri23R-containing fusion proteins are colored brown (see pictograms at the left for illustration); the main chains of interacting residues from the binding interfaces are marked by red and green color, respectively.
Residues involved in the putative binding interfaces and intermolecular contacts between them were evaluated using PDBePISA. The analysis of the intermolecular salt bridges and hydrogen bonds in the representative structures of the best clusters for the Gal4BD-MBP/Gal4AD-Bri23R and Gal4AD-MBP/Gal4BD-Bri23R complexes are shown in  Table 2 for description) obtained for the best 30 HADDOCK-generated complexes are shown as ribbon-mesh pictograms generated by PDBePISA. MBP-containing fusion proteins are colored gold, while Bri23Rcontaining fusion proteins are colored brown (see pictograms at the left for illustration); the main chains of interacting residues from the binding interfaces are marked by red and green color, respectively.  (Supplementary Tables S2 and S3). As we mentioned above, in addition to the evaluation of the putative complexes between MBP-and Bri23R-containing fusion proteins interacting in the Y2HS, the same approach, in which AlphaFold2-based modeling was complemented with HADDOCK- Three main clusters were identified from the 30 lowest-energy structures generated for either the Gal4BD-MBP/Gal4AD-Bri23R or Gal4AD-MBP/Gal4BD-Bri23R complex ( Table 2). The average intermolecular energies of these clusters were higher in the first complex (from 360.1 to 313.8 kcal/mol) compared to those in the second complex (from 307.4 to 152.9 kcal/mol). However, the best models for the Gal4AD-MBP/Gal4BD-Bri23R complex were compatible with those for the Gal4BD-MBP/Gal4AD-Bri23R complex, while two other less energy-efficient clusters were similar to each other ( Figure 3B). The major contribution in the intermolecular energies was from the electrostatic component, which indicates a significant impact of polar interactions in the formation of the complexes.
Residues involved in the putative binding interfaces and intermolecular contacts between them were evaluated using PDBePISA. The analysis of the intermolecular salt bridges and hydrogen bonds in the representative structures of the best clusters for the Gal4BD-MBP/Gal4AD-Bri23R and Gal4AD-MBP/Gal4BD-Bri23R complexes are shown in Supplementary Tables S2 and S3, respectively. The binding interfaces and residues participating in either polar or hydrophobic contacts are marked in Figure 3 and colormapped in the sequences of the interacting proteins in Supplementary Figures S4 and S5, respectively. In accordance with enhanced intermolecular energies, a higher number of interfacial residues involved in the intermolecular polar contacts were detected in the best models of the Gal4BD-MBP/Gal4AD-Bri23R complex compared to those of the Gal4AD-MBP/Gal4BD-Bri23R complex (Supplementary Tables S2 and S3). These residues formed from six to twelve polar contacts between the Gal4BD-MBP and Gal4AD-Bri23R molecules, the majority of which were between MBP and Bri23R components of the fusion proteins (Supplementary Table S2). Two to four polar contacts were found in the representative structures of the Gal4AD-MBP/Gal4BD-Bri23R complex (Supplementary Table S3). At the same time, the Van der Waals energy of the intermolecular interactions was similar for both Gal4BD-MBP/Gal4AD-Bri23R and Gal4AD-MBP/Gal4BD-Bri23R complexes ( Table 2).
The location of the interacting residues in the binding interfaces varied between different clusters of each complex, though the residues from Bri23R and at least two αhelices (α1 and α2) of MBP, as well as their nearest surroundings, were always involved ( Figure 3, Supplementary Figures S4 and S5). In addition to the intermolecular contacts between MBP and Bri23R, other components of the fusion proteins including Gal4AD, Gal4BD and connecting linkers were involved in the complex formation, which was in accordance with activation of the Gal4-induced metabolic pathways in the yeast cells expressing complementing pairs of fusion proteins: either Gal4BD-MBP and Gal4AD-Bri23R or Gal4AD-MBP and Gal4BD-Bri23R (Supplementary Tables S2 and S3).
As we mentioned above, in addition to the evaluation of the putative complexes between MBP-and Bri23R-containing fusion proteins interacting in the Y2HS, the same approach, in which AlphaFold2-based modeling was complemented with HADDOCK-based docking and PDBePISA-based structural analysis was applied to evaluate a suggested interaction between MBP and Bri23R in their native forms. Since the predominant participation of the two α-helices of MBP in the interaction with Bri23R was detected during docking of the fusion proteins; the modeling of the MBP/Bri23R complex was carried out in two variants, specifying either the entire solvent accessible surface of MBP as AIR or only amino acid residues from the alpha-helices. Two main clusters were identified among the 30 lowestenergy structures in either calculation (Table 2 and Figure 4). More energy-efficient clusters were obtained in the first case; the major contribution in the increasing intermolecular energies was from the electrostatic component, which was associated with an enhanced amount of the interfacial polar contacts (Supplementary Table S4 and Figure S6). All four clusters were less energy-efficient than either the Gal4BD-MBP/Gal4AD-Bri23R complexes or the best cluster of the Gal4AD-MBP/Gal4BD-Bri23R complex; however, the binding energies were compatible with or prevailed over those of two less energy-efficient clusters of the Gal4AD-MBP/Gal4BD-Bri23R complex. Given the small size of the Bri23R surface available for intermolecular interactions, such binding energies were rather satisfactory. Moreover, PDBePISA-based analysis showed that cluster # 30 has not only the highest solvation free energy gain from the hydrophobic interactions but the lowest p-value, indicating that the suggested interface is supposed to be interaction-specific (Supplementary Table S4).
The analysis of a pattern of the binding interfaces in the MBP surface revealed that similar to the fusion proteins and independently of using or not using the α-helical part of MBP as AIR, at least two α-helices of MBP are predominantly involved in the intermolecular interaction ( Figure 4 and Supplementary Figure S6). Further analysis of the four clusters generated for the MBP/Bri23R complex shows that in Bri23R, both free termini of the β-hairpin and the loop connecting the two β-strands are solvent accessible and involved in the intermolecular interactions (Figure 4 and Supplementary Figure S6). Moreover, according to Figure 4, the interaction with Bri23R causes some structural ordering of MBP, which twists around the β-hairpin so that the size of the internal space between the α-helices is adjusted to the size of the β-hairpin.
Thus, the computational approach, which included AlphaFold2-based modeling and HADDOCK-based docking complemented the results obtained in the yeast two-hybrid system and allowed not only to assume the formation of a complex between MBP and the C-terminal β-hairpin of Bri2, but also the predominant participation of the alpha-helical regions of MBP in the interaction.
It should be noted that β-hairpins are widespread in the protein kingdom and are involved in many protein-protein interactions. This motif is also popular in synthetically designed proteins and peptides due to its stability and adaptability to broad functions [38]. The twisted shape, open edges for hydrogen bonding, and the hydrophobic surfaces of β-hairpin structures provide a unique ability to fold in a variety of ways. However, the relationship between the sequence of β-hairpins and the way they are assembled is far from understood [39]. β-hairpins are attracting great attention as a regulator of amyloidogenesis, which either plays a role in amyloid assembly as key building blocks of amyloid fibrils or preventing the formation of toxic oligomers [40]. It should also be noted that the Simons group has previously shown that MBP forms a barrier in the cytoplasmic space that allows size selection of proteins, thus preventing the diffusion of most soluble and membrane proteins into the myelin sheath [41]. They also showed that the binding of MBP to the inner membrane mediates the formation of a dense protein network, which is similar to the network formed by amyloid fibrils. These amyloid-like interactions provide the molecular basis for protein extrusion and myelin membrane zipping [42]. We suggested that amyloid-like structures associated with MBP accumulation in the inner leaflet of the membrane bilayer can serve as a seed that attracts native Bri23 peptide and initiates its interaction with MBP. Thus, the prevalence and physiological significance of the observed interaction between MBP and the C-terminal β-hairpin of Bri2 should be carefully studied in vivo and in vitro.

Conclusions
Two subjects of this study, myelin basic protein (MBP) and integral transmembrane protein II associated with familial British and Danish dementia (ITM2B/Bri2 or Bri2), are involved in neuronal development and degeneration. However, the disordered structure of MBP and the multistep complex proteolytic processing of Bri2 complicate their study and, as a result, the exact physiological and pathological roles of these proteins are far from understood. Here, the interaction of MBP with the C-terminal peptide of Bri2 was discovered using a yeast two-hybrid system that allows the interaction of proteins in their natural environment to be detected. The discovered peptide (Bri23R) mimicked the native one that appeared during Bri2 maturation and had physiological activity as a modulator of amyloid precursor protein processing. In this regard, computational methods of structural biology, including the artificial intelligence system AlphaFold2 and high ambiguity driven protein-protein docking (HADDOCK), were used to obtain a mechanistic explanation of the detected protein-protein interaction and to elucidate the possible structure of the MBP complex with the Rri23R peptide. As expected, MBP was mostly unstructured, although it has well-defined α-helical regions, while Bri23R forms a stable β-hairpin. Simulation of the interaction between MBP and Bri23R showed that MBP twists around the peptide, while the internal space between MBP α-helices is adjusted to the size of the β-hairpin. It is known that β-hairpins are ubiquitous protein structural motives found in molecules throughout the tree of life, which can serve as either inductors or inhibitors of amyloidogenesis. In this regard, since the association of MBP to the inner leaflet of the membrane bilayer shares features with amyloid fibril formation and native Bri23 is known to inhibit amyloid aggregation, the found MBP-Bri23R interaction indicates that Bri23 may serve as a peptide chaperon for MBP, thus participating in myelin membrane assembly.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/cryst12020197/s1, Figure S1. Validation assay in yeasts to confirm MBP-Bri interaction. Figure S2. The names and sequences of the fusion proteins, whose 3D models were used in the study. Figure S3. Quality evaluation of the AlphaFold2-built 3D models of the fusion proteins presented as a per-residue confidence score (pLDDT). Figure S4. Sequences of the fusion proteins with color mapping of binding interfaces on the representative structures of the lowest energy clusters were obtained for the best 30 HADDOCK-generated Gal4BD-MBP/Gal4AD-Bri23R complexes. Figure S5. Sequences of the fusion proteins with color mapping of binding interfaces on the representative structures of the lowest energy clusters were obtained for the best 30 HADDOCK-generated Gal4AD-MBP/Gal4BD-Bri23R complexes. Figure S6. Sequences of MBP and Bri23R with color mapping of binding interfaces on the representative structures of the lowest energy clusters were obtained for the best 30 HADDOCK-generated MBP/Bri23R complexes. Table S1. Unique clones were identified by means of a two-hybrid yeast system. Table S2. Polar interactions in the intermolecular interfaces of the representative structures of the lowest energy clusters were obtained for the best 30 HADDOCKgenerated Gal4BD-MBP/Gal4AD-Bri23R complexes. Table S3. Polar interactions in the intermolecular interfaces of the representative structures of the lowest energy clusters were obtained for the best 30 HADDOCK-generated Gal4AD-MBP/Gal4BD-Bri23R complexes. Table S4. The solvation free energy and polar interactions of the intermolecular interfaces of the representative structures of the lowest energy clusters were obtained for the best 30 HADDOCK-generated MBP/Bri23R complexes. Supplementary archive. PDB files of HADDOCK-created models.