Sequence environment of mutation affects stability and folding in collagen model peptides of osteogenesis imperfecta

Designing collagens to shed light on the multi-scale structure-function mapping of matrix disorders

Collagens are the most abundant structural proteins in the extracellular matrix of animals and play crucial roles in maintaining the structural integrity and mechanical properties of tissues and organs while mediating important biological processes. Fibrillar collagens have a unique triple helix structure with a characteristic repeating sequence of (Gly-X-Y)n. Variations within the repetitive sequence can cause misfolding of the triple helix, resulting in heritable connective tissue disorders. The most common variations are single-point missense mutations that lead to the substitution of a glycine residue with a bulkier amino acid (Gly → X). In this review, we will first discuss the importance of collagen's triple helix structure and how single Gly substitutions can impact its folding, structure, secretion, assembly into higher-order structures, and biological functions. We will review the role of "designer collagens," i.e., synthetic collagen-mimetic peptides and recombinant bacterial collagen as model systems to include Gly → X substitutions observed in collagen disorders and investigate their impact on structure and function utilizing in vitro studies. Lastly, we will explore how computational modeling of collagen peptides, especially molecular and steered molecular dynamics, has been instrumental in probing the effects of Gly substitutions on structure, receptor binding, and mechanical stability across multiple length scales.

Association between N-Terminal Pyrenes Stabilizes the Collagen Triple Helix

Collagen model peptides featuring the fluorophore pyrene at their N-termini have been synthesized, and their thermal denaturation has been examined using circular dichroism (CD) and fluorescence spectroscopies. Flanking the (Pro-Hyp-Gly)7 core of the peptide monomers at positions 1 and/or 23 in the primary sequence, Lys residues were introduced to ensure water solubility. Triple helices derived from such peptides show a broad excimer emission at ∼480 nm, indicative of interaction between the pyrene units. CD experiments show that the fluorophores enhance helix stability primarily through entropic effects. Unfolding temperatures (Tm) increase by up to 7 °C for systems with N-terminal lysine residues and by up to 21 °C for systems in which the first-position Lys is replaced by Ala. Tm values derived from fluorescence measurements (at 50 μM) typically lie within ∼1 °C of those obtained using CD (at 200 μM). Computational modeling in a water continuum using B3LYP-GD3 and M06-2X functionals predicts that face-to-face association of fluorophores can occur while H-bonding within the [(POG)n]3 assembly is retained. Such parallel stacking is consistent with hydrophobically driven stabilization. Labeling collagen peptides with pyrene is a synthetically simple way to promote triple helicity while providing a means to obtain Tm data on relatively dilute samples.

A solution structure analysis reveals a bent collagen triple helix in the complement activation recognition molecule mannan-binding lectin

Collagen triple helices are critical in the function of mannan-binding lectin (MBL), an oligomeric recognition molecule in complement activation. The MBL collagen regions form complexes with the serine proteases MASP-1 and MASP-2 in order to activate complement, and mutations lead to common immunodeficiencies. To evaluate their structure-function properties, we studied the solution structures of four MBL-like collagen peptides. The thermal stability of the MBL collagen region was much reduced by the presence of a GQG interruption in the typical (X-Y-Gly)n repeat compared to controls. Experimental solution structural data were collected using analytical ultracentrifugation and small angle X-ray and neutron scattering. As controls, we included two standard Pro-Hyp-Gly collagen peptides (POG)_10-13, as well as three more peptides with diverse (X-Y-Gly)n sequences that represented other collagen features. These data were quantitatively compared with atomistic linear collagen models derived from crystal structures and 12,000 conformations obtained from molecular dynamics simulations. All four MBL peptides were bent to varying degrees up to 85^o in the best-fit molecular dynamics models. The best-fit benchmark peptides (POG)_n were more linear but exhibited a degree of conformational flexibility. The remaining three peptides showed mostly linear solution structures. In conclusion, the collagen helix is not strictly linear, the degree of flexibility in the triple helix depends on its sequence, and the triple helix with the GQG interruption showed a pronounced bend. The bend in MBL GQG peptides resembles the bend in the collagen of complement C1q and may be key for lectin pathway activation.

Recent Advances in Collagen Mimetic Peptide Structure and Design

