The p53‐Dependent Expression of Frataxin Controls 5‐Aminolevulinic Acid‐Induced Accumulation of Protoporphyrin IX and Photo‐Damage in Cancerous Cells

Mitochondrial frataxin is involved in various functions such as iron homeostasis, iron–sulfur cluster biogenesis, the protection from oxidative stress and apoptosis and acts as a tumor suppressor protein. We now show that the expression of frataxin is stimulated in a p53‐dependent manner and prove that frataxin is a direct p53 target gene by showing that the p53‐responsive element in the promoter of the mouse frataxin gene is bound by p53. The bacterial expression of human frataxin stimulated maturation of human ferrochelatase, which catalyzes the insertion of iron into protoporphyrin at the last step of heme biosynthesis. Overexpression of frataxin in human cancer A431 and HeLa cells lowered 5‐aminolevulinic acid(ALA)‐induced accumulation of protoporphyrin and induced resistance to ALA‐induced photo‐damage, whereas p53 silencing with siRNA in non tumor HEK293T cells down‐regulated the expression of frataxin and increased the accumulation of protoporphyrin. Thus, the decrease of the expression of frataxin unregulated by p53 in tumor cells enhances ALA‐induced photo‐damage, by down‐regulation of mitochondrial functions.     

Iron-sulfur cluster biosynthesis. Characterization of frataxin as an iron donor for assembly of [2Fe-2S] clusters in ISU-type proteins

