Structure of a 14-3-3ε:FOXO3apS253 Phosphopeptide Complex Reveals 14-3-3 Isoform-Specific Binding of Forkhead Box Class O Transcription Factor (FOXO) Phosphoproteins

The transcriptional activity of Forkhead Box O3 (FOXO3a) is inactivated by AKT-mediated phosphorylation on Serine 253 (S253), which enables FOXO3a binding to 14-3-3. Phosphorylated FOXO3a binding to 14-3-3 facilitates the nuclear exclusion of FOXO3a, causing cancer cell proliferation. The FOXO3a/14-3-3 interaction has, therefore, emerged as an important therapeutic target. Here, we report a comprehensive analysis using fluorescence polarization, isothermal titration calorimetry, small-angle X-ray scattering, X-ray crystallography, and molecular dynamics simulations to gain molecular-level insights into the interaction of FOXO3apS253 phosphopeptide with 14-3-3ε. A high-resolution structure of the fluorophore-labeled FOXO3apS253:14-3-3ε complex revealed a distinct mode of interaction compared to other 14-3-3 phosphopeptide complexes. FOXO3apS253 phosphopeptide showed significant structural difference in the positions of the -3 and -4 Arg residues relative to pSer, compared to that of a similar phosphopeptide, FOXO1pS256 bound to 14-3-3σ. Moreover, molecular dynamics studies show that the significant structural changes and molecular interactions noticed in the crystal structure of FOXO3apS253:14-3-3ε are preserved over the course of the simulation. Thus, this study reveals structural differences between the binding to 14-3-3 isoforms of FOXO1pS256 versus FOXO3apS253, providing a framework for the rational design of isoform-specific FOXO/14-3-3 protein–protein interaction inhibitors for therapy.


INTRODUCTION
The Forkhead box (Fox) class O transcription factor (FOXO) family contains four members: FOXO1, FOXO3a, FOXO4, and FOXO6, sharing a high homology between each other. 1,2 Among these FOXO proteins, FOXO3a is involved in cellular apoptosis, proliferation, cell cycle progression, DNA damage, and tumorigenesis 3,4 and is widely implicated in breast, liver, colon, prostate, bladder, and nasopharyngeal cancer. 5 FOXO3a is a multidomain protein made up of an N-terminal DNA binding forkhead domain (FKH), nuclear localization signal (NLS), nuclear export signal (NES), and C-terminal transactivation domain (TAD) ( Figure 1A). The transcriptional activity of FOXO3a is regulated by various post-translation modifications such as phosphorylation, acetylation, ubiquitination, and methylation. 1,6 The phosphatidylinositol 3-kinase (PI3K) and protein kinase B (AKT) pathways induce phosphorylation at three key residues in FOXO3a (T24, S253, and S315), which contain an AKT consensus phosphorylation motif (RxRxxS/T, where x is any residue) in all human FOXO proteins; however, FOXO6 lacks the Cterminal phosphorylation site (corresponding to S315 of FOXO3a) ( Figure 1A). 7,2 The S253 phosphorylation is critical for AKT-mediated nuclear exclusion of FOXO3a. 8,9 AKT-mediated phosphorylations enable FOXO binding to the dimeric form of the adapter protein 14-3-3 due to the overlap between the recognition motifs of AKT and 14-3-3 (RxxpS/TxP, where p is phosphate and x is any residue). This results in the translocation of FOXO from the nucleus, which subsequently inhibits the transcriptional activity of FOXO. 7 Moreover, the basic residues in the S253 phosphorylation motif, located at the end of the DNA binding forkhead domain, are important for DNA binding activity. 10 Binding of 14-3-3 to phosphorylated FOXO3a is suggested to mask the nuclear localization sequence or DNA binding site and, thereby, inactivate the FOXO3a transcription activity. 11,12 Modulation of the FOXO/14-3-3 association is a potential anticancer strategy because FOXO is typically inactivated in human cancers by phosphorylation rather than mutation, and its reactivation inhibits cancer outgrowth. 5,13−15 14-3-3 proteins are a family of evolutionarily conserved modulator proteins that bind to specific Ser/Thr-phosphorylated motifs on target proteins. In humans, there are seven 14-3-3 isoforms designated by Greek letters (β, γ, ε, η, σ, τ, and ζ), which are known to form homo-and heterodimers. 16,17 Each 14-3-3 protomer harbors the phosphopeptide binding groove. The seven 14-3-3 isoforms share 65−88% of sequence identity between them, but the phosphopeptide binding regions are highly conserved across the different isoforms. 17 The highly conserved phosphomotif binding groove, as well as overlapping and redundant functions, made it difficult to exploit 14-3-3s for therapeutic intervention. Despite their similarities, the isoformspecific functions and isoform-specific binding to a particular phosphomotif have been established in many instances. 18−23 In a recent study, 13 a biodiversity peptide library screening identified the hit peptide probe 9J10 that modulates the FOXO3a signaling pathway by interacting with the 14-3-3 protein family. The pull-down and mass spectrometry analysis revealed 14-3-3ε is one among the abundant 14-3-3 proteins interacting with 9J10. Moreover, there is a sharp reduction in the interaction between 14-3-3ε and FOXO3a in cells expressing the 9J10 probe signifying FOXO3a interaction with 14-3-3ε. In this study, the binding of the FOXO3a pS253 phosphopeptide motif ( 248 RRRAVpSMDNSN 258 ) to 14-3-3ε was measured by fluorescence polarization (FP) assay and isothermal titration calorimetry (ITC) and further structurally characterized using small-angle X-ray scattering (SAXS), X-ray crystallography, and molecular dynamics (MD) simulations. The FOXO3a pS253 :14-3-3ε complex structure solved at 1.85 Å resolution revealed a distinct structural feature compared with the highly similar FOXO1 p S 2 5 6 phosphopeptide ( 251 RRRAApSMDNNSK 262 ) bound to 14-3-3σ. This study reports the structural evidence of 14-3-3 isoform-specific binding of highly similar phosphopeptide motifs and provides a framework for the rational design of inhibitors with enhanced isoform selectivity against FOXO/14-3-3 protein−protein interactions.

