Free energy simulation to investigate the effect of amino acid sequence environment on the severity of osteogenesis imperfecta by glycine mutations in collagen

Kyung Hoon Lee, Mark M.Banaszak Holl

Research output: Contribution to journalArticleResearchpeer-review

2 Citations (Scopus)


Molecular dynamics simulations were carried out to calculate the free energy change difference of two collagen-like peptide models for Glya→Ser mutations causing two different osteogenesis imperfecta phenotypes. These simulations were performed to investigate the impact of local amino acid sequence environment adjacent to a mutation site on the stability of the collagen. The average free energy differences for a Glya→Ser mutant relative to a wild type are 3.4 kcal/mol and 8.2 kcal/mol for a nonlethal site and a lethal site, respectively. The free energy change differences of mutant containing two Ser residues relative to the wild type at the nonlethal and lethal mutation sites are 4.6 and 9.8 kcal/mol, respectively. Although electrostatic interactions stabilize mutants containing one or two Ser residues at both mutation sites, van der Waals interactions are of sufficient magnitude to cause a net destabilization. The presence of Gln and Arg near the mutation site, which contain large and polar side chains, provide more destabilization than amino acids containing small and nonpolar side chains.

Original languageEnglish
Pages (from-to)401-409
Number of pages9
Issue number6
Publication statusPublished - 1 Jun 2011
Externally publishedYes


  • collagen
  • free energy simulation
  • osteogenesis imperfecta
  • stochastic boundary molecular dynamics simulation
  • thermodynamic integration

Cite this