ISU (eukaryotes) and IscU (prokaryotes) are a homologous family of proteins that appear to provide a platform for assembly of [2Fe-2S] centers prior to delivery to an apo target protein. The intermediate [2Fe-2S] ISU-bound cluster is formed by delivery of iron and sulfur to the apo ISU, with the latter delivered through an IscS-mediated reaction. The identity of the iron donor has thus far not been established. In this paper we demonstrate human frataxin to bind from six to seven iron ions. Iron binding to frataxin has been quantitated by iron-dependent fluorescence measurements [K(D)(Fe(3+)) approximately 11.7 microM; (K(D)(Fe(2+)) approximately 55.0 microM] and isothermal titration calorimetry (ITC) [K(D)(Fe(3+)) approximately 10.2 microM]. Enthalpies and entropies for ferric ion binding were determined from calorimetric measurements. Both fluorescence (K(D) 0.45 microM) and ITC measurements (K(D) 0.15 microM) demonstrate holo frataxin to form a complex with ISU with sub-micromolar binding affinities. Significantly, apo frataxin does not bind to ISU, suggesting an important role for iron in cross-linking the two proteins and/or stabilizing the structure of frataxin that is recognized by ISU. Holo frataxin is also shown to mediate the transfer of iron from holo frataxin to nucleation sites for [2Fe-2S] cluster formation on ISU. We have demonstrated elsewhere [J. Am. Chem. Soc. 2002, 124, 8774-8775] that this iron-bound form of ISU is viable for assembly of holo ISU, either by subsequent addition of sulfide or by NifS-mediated sulfur delivery. Provision of holo frataxin and inorganic sulfide is sufficient for cluster assembly in up to 70% yield. With NifS as a sulfur donor, yields in excess of 70% of holo ISU were obtained. Both UV-vis and CD spectroscopic characteristics were found to be consistent with those of previously characterized ISU proteins. The time course for cluster assembly was monitored from the 456 nm absorbance of holo ISU formed during the [2Fe-2S] cluster assembly reaction. A kinetic rate constant k(obs) approximately 0.075 min(-)(1) was determined with 100 microM ISU, 2.4 mM Na(2)S, and 40 microM holo frataxin in 50 mM Tris-HCl (pH 7.5) with 4.3 mM DTT. Similar rates were obtained for NifS-mediated sulfur delivery, consistent with iron release from frataxin as a rate-limiting step in the cluster assembly reaction.     

Structural Characterization of Missense Mutations Using High Resolution Mass Spectrometry: A Case Study of the Parkinson’s-Related Protein,DJ-1

Missense mutations that lead to the expression of mutant proteins carrying single amino acid substitutions are the cause of numerous diseases. Unlike gene lesions, insertions, deletions, nonsense mutations, or modified RNA splicing, which affect the length of a polypeptide, or determine whether a polypeptide is translated at all, missense mutations exert more subtle effects on protein structure, which are often difficult to evaluate. Here, we took advantage of the spectral resolution afforded by the EMR Orbitrap platform, to generate a mass spectrometry-based approach relying on simultaneous measurements of the wild-type protein and the missense variants. This approach not only considerably shortens the analysis time due to the concurrent acquisition but, more importantly, enables direct comparisons between the wild-type protein and the variants, allowing identification of even subtle structural changes. We demonstrate our approach using the Parkinson’s-associated protein, DJ-1. Together with the wild-type protein, we examined two missense mutants, DJ-1_A104T and DJ-1_D149A, which lead to early-onset familial Parkinson’s disease. Gas-phase, thermal, and chemical stability assays indicate clear alterations in the conformational stability of the two mutants: the structural stability of DJ-1_D149A is reduced, whereas that of DJ-1_A104T is enhanced. Overall, we anticipate that the methodology presented here will be applicable to numerous other missense mutants, promoting the structural investigations of multiple variants of the same protein.

Open image in new window

Graphical Abstract ?

     

In silico analysis of nonsynonymous single nucleotide polymorphisms of the human adiponectin receptor 2 (ADIPOR2) gene

Polymorphisms of the ADIPOR2 gene are frequently linked to a higher risk of developing diseases including obesity, type 2 diabetes and cardiovascular diseases. Though mutations of the ADIPOR2 gene are detrimental, there is a lack of comprehensive in silico analyses of the functional and structural impacts at the protein level. Considering the involvement of ADIPOR2 in glucose uptake and fatty acid oxidation, an in silico functional analysis was conducted to explore the possible association between genetic mutations and phenotypic variations. A genomic analysis of 82 nonsynonymous SNPs in ADIPOR2 was initiated using SIFT followed by the SNAP2, nsSNPAnalyzer, PolyPhen-2, SNPs&GO, FATHMM and PROVEAN servers. A total of 10 mutations (R126W, L160Q, L195P, F201S, L235R, L235P, L256R, Y328H, E334K and Q349H) were predicted to have deleterious effects on the ADIPOR2 protein and were therefore selected for further analysis. Theoretical models of the variants were generated by comparative modeling via MODELLER 9.16. A protein structural analysis of these amino acid variants was performed using SNPeffect, I-Mutant, ConSurf, Swiss-PDB Viewer and NetSurfP to explore their solvent accessibility, molecular dynamics and energy minimization calculations. In addition, FTSite was used to predict the ligand binding sites, while NetGlycate, NetPhos2.0, UbPerd and SUMOplot were used to predict post-translational modification sites. All of the variants showed increased free energy, though F201S exhibited the highest energy increase. The root mean square deviation values of the modeled mutants strongly indicated likely pathogenicity. Remarkably, three binding sites were detected on ADIPOR2, and two mutations at positions 328 and 201 were found in the first and second binding pockets, respectively. Interestingly, no mutations were found at the post-translational modification sites. These genetic variants can provide a better understanding of the wide range of disease susceptibility associated with ADIPOR2 and aid the development of new molecular diagnostic markers for these diseases. The findings may also facilitate the development of novel therapeutic elements for associated diseases.     

Competitive adsorption of bacteriophage T4 lysozyme stability variants at hydrophilic glass surfaces

The competitive adsorption behavior exhibited by the wild-type T4 lysozyme and two of its structural stability variants was studied by 125I radioisotope labeling. The mutant lysozymes were produced by substitution of the isoleucine residue at position 3 in the wild type with a tryptophan residue, resulting in a protein with lower structural stability, or with a cysteine residue, resulting in a protein with higher structural stability. Adsorption kinetics were recorded for binary protein mixtures in contact with a clean glass surface, in which one variant had been radiolabeled and the other had not. All pair permutations were tested. The kinetic data show that in instances in which exchange reactions between adsorbed protein and dissolved protein occur, they occur such that more stable variants are removed from the surface by less stable variants. The less stable proteins thus exhibited an advantage in competitive adsorption over the more stable proteins, in these tests.     

Analysis of the relationships between evolvability,thermodynamics, and the functions of intrinsically disordered proteins/regions

The evolvability of proteins is not only restricted by functional and structural importance, but also by other factors such as gene duplication, protein stability, and an organism's robustness. Recently, intrinsically disordered proteins (IDPs)/regions (IDRs) have been suggested to play a role in facilitating protein evolution. However, the mechanisms by which this occurs remain largely unknown. To address this, we have systematically analyzed the relationship between the evolvability, stability, and function of IDPs/IDRs. Evolutionary analysis shows that more recently emerged IDRs have higher evolutionary rates with more functional constraints relaxed (or experiencing more positive selection), and that this may have caused accelerated evolution in the flanking regions and in the whole protein. A systematic analysis of observed stability changes due to single amino acid mutations in IDRs and ordered regions shows that while most mutations induce a destabilizing effect in proteins, mutations in IDRs cause smaller stability changes than in ordered regions. The weaker impact of mutations in IDRs on protein stability may have advantages for protein evolvability in the gain of new functions. Interestingly, however, an analysis of functional motifs in the PROSITE and ELM databases showed that motifs in IDRs are more conserved, characterized by smaller entropy and lower evolutionary rate, than in ordered regions. This apparently opposing evolutionary effect may be partly due to the flexible nature of motifs in IDRs, which require some key amino acid residues to engage in tighter interactions with other molecules. Our study suggests that the unique conformational and thermodynamic characteristics of IDPs/IDRs play an important role in the evolvability of proteins to gain new functions.     

Formation of protein-coated iron minerals

The ability of iron to cycle between Fe(2+) and Fe(3+) forms has led to the evolution, in different forms, of several iron-containing protein cofactors that are essential for a wide variety of cellular processes, to the extent that virtually all cells require iron for survival and prosperity. The redox properties of iron, however, also mean that this metal is potentially highly toxic and this, coupled with the extreme insolubility of Fe(3+), presents the cell with the significant problem of how to maintain this essential metal in a safe and bioavailable form. This has been overcome through the evolution of proteins capable of reversibly storing iron in the form of a Fe(3+) mineral. For several decades the ferritins have been synonymous with the function of iron storage. Within this family are subfamilies of mammalian, plant and bacterial ferritins which are all composed of 24 subunits assembled to form an essentially spherical protein with a central cavity in which the mineral is laid down. In the past few years it has become clear that other proteins, belonging to the family of DNA-binding proteins from starved cells (the Dps family), which are oligomers of 12 subunits, and to the frataxin family, which may contain up to 48 subunits, are also able to lay down a Fe(3+) mineral core. Here we present an overview of the formation of protein-coated iron minerals, with particular emphasis on the structures of the protein coats and the mechanisms by which they promote core formation. We show on the one hand that significant mechanistic similarities exist between structurally dissimilar proteins, while on the other that relatively small structural differences between otherwise similar proteins result in quite dramatic mechanistic differences.     

The emulsifying properties of β-lactoglobulin genetic variants A, B and C

The emulsifying properties of three genetic variants of β-lactoglobulin (β-lac) (the A, B and C variants) are investigated as a function of protein concentration. Differences in the emulsifying properties and emulsion stability are explained in view of the known differences in physico-chemical and structural/conformational properties of the β-lac variants. β-lac A forms the finest emulsion droplets, and β-lac C the largest droplets. The order of decreasing emulsifying ability (A>B>C) can be explained in terms of differences in the molecular structure, and conformational stability of the variant proteins. The creaming stability, when compared at the same particle size, is greatest for β-lac C, with β-lac A and B having a similar and lower stability. The differences in creaming stability may arise from a higher surface coverage for the β-lac C droplets at an equivalent particle size. The storage stability is lower for β-lac A than for β-lac B and C, which both show a similar behaviour. Storage stability differences are discussed in terms of differences in molecular structure, conformational stability, interfacial viscosity and surface coverage for the three variants.     

Folding and self-association of atTic20 in lipid membranes: implications for understanding protein transport across the inner envelope membrane of chloroplasts

Background

We have previously shown that the P. gingivalis HmuY hemophore-like protein binds heme and scavenges heme from host hemoproteins to further deliver it to the cognate heme receptor HmuR. The aim of this study was to characterize structural features of HmuY variants in the presence and absence of heme with respect to roles of tryptophan residues in conformational stability.

Results

HmuY possesses tryptophan residues at positions 51 and 73, which are conserved in HmuY homologs present in a variety of bacteria, and a tryptophan residue at position 161, which has been found only in HmuY identified in P. gingivalis strains. We expressed and purified the wildtype HmuY and its protein variants with single tryptophan residues replaced by alanine or tyrosine residues. All HmuY variants were subjected to thermal denaturation and fluorescence spectroscopy analyses. Replacement of the most buried W161 only moderately affects protein stability. The most profound effect of the lack of a large hydrophobic side chain in respect to thermal stability is observed for W73. Also replacement of the W51 exposed on the surface results in the greatest loss of protein stability and even the large aromatic side chain of a tyrosine residue has little potential to substitute this tryptophan residue. Heme binding leads to different exposure of the tryptophan residue at position 51 to the surface of the protein. Differences in structural stability of HmuY variants suggest the change of the tertiary structure of the protein upon heme binding.

Conclusions

Here we demonstrate differential roles of tryptophan residues in the protein conformational stability. We also propose different conformations of apo- and holoHmuY caused by tertiary changes which allow heme binding to the protein.     

Deciphering the Role of Filamin B Calponin-Homology Domain in Causing the Larsen Syndrome,Boomerang Dysplasia,and Atelosteogenesis Type I Spectrum Disorders via a Computational Approach

Filamins (FLN) are a family of actin-binding proteins involved in regulating the cytoskeleton and signaling phenomenon by developing a network with F-actin and FLN-binding partners. The FLN family comprises three conserved isoforms in mammals: FLNA, FLNB, and FLNC. FLNB is a multidomain monomer protein with domains containing an actin-binding N-terminal domain (ABD 1–242), encompassing two calponin-homology domains (assigned CH1 and CH2). Primary variants in FLNB mostly occur in the domain (CH2) and surrounding the hinge-1 region. The four autosomal dominant disorders that are associated with FLNB variants are Larsen syndrome, atelosteogenesis type I (AOI), atelosteogenesis type III (AOIII), and boomerang dysplasia (BD). Despite the intense clustering of FLNB variants contributing to the LS-AO-BD disorders, the genotype-phenotype correlation is still enigmatic. In silico prediction tools and molecular dynamics simulation (MDS) approaches have offered the potential for variant classification and pathogenicity predictions. We retrieved 285 FLNB missense variants from the UniProt, ClinVar, and HGMD databases in the current study. Of these, five and 39 variants were located in the CH1 and CH2 domains, respectively. These variants were subjected to various pathogenicity and stability prediction tools, evolutionary and conservation analyses, and biophysical and physicochemical properties analyses. Molecular dynamics simulation (MDS) was performed on the three candidate variants in the CH2 domain (W148R, F161C, and L171R) that were predicted to be the most pathogenic. The MDS analysis results showed that these three variants are highly compact compared to the native protein, suggesting that they could affect the protein on the structural and functional levels. The computational approach demonstrates the differences between the FLNB mutants and the wild type in a structural and functional context. Our findings expand our knowledge on the genotype-phenotype correlation in FLNB-related LS-AO-BD disorders on the molecular level, which may pave the way for optimizing drug therapy by integrating precision medicine.     

Rapid Degradation of SARS-CoV-2 Spike S Protein by A Specific Serine Protease

The S protein of SARS-CoV-2 is a crucial structural and functional component for virus entry. Due to the constant mutation of the virus, there are very limited ways to prevent and control COVID-19. This experiment used a macroscopic SDS-PAGE method and proved that the S protein of wild-type SARS-CoV-2 virus, especially the S1 subunit, is very sensitive to alkaline serine protease with acidic pI (ASPNJ), NJ represents Neanthes japonica (Izuka) from which ASP is purified). ASPNJ cleaves proteins when the carbonyl group of the peptide bond is contributed by arginine or lysine. ASPNJ can degrade the S protein very quickly and effectively in vitro with relative selectivity. It can be inferred that the S, S1 and RBD of SARS-CoV-2 variants can also be easily degraded by ASPNJ. This rapid and strong degradation of the S protein by ASPNJ may become a potential new treatment strategy.     