Collagen mimetic peptides (CMPs) fold into a polyproline type II triple helix, allowing the study of the structure and function (or misfunction) of the collagen family of proteins. This Perspective will focus on recent developments in the use of CMPs toward understanding the structure and controlling the stability of the triple helix. Triple helix assembly is influenced by various factors, including the single amino acid propensity for the triple helix fold, pairwise interactions between these amino acids, and long-range effects observed across the helix, such as bend, twist, and fraying. Important progress in creating a comprehensive and predictive understanding of these factors for peptides with exclusively natural amino acids has been made. In contrast, several groups have successfully developed unnatural amino acids that are engineered to stabilize the triple helical structure. A third approach to controlling the triple helical structure includes covalent cross-linking of the triple helix to stabilize the assembly, which eliminates the problematic equilibrium of unfolding into monomers and enforces compositional control. Advances in all these areas have resulted in significant improvements to our understanding and control of this important class of protein, allowing for the design and application of more chemically complex and well-controlled collagen mimetic biomaterials.

Osteogenesis Imperfecta: Search for Mutations in Patients from the Republic of Bashkortostan (Russia)

Osteogenesis imperfecta (OI) is an inherited disease of bone characterized by increased bone fragility. Here, we report the results of the molecular architecture of osteogenesis imperfecta research in patients from Bashkortostan Republic, Russia. In total, 16 mutations in COL1A1, 11 mutations in COL1A2, and 1 mutation in P3H1 and IFIMT5 genes were found in isolated states; 11 of them were not previously reported in literature. We found mutations in CLCN7, ALOX12B, PLEKHM1, ERCC4, ARSB, PTH1R, and TGFB1 that were not associated with OI pathogenesis in patients with increased bone fragility. Additionally, we found combined mutations (c.2869C>T, p. Gln957* in COL1A1 and c.1197+5G>A in COL1A2; c.579delT, p. Gly194fs in COL1A1 and c.1197+5G>A in COL1A2; c.2971G>C, p. Gly991Arg in COL1A2 and c.212G>C, p.Ser71Thr in FGF23; c.-14C>T in IFITM5 and c.1903C>T, p. Arg635* in LAMB3) in 4 patients with typical OI clinic phenotypes.

Terminal repeats impact collagen triple-helix stability through hydrogen bonding

Nearly 30% of human proteins have tandem repeating sequences. Structural understanding of the terminal repeats is well-established for many repeat proteins with the common α-helix and β-sheet foldings. By contrast, the sequence–structure interplay of the terminal repeats of the collagen triple-helix remains to be fully explored. As the most abundant human repeat protein and the most prevalent structural component of the extracellular matrix, collagen features a hallmark triple-helix formed by three supercoiled polypeptide chains of long repeating sequences of the Gly–X–Y triplets. Here, with CD characterization of 28 collagen-mimetic peptides (CMPs) featuring various terminal motifs, as well as DSC measurements, crystal structure analysis, and computational simulations, we show that CMPs only differing in terminal repeat may have distinct end structures and stabilities. We reveal that the cross-chain hydrogen bonding mediated by the terminal repeat is key to maintaining the triple-helix's end structure, and that disruption of it with a single amide to carboxylate substitution can lead to destabilization as drastic as 19 °C. We further demonstrate that the terminal repeat also impacts how strong the CMP strands form hybrid triple-helices with unfolded natural collagen chains in tissue. Our findings provide a spatial profile of hydrogen bonding within the CMP triple-helix, marking a critical guideline for future crystallographic or NMR studies of collagen, and algorithms for predicting triple-helix stability, as well as peptide-based collagen assemblies and materials. This study will also inspire new understanding of the sequence–structure relationship of many other complex structural proteins with repeating sequences.

Collagen mimetic peptides (CMPs) only differing in terminal repeat have distinct stabilities and end structures due to a spatial hydrogen bonding profile that is useful for future crystallography, algorithm prediction, and materials of collagen.

Three Patient Kindred with a Novel Phenotype of Osteogenesis Imperfecta due to a COL1A1 Variant