Binding of FOXO3a Phosphopeptide to 14-3-3ε.
Fluorescence polarization assay was used to quantify the affinity of the tetramethylrhodamine (TAMRA)-labeled FOXO3a pS253 phosphopeptide (TAMRA-248 RRRAVpSer-MDNSN 258 ) to the wild-type (WT) 14-3-3ε protein. The phosphorylated FOXO3a pS253 peptide binds with a K d of 5.4 ± 0.7 μM ( Figure 1B). The affinity of the unlabeled FOXO3a pS253 phosphopeptide ( 248 RRRAVpSerMDNSN 258 ) to the 14-3-3ε protein was measured using isothermal titration calorimetry (ITC) ( Figure 1C), resulting in a comparable K d of 6.7 ± 1.8 μM, validating no significant contribution of TAMRA on the peptide binding affinity. The calculated binding enthalpy ΔH, the K d , the stoichiometry, n, Gibb's free energy, ΔG, and entropy, ΔS, were reported from the ITC experiment (Table S1). Our K d determined by FP and ITC is consistent with those of previously reported 14-3-3ε− phosphopeptide interactions. 13,22 2.2. SAXS Analysis of the 14-3-3ε:FOXO3a p253 Complex. The SAXS patterns of WT 14-3-3ε and its complex with a 6-fold molar excess of FOXO3a pS253 phosphopeptide (TAMRA: 248 RRRAVpSMDNSN 258 ) were measured in order to determine the solution size and shape ( Figure 1C). Inspection of the low angle of the Guinier plots reveals good data quality and no protein aggregation ( Figure 1D, inset). 14-3-3ε has a radius of gyration (R g ) of 32.23 ± 0.37 Å and a maximum dimension (D max ) of 99 ± 2 Å ( Figure S1). The R g values calculated from the real-space P(r) and the Guinier  Table 1). The estimated R g , D max , and molecular mass of the full length indicate that 14-3-3ε is dimeric at the concentrations used. The 14-3-3ε contains a flexible C-terminus, which is missing in the available crystal structures, which could be the cause of the increase in the solution R g and D max values compared to those of the dimeric crystal structure of the 14-3-3ΔC protein (R g value of 28.9 Å and D max of 85.1 Å). 17 The ab initio solution shape of 14-3-3ε was reconstructed and had a good fit to the experimental data with a discrepancy of χ 2 = 1. Ten independent reconstructions produced a similar envelope with a normalized spatial discrepancy (NSD) of 1.3 ± 0.07, and the average structure is shown in Figure 1E. The known atomic structure of 14-3-3εΔC (PDB ID: 2BR9) 17 was well positioned inside the lowresolution solution model, demonstrating a match between the solution and crystal structure. The SAXS derived parameters R g , D max , and molecular mass of 14-3-3ε mixed with a 1:6 molar ratio of FOXO3a pS253 reveal no significant oligomerization and conformational changes in 14-3-3ε upon addition of the phosphopeptide (Table 1).