Molecular Docking and Molecular Dynamics Study on the Effect of ERCC1 Deleterious Polymorphisms in ERCC1-XPF Heterodimer

Excision repair cross complementation group 1 (ERCC1) is an important protein in the nucleotide excision repair (NER) pathway, which is responsible for removing DNA adducts induced by platinum based compounds. The heterodimer ERCC1-XPF is one of two endonucleases required for NER. Genetic variations or polymorphisms in ERCC1 gene alter DNA repair capacity. Reduced DNA repair (NER) capacity may result in tumors and enhances cisplatin chemotherapy in cancer patients, which functions by causing DNA damage. Therefore, ERCC1 variants have the potential to be used as a strong candidate biomarker in cancer treatments. In this study we identified five variants V116M, R156Q, A199T, S267P, and R322C of ERCC1 gene as highly deleterious. Further structural and functional analysis has been conducted for ERCC1 protein in the presence of three variants V116M, R156Q, and A199T. Occurrence of theses variations adversely affected the regular interaction between ERCC1 and XPF protein. Analysis of 20 ns molecular dynamics simulation trajectories reveals that the predicted deleterious variants altered the ERCC1-XPF complex stability, flexibility, and surface area. Notably, the number of hydrogen bonds in ERCC1-XPF mutant complexes decreased in the molecular dynamic simulation periods. Overall, this study explores the link between the ERCC1 deleterious variants and cisplatin chemotherapy for various cancers with the help of molecular docking and molecular dynamic approaches.     

