4.1. Physicochemical Profiling and Screening of Peptide Candidates
Physicochemical characterization using ProtParam, ProtScale, and solubility prediction tools confirmed that all 30 designed peptides exhibit favorable stability profiles (instability indices < 40), supporting their potential for biological applications (
Table 1). Systematic screening identified distinct functional classes among the candidates. The peptides displayed a broad range of positive charges (21 – 43 residues) and predominantly hydrophilic character (GRAVY values: –0.368 to –1.328), essential for electrostatic interactions with nucleic acids. All candidates exceeded the minimum solubility threshold (> 0.45), with values ranging from 0.647 to 0.795, ensuring adequate solubility for delivery applications. Predicted half-lives were uniformly favorable across biological systems, indicating robust metabolic stability.
| Peptide_# | Aliphatic Index | Instability Index | GRAVY | Half-Life (h) (Mammalian Reticulocytes, in vitro), (Yeast, in vivo), (E. Coli, in vivo) | Number of Positively Charged Residues | Number of Amino Acids | Molecular Weight (kDa) | pI | Predicted Scaled Solubility (> 0.45) | Hydrophobicity Scale |
|---|
| Peptide_1 | 51.04 | 21.57 | –0.783 | 30, > 20, > 10 | 21 | 96 | 10.33 | 10.96 | 0.718 | –0.74 |
| Peptide_2 | 55.78 | 19.04 | –0.498 | 30, > 20, > 10 | 21 | 116 | 12.25 | 10.96 | 0.694 | –0.42 |
| Peptide_3 | 66.72 | 22.54 | –0.368 | 30, > 20, > 10 | 21 | 116 | 12.39 | 10.96 | 0.647 | –0.29 |
| Peptide_4 | 51.28 | 30.08 | –1.083 | 30, > 20, > 10 | 35 | 109 | 12.13 | 12.17 | 0.733 | –1.10 |
| Peptide_5 | 42.61 | 32.11 | –1.277 | 30, > 20, > 10 | 35 | 115 | 13.04 | 11.96 | 0.721 | –1.28 |
| Peptide_6 | 42.61 | 33.18 | –1.277 | 30, > 20, > 10 | 35 | 115 | 13.04 | 11.96 | 0.721 | –1.27 |
| Peptide_7 | 45.66 | 25.17 | –0.897 | 30, > 20, > 10 | 35 | 129 | 13.92 | 12.17 | 0.731 | –0.89 |
| Peptide_8 | 53.24 | 36.6 | –0.840 | 30, > 20, > 10 | 29 | 105 | 11.24 | 12.16 | 0.738 | –0.84 |
| Peptide_9 | 59.47 | 11.99 | –0.713 | 30, > 20, > 10 | 25 | 94 | 10.14 | 11.29 | 0.736 | –0.70 |
| Peptide_10 | 53.24 | 36.6 | –0.840 | 30, > 20, > 10 | 29 | 105 | 11.24 | 12.16 | 0.694 | –0.84 |
| Peptide_11 | 51.28 | 30.08 | –1.083 | 30, > 20, > 10 | 35 | 109 | 12.13 | 12.17 | 0.695 | –1.10 |
| Peptide_12 | 45.66 | 25.17 | –0.897 | 30, > 20, > 10 | 35 | 129 | 13.92 | 12.17 | 0.738 | –0.89 |
| Peptide_13 | 55.50 | 26.48 | –0.780 | 30, > 20, > 10 | 35 | 129 | 14.06 | 12.17 | 0.738 | –0.76 |
| Peptide_14 | 55.50 | 28.32 | –0.780 | 30, > 20, > 10 | 35 | 129 | 14.06 | 12.17 | 0.738 | –0.76 |
| Peptide_15 | 50.52 | 19.46 | –0.782 | 30, > 20, > 10 | 21 | 97 | 10.43 | 10.96 | 0.721 | –0.74 |
| Peptide_16 | 49.00 | 15.41 | –0.959 | 30, > 20, > 10 | 25 | 100 | 11.05 | 11.04 | 0.731 | –0.93 |
| Peptide_17 | 61.70 | 19.48 | –0.808 | 30, > 20, > 10 | 25 | 100 | 11.18 | 11.04 | 0.730 | –0.77 |
| Peptide_18 | 42.61 | 32.11 | –1.277 | 30, > 20, > 10 | 35 | 115 | 13.04 | 11.96 | 0.730 | –1.28 |
| Peptide_19 | 45.29 | 36.36 | –1.328 | 30, > 20, > 10 | 35 | 121 | 14.08 | 11.75 | 0.695 | –1.29 |
| Peptide_20 | 50.00 | 32.18 | –1.016 | 30, > 20, > 10 | 35 | 141 | 16.01 | 11.75 | 0.695 | –0.96 |
| Peptide_21 | 47.93 | 29.29 | –0.959 | 30, > 20, > 10 | 35 | 135 | 14.96 | 11.96 | 0.795 | –0.93 |
| Peptide_22 | 42.61 | 33.18 | –1.277 | 30, > 20, > 10 | 35 | 115 | 13.04 | 11.96 | 0.795 | –1.27 |
| Peptide_23 | 44.14 | 39.46 | –1.055 | 30, > 20, > 10 | 29 | 111 | 12.14 | 11.94 | 0.795 | –1.03 |
| Peptide_24 | 43.75 | 39.7 | –1.049 | 30, > 20, > 10 | 29 | 112 | 12.20 | 11.94 | 0.730 | –1.02 |
| Peptide_25 | 43.75 | 37.41 | –1.060 | 30, > 20, > 10 | 29 | 112 | 12.29 | 11.94 | 0.795 | –1.03 |
| Peptide_26 | 43.75 | 38.67 | –1.060 | 30, > 20, > 10 | 29 | 112 | 12.29 | 11.94 | 0.795 | –1.03 |
| Peptide_27 | 39.85 | 16.9 | –1.098 | 30, > 20, > 10 | 37 | 133 | 14.25 | 11.21 | 0.795 | –1.08 |
| Peptide_28 | 39.85 | 17.96 | –1.098 | 30, > 20, > 10 | 37 | 133 | 14.25 | 11.21 | 0.795 | –1.08 |
| Peptide_29 | 38.97 | 15.19 | –1.282 | 30, > 20, > 10 | 43 | 136 | 15.00 | 11.28 | 0.721 | –1.27 |
| Peptide_30 | 38.97 | 17.65 | –1.282 | 30, > 20, > 10 | 43 | 136 | 15.00 | 11.28 | 0.721 | –1.27 |
Through multi-parameter optimization, seven peptides (Peptides 9, 15, 16, 21, 24, 27, and 28) were selected as lead candidates:
Peptide_9 exhibited the highest stability (Instability Index: 11.99), optimal hydrophilicity (GRAVY: -0.736), moderate positive charge (25 residues), and a compact structure (94 aa, 10.14 kDa), suggesting efficient membrane translocation and structural integrity.
Peptide_15 combined high stability (Instability Index: 19.46), moderate size (97 aa, 10.43 kDa), moderate positive charge (21 residues), and favorable physicochemical properties (GRAVY: –0.782), advantageous for crossing biological barriers.
Peptide_16 showed excellent stability (Instability Index: 15.41), balanced hydrophobicity (GRAVY: –0.959), moderate positive charge (25 residues), and moderate size (100 aa, 11.05 kDa), favoring cellular uptake without aggregation.
Peptide_21 (135 aa, 14.9 kDa) demonstrated the highest solubility (0.795), substantial positive charge (35 residues), and acceptable stability (Instability Index: 29.29), supporting strong DNA-binding in aqueous formulations.
Peptide_24, despite a higher instability index (39.70), was included due to its exceptional solubility (0.730), moderate charge (29 residues), balanced hydrophobicity (GRAVY: –1.049), and optimal size (112 aa, 12.20 kDa), making it suitable for formulations requiring high solubility.
Peptide_27 (133 aa, 14.3 kDa) was chosen for enhanced DNA-binding capacity, with 37 positively charged residues, good stability (Instability Index: 16.90), and balanced hydrophilicity (GRAVY: –1.098).
Peptide_28 (133 aa, 14.3 kDa) mirrored Peptide_27’s charge distribution (37 positive residues) and stability (Instability Index: 17.96), serving as a structural comparator.
Selection criteria emphasized balanced functional attributes: Positive charge density (25 – 37 residues), hydrophobicity within the CPP-preferred range (GRAVY: –0.7 to –1.1), molecular weight (10 – 15 kDa), structural stability (instability index < 30, with exceptions for compensatory traits), and high solubility (> 0.69). Peptides with extreme hydrophilicity (GRAVY < –1.25) were deprioritized due to potential membrane penetration limitations. The aliphatic indices (38.97 – 66.72) and pI values (10.96 – 12.17) of selected candidates further supported their thermostability and cationic nature, essential for DNA binding and cellular interactions. This multi-parameter strategy ensured the selection of candidates with the highest potential for success in downstream functional assays.
4.2. Structural and Biophysical Profiling of Lead Candidates
Hydrophobicity analysis using Kyte-Doolittle profiling revealed distinct amphipathic patterns across the seven selected peptides, all exhibiting the alternating hydrophobic-hydrophilic domains characteristic of effective CPPs (
Table 1 and Supplementary Table S2). Three candidates emerged with particularly favorable profiles:
Peptide_9 demonstrated optimal structural characteristics with harmonized amphipathicity, featuring moderate hydrophobic peaks (+2.0 to +2.5) in terminal regions balanced by extensive hydrophilic segments. This pattern, combined with its exceptional stability (Instability Index: 11.99) and compact architecture (94 aa), suggests ideal properties for membrane translocation while maintaining structural integrity.
Peptide_27 displayed stronger hydrophobic character (+2.5 to +3.0 peaks) coupled with high positive charge density (37 aa), indicating enhanced DNA-binding capacity, though with potential aggregation concerns that would require formulation optimization.
Peptide_24 showed exceptionally balanced hydropathy with moderate hydrophobic domains (+1.5 to +2.0) distributed throughout its sequence, complemented by outstanding solubility (0.730) and optimal molecular size (112 aa). This profile suggests low aggregation propensity while retaining sufficient hydrophobicity for membrane interaction, making it particularly suitable for stable formulation development.
Comparative analysis identified specific limitations among the remaining candidates. Peptides_15 and _16 were characterized by predominant hydrophilicity, which may reduce their membrane penetration capability. Peptide_21 was found to exhibit a complex hydropathy profile with partially exposed hydrophilic surfaces. Despite its moderate-to-large size (135 aa, ~14.9 kDa), previous studies have demonstrated that CPPs of comparable dimensions can effectively traverse cellular membranes when amphipathicity and charge are optimally distributed (
23,
24). Accordingly, peptide size was considered a minor, non-limiting factor for cellular uptake. Structural modeling using I-TASSER confirmed that Peptide_9 adopts a dual α-helical architecture conducive to both membrane interaction and DNA binding. In contrast, Peptide_24 was shown to exhibit balanced amphipathicity and favorable solvent accessibility, albeit with slightly reduced helicity.
Based on integrated analysis of hydrophobicity, structural stability, charge distribution, and molecular architecture, three candidates were prioritized for functional characterization:
Peptide_9 as the lead candidate for optimal balanced properties.
Peptide_27 for high-capacity DNA binding applications.
Peptide_24 for formulations requiring exceptional stability and solubility.
Further energy and charge profile analyses of Peptide_9 and Peptide_24 were performed using the Protein-Sol server. It was revealed that Peptide_24 exhibits a higher potential for effective endosomal escape and enhanced gene delivery (Supplementary Figure S1). Endosomal escape peptides are typically characterized by pH-sensitive behavior, whereby structural stability and limited membrane interaction are maintained under neutral pH conditions, such as those found in extracellular environments. In contrast, significant conformational changes are induced under acidic pH conditions representative of the endosomal compartment (approximately pH 5.0 - 6.5). Energy profile analysis, expressed as joules per amino acid, demonstrated that within the acidic pH range (pH 5.0 - 6.0) and across various ionic strengths, Peptide_24 undergoes a markedly greater increase in energy compared to Peptide_9. This increase was indicative of enhanced structural instability, a stronger propensity for conformational rearrangement, and increased membrane interaction. Such properties were shown to facilitate the transition of the peptide into an active amphipathic helical conformation, a critical feature for endosomal membrane disruption and efficient cargo release.
Additionally, charge profile analysis further supports Peptide_24's superior performance. At acidic pH, Peptide_24 maintains a slightly higher net positive charge compared to Peptide_9. This higher charge density is essential for forming stable complexes with negatively charged genetic cargo as well as for stronger electrostatic interactions with endosomal lipid membranes. Furthermore, this higher charge enhances the peptide’s ability to engage in mechanisms such as the "proton sponge" effect, which aids in endosomal escape (
25,
26). The results from both energy and charge profiles strongly position Peptide_24 as the candidate for efficient endosomal escape and successful gene delivery.
4.3. Evaluation of mRNA Expression for Peptide_24
The mRNA for Peptide_24 [318 nt, 106 (+ 6 His-tag) aa] was analyzed for its thermodynamic and translational characteristics (
Figures 1 and
2). Its composition is 102 A, 59 C, 91 G, and 68 U, yielding a 46.5% G + C content, providing a good balance for stability and efficient expression. The secondary structure shows a stable monomeric fold (MFE = –75.9 kcal/mol) to protect the mRNA from degradation and extend its lifespan. Analysis of dimer formation (ΔGmonomer ≈ –76 vs. ΔGdimer ≈ –30 kcal/mol) confirmed that the single-stranded form is thermodynamically favored, reducing the risk of molecule aggregation. A localized constraint was identified near the start codon: The AUG codon is partially base-paired within a stem-loop, which slightly hinders ribosome access. However, analysis of structural ensembles showed the molecule is flexible (MFE frequency = 0.01%, diversity = 80.04), meaning it constantly shifts between different shapes. This dynamic behavior allows the AUG codon to be temporarily exposed for translation initiation as seen in the expression profile (
Figure 4A).
Structural and functional dynamics of the Peptide-24 mRNA; the minimum free energy illustrates a stable monomeric fold (MFE = -75.9 kcal/mol), contributing to mRNA longevity. A key structural constraint is observed near the start codon (AUG) (highlighted in red), where partial base-pairing within a stem-loop moderately limits ribosome accessibility.
Structural and functional analysis of the Peptide-24 mRNA reveals secondary-structure features that contribute to its stability. The centroid secondary structure shows that a notable structural constraint forms near the start codon (AUG) (highlighted in red), where partial base-pairing within a stem–loop moderately reduces ribosomal accessibility and may influence translation initiation.
HEX docking predicted poses 1 and 2 of Peptide_24 with a double-stranded DNA fragment (extracted from PDB: 1AIO) and detailed view of the Peptide_24 sequence and motifs.
This mild, inherent limitation on translation initiation aligns with biosafety objectives. While Peptide_24's cationic density is crucial for interacting with its cargo, it also carries a risk of cytotoxicity through non-specific binding to cell membranes. Therefore, a moderated initiation rate helps prevent excessive protein production and subsequent cellular stress, acting as a built-in safety mechanism. This controlled expression level is intentionally maintained in our design to provide a safety margin during initial validation tests.
4.4. Docking Analysis of Peptide-DNA Interactions
To assess the DNA-binding potential of the designed peptides, molecular docking was performed using HEX. HEX is a fast rigid-body docking program that employs spherical polar Fourier correlations to evaluate both shape complementarity and electrostatic interactions between macromolecules. This algorithm allows efficient sampling of rotational and translational degrees of freedom, making it suitable for modeling peptide–nucleic acid interactions. Prior to docking, a representative double-strand DNA fragment was extracted from PDB entry 1AIO, and its structure was energy-minimized to remove steric clashes and optimize backbone geometry.
Docking of Peptide_24 against the optimized DNA fragment revealed two energetically favorable binding modes, with total docking energies (E_total) of −574.3 (Pose 1) and −545.0 (Pose 2), respectively (
Figure 3 and Supplementary Figure S2). These two docking poses correspond to distinct DNA-interaction regions, both localized in segments of the peptide enriched in positively charged residues (Lys and Arg). Structural inspection of the complexes indicated that these cationic residues are oriented toward the negatively charged phosphate backbone of DNA, resulting in strong hydrogen and electrostatic attractions and good spatial complementarity.
Biochemical characterization of purified peptide and its DNA-binding activity; A, sodium dodecyl sulfate gel electrophoresis (SDS-PAGE) analysis of purified peptide; samples were heat-denatured and separated on a 15% polyacrylamide gel. Lanes (left to right): Elution 1-3, protein; B, gel retardation assay evaluating peptide-DNA complex formation with increasing peptide concentrations [Lane 1: Plasmid (P; 70 ng); MWM: 1Kb DNA molecular weight marker; 2: P+ 5 μg peptide; 3: P+ 10 μg peptide; 4: P+ 15 μg peptide; 5: P+ 20 μg peptide; 6: P+ 30 μg peptide; 7: P+ 40 μg peptide; 8: P+ 50 μg peptide]; progressive retardation of plasmid migration indicates effective peptide-DNA complex formation, confirming the DNA-binding capacity of the purified peptide.
The pose with the lower energy score (−574.3) exhibited tighter packing and more extensive electrostatic interfaces, suggesting it represents the primary binding mode, whereas the second pose (−545.0) may reflect an alternative register of interaction. In both cases, the interaction pattern is consistent with non-sequence-specific DNA binding mediated by hydrogen and electrostatic interactions. This mechanism is commonly observed in cell-penetrating and DNA-condensing peptides. Thus, docking results support the ability of Peptide_24 to form stable peptide–DNA complexes through multiple positively charged regions, which is essential for efficient gene delivery.