Osteogenesis imperfecta (OI) is characterized by fractures and progressive bone deformities. Fracture rates peak during the toddler and adolescent years and decline during adulthood but do not stop entirely. We describe a kindred, the affected members of which were the mother and two sons, who presented with an apparently unique phenotype of OI. Our patients demonstrated a pattern of prenatal bone deformities followed by multiple, nontraumatic long bone fractures within the first two years of life and then an absence of nontraumatic fractures thereafter. No extra-skeletal manifestations have been noted to date. The mother did not receive bisphosphonate therapy but had no nontraumatic fractures after the age of five months. Intravenous bisphosphonate therapy was started for both sons within two months of birth, with the most recent infusions at age 18 months and 28 months in Patients 2 and 3, respectively. Two patients harbored a variant of uncertain significance in the COL1A1 gene. This heterozygous variant, c.3548C>T; p.(Pro1183Leu), is listed in the OI Variant Database as affecting only one other individual with osteopenia. We describe three family members with a unique presenting phenotype of OI, characterized by cessation of nontraumatic fractures after the first two years of life.

Histopathology of osteogenesis imperfecta bone. Supramolecular assessment of cells and matrices in the context of woven and lamellar bone formation using light, polarization and ultrastructural microscopy

Diaphyseal long bone cortical tissue from 30 patients with lethal perinatal Sillence II and progressively deforming Sillence III osteogenesis imperfecta (OI) has been studied at multiple levels of structural resolution. Interpretation in the context of woven to lamellar bone formation by mesenchymal osteoblasts (MOBLs) and surface osteoblasts (SOBLs) respectively demonstrates lamellar on woven bone synthesis as an obligate self-assembly mechanism and bone synthesis following the normal developmental pattern but showing variable delay in maturation caused by structurally abnormal or insufficient amounts of collagen matrix. The more severe the variant of OI is, the greater the persistence of woven bone and the more immature the structural pattern; the pattern shifts to a structurally stronger lamellar arrangement once a threshold accumulation for an adequate scaffold of woven bone has been reached. Woven bone alone characterizes lethal perinatal variants; variable amounts of woven and lamellar bone occur in progressively deforming variants; and lamellar bone increasingly forms rudimentary and then partially compacted osteons not reaching full compaction. At differing levels of microscopic resolution: lamellar bone is characterized by short, obliquely oriented lamellae with a mosaic appearance in progressively deforming forms; polarization defines tissue conformations and localizes initiation of lamellar formation; ultrastructure of bone forming cells shows markedly dilated rough endoplasmic reticulum (RER) and prominent Golgi bodies with disorganized cisternae and swollen dispersed tubules and vesicles, structural indications of storage disorder/stress responses and mitochondrial swelling in cells with massively dilated RER indicating apoptosis; ultrastructural matrix assessments in woven bone show randomly oriented individual fibrils but also short pericellular bundles of parallel oriented fibrils positioned obliquely and oriented randomly to one another and in lamellar bone show unidirectional fibrils that deviate at slight angles to adjacent bundles and obliquely oriented fibril groups consistent with twisted plywood fibril organization. Histomorphometric indices, designed specifically to document woven and lamellar conformations in normal and OI bone, establish ratios for: i) cell area/total area X 100 indicating the percentage of an area occupied by cells (cellularity index) and ii) total area/number of cells (pericellular matrix domains). Woven bone is more cellular than lamellar bone and OI bone is more cellular than normal bone, but these findings occur in a highly specific fashion with values (high to low) encompassing OI woven, normal woven, OI lamellar and normal lamellar conformations. Conversely, for the total area/number of cells ratio, pericellular matrix accumulations in OI woven are smallest and normal lamellar largest. Since genotype-phenotype correlation is not definitive, interposing histologic/structural analysis allowing for a genotype-histopathologic-phenotype correlation will greatly enhance understanding and clinical management of OI.

Exploring the potential of the recombinant human collagens for biomedical and clinical applications: a short review