Redesign the α/β Fold to Enhance the Stability of Mannanase Man23 from <Emphasis Type="Italic">Bacillus subtilis</Emphasis>

In this work, we engineered the α/β fold of mannanase Man23 based on its molecular structure analysis to obtain more stable variants. By introducing 31 single-site mutations in the α/β fold and shuffling them, the incorporation of four mutations (K178R, K207R, N340R, and S354R) displayed a good balance between high activity and stability at higher temperature and broader pH. This quartet variant was characterized by an almost threefold increased activity and a sevenfold increased stability compared to native mannanase Man23. Our results suggest that such work is safe to increase our target protein stability with no loss of activity.     

Disease-Causing Mutation in Extracellular and Intracellular Domain of FGFR1 Protein: Computational Approach

In-depth computationally based structural analysis of human fibroblast growth factor type 1 (FGFR1) protein carrying disease-causing mutation was performed in this study. Gain or loss of function due to missense mutations in FGFR1 is responsible for a variety of disorders including Kallmann syndrome, Apert syndrome, Pfeiffer syndrome, Crouzon syndrome, etc. The mutant model of the human FGFR1 protein was subjected to various in silico analysis, and most deleterious SNPs were screened out. Furthermore, docking and long molecular dynamics simulations were carried out with an intention of studying the possible impact of these mutations on the protein structure and hence its function. Analysis of various structural properties—especially of those of the functionally important regions: the extracellular immunoglobulin domain and intracellular Tyrosine kinase domain—gave some insights into the possible structural characteristics of the disease mutant and the wild-type forms of the protein. In a nutshell, compared to the wild-type protein, the mutant structures V273M and S685F are associated with significant changes, and the functionally important regions seem to adopt such structures that are not conducive for the wild-type-like functionality.     