Overall
Structure of the 14-3-3ε and FOXO3a pS253 Phosphopeptide Complex. The crystal of TAMRA-labeled FOXO3a pS253 phosphopeptide bound to 14-3-3εΔC (aa 1-232) diffracted up to 1.85 Å resolution and belongs to the C222 1 space group. The complex structure was solved, and the model was built and refined to a final R factor and R free of 16.55% and 18.88%, respectively. The data collection, processing, and refinement statistics are given in Table 2. The asymmetric unit contains one monomer, and the biological dimer of 14-3-3ε is formed with a symmetry-related molecule generated by the crystallographic 2-fold axis. Dimeric 14-3-3ε displays the characteristic W-like shape with each of the monomers harboring one FOXO3a pS253 phosphopeptide ( Figure 2A). We found an interpretable electron density for 231 residues of 14-3-3ε and 8 out of 11 residues of FOXO3a pS253 phosphopeptide ( Figure 2B). In addition, we observed a clear electron density for the TAMRA fluorophore that is attached to the N-terminus of the FOXO3a pS253 phosphopeptide for estimating the binding affinity of the peptide using the FP assay ( Figure 2B). The monomeric 14-3-3ε structure is composed of nine helices (αA-αI) with N-and C-terminal domains spanning between the αA-αF and αG-αI helices, respectively ( Figure 2A). The FOXO3a pS253 phosphopeptide motif ( 248 RRRAVpSerMDNSN 258 ) binds in the 14-3-3ε canonical peptide binding groove between the N-and Cterminal domains (Figure 2A). The phosphorylated residue FOXO3a pS253 interacts with the conserved binding site residues Lys50, Arg57, Arg130, and Tyr131 of 14-3-3ε. Additionally, salt bridge interactions were observed for FOXO3a R249 and FOXO3a D255 with Glu183 and Lys50 of 14-3-3ε, respectively ( Figure 2C and Table S2). The main chain atoms of FOXO3a V252 and FOXO3a M254 are stabilized by hydrogen bonded interactions with Asn227 and Asn176, Lys50 of 14-3-3ε, respectively (Table S1). However, the hydrophobic side chains of FOXO3a V252 and FOXO3a M254 are surrounded by the hydrophobic residues Leu230, Leu175, Leu223, Ile220, and Leu219 ( Figure 2D). Interestingly, the FOXO3a R249 side chain is pointed toward FOXO3a pS253 and makes an intramolecular salt bridge interaction ( Figure 2D). The TAMRA moiety of the FOXO3a pS253 peptide does not make any significant interaction with its 14-3-3ε or FOXO3a pS253 peptide residues in the asymmetric unit. It instead makes extensive interactions with TAMRA, Gln68 of 14-3-3ε, and the FOXO3a R248 moiety from the neighboring symmetry-related molecule in the crystal lattice ( Figure S2). However, the SAXS studies demonstrated that the TAMRA-labeled FOXO3a peptide does not affect the overall conformation and the oligomerization of 14-3-3ε in solution, thus precluding any TAMRA induced artifacts (Table 1).