Natural animal collagen and its recombinant collagen are favourable replacements in human tissue engineering due to their remarkable biomedical property. However, this exploitation is largely restricted due to the potential of immunogenicity and virus contamination. Exploring new ways to produce human collagen is fundamental to its biomedical and clinical application. All human fibrillar collagen molecules have three polypeptide chains constructed from a repeating Gly-Xaa-Yaa triplet, where Xaa and Yaa represent one random amino acid. Using cDNA techniques to modify several repeat sequences of the cDNA fragment, a novel human collagen, named recombinant human-like collagen (rHLC), with low immunogenicity and little risk from hidden virus can be engineered and notably tailored to specific applications. Human-like collagen (HLC) was initially used as a coating to modify the tissue engineering scaffold, and then used as the scaffold after cross-link agents were added to increase its mechanical strength. Due to its good biocompatibility, low immunogenicity, stabilised property, and the ability of mass production, HLC has been widely used in skin injury treatments, vascular scaffolds engineering, cartilage, bone defect repair, skincare, haemostatic sponge, and drug delivery, including coating with medical nanoparticles. In this review, we symmetrically reviewed the development, recent advances in design and application of HLC, and other recombinant human collagen-based biomedicine potentials. At the end, future improvements are also discussed.

Triple-Helix-Stabilizing Effects in Collagen Model Peptides Containing PPII-Helix-Preorganized Diproline Modules

Collagen model peptides (CMPs) serve as tools for understanding stability and function of the collagen triple helix and have a potential for biomedical applications. In the past, interstrand cross-linking or conformational preconditioning of proline units through stereoelectronic effects have been utilized in the design of stabilized CMPs. To further study the effects determining collagen triple helix stability we investigated a series of CMPs containing synthetic diproline-mimicking modules (ProMs), which were preorganized in a PPII-helix-type conformation by a functionalizable intrastrand C2 bridge. Results of CD-based denaturation studies were correlated with calculated (DFT) conformational preferences of the ProM units, revealing that the relative helix stability is mainly governed by an interplay of main-chain preorganization, ring-flip preference, adaptability, and steric effects. Triple helix integrity was proven by crystal structure analysis and binding to HSP47.

Corrigendum: Genotype-Phenotype Association Analysis Reveals New Pathogenic Factors for Osteogenesis Imperfecta Disease

[This corrects the article DOI: 10.3389/fphar.2019.01200.].

Adverse effects of Alport syndrome-related Gly missense mutations on collagen type IV: Insights from molecular simulations and experiments

Patients with Alport syndrome (AS) exhibit blood and elevated protein levels in their urine, inflamed kidneys, and many other abnormalities. AS is attributed to mutations in type IV collagen genes, particularly glycine missense mutations in the collagenous domain of COL4A5 that disrupt common structural motifs in collagen from the repeat (Gly-Xaa-Yaa)_n amino acid sequence. To characterize and elucidate the molecular mechanisms underlying how AS-related mutations perturb the structure and function of type IV collagen, experimental studies and molecular simulations were integrated to investigate the structure, stability, protease sensitivity, and integrin binding affinity of collagen-like proteins containing amino acid sequences from the α5(IV) chain and AS-related Gly missense mutations. We show adverse effects where (i) three AS-related Gly missense mutations significantly reduced the structural stability of the collagen in terms of decreased melting temperatures and calorimetric enthalpies, in conjunction with a collective drop in the external work needed to unfold the peptides containing mutation sequences; (ii) due to local unwinding around the sites of mutations, these triple helical peptides were also degraded more rapidly by trypsin and chymotrypsin, as these enzymes could access the collagenous triple helix more easily and increase the number of contacts; (iii) the mutations further abolished the ability of the recombinant collagens to bind to integrins and greatly reduced the binding affinities between collagen and integrins, thus preventing cells from adhering to these mutants. Our unified experimental and computational approach provided underlying insights needed to guide potential therapies for AS that ameliorate the adverse effects from AS disease onset and progression.

Genotype-Phenotype Association Analysis Reveals New Pathogenic Factors for Osteogenesis Imperfecta Disease