Adaptive Modelling of Mutated FMO3 Enzyme Could Unveil Unexplored Scenarios Linking Variant Haplotypes to TMAU Phenotypes

Background: Trimethylaminuria (TMAU) is a rare genetic disease characterized by the accumulation of trimethylamine (TMA) and its subsequent excretion trough main body fluids, determining the characteristic fish odour in affected patients. We realized an experimental study to investigate the role of several coding variants in the causative gene FMO3, that were only considered as polymorphic or benign, even if the available literature on them did not functionally explain their ineffectiveness on the encoded enzyme. Methods: Mutational analysis of 26 TMAU patients was realized by Sanger sequencing. Detected variants were, subsequently, deeply statistically and in silico characterized to determine their possible effects on the enzyme activity. To achieve this goal, a docking prediction for TMA/FMO3 and an unbinding pathway study were performed. Finally, a TMAO/TMA urine quantification by 1H-NMR spectroscopy was performed to support modelling results. Results: The FMO3 screening of all patients highlighted the presence of 17 variants distributed in 26 different haplotypes. Both non-sense and missense considered variants might impair the enzymatic kinetics of FMO3, probably reducing the interaction time between the protein catalytic site and TMA, or losing the wild-type binding site. Conclusions: Even if further functional assays will confirm our predictive results, considering the possible role of FMO3 variants with still uncertain effects, might be a relevant step towards the detection of novel scenarios in TMAU etiopathogenesis.     

