Selection of the peptide sequences followed by the docking studies
The TNFα structure was selected for peptide selection and docking studies as shown in Fig 1A, B.A single monomer subunit of the trimer was used for further studies. The potential binding site was extracted, using Acsite software (Fig 1C). This cavity was used in the docking studies of the peptides as described below. After docking was done with two different docking algorithms namely, Flex-X [5] and Genetic Optimisation for Ligand Docking (GOLD) [6], the interactive maps of the docked molecules were made in LIGPLOT [7] and the results were then analyzed to find the peptides that could bind with higher affinity to the TNFα and nine of the peptides were docked with varied binding energies. To confirm the results, docking was repeated with increased stringency by using a reduced cavity size as a receptor, which resulted in seven matching peptides which were again modeled with TNFα (Protein Data Bank, PDB ID: 1TNF) as a template using MODELLER software [8]. The scoring function, X-score was used to compare the binding energies [9]. The peptides were selected after visualizing them in Insight II. The peptides docked in a similar fashion, with binding energy within a cutoff of ±2.2 kcal/mol,and similar interacting residues of the receptor protein were shortlisted for solid-phase peptide synthesis. The selected hydrophobic peptide “PIYLGGVFQ” (after in vitro assays) as shown in Fig 1D was docked with the TNFα using Flex-X and GOLD (Fig 2A-D), and the interaction maps obtained from LIGPLOT(Fig 2E).
The docking of peptides was also done through the Docking INCrementally (DINC) server [10] using peptide and TNF-alpha (PDB ID: 2AZ5) [11] as shown in Fig 2F. The peptide and protein structures were prepared in the Chimera program. For this molecular docking, the grid center X= -19.2, Y=74.5, Z= 33.8 Å, and box size 24X24X24 were selected around the co-crystalized ligand of the TNFα.
Molecular dynamics (MD) simulation analysis
To validate the stability of the peptide inside the active site of TNFα, molecular dynamic simulation studies were done, using Assisted Model Building with Energy Refinement (AMBER)-12 software [12]. It provided information about the important binding interactions and free energy of the ligand, as well as the contribution of each residue to total binding interaction energies (Table 1). Protein structure (PDB ID: 2AZ5) has four chains: A, B, C, and D, where chain A-B isidentical to chain C-D of the structure. Each A-B and C-D chains have the same bound co-crystallized ligand (307) in PDB (named as the standard, STD). As per some reports, A and B chains were used from PDB: 2AZ5 for computational studies. Therefore A-B chain was used in to dock the peptide into the active site of TNFα [13]–[15]. Similarly, inMD simulation studies, chain A-B from the above structure was used as the peptide was docked into the active site of TNFα [11]. Since it is a structure obtained by X-ray diffraction studies of ligands bound to TNFα, STD was used as a reference molecule. This facilitated our understanding of the peptide binding mechanism (H-bonding, stability, and energy score).
Table 1. The binding energy between TNFα and ligands (Std and Peptide) interactions, as derived from the Molecular Dynamics studies.
Energy Component (in kcal/mol)
|
STD
|
Peptide
|
EVDWAALS
|
-47.66
|
-68.88
|
EEL
|
-5.60
|
-86.73
|
EGB
|
15.80
|
127.96
|
ESURF
|
-5.85
|
-10.57
|
ΔGgas
|
-53.25
|
155.60
|
ΔGsolv
|
9.95
|
117.39
|
ΔGbind
|
-43.31
|
-38.21
|
Initially, the molecular docking studies were performed to place the STD and peptide molecules into the active site of TNFα, followed by the MD simulation. For the preparation of ligands and protein, the initial parameters and topology files (frcmod, prmtop, inpcrd, etc.) for the ligands (STD and peptide), protein, and protein-ligand complexes (TNFα-STD and TNFα-peptide) were generated using Antechamber [16] and tleap program of AMBER software by implementing, General Amber Force Field (GAFF) [17] and Amber ff99SB force field [18]. An explicit TIP3P water model was used to solvate the systems [19], and the solvation box was extended to 20 Å in all directions of the cubic box (solute forming). The systems generated were minimized followed by gradual heating from 0 to 300K, and density equilibration was done under the NPT ensemble. Then a constant pressure equilibration of 1 atm pressure (pressure relaxation time of 2.0 ps) for 1 ns at 300 K was applied. Subsequently, a 20 ns production run was done under NPT ensemble with a non-bonded interaction cut-off distance (12 Å). So, long-range electrostatic interactions were subjected to Particle-Mesh Ewald (PME) method [20], then the bulk effect simulation was performed by allowing periodic boundary conditions. The relative binding free energy was calculated for TNFα-ligand complex formation using the molecular mechanics-generalized born surface area (MM-GBSA) method [21] on the last four ns trajectory acquired from MD simulations of the protein-ligand complex formation to ensure conformational sampling and to obtain reliable binding free energy values. The Root Mean Square Deviation (RMSD) and B-factor analysis were performed with the CPPTRAJ program in AMBER [22], and the hydrogen-bond occupancy analysis using the visual molecular dynamics (VMD) package [23].
After placing the known standard (STD) and our experimental peptide into the TNFα active site, the binary complexes (TNFα-STD and TNFα-peptide) formed after docking studies were subjected to MD simulation analysis for 20 ns. This was performed for the evaluation of the binding affinity of various ligands with the TNFα under multiple dynamic conditions. The study of the whole-system RMSD, ligand RMSD, protein-backbone RMSD, and B-factor/atomic fluctuations wasdone. The entire protein, and the protein backbone RMSD data, indicated that both binary systems were stable during the MD simulation study (Fig. 3A and B). The protein backbone and the whole protein were completely stable in binary complexes during the 20 ns simulation runs. The ligands were also found stable during the simulations especially after the 5 ns runs. The B-factor/atomic fluctuation data suggested that the fluctuations of atoms of the protein were very less during the simulations and most of the atoms/residues involved in atomic displacement were part of the loop region only (Fig. 3C and D).
The RMSD deviations for the receptor in both the binary complexes were ~0.75 Å and between the binary complexes was ~0.91 Å which are within the range of acceptability (0 ≤ 2 Å) and confirmed protein structure stability and shape in the dynamic state (Fig. 3A). The important active site residues (IleA58, LeuA57, TyrA59, SerA60, LysA98, ProB117, IleB118, TyrB119, GlyB121, GlyB122, TyrB151, etc.) were found to be identical in the binary complexes. The major RMSD deviation was observed for SerA86, ProB70, ArgB103, and minor RMSD deviation was noticed for GluA23, SerA71, CysA101, GluB23, SerB86, etc. The amino acid residues which were deviated in RMSD analysis were part of loop regions and found to be far from the core region of the active site. It was observed that no major RMSD deviations of ligands structure, either before or after the MD simulation studies were found in both binary complexes and were found to be less than 2 Å. This minor deviation shows that the ligands remained stable and consistently bound inside the active site of TNFα throughout the simulation runs.
It is well known that the active site of TNFα is completely hydrophobic and crystalized ligand (STD) also does not possess any hydrogen bond interactions with the key active site residues [11]. A similar observation was noticed after the 20 ns MD simulation runs with STD. The designed peptide possesses various amino acids in its structure and due to the availability of various H-bond donors as well as acceptor centers, multiple H-bond interactions were observed during the molecular docking analysis (H-bond interactions with Lys98B, Pro117B, Tyr119B, and Gly121B) and similarly during the MD simulation. The H-bond occupancy results suggested that the H-bond interactions between protein and peptide were very weak. The per-residue decomposition analysis data suggested that Tyr119A, Leu57B, Tyr59B, Tyr119B, and Leu120B have a major contribution to the binding of STD in the active site of TNFα (Fig. 3D). While Tyr151A shows a slightly negative impact on STD binding. The residues Tyr59A, Tyr119A, Tyr151A, Leu57B, Tyr59B, Tyr119B, Leu120B, and Gly122B also show a greater contribution to the binding of peptides in the active site of TNFα. Figure 3E shows the binding poses of the STD with TNFα (PDB: 2AZ5) while the peptide docked into the active site has been depicted in Fig 3F.
The ΔG bind or the binding free energy contribution reported inhibitor (307/STD) demonstrated a strong binding affinity (-43.31 kcal/mol) while the designed peptide exhibits a comparable but lower binding affinity (-38.21 kcal/mol) for the TNFα protein.
Synthesis and purification of peptide inhibitor for human TNFα protein
The synthesis was done by solid-phase peptide synthesis (SPPS) process using the Fmoc strategy. The nature of the peptides essential for their purification using the peptide calculator (http://www.innovagen.se/ custom-peptide-synthesis/ peptide-property-calculator/ peptide-property-calculator.asp) showed that the peptide is highly hydrophobic (Fig 1D). After SPPS, they were cleaved and desalted on the LH-20 column. Then the peptides were dissolved in methanol and concentrated with rotavapor followed by lyophilization. The yield of the peptide was: between 60 and 80 mg. The peptides were purified using reversed-phase HPLC (Fig S2C). Briefly, 25 μl of the peptide with appropriate dilution was injected and resolved with the mobile phase (methanol). The identity and molecular weights of the peptides were confirmed by mass spectrometry. The micromolar quantity of the peptide was dissolved in 0.1% TFA. The solution was mixed with an equal amount of matrix 4-hydroxy cyanocinnamic acid (HCCA), prepared in acetonitrile and 0.1% TFA (1:2). The matrix-analyte sample (1 µl) was loaded on the steel plate and left for crystallization for 5 min. The compound was then analyzed by MALDI-TOF/TOF mass spectrometer. The mass of all the peptides matched the calculated mass and the peak was clean and single (Fig.S2B).
Cytotoxicity assay
The Wehi-164 murine fibrosarcoma cell line undergoes apoptosis by TNFα stimulation in the presence of actinomycin D (AcD), which makes it suitable for suchan experiment [25]. On the other hand, immune cells such as J774 macrophages, Jurkat cells, etc., show cell proliferation in response to TNF-α. So, in this study, 70-80% confluent Wehi-164 cells,with more than 90% viability were used.In 96-well plates, 5x104 cells per wellwere taken in the presence or absence of AcD, 0.2μg/ml, TNFα (100 ng/ml), peptide,andsuramin (100 and 200 μM), premixed at 4°C for 1h. After 20h incubation at 37°C, 25 µl of MTT solution (5 mg/ml in 1X PBS) was added and incubated for 4 h. Then isopropanol (with 0.04 M HCl) was used to dissolve the formazan before measuring the OD at 570 nm with background subtraction at 650 nm (Molecular Devices, USA) to determine the TNFα-induced cytotoxicity.
Preparation of nuclear extract
The cell nuclear protein was extracted using the method of Schreiber and coworkers [26] with a minor modification. Briefly, A549 cells at a concentration of 106per well were incubated with or without TNFα (100 ng/ml) and different concentrations of peptide or suramin,for 45 min followed by washing with PBS. These treated cells were then resuspended in the cell lysis buffer [10 mM KCl, 10 mM HEPES (pH 7.5), 0.1 mM EDTA, 0.5% NP40, 1 mM DTT, and 0.5 mM PMSF with mammalian protease inhibitor or PI (Sigma)]. The cells were kept in ice for 25-30 minutes with intermittent mixing, then vortexed followed by centrifugation at 12,000 g for 10 min at 4°C. The pellet containing the nuclei was washed with the cell lysis buffer before resuspending in the nuclear extraction buffer (0.4 M NaCl, 20 mM HEPES with pH 7.5, 1 mM EDTA, 1 mM PMSF, 1 mM DTT with PIsolution) and incubated in ice for half an hour. The nuclear protein extract was then collected as the supernatant after centrifugation (12,000 g at 4°C for 15 min).
Electrophoretic Mobility Shift Assay (EMSA)
A dsDNA probe for NFκB (Promega) was used for the gel shift assay after labeling with radioactive ATP [γ-32P] by T4 polynucleotide kinase using the manufacturer’s instructions. The reaction mixture contained 2 μl of 5x binding buffer (50 mM Tris–HCl (pH 7.5), 20% glycerol, 250 mM NaCl, 5 mM MgCl2, 2.5 mM DTT, 2.5 mM EDTA, and 0.25 mg/ml poly-dI-dC) and 2 μg of nuclear extract and nuclease-free water. Aprobe (1 μl, 106cpm) was added to initiate the reaction for 30 min at room temperature. The protein binding specificity with the DNA was assessed by competition reactions, where 20-fold molar excess of oligonucleotides (unlabeled) was added to each tubeprior to addinga radio-labeled probe. All the prepared samples were subjected to native gel electrophoresis (4% polyacrylamide gel). The radiographic gels were scanned (Fuji FLA-2000 Phosphor imager) for further analysis.
Western blot analysis
The nuclear protein extracts from the TNF-stimulated A549 cells were subjected to SDS-PAGE and transferred to a nitrocellulose membrane (Millipore, USA) using the transfer buffer (192 mM glycine, 25 mM Tris, 20% methanol) for 20 min at 20 V semi-dry transfer apparatus (BioRad, USA). Blocking was done with 1% BSA for 1 hour at room temperature. After three times (5 min each) washing with 1X TBST, the membrane was incubated overnight with anti-NFκB (p65) and anti-β actin antibodies at 4ºC. It was then probed with anti-rabbit IgG HRP for 1 h after the TBST wash. The blot was washed again and developed using 4-chloronaphthol (Sigma). The band intensities were measured and the densitometric analysis was performed with the Alpha Digi Doc software tool (Alpha-Innotech Corporation, USA).
Immunocytochemical assay
Inhibition of the TNFα-induced NFκB activation resulting in NFκBRelA/p65 nuclear translocation in A549cells by the peptide was visualized by immunocytochemistry assays.The A549 is a human lung epithelial cell line that responds well to the stimulation by TNF-α resulting in the NFκB activation [25].The cells were seeded at a concentration of 2.5 x 104 cells/well and cultured overnight before incubating for 45 min with TNFα (100 ng/ml) with and without the peptide (premixed for 1 hour). The cells were then washed with isotonic buffer PBS followed by fixing and blocking (PBS with 0.05% sodium azide and 1% BSA) for 30 min and the cell permeabilization was done with Triton x-100 (1%) for 30 min. For the nuclear staining, 10 µl of propidium iodide (0.5mg/ml) was mixed in the well and incubated at 4°C for 30 min in dark.
The cell staining for p65 was done using the supplier’s protocol (Santa Cruz, USA). Briefly, cells were incubated in theprimary antibody against NFκB p65 (Santa Cruz, USA, sc-372)diluted in blocking buffer (overnight, 4°C), followed byFITC-conjugated secondary anti-rabbit IgG (1:200) was added to the cell and incubated for 1 hour. Cells were washed thrice with PBS and scanned under a fluorescence microscope (Nikon, Japan).
Inhibition of cell surface binding of TNF-α
The U937 cells were plated at a density of 2×104 cells per well (96-well plate) and stimulated with human TNFα (100 ng/ml) with and without the anti-TNFα peptide for 1 hour at room temperature. Then the cells were treated with the primary and FITC-labeled secondary antibodies as described above. The TNF-α binding to receptors was analyzed by fluorescence microscope and image processing was done with ImageJ (NIH).
Analysis of LPS-induced cell surface TNFα by flow cytometry
The peptide decreases the TNFα signaling through NFκB and therefore it should decrease surface TNFα as well. As an indirect approach to evaluate the effect of the peptide to quantitate surface TNFα with and without the peptide, the expression of surface TNFα was measured by flow cytometry. For this experiment, U937, a pro-monocytic, human myeloid leukemia cell line, was used because it expresses the receptors of TNF-α (TNFR1 and TNFR2) in a substantial number [24]. For all flow cytometry experiments, LPS (1 μg/ml) was dissolved in RPMI. The peptide was dissolved in DMSO (1% v/v) before diluting with RPMI. The U937 cells were transferred into microcentrifuge tubes at a concentration of 4×106/ml. The cells were treated with 1 μg/ml LPS with or without the peptide (200 μM) for 2 h. Cells were then resuspended in 200 µl FACS buffer. All the cells were labeled with monoclonal anti-human TNFα antibody (Peprotech, USA) and subsequently washed with 1X PBS followed by treatment with FITC-conjugated secondary anti-mouse antibody. After the wash flow cytometry was performed. A total of 50,000 events were collected using a FACS Calibur flow cytometry system (BD Biosciences, USA), and data were analyzed using the FlowJo program.
Animal experiments
Animal experiments were performed according to the institutional guidelines and the protocol was approved by the animal ethics committee of the National Institute of Immunology (NII), New Delhi. A total of 20 female DBA1/J, 6 to 8 weeks old mice were procured from the inbred facility of NII. Mice were housed in up to 6 mice per cage in a room maintained at 23 ± 2°C with 50 ± 10% humidity and 12-h light/12-h dark cycles. The animals were allowed free access to tap water and regular rodent chow.
Induction of CIA and treatment
Induction of CIA was done in accordance with our previous study [27], with some minor modifications. Briefly, the mice were acclimatized for 7 days and divided into four groups, healthy, CIA, CIA+peptide (5 mg/kg), and CIA+peptide(10 mg/kg). Thereafter, all mice except the healthy group were immunized by intradermal injections of an emulsion containing 100 μg of immunization grade bovine type II collagen (CII) in the Freund's complete adjuvant at the base of the tail. On day 21, a booster containing 100 μg CII emulsified with Freund's incomplete adjuvant was administered and observed for disease development. Treatment was initiated after 12 days of the booster dose when the disease severity of all groups was observed to be the maximum. The peptide was dissolved in DMSO (≤ 1% v/v) followed by dilution in 1X sterile PBS beforeadministering intraperitoneally (i.p) in the mice, three times a week. The healthy group was treated with a vehicle (1% DMSO in PBS (v/v)) alone. A schematic presentation of the disease induction and treatment protocol is shown in Fig 5A. All the groups of mice were regularly scored for disease indicators (inflammation, redness, severity, etc.) according to their severity until termination. All the reagents used for the generation of the CIA models were obtained from Chondrex (USA).
Clinical evaluation of the disease
To assess the disease characteristics, all the mice were scored for disease extent at 3-4 day intervals for ~20 days following the booster immunization as described in our previous study [28]. Both paws and joints of fore and hind limbs were scored as follows: 0: Normal; 1: redness and swelling of the paw or one digit; 2: swelling and redness of ankle and wrist; 3: severe redness and swelling of entire paw including digits and 4: severe arthritis of the limb including entire paw and digits and expressed as the mean of all 4 limbs of a mouse.
Micro-CT scanning and image processing
Micro-CT scanning of hind-limbs from healthy control, CIA, CIA+peptide (5mg/kg), and CIA+peptide (10mg/kg) groups were performed at 70kVp, 114µA by exposing the bones to radiation for 300ms. The medium used for scanning was formalin. Scanning at 10µm resolution was performed using a 0.1 aluminum filter to reduce the beam hardening effect. Once the scan was completed, the in-built reconstruction algorithm was used to develop the cross-sectional images. The range of global threshold, 654.2-1239.7 mg HA/ccm, was used to develop the 3D geometry of the bones.
Statistical analyses
The statistical data analyses were performed using GraphPad Prism software. Three experimental values(n) are given as mean ± standard deviation (S.D.) in the cell culture assays. The student’st-test was employed to calculate the differences from the respective controls for each paired experiment. The data from animal experiments have been expressed as mean ± SD (n = 5 or 6). As applicable, p-values were calculated both using unpaired t-tests, and one-way or two-way ANOVA tests; andp ≤ 0.05 was significant.