Osteogenesis imperfecta (OI), mainly caused by structural abnormalities of type I collagen, is a hereditary rare disease characterized by increased bone fragility and reduced bone mass. Clinical manifestations of OI mostly include multiple repeated bone fractures, thin skin, blue sclera, hearing loss, cardiovascular and pulmonary system abnormalities, triangular face, dentinogenesis imperfecta (DI), and walking with assistance. Currently, 20 causative genes with 18 subtypes have been identified for OI, of them, variations in COL1A1 and COL1A2 have been demonstrated to be major causative factors to OI. However, the complexity of the bone formation process indicates that there are potential new pathogenic genes associated with OI. To comprehensively explore the underlying mechanism of OI, we conducted association analysis between genotypes and phenotypes of OI diseases and found that mutations in COL1A1 and COL1A2 contributed to a large proportion of the disease phenotypes. We categorized the clinical phenotypes and the genotypes based on the variation types for those 155 OI patients collected from literature, and association study revealed that three phenotypes (bone deformity, DI, walking with assistance) were enriched in two variation types (the Gly-substitution missense and groups of frameshift, nonsense, and splicing variations). We also identified four novel variations (c.G3290A (p.G1097D), c.G3289C (p.G1097R), c.G3289A (p.G1097S), c.G3281A (p.G1094D)) in gene COL1A1 and two novel variations (c.G2332T (p.G778C), c.G2341T (p.G781C)) in gene COL1A2, which could potentially contribute to the disease. In addition, we identified several new potential pathogenic genes (ADAMTS2, COL5A2, COL8A1) based on the integration of protein–protein interaction and pathway enrichment analysis. Our study provides new insights into the association between genotypes and phenotypes of OI and novel information for dissecting the underlying mechanism of the disease.

Molecular underpinnings of integrin binding to collagen-mimetic peptides containing vascular Ehlers-Danlos syndrome-associated substitutions

Collagens carry out critical extracellular matrix (ECM) functions by interacting with numerous cell receptors and ECM components. Single glycine substitutions in collagen III, which predominates in vascular walls, result in vascular Ehlers-Danlos syndrome (vEDS), leading to arterial, uterine, and intestinal rupture and an average life expectancy of <50 years. Collagen interactions with integrin α₂β₁ are vital for platelet adhesion and activation; however, how these interactions are impacted by vEDS-associated mutations and by specific amino acid substitutions is unclear. Here, we designed collagen-mimetic peptides (CMPs) with previously reported Gly → Xaa (Xaa = Ala, Arg, or Val) vEDS substitutions within a high-affinity integrin α₂β₁-binding motif, GROGER. We used these peptides to investigate, at atomic-level resolution, how these amino acid substitutions affect the collagen III-integrin α₂β₁ interaction. Using a multitiered approach combining biological adhesion assays, CD, NMR, and molecular dynamics (MD) simulations, we found that these substitutions differentially impede human mesenchymal stem cell spreading and integrin α₂-inserted (α₂I) domain binding to the CMPs and were associated with triple-helix destabilization. Although an Ala substitution locally destabilized hydrogen bonding and enhanced mobility, it did not significantly reduce the CMP-integrin interactions. MD simulations suggested that bulkier Gly → Xaa substitutions differentially disrupt the CMP-α₂I interaction. The Gly → Arg substitution destabilized CMP-α₂I side-chain interactions, and the Gly → Val change broke the essential Mg²⁺ coordination. The relationship between the loss of functional binding and the type of vEDS substitution provides a foundation for developing potential therapies for managing collagen disorders.

Effects of flexibility of the α2 chain of type I collagen on collagenase cleavage

Cleavage of collagen by collagenases such as matrix metalloproteinase 1 (MMP-1) is a key step in development, tissue remodeling, and tumor proliferation. The abundant heterotrimeric type I collagen composed of two α1(I) chains and one α2(I) chain is efficiently cleaved by MMP-1 at a unique site in the triple helix, a process which may be initiated by local unfolding within the peptide chains. Atypical homotrimers of the α1(I) chain, found in embryonic and cancer tissues, are very resistant to MMP cleavage. To investigate MMP-1 cleavage, recombinant homotrimers were constructed with sequences from the MMP cleavage regions of human collagen chains inserted into a host bacterial collagen protein system. All triple-helical constructs were cleaved by MMP-1, with α2(I) homotrimers cleaved efficiently at a rate similar to that seen for α1(II) and α1(III) homotrimers, while α1(I) homotrimers were cleaved at a much slower rate. The introduction of destabilizing Gly to Ser mutations within the human collagenase susceptible region of the α2(I) chain did not interfere with MMP-1 cleavage. Molecular dynamics simulations indicated a greater degree of transient hydrogen bond breaking in α2(I) homotrimers compared with α1(I) homotrimers at the MMP-1 cleavage site, and showed an extensive disruption of hydrogen bonding in the presence of a Gly to Ser mutation, consistent with chymotrypsin digestion results. This study indicates that α2(I) homotrimers are susceptible to MMP-1, proves that the presence of an α1(I) chain is not a requirement for α2(I) cleavage, and supports the importance of local unfolding of α2(I) in collagenase cleavage.