Rational Design of Practically Important Enzymes

Majority of native enzymes are poorly applicable for practical usage: that is why different methods of enzyme modification are used to obtain the biocatalysts with appropriate characteristics. Development of genome sequencing and various modern approaches in protein engineering allow one to identify protein of interest and to improve the enzyme properties for a particular process. This review describes the results on development of novel biocatalysts based on bioinformatics and rational design. New genes encoding formate dehydrogenase (FDH) from bacterium Staphylococcus aureus, yeasts Ogataea parapolymorpha and Saccharomyces cerevisiae and moss Physcomitrella patens (SauFDH, OpaFDH, SceFDH and PpaFDH, respectively), have been cloned. New FDHs were produced in the active form and characterized. SauFDH was shown to have at least 2-fold higher catalytic constant than other known FDHs. OpaFDH has catalytic parameters as good as those for soy FDH mutant forms, and in addition, is more thermostable. Apo- and holo-forms of SauFDH have been crystallized. Mutation of two Cys residues in Pseudomonas sp.101 enzyme (PseFDH) yields enzyme preparations with improved kinetic parameters and enhanced thermal and chemical stability. New generation of PseFDH preparations with the coenzyme specificity changed from NAD⁺ to NADP⁺ have been obtained. The effect of ionic liquids on the catalytic properties and thermal stability of six wild-type recombinant FDHs, and a number of their mutants, have been studied. In case of D-amino acid oxidase (DAAO), single-point mutations have been combined to create multi-point mutants. The introduced amino acid replacements have been shown to exert an additive effect, improving both kinetic parameters and increasing thermal and chemical stability. DAAO genes from Hansenula polymorpha yeast have been cloned. α-Amino acid ester hydrolase (AEH) gene has been cloned and expressed in the active form in E. coli. Structural modeling has been performed and the effectiveness in amino beta-lactams synthesis studied. The structure of a single-chain penicillin acylase from Alcaligenes faecalis (scAfPA) has been modeled and two variants of scAfPA gene was generated by PCR. Both variants have been expressed in E. coli, isolated and characterized. Catalytic properties of scAfPA were slightly better than those of its natural heterodimer.     