Structural
Differences between the FOX-O3a pS253 :14-3-3ε and FOXO1 pS253 :14-3-3σ Complexes. Three distinct phosphopeptide binding motifs for 14-3-3 have been identified and defined as mode 1 (RSxpS/T-xP; PDB ID: 1QJB), mode 2 (RxxxpS/T-xP; PDB ID: 1QJA) and mode 3 (pS/pTX 1−2 -COOH), wherein pS and pT are phosphorylated Ser and Thr, respectively, and X is any amino acid. 24,25 In mode I and mode 2, proline is preferred at position +2 with respect to the pSer or pThr, and arginine, at −3 and −4 in mode 1 and mode 2, respectively. In mode III, the pSer or pThr is the penultimate residue in the C-terminal tail of a binding partner. In addition, the recently solved structure of FOXO1 pS256 bound 14-3-3σ ( 252 RRRAApSMD 258 ; PDB ID:  6QZS), 26 which is highly similar to FOXO3a pS253 , was available for structural comparison. In the FOXO1 pS256 :14-3-3σ complex structure, the residues from Arg 252 to Asp 258 spanning from the −4 to +2 positions relative to pSer are built into the electron density out of the 12 amino acids of the FOXO1 pS256 peptide ( 251 RRRAApSMDNNSK 262 ). Here in the FOXO3a pS253 :14-3-3ε structure, the residues from Arg 248 to Asp 255 spanning from the −5 to +2 positions relative to pSer are built into the electron density from the 11 amino acid FOXO3a pS253 peptide ( 248 RRRAVpSerMDNSN 258 ). Therefore, the 14-3-3 binding motifs in both the FOXO3a pS253 and FOXO1 pS256 peptides are highly similar in length and sequence, except there is one amino acid where FOXO3a V252 is replaced by FOXO1 A255 at the −1 position relative to pSer. The FOXO3a pS253 and FOXO1 pS256 motifs are not positioned at the C-terminus of the full-length proteins and have the canonical binding sequence (RRRAV/ApSMD), thereby representing either mode 1 or 2 in this case. The mode 1 and mode 2 motifs have invariant Pro at the +2 position with respect to pSer, while FOXO3a/ FOXO1 has negatively charged Asp at the +2 position with respect to pSer. The presence of a Pro at the +2 position is known to introduce a kink in the mode 1 and 2 motifs, which leads the peptide chain to exit the central binding channel and cross the 14-3-3 monomer in the direction of the adjacent monomer ( Figure 3A and B). Interestingly, the FOXO3a/ FOXO1 phosphopeptide Asp at the +2 position makes the kink, which enables the peptide chain to exit from the binding groove and is more similar to that of the mode 2 peptide motif than the mode 1 motif (Figure 3). A significant difference is observed in the side chain conformation +2 Asp between the FOXO3a pS253 and FOXO1 pS256 motifs. The negatively charged carboxylic acid group of the +2 Asp is pulled toward the positively charged phosphopeptide binding pocket, makes a salt bridge interaction with the phosphopeptide binding residue, Lys50, and is the cause for the kink to exit the central binding channel in the FOXO3a pS253 peptide ( Figure 3C). In contrast, the +2 Asp is pointed outward and solvent-exposed in the FOXO1 pS256 peptide ( Figure 3C). Another interesting feature in the FOXO3a pS253 peptide structure is the salt bridge interaction between the pSer phosphate oxyanion and Lys50, which is a highly conserved residue across the different isoforms of 14-3-3. Lys50 interacts with the pSer phosphate oxyanion in the mode 1 and mode 2 motifs, but this interaction is missing in the so far known structures of 14-3-3ε and FOXO1 ( Figure S2). The Lys50 side chain is either disordered (PDB ID: 2BR9) 17 or oriented outward (PDB ID: 6EIH) 27 or occluded (PDB ID: 3UAL) 22 to prevent the interaction with the pSer phosphate oxyanion ( Figure S3). However, the carboxylic acid group of Asp +2 stabilizes the Lys50 position similar to the mode 1 and mode 2 motifs and preserves the Lys50 interaction with the phosphopeptide motif. Moreover, a subtle difference is observed in the side chain conformation of the +1 Met position between the FOXO3a pS253 and FOXO1 pS256 peptides ( Figure 3C). The FOXO3a Met side chain is well surrounded by the hydrophobic residues Gly172, Leu175, Ile220, and Leu223.  A second interesting aspect of the FOXO3a/FOXO1 phosphopeptide is the position of arginines, as both the phosphopeptides contain consecutive Arg's at the −3, −4, and −5 positions with respect to pSer, whereas the mode 1 and mode 2 peptides have Arg at either −3 or −4 with respect to pSer and have been demonstrated to be critical for binding. 28 In the mode 1 structure, −3 Arg enables a kink in the phosphopetide to exit the central pocket in the direction of the adjacent 14-3-3 monomer, but there is no clear electron density to locate its side chain ( Figure 3A). Similarly, −3 Arg introduces a kink in the FOXO1 peptide to exit the central binding groove while the peptide chain returns into the pocket at −4 Arg and there is no clear side chain electron density for both Arg's at the −3 and −4 positions with respect to pSer ( Figure 3C and E). 26 In contrast, in the mode 2 peptide, a large kink is introduced by −4 Arg, and its side chain folds back to point toward the pSer forming an intramolecular salt bridge between the guanidine group of Arg and the pSer phosphate oxyanion ( Figure 3B). In the FOXO3a phosphopeptide, the backbone positions of the −3 and −4 Arg residues of FOXO3a are significantly different compared to those from the mode 1, mode 2, and FOXO1 peptides ( Figure 3A−C). The difference is mainly due to the fact that the mode 2 peptide has the aromatic residue Tyr at the pSer −2 position ( 3 RLYHpSLPA 10 ) instead of Ala at pSer −2 in FOXO3a pS253 . The side chain of −2 Tyr occupies the central binding groove, and −3 Arg exits the central binding groove ( Figure 3B). However, the FOXO3a −3 Arg continues to lie in the central binding groove, and only −4 Arg makes a large kink and its side chain points toward the pSer to form an intramolecular interaction.
In both the mode 2 and FOXO3a structures, the −4 Arg guanidino group is located at the same position, and Arg is further stabilized by an additional salt bridge between the δ NH of the arginine guanidino group and the carboxyl side chain of the highly conserved Glu180 ( Figure 3B). We observed a similar intramolecular interaction between −2 Arg and the pSer of 9J10 phosphopeptide (LNTRPGRRRNpSN) in our recently reported structure ( Figure 3D; PDB ID: 7C8E). 13 These results show that Arg at the −2, −3, or −4 position brings significant structural differences in the binding mode of phosphopeptides.