Collagen Gly missense mutations: Effect of residue identity on collagen structure and integrin binding

Gly missense mutations in type I collagen, which replace a conserved Gly in the repeating (Gly-Xaa-Yaa)_n sequence with a larger residue, are known to cause Osteogenesis Imperfecta (OI). The clinical consequences of such mutations range from mild to lethal, with more serious clinical severity associated with larger Gly replacement residues. Here, we investigate the influence of the identity of the residue replacing Gly within and adjacent to the integrin binding ⁵⁰²GFPGER⁵⁰⁷ sequence on triple-helix structure, stability and integrin binding using a recombinant bacterial collagen system. Recombinant collagens were constructed with Gly substituted by Ala, Ser or Val at four positions within the integrin binding region. All constructs formed a stable triple-helix structure with a small decrease in melting temperature. Trypsin was used to probe local disruption of the triple helix, and Gly to Val replacements made the triple helix trypsin sensitive at three of the four sites. Any mutation at Gly505, eliminated integrin binding, while decreased integrin binding affinity was observed in the replacement of Gly residues at Gly502 following the order Val > Ser > Ala. Molecular dynamics simulations indicated that all Gly replacements led to transient disruption of triple-helix interchain hydrogen bonds in the region of the Gly replacement. These computational and experimental results lend insight into the complex molecular basis of the varying clinical severity of OI.

The Proline/Glycine-Rich Region of the Biofilm Adhesion Protein Aap Forms an Extended Stalk that Resists Compaction

Staphylococcus epidermidis is one of the primary bacterial species responsible for healthcare-associated infections. The most significant virulence factor for S. epidermidis is its ability to form a biofilm, which renders the bacteria highly resistant to host immune responses and antibiotic action. Intercellular adhesion within the biofilm is mediated by the accumulation-associated protein (Aap), a cell wall-anchored protein that self-assembles in a zinc-dependent manner. The C-terminal portion of Aap contains a 135-aa-long, proline/glycine-rich region (PGR) that has not yet been characterized. The region contains a set of 18 nearly identical AEPGKP repeats. Analysis of the PGR using biophysical techniques demonstrated the region is a highly extended, intrinsically disordered polypeptide with unusually high polyproline type II helix propensity. In contrast to many intrinsically disordered polypeptides, there was a minimal temperature dependence of the global conformational state of PGR in solution as measured by analytical ultracentrifugation and dynamic light scattering. Furthermore, PGR was resistant to conformational collapse or α-helix formation upon the addition of the osmolyte trimethylamine N-oxide or the cosolvent 2,2,2-trifluoroethanol. Collectively, these results suggest PGR functions as a resilient, extended stalk that projects the rest of Aap outward from the bacterial cell wall, promoting intercellular adhesion between cells in the biofilm. This work sheds light on regions of low complexity often found near the attachment point of bacterial cell wall-anchored proteins.

CD and NMR investigation of collagen peptides mimicking a pathological Gly-Ser mutation and a natural interruption in a similar highly charged sequence context

Even a single Gly substitution in the triple helix domain of collagen leads to pathological conditions while natural interruptions are suggested to play important functional roles. Two peptides-one mimicking a pathological Gly-Ser substitution (ERSEQ) and the other one modeling a similar natural interruption sequence (DRSER)-are designed to facilitate the comparison for elucidating the molecular basis of their different biological roles. CD and NMR investigation of peptide ERSEQ indicates a reduction of the thermal stability and disruption of hydrogen bonding at the Ser mutation site, providing a structural basis of the OI disease resulting from the Gly-Ser mutation in the highly charged RGE environment. Both CD and NMR real-time folding results indicate that peptide ERSEQ displays a comparatively slower folding rate than peptide DRSER, suggesting that the Gly-Ser mutation may lead to a larger interference in folding than the natural interruption in a similar RSE context. Our studies suggest that unlike the rigid GPO environment, the abundant R(K)GE(D) motif may provide a more flexible sequence environment that better accommodates mutations as well as interruptions, while the electrostatic interactions contribute to its stability. These results shed insight into the molecular features of the highly charged motif and may aid the design of collagen biomimetic peptides containing important biological sites.

