Modulation of biomineralization morphology by phosphorylated collagen peptides: insights into osteogenesis imperfecta pathophysiology

Osteogenesis imperfecta (OI) is a hereditary skeletal disorder characterized by bone fragility and deformities, primarily attributed to defects in type I collagen, the most abundant structural protein in humans. Multiple phosphorylation sites have been detected within collagen, suggesting that phosphorylation may influence mineralization processes, thereby impacting the development of OI. In this study, we investigated the modulation of biomineralization morphology by phosphorylated collagen peptides mimicking Gly-Ser mutations in osteogenesis imperfecta. A series of collagen peptide sequences, including GPO13S, GPO13pS, GPO12S, GPO12pS, GPO11S, and GPO11pS, were synthesized to explore the role of phosphorylation in peptide stability and its templating effect on biomineralization. The CD results indicated that the phosphorylation of Gly-pSer mutants reduces the stability of collagen peptides. SEM images revealed that phosphorylated peptides acted as templates, guiding the morphology of calcium carbonate into either olive-like or spherical structures, depending on their conformational state of the peptides. Non-phosphorylated peptides maintained a calcite crystal structure. The XRD patterns predominantly exhibited peaks associated with calcite and vaterite for GPO13pS-CaCO3, GPO12pS-CaCO3, and GPO11pS-CaCO3, and peaks associated with calcite for GPO13S-CaCO3, GPO12S-CaCO3, and GPO11S-CaCO3, indicating a transformation of mesocrystals influenced by peptide phosphorylation. Our findings elucidate the crucial role of phosphorylated collagen peptides in mediating biomineralization morphology and polymorph selection, offering insights into the complex pathophysiology of OI.

Effect of alanine versus serine at position 88 of human transthyretin mutants on the protein stability

Human transthyretin (TTR) is a homo-tetrameric plasma protein associated with a high percentage of β-sheet forming amyloid fibrils. It accumulates in tissues or extracellular matrices to cause amyloid diseases. Free energy simulations with thermodynamic integration based on all-atom molecular dynamics simulations have been carried out to analyze the effects of the His88 → Ala and Ser mutations on the stability of human TTR. The calculated free energy change differences (ΔΔG) caused by the His88 → Ala and His88 → Ser mutations are -1.84 ± 0.86 and 7.56 ± 0.55 kcal/mol, respectively, which are in excellent agreement with prior reported experimental values. The simulation results show that the H88A mutant is more stable than the wild type, whereas the H88S mutant is less stable than the wild type. The free energy component analysis shows that the contribution to the free energy change difference (ΔΔG) for the His88 → Ala and His88 → Ser mutations mainly arise from electrostatic and van der Waals interactions, respectively. The electrostatic term stabilizes the H88A mutant more than the wild type, but the van der Waals interaction destabilizes the H88S mutant relative to the wild type. Individual residue contributions to the free energy change show neighboring residues exert stabilizing and destabilizing influence on the mutants. The implications of the simulation results for understanding the stabilizing and destabilizing effect and its contribution to protein stability are discussed.

Free energy simulations to study mutational effect of a conserved residue, Trp24, on stability of human serum retinol-binding protein

Human serum retinol-binding protein (RBP) is a plasma transport protein for vitamin A. RBP is a prime subclass of lipocalins, which bind nonpolar ligands within a β-barrel. To understand the role of Trp 24, one of the highly conserved residues in RBP, free energy simulations have been carried out to understand the effects of the mutations from Trp at position 24 to Leu, Phe, and Tyr in the apo-RBP on its thermal stability. We examine various unfolded systems to study the dependence of the free energy differences on the denatured structure. Our calculated free energy difference values for the three mutations are in excellent agreement with the experimental values when the initial coordinates of the seven-residue peptide segments truncated from the crystal structure are used for the denatured systems. Our free energy change differences for the Trp→Leu, Trp→Phe, and Trp→Tyr mutations are 2.50 ± 0.69, 2.58 ± 0.50, and 2.49 ± 0.48 kcal/mol, respectively, when the native-like seven-residue peptides are used as models for the denatured systems. The main contributions to the free energy change differences for the Trp24→Leu and Trp24→Phe mutations are mainly from van der Waals and covalent interactions, respectively. Electrostatic, van der Waals and covalent terms equally contribute to the free energy change difference for the Trp24→Tyr mutation. The free energy simulation helps understand the detailed microscopic mechanism of the stability of the RBP mutants relative to the wild type and the role of the highly conserved residue, Trp24, of the human RBP.Communicated by Ramaswamy H. Sarma.

Simulation analysis of selective alanine mutation effect on stability of human prion protein

Prion diseases are neurodegenerative disorders caused by spongiform degeneration of the brain. Understanding the fundamental mechanism of prion protein aggregation caused by mutations is very crucial to resolve the pathology of prion diseases. To help understand the roles of individual residues on the stability of the human prion protein, the computational method of free energy simulations based on atomistic molecular dynamics trajectories is applied to Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations of the human prion protein. The simulations show that all three alanine mutations destabilize the human prion protein. The calculated free energy change differences, ΔΔG, for the Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations are in good agreement with the experimental values. The significant destabilizing effects on the mutants relative to the wild-type protein arise from van der Waals terms. Furthermore, our free energy decomposition analysis shows that the major contribution to destabilizing the V180A and V209A mutants relative to the wild-type protein is originated from van der Waals interactions from residues near the mutation sites. In contrast, the contribution to destabilizing the F175A mutant is mainly caused by van der Waals interactions from residues near and far away from the mutation site. Our results show that the free energy simulation with a thermodynamic integration approach for selected alanine scanning mutations is beneficial for understanding the detailed mechanism of human prion protein destabilization, specific residues' role, and the hydrophobic effect on protein stability.Communicated by Ramaswamy H. Sarma.

Free energy simulations to understand the effect of Met → Ala mutations at positions 205, 206 and 213 on stability of human prion protein

Prion diseases are a family of infectious amyloid diseases affecting human and animals. Prion propagation in transmissible spongiform encephalopathies is associated with the unfolding and conversion of normal cellular prion protein into its pathogenic scrapie form. Understanding the fundamentals of prion protein aggregation caused by mutations is crucial to unravel the pathology of prion diseases. To help understand the contributions of individual residues to the stability of the human prion protein, we have carried out free energy simulations based on atomistic molecular dynamics trajectories. We focus on Met → Ala mutations at positions 205, 206 and 213, which are mostly buried residues located on helix 3 of the protein. The simulations predicted that all three mutations destabilize the prion protein. Changes in unfolding free energy upon mutation, ∆∆G, are 3.10 ± 0.79, 2.00 ± 0.26 and 3.06 ± 0.66 kcal/mol for M205A, M206A and M213A, respectively, in excellent agreement with the corresponding experimental values of 3.09 ± 0.28, 1.50 ± 0.34 and 3.12 ± 0.27 kcal/mol [T. Hart et al. (2009) PNAS 106, 5651-5656]. Component analysis indicates that the major contributions to the loss of protein stability arise from van der Waals interactions for the M205A and M206A mutations, and from van der Waals and covalent energy terms for M213A. Interestingly, while free energy contributions from a majority of residues neighboring the mutation sites tend to stabilize the wild type, there are a few residues stabilizing the mutant side chains. Our results show that this approach to free energy calculation can be very useful for understanding the detailed mechanism of human prion protein stability.

A perspective on structural and computational work on collagen

Collagen is the single most abundant protein in the extracellular matrix in the animal kingdom, with remarkable structural and functional diversity and regarded one of the most useful biomaterials.