CONCLUSIONS
In summary, we have characterized the binding mechanism of the FOXO3a pS253 phosphopeptide to 14-3-3ε using FP, ITC, SAXS, X-ray crystallography, and MD simulation studies. The cocrystal structure demonstrated that the TAMRA-labeled FOXO3a pS253 phosphopeptide binds to 14-3-3ε in a distinct manner compared to that of the highly homologous FOXO1 pS256 peptide binding to 14-3-3σ. Further, the FOXO3a and FOXO1 phosphopeptides resemble the mode 2 and mode 1 peptide motifs for binding to 14-3-3, respectively. Though closely resembling mode 2, FOXO3a showed a significant difference in the backbone positions of the −3 and −4 Arg residues compared to those from a canonical mode 2 peptide. Moreover, MD simulation studies confirmed that the structural changes noticed in the crystal structure of FOXO3a pS253 :14-3-3ε and the observed molecular interactions are preserved over the course of the simulation. Together, our data illustrate that the highly similar FOXO phosphopeptides adopt a distinct binding interaction with the conserved phosphomotif binding groove in 14-3-3s. Previous studies strongly support the notion of the 14-3-3 isoform-specific binding of phosphopeptide. 18−23 Interestingly, a recent study suggests the hierarchized profile in binding affinities of 14-3-3 isoforms to the phosphopeptide. 29 Probing the FOXO3a−14-3-3 interaction has potential therapeutic significance because FOXO3a is frequently inactivated in many cancers by Aktmediated phosphorylation and subsequent 14-3-3 binding. 5,30 The observed isoform-specific binding lays the groundwork for developing protein−protein interaction inhibitors with enhanced isoform selectivity against FOXO:14-3-3 complexes.