Local amino acid sequence patterns dominate the heterogeneous phenotype for the collagen connective tissue disease Osteogenesis Imperfecta resulting from Gly mutations

Osteogenesis Imperfecta (OI), a hereditary connective tissue disease in collagen that arises from a single Gly → X mutation in the collagen chain, varies widely in phenotype from perinatal lethal to mild. It is unclear why there is such a large variation in the severity of the disease considering the repeating (Gly-X-Y)n sequence and the uniform rod-like structure of collagen. We systematically evaluate the effect of local (Gly-X-Y)n sequence around the mutation site on OI phenotype using integrated bio-statistical approaches, including odds ratio analysis and decision tree modeling. We show that different Gly → X mutations have different local sequence patterns that are correlated with lethal and nonlethal phenotypes providing a mechanism for understanding the sensitivity of local context in defining lethal and non-lethal OI. A number of important trends about which factors are related to OI phenotypes are revealed by the bio-statistical analyses; most striking is the complementary relationship between the placement of Pro residues and small residues and their correlation to OI phenotype. When Pro is present or small flexible residues are absent nearby a mutation site, the OI case tends to be lethal; when Pro is present or small flexible residues are absent further away from the mutation site, the OI case tends to be nonlethal. The analysis also reveals the dominant role of local sequence around mutation sites in the Major Ligand Binding Regions that are primarily responsible for collagen binding to its receptors and shows that non-lethal mutations are highly predicted by local sequence considerations alone whereas lethal mutations are not as easily predicted and may be a result of more complex interactions. Understanding the sequence determinants of OI mutations will enhance genetic counseling and help establish which steps in the collagen hierarchy to target for drug therapy.

Osteogenesis imperfecta missense mutations in collagen: structural consequences of a glycine to alanine replacement at a highly charged site

Glycine is required as every third residue in the collagen triple helix, and a missense mutation leading to the replacement of even one Gly in the repeating (Gly-Xaa-Yaa)(n) sequence with a larger residue leads to a pathological condition. Gly to Ala missense mutations are highly underrepresented in osteogenesis imperfecta (OI) and other collagen diseases, suggesting that the smallest replacement residue, Ala, might cause the least structural perturbation and mildest clinical consequences. The relatively small number of Gly to Ala mutation sites that do lead to OI must have some unusual features, such as greater structural disruption because of local sequence environment or location at a biologically important site. Here, peptides are used to model a severe OI case in which a Gly to Ala mutation is found within a highly stabilizing Lys-Gly-Asp sequence environment. Nuclear magnetic resonance, circular dichroism, and differential scanning calorimetry studies indicate this Gly to Ala replacement leads to a substantial loss of triple-helix stability and nonequivalence of the Ala residues in the three chains such that only one of the three Ala residues is capable of forming a good backbone hydrogen bond. Examination of reported OI Gly to Ala mutations suggests their preferential location at known collagen binding sites, and we propose that structural defects caused by Ala replacements may lead to pathology when they interfere with interactions.

Osteogenesis imperfecta model peptides: incorporation of residues replacing Gly within a triple helix achieved by renucleation and local flexibility

Missense mutations, which replace one Gly with a larger residue in the repeating sequence of the type I collagen triple helix, lead to the hereditary bone disorder osteogenesis imperfecta (OI). Previous studies suggest that these mutations may interfere with triple-helix folding. NMR was used to investigate triple-helix formation in a series of model peptides where the residue replacing Gly, as well as the local sequence environment, was varied. NMR measurement of translational diffusion coefficients allowed the identification of partially folded species. When Gly was replaced by Ala, the Ala residue was incorporated into a fully folded triple helix, whereas replacement of Gly by Ser or Arg resulted in the presence of some partially folded species, suggesting a folding barrier. Increasing the triple-helix stability of the sequence N-terminal to a Gly-to-Ser replacement allowed complete triple-helix folding, whereas with the substitution of Arg, with its large side chain, the peptide achieved full folding only after flexible residues were introduced N-terminal to the mutation site. These studies shed light on the factors important for accommodation of Gly mutations within the triple helix and may relate to the varying severity of OI.