Effects of a chemical chaperone on genetic mutations in ��-galactosidase A in Korean patients with Fabry disease

Fabry disease is an X-linked inborn error of glycosphingolipid catabolism that results from mutations in the gene encoding the α-galactosidase A (GLA) enzyme. We have identified 15 distinct mutations in the GLA gene in 13 unrelated patients with classic Fabry disease and 2 unrelated patients with atypical Fabry disease. Two of the identified mutations were novel (i.e., the D231G missense mutation and the L268delfsX1 deletion mutation). This study evaluated the effects of the chemical chaperones 1-deoxygalactonojirimycin (DGJ) on the function of GLA in vitro, in cells containing missense mutations in the GLA gene. Nine missense and a nonsense mutations, including one novel mutation were cloned into mammalian expression vectors. After transient expression in COS-7 cells, GLA enzyme activity and protein expression were analyzed using fluorescence spectrophotometry and Western blot analysis, respectively. DGJ enhanced GLA enzyme activity in the M42V, I91T, R112C and F113L mutants. Interestingly, the I91T and F113L mutations are associated with the atypical form of Fabry disease. However, DGJ treatment did not have any significant effect on the GLA enzyme activity and protein expression of other mutants, including C142W, D231G, D266N, and S297F. Of note, GLA enzyme activity was not detected in the novel mutant (i.e., D231G), although protein expression was similar to the wild type. In the absence of DGJ, the E66Q mutant had wild-type levels of GLA protein expression and approximately 40% GLA activity, indicating that E66Q is either a mild mutation or a functional single nucleotide polymorphism (SNP). Thus, the results of this study suggest that the chemical chaperone DGJ enhances GLA enzyme activity and protein expression in milder mutations associated with the atypical form of Fabry disease.     

Are Pro<Superscript>8</Superscript>/Pro<Superscript>18</Superscript> really critical for functional dynamic behavior of human endostatin N-terminal peptide? A comparative molecular dynamics study

Endostatin which is derived from the non-collagenous domain 1 of collagen XVIII and is a recently identified broad spectrum anti-angiogenesis agent, inhibits 65 different tumor types. The N-terminal fragment of endostatin protein (ES) has the same antitumor, antimigration and antipermeability effects as the entire protein. In the current study, we modeled two mutant variants of ES with two mutation sites (M1-ES (Pro⁸ → Ala) and M2-ES (Pro¹⁸ → Ala)) and tried to understand proline’s effect on the peptide structure/stability by introducing P⁸A/P¹⁸A mutations, and then in order to gain functional insight into mutation caused by amino acid substitution to the peptide structure/function, these effects were predicted using computational tools. From the RMSD analyses, it can be concluded that dynamic behavior of wild-type and mutant structures was not significantly different from each other and all systems reached equilibrium. The RMSF analysis also revealed that the M2-ES has smaller overall flexibility than the WT-ES and M1-ES structures. The radius of gyration analysis then confirmed the structure of M2-ES compared to wild-type and M1 variant becomes more compact during simulation of our systems. Finally, molecular dynamics simulation analysis shows that replacement of Pro residue with Ala is able to induce a distinct β-sheet in both mutant structures. Indeed, the docking analysis shows the WT-ES and M2-ES bind to the same region of α_vβ₃ integrin, suggesting similar interaction pattern with a relatively equal binding energy into this receptor. Our results speculated that the P⁸A/P¹⁸A replacements confer no improvement (or no tangible weakness) in the peptide biological activity although is able to change structural conformation of N-terminal fragment of human endostatin protein.