Cloning, Expression, and Purification.
The WT 14-3-3ε and 14-3-3εΔC (aa 1-232) genes were PCR amplified and then inserted into the pET28a vector (Novagen 69864-3) containing an N-terminal 6xHis tag with a thrombin cleavage site, using the BamHI and XhoI restriction enzyme sites. The plasmid was later expressed in E. coli BL21 DE3 cells. The culture was grown in Luria broth media containing kanamycin (50 μg/mL) at 37°C. The cells were induced with 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) once the O.D.600 nm reached 0.6 and continuously grown for 4 h at 37°C. The cells were harvested by centrifugation and resuspended in the buffer containing 50 mM Tris pH 7.5, 300 mM NaCl, 1 mM DTT, 0.1 mM PMSF, and 1 protease inhibitor tablet followed by lysis by sonication. The lysate was centrifuged, and the supernatant was passed through a 5 mL Ni-NTA IMAC column, followed by washes with 50 mM Tris pH 7.5, 300 mM NaCl, 1 mM DTT, and 0−20 mM imidazole. The bound protein was eluted using 50 mM Tris pH 7.5, 300 mM NaCl, 1 mM DTT, and 500 mM imidazole. The affinity purified protein was then injected into a S75 Superdex size exclusion column. For crystallization, the 14-3-3εΔC protein was subjected to thrombin digestion to remove the 6xHis-tag prior to purifying using a size exclusion column.

Fluorescence Polarization (FP).
FP assay was performed using black 96 well plates (LVL Technologies). In direct binding FP experiments, 15 nM TAMRA-conjugated synthetic FOXO3a pS253 phosphopeptide (Designer Bioscience Ltd., UK) was mixed with increasing concentrations of WT 14-3-3ε protein (from 0.002 μM to 63 μM) in triplicate. The assay volume was 30 μL per reaction with 3× working concentrations of protein and peptide, 10 μL each prepared in assay buffer (10 mM HEPES pH 7.5, 150 mM NaCl, 0.5 mM DTT, and 0.005% Tween 20). The plate was incubated for 10 min at room temperature, and relative fluorescence was measured with the excitation wavelength of 530 nm and an emission wavelength of 610 nm using the TECAN infinite M1000 Pro microplate reader (TECAN). The degree of polarization ratio of fluorescence intensities parallel and perpendicular was calculated by the reader using iControl software and expressed in millipolarization units (mP). The FP values were plotted against the log of the 14-3-3ε protein concentrations, and the binding affinity, K d was determined from the generated sigmoidal curve as analyzed in GraphPad Prism 6.

Isothermal Titration Calorimetry.
Wild-type 14-3-3ε was dialyzed in a buffer containing 20 mM HEPES, pH 7.5, 150 mM NaCl, and 1 mM DTT. Unlabeled FOXO3a pS253 phosphopeptide (GaloreTx Pharmaceuticals Pvt. Ltd., India) was dissolved in the same buffer as mentioned above. ITC experiments were performed using MicroCal ITC 200 at 25°C in high-feedback mode with stirring speed 800 rpm and filter period time 5 s. 50 μM FOXO3a pS253 phosphopeptide was titrated with 20 injections of 590 μM WT 14-3-3ε, with each 2 μL injection with a 150 s interval between each injection. A control experiment was conducted with identical buffer conditions and peptide in the cell replaced by buffer, and this control experimental result was used to subtract to account for the heat of dilution. The data was analyzed by origin 7 software, and the results were fit by using a one set of sites model.

