The C-terminal region of DKK1 cysteine-rich domain is crucial for interaction with LRP6 co-receptors:
The crystal structure of the LRP6 receptor with DKK1 protein (PDB ID: 3S8V) was retrieved from Protein Data Bank (https://www.rcsb.org/). The interface residues of the complex were visualized and interpreted using UCSF chimera. The FASTA sequence of LRP6 and DKK1 were retrieved and interface residues were identified from the software, after analyzing both a small stretch of the DKK1 C-terminal region amino acids (16 residues: Asp673 to Arg688) was found to be interacting majorly with LRP6 co-receptor (Table 1 and inner figure). Further, docking was proceeded with the LRP6 co-receptor with the 16 amino acid chain from the DKK1 protein using the Pydock server. A blind docking was performed without specifying the binding site of the Ligand during the docking process.
The results obtained were clustered and examined by comparing the binding energies and different docked conformations of the peptide generated at the end. Anticipatedly, the peptides were found to have interacted within the LRP6 co-receptor accurately at the exact place of the DKK1 protein binds (Fig. 2). It is well understood that the 16 amino acid sequences can hinder the binding of LRP6 with DKK1 protein independently. Then the peptide interface with the receptor was determined using SPPIDER and the interacting residues were recognized.
Table 1. Interacting residues at the interface of LRP6 and DKK1
LRP6 co-receptor
|
DKK1
|
R639 E663 D673 Y677 K684 I681 R688 Y706 E708 R751 W767 G768 G769 K770 P771 N789 V790 G791 R792 L810 D811 T812 N813 L814 D830 D831 P833 H834
|
Q184 E185 H204 F205 W206 S207 I209 G217 Q218 V219 T221 K222 R224 R225 S228 L231 E232 I233 F234 R236 Q248 S258 R259 L260
|
Computational Alanine Scanning Analysis Of The Peptide:
After analyzing the results of 16aa inhibitory peptide’s specificity and accuracy of the binding site and binding mode to the target protein, we aimed to get some combinations of the 16aa peptide and decided to select alanine scanning for this purpose. Alanine scanning of the 16aa inhibitory peptide was performed to identify the best alanine mutated inhibitory peptides, which interacted with the critical residues. These were implied in the association of LRP6 with the DKK1 protein, which thereby could allow us to shorten the non-essential portion of the peptide. The amino acids were mutated to alanine consecutively one by one using BIOVIA Discovery Studio 2020 software. Alanine substituted peptides were then modelled and minimized using the UCSF Chimera program. Subsequently, docking studies were performed for each mutated peptide with the LRP6 receptor, and the binding energies of the complex and their three-dimensional conformations of binding interfaces were obtained.
It can be observed from the values presented in Table 2, that the mutating residues Asp673 to Arg688 with alanine have no significant effect on the complex steadiness in terms of binding energy and also on structure confirmation. Moreover, we also calculated the molecular mechanistic values, such as electrostatics, desolvation, and VdW (Van der Waals) for each peptide. VdW is generally a weak force, which enhances the attraction and binding efficiency of a ligand to its receptor, and also it is concerned in the maintenance and stabilization of a drug-receptor complex.
Table 2. Results of alanine scanning and molecular docking of inhibitory and mutated peptides with LRP6
Code
|
Conf
|
Electrostatics
|
Desolvation
|
VdW
|
Total
|
Rank
|
16aa
|
DNRIYWTDISLKTISR
|
-11.990
|
-14.493
|
21.969
|
-24.286
|
6
|
16aa1
|
ANRIYWTDISLKTISR
|
-28.715
|
-2.747
|
41.939
|
-27.268
|
2
|
16aa2
|
DARIYWTDISLKTISR
|
-15.812
|
-7.565
|
-0.283
|
-23.405
|
11
|
16aa3
|
DNAIYWTDISLKTISR
|
-13.553
|
-14.045
|
48.960
|
-22.703
|
15
|
16aa4
|
DNRAYWTDISLKTISR
|
-17.432
|
-7.877
|
22.974
|
-23.011
|
14
|
16aa5
|
DNRIAWTDISLKTISR
|
-9.744
|
-18.103
|
57.904
|
-22.056
|
17
|
16aa6
|
DNRIYATDISLKTISR
|
-19.637
|
-6.813
|
39.339
|
-22.516
|
16
|
16aa7
|
DNRIYWADISLKTISR
|
-13.829
|
-12.985
|
17.198
|
-25.093
|
4
|
16aa8
|
DNRIYWTAISLKTISR
|
-12.543
|
-17.823
|
30.949
|
-27.271
|
1
|
16aa9
|
DNRIYWTDASLKTISR
|
-15.562
|
-13.317
|
30.765
|
-25.803
|
3
|
16aa10
|
DNRIYWTDIALKTISR
|
-8.805
|
-21.140
|
65.568
|
-23.388
|
12
|
16aa11
|
DNRIYWTDISAKTISR
|
-11.126
|
-15.316
|
31.547
|
-23.288
|
13
|
16aa12
|
DNRIYWTDISLATISR
|
-19.478
|
-10.312
|
60.557
|
-23.734
|
10
|
16aa13
|
DNRIYWTDISLKAISR
|
-16.639
|
-7.245
|
-2.822
|
-24.166
|
7
|
16aa14
|
DNRIYWTDISLKTASR
|
-22.078
|
-1.807
|
-0.394
|
-23.924
|
9
|
16aa15
|
DNRIYWTDISLKTIAR
|
-12.001
|
-14.568
|
21.417
|
-24.427
|
5
|
16aa16
|
DNRIYWTDISLKTISA
|
-15.186
|
-9.681
|
8.077
|
-24.060
|
8
|
Different 16 amino acid length short-peptides (Asp673 to Arg688) derived from alanine scanning was further docked and analyzed to validate its binding efficiency and binding site interactions to LRP6. We also ensured whether all these peptides, including the unmutated 16aa inhibitory peptide, were also binding to the target site of DKK1, which was highly involved in LRP6 contact. We utilized three docking platforms pyDock, Hawkdock, and Cluspro 2.0 to validate extensively and to confirm the identical output results of our docking studies. The results presented in Fig. 3 show that the best poses retrieved from each software resulted in conformations very close to the original 16 amino acid peptide, thereby validating its effectiveness as a prevailing DKK1 inhibitor. Figure 4 illustrates the possible mechanism of action of the designed and screened peptides on the pre-associated LRP6 and DKK1. Figure 4A displays the complex of 16aa8-LRP6 (Cyan & Orange) and 16aa8-DKK1 (Pink & Green). Figure 4B is the sandwich view of the same complex. This view gives a clear idea about how the peptide is blocking the association of DKK1 to LRP6 and also gives the vision to recognize the dual effect and binding efficiency of the 16aa8 peptide to inhibit both Lrp6 and DKK1 at their target site. The same tactic was further verified also to all the screened mutated peptides and unmutated inhibitory peptides and interestingly, all peptides were shown the same kind of effect. Figure 4C is the docking result of 16aa8 shows the details of interaction with the key residues of DKK1 identified earlier.
Further, the initially docked complexes of 16aa and screened mutated peptides (16aa1, 16aa6-8) with target protein LRP6 (Chain A from PDB ID: 3S8V) were subjected to MD Simulation studies for analyzing the stability in terms of RMSD and molecular interactions during the 50ns period. The stability of the protein-peptide complexes was analyzed using the RMSD of backbone atoms. The calculated average RMSD of these complexes were in between the 0.2 to 0.6 nm range (Fig. 5). At the end of the 50 ns simulation, all the protein-peptide complexes attained their stable conformations. The calculated RMSF revealed the occurrence of slight fluctuations in the loop region. The active site residues were stable at 0.3 nm deviation.
Physiochemical Properties Of The Peptide:
Several tools such as ExPASy (ProtParam) and ToxinPred were used to identify the half-life and toxicity of the peptide. Half-life is defined as calculating the time taken for half of the quantity of a particular drug or any compounds in a cell to get vanished after its process. ProtParam relies on the "N-end rule", which relates to the half-life of a protein to the individuality of its N-terminal residue. Considering the N-terminal of the sequence considered is F (Phe), the estimated half-life is 1.1 to 4.4 hours (mammalian reticulocytes, in vitro) for all peptides (Table 3). The instability index (II) is computed to be -0.832 to 0.23 this classifies the peptide 16aa7 as unstable and remaining are stable. The physical stability of peptides influenced by both intrinsic and external factors. Prediction of peptide stability suggests the precautionary steps be taken to make the therapeutically potent unstable peptides to stable one through certain biochemical analysis. ToxinPred was used to predict the toxicity of peptides with desired toxicity by mutating the minimum number of amino acids. ‘ToxinPred’ tool allows submitting query peptide in FASTA format and to optimize the peptide sequence to get maximum/minimum/desired toxicity based upon the Quantitative Matrix-based position-specific scores by comparing with the toxic/non-toxic data set. The predicted ‘ToxinPred’ toxicity results of the peptides indicated that it is non-toxic compared to the mutated peptides in silico.
Table 3. Results of physicochemical, toxicity, and ADMET properties of the 16aa inhibitory and 16aa1, 16aa7-9 peptides
Name
|
MlogP
|
S + logD
|
Molecular weight
|
Instability index*
|
Half-life*
|
Toxicity **
|
16aa
|
-10.251
|
-5.818
|
1967.268
|
0.23
|
1.1 hr
|
Non-toxin
|
16aa1
|
-9.484
|
-5.544
|
1923.258
|
0.23
|
4.4 hrs
|
Non-toxin
|
16aa7
|
-9.757
|
-5.877
|
1937.242
|
-5.08
|
1.1hr
|
Non-toxin
|
16aa8
|
-9.484
|
-5.64
|
1923.258
|
0.23
|
1.1 hr
|
Non-toxin
|
16aa9
|
10.627
|
-6.41
|
1925.187
|
0.23
|
1.1hr
|
Non-toxin
|
MlogP. Moriguchi estimation of LogP; S + logD – Octanol water partition co-efficient; *results from Protparam; **Results from Toxinpred.
Further, we predicted the fraction unbound and total clearance of the unmutated and mutated peptides screened using pkCSM. In drug development, the fraction unbound in plasma offers a better perceptive of the pharmacokinetic properties of the peptides and helps in choosing the best one. The unbound fractions are active and bind to proteins to make drug-protein complex. The resulted unbound fraction ranges from 0.371 to 0.373 for the top five compounds. The total clearance is examined as the proportionally constant (CLtot) and it is the mixture of renal and hepatic clearance. The values range from − 0.42 to -0.898 Log(ml/min/kg) indicates the adaptive nature of these compounds' physiological environment while carrying out therapeutic validations. Then the maximum tolerated dose (MRTD) was calculated which provides the toxic dose threshold in humans. The values obtained were 0.438 for all the five peptides and this indicates that the compounds are non-toxic to humans.