SAXS Data Collection and
Analysis. SAXS data of WT 14-3-3ε were measured at 1.6 mg/mL in a buffer containing 50 mM Tris/HCl, pH 7.5, and 200 mM NaCl using a BIOSAXS-1000 Rigaku instrument with a microfocus X-ray generator (1.5418 Å wavelength). The measurement with FOXO3a phosphopeptide was performed with a 6-fold excess of the peptide. For each protein sample, data were collected for a corresponding buffer under identical experimental conditions, providing a background scattering to subtract from the scattering of the sample. All the data processing steps were performed using PRIMUS. 31 The low-q region of the Guinier plot was used to analyze for protein aggregation. The forward scattering I(0) and the radius of gyration, R g , were computed using the Guinier approximation assuming that at very small angles (q < 1.3/R g ) the intensity is represented as I(q) = I(0) exp(−(qR g )2/3). The distance distribution function P(r) of the maximum particle dimension, D max , and the radius of gyration, R g , were also computed from the extended scattering patterns using the indirect transform package GNOM. 32 GASBOR 33 was used to build ab initio low-resolution models considering low angle data (q < 2 nm −1 ). The ab initio solution shape was obtained by superposition of 10 independent model reconstructions with the program SUPCOMB 34 and then averaged using DAMAVER. 35 The ab initio model was superposed with the X-ray structure using SUPCOMB. 4.5. Crystallization and Structure Determination. The purified, 5 mg/mL 14-3-3εΔC protein was incubated with the FOXO3a pS253 peptide (TAMRA-RRRAVpSMDNSN) at a ratio of 1:1 (0.15 mM:0.15 mM) supplemented with 20 mM HEPES, pH 7.5, 0.05 mM EDTA, and 1 mM MgCl 2 . The protein mixture was screened using Hampton and Jena Bioscience (Basic and Classic) condition kits using sitting drop at 16°C. Diffracting quality crystals were later obtained using the hanging drop vapor diffusion method at 16°C, in conditions containing 30% MPD, 10% PEG 4000, 0.1 M imidazole HCl, pH 8.0. These crystals were flash frozen in liquid nitrogen and shipped to the PXBL-21 beamline of the Indus-2 synchrotron facility (RRCAT, Indore, India). 36 The diffraction data were collected up to 1.85 Å at 100 K using an X-ray beam with a wavelength of 0.97 Å and a MARCCD 225 mm detector. The diffraction images were processed and integrated using iMOSFLM 37 and later scaled using the Scala program 38 (Evans, 2006) in the CCP4 suite. 39,40 The initial structure solution of the protein was obtained by the MR method using the Phaser MR module, 41 with the monomer unit of 14-3-3εΔC as the starting model (PDB ID: 2BR9). 17 Further model building was done in Coot 42 followed by refinement using the REFMAC5 program. 43 The atomic coordinate and structure factor have been deposited in the RCSB PDB with accession code 7V9B. 4.6. Molecular Dynamics (MD) Simulations. MD simulations lasting 100 ns were performed for complexes using the Amber99SB force field of the Amber18 package. 44 Hydrogens were added, and the protonation state of ionizable groups was chosen to pH 7.0. Missing residues were added, and the terminal residues were patched with acetyl at the Nterminal and an N-methyl group at the C-terminal of the peptide and protein, as they were not the natural termini. The parameters for pSer were obtained as described in ref 45. Energy minimization and MD simulations were carried out with PMEMD MD. Minimization was performed by the steepest descent method for the first 3000 steps followed by 2000 steps of the conjugated gradient method. First, the two ACS Omega http://pubs.acs.org/journal/acsodf Article peptides within the complex were minimized with the restraints of a harmonic force constant of 100.0 kcal/(mol· Å 2 ) applied to all protein heavy atoms. The protein and peptide were then minimized with no restraints. Each minimized complex was inserted in a water box of TIP3P water that extended 10 Å from the complex and neutralized by adding counterions. Water molecules alone were minimized while restraining all protein and peptide heavy atoms. Finally, all restraints were removed, and the total system (protein, peptide, and solvent) was minimized. After energy minimization, the system was slowly heated from 0 to 300.0 K over 100 ps in the NPT ensemble under 1 atm of pressure. Equilibration was done in three steps of 100 ps each starting with 10.0 kcal/ (mol·Å 2 ) harmonic restraint weight followed by 1.0 kcal/(mol· Å 2 ) on all heavy atoms and ending with a final equilibration with no restraints. Production runs of 100 ns were carried out without any restraint at a temperature of 300.0 K and a pressure of 1 atm using a time step of 2 fs. A 12.0 Å cutoff was used as nonbonded cutoff, and bond lengths involving bonds to hydrogen atoms were constrained using the SHAKE algorithm. 46 The coordinates were stored every 20 ps. RMSD/RMSF was measured by utilizing the Amber analyses tool PTRAJ. The hydrogen-bonding analysis tool incorporated in the VMD software 47 was utilized to calculate the hydrogen bond interaction analyses of the MD trajectories.
■ ASSOCIATED CONTENT