Predicting the clinical lethality of osteogenesis imperfecta from collagen glycine mutations

Osteogenesis imperfecta (OI), or brittle bone disease, often results from missense mutation of one of the conserved glycine residues present in the repeating Gly-X-Y sequence characterizing the triple-helical region of type I collagen. A composite model was developed for predicting the clinical lethality resulting from glycine mutations in the alpha1 chain of type I collagen. The lethality of mutations in which bulky amino acids are substituted for glycine is predicted by their position relative to the N-terminal end of the triple helix. The effect of a Gly --> Ser mutation is modeled by the relative thermostability of the Gly-X-Y triplet on the carboxy side of the triplet containing the substitution. This model also predicts the lethality of Gly --> Ser and Gly --> Cys mutations in the alpha2 chain of type I collagen. The model was validated with an independent test set of six novel Gly --> Ser mutations. The hypothesis derived from the model of an asymmetric interaction between a Gly --> Ser mutation and its neighboring residues was tested experimentally using collagen-like peptides. Consistent with the prediction, a significant decrease in stability, calorimetric enthalpy, and folding time was observed for a peptide with a low-stability triplet C-terminal to the mutation compared to a similar peptide with the low-stability triplet on the N-terminal side. The computational and experimental results together relate the position-specific effects of Gly --> Ser mutations to the local structural stability of collagen and lend insight into the etiology of OI.     

Detection of mutations in human type I collagen mRNA in osteogenesis imperfecta by indirect RNase protection

A method for detecting a wide variety of mutations within type I collagen has been developed and evaluated on a series of patients with osteogenesis imperfecta. RNA, extracted from the nuclear and cytoplasmic compartment of cultured fibroblasts from affected individuals, is hybridized with antisense single-stranded cDNA to the alpha 1(I) mRNA. The hybrid is digested with RNase A and T1 under varying degrees of stringency. The resistant RNA bands are separated by electrophoresis in agarose, transferred to nitrocellulose, and hybridized with antisense cRNA colinear with the protecting probe. This approach is capable of detecting previously defined mutations such as 252-base pair deletion and a 1-base pair mismatch within the alpha 1(I) mRNA. The method appears to be particularly useful in detecting abnormalities of RNA processing that behave as an insert or deletion within the mature mRNA. The procedure should be generally applicable for the identification and localization of any mutation within an entire gene if the gene is expressed as an RNA and a complete cDNA for the mRNA is available.     

Effect of osteogenesis imperfecta mutations on free energy of collagen model peptides: a molecular dynamics simulation

We have carried out stochastic boundary molecular dynamics simulations to estimate free energy changes for substitutions of Gly with Val, Arg and Trp residues in a collagen-like peptide. The relative free energy change differences of mutants containing a Val, an Arg and a Trp relative to the wild type are 5.7, 8.1 and 9.5 kcal/mol, respectively. The corresponding free energy change differences of mutants containing two mutated residues are on average 7.6, 10.5 and 14.7 kcal/mol, respectively. We show that the free energy change differences are correlated with the severity of OI from statistical analysis and mechanical properties of the individual tropocollagen molecules. This simulation result indicates an atomistic-level mechanistic understanding of the effect of OI mutations in terms of stability of the mutants relative to the wild type, which could eventually provide a new strategy for diagnosis and treatment of the disease.     

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

Molecular dynamics simulations were carried out to calculate the free energy change difference of two collagen-like peptide models for Gly --> 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 Gly --> 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.     

Contribution of salt bridges near the surface of a protein to the conformational stability

Salt bridges play important roles in the conformational stability of proteins. However, the effect of a surface salt bridge on the stability remains controversial even today; some reports have shown little contribution of a surface salt bridge to stability, whereas others have shown a favorable contribution. In this study, to elucidate the net contribution of a surface salt bridge to the conformational stability of a protein, systematic mutant human lysozymes, containing one Glu to Gln (E7Q) and five Asp to Asn mutations (D18N, D49N, D67N, D102N, and D120N) at residues where a salt bridge is formed near the surface in the wild-type structure, were examined. The thermodynamic parameters for denaturation between pH 2.0 and 4.8 were determined by use of a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. The denaturation Gibbs energy (DeltaG) of all mutant proteins was lower than that of the wild-type protein at pH 4, whereas there was little difference between them near pH 2. This is caused by the fact that the Glu and Asp residues are ionized at pH 4 but protonated at pH 2, indicating a favorable contribution of salt bridges to the wild-type structure at pH 4. Each contribution was not equivalent, but we found that the contributions correlate with the solvent inaccessibility of the salt bridges; the salt bridge contribution was small when 100% accessible, while it was about 9 kJ/mol if 100% inaccessible. This conclusion indicates how to reconcile a number of conflicting reports about role of surface salt bridges in protein stability. Furthermore, the effect of salts on surface salt bridges was also examined. In the presence of 0.2 M KCl, the stability at pH 4 decreased, and the differences in stability between the wild-type and mutant proteins were smaller than those in the absence of salts, indicating the compensation to the contribution of salt bridges with salts. Salt bridges with more than 50% accessibility did not contribute to the stability in the presence of 0.2 M KCl.     

The role of interchain disulphide bridges in the conformational stability of human immunoglobulin G1 subclass. Hydrogen-deuterium exchange studies.

The hydrogen-deuterium exchange data of human immunoglobulin G1 (IgG1) are interpreted by assuming fast fluctuations of the protein conformation, through which the peptide groups become exposed to the solvent. The probability of solvent exposure of peptide hydrogens reflects a rather loose conformation for native IgG in comparison with other globular proteins. The probability of solvent exposure is greater than 10(-3) for half of the peptide groups, which shows that the conformational transitions by which these groups are exposed to the solvent are accompanied by changes in standard free energy less than 17 kJ/mol (4 kcal/mol). In the range of pH 6.2-8.45, at 25 degrees C no gross conformational changes are reflected in the hydrogen-deuterium exchange behaviour of the native, the reduced-nonalkylated-reassociated and the reduced-S-alkylated-reassociated IgG1. No difference could be detected in the conformational stability of the native and reoxidised reassociated IgG1 proteins. The lack of inter-subunit disulphide bridges in S-alkylated-reassociated molecules results in an increased conformational motility. This destabilization of protein conformation affects about 90% of the peptide groups covered by the measurements, and corresponds to changes in standard free energy of 8 kJ/mol on the average.     

The severity of osteogenesis imperfecta: a comparison to the relative free energy differences of collagen model peptides

Molecular dynamics simulations were carried out to calculate free energy differences between the folded and unfolded states of wild type and mutant collagen model peptides. The calculated stability of the collagen models was compared with the severity of osteogenesis imperfecta. Free energy differences of Gly → Xaa (Xaa: Ser, Cys, Glu, and Asp) mutations between the wild type and the mutants at position 15 of the model peptide were 3.8, 4.2, 5.6, and 8.8 kcal/mol, respectively. The corresponding free energy differences of a second Gly mutation at the same position in different chains were, on average, 1.3, 1.5, 2.9, and 5.4 kcal/mol, respectively. Free energy simulations were also performed to estimate the relative stability between an oxidized form and a reduced form of the mutants containing two Cys residues, which indicated that the mutant of the collagen-like peptide containing an intramolecular disulfide bond was more stable than the mutant containing one Cys residue but less stable than the wild type. The calculated free energy differences between an oxidized and a reduced form of the mutants containing two Cys residues are 0.8 and 2.6 kcal/mol for the disulfide bonds between Chains A and B and between Chains A and C, respectively.     

Structural studies of lysyl-tRNA synthetase: conformational changes induced by substrate binding

Lysyl-tRNA synthetase is a member of the class II aminoacyl-tRNA synthetases and catalyses the specific aminoacylation of tRNA(Lys). The crystal structure of the constitutive lysyl-tRNA synthetase (LysS) from Escherichia coli has been determined to 2.7 A resolution in the unliganded form and in a complex with the lysine substrate. A comparison between the unliganded and lysine-bound structures reveals major conformational changes upon lysine binding. The lysine substrate is involved in a network of hydrogen bonds. Two of these interactions, one between the alpha-amino group and the carbonyl oxygen of Gly 216 and the other between the carboxylate group and the side chain of Arg 262, trigger a subtle and complicated reorganization of the active site, involving the ordering of two loops (residues 215-217 and 444-455), a change in conformation of residues 393-409, and a rotation of a 4-helix bundle domain (located between motif 2 and 3) by 10 degrees. The result of these changes is a closing up of the active site upon lysine binding.     

Spectral studies of conformational changes induced in beta-glucuronidase by saccharo-1,4-lactone.

Ultraviolet difference spectroscopy has been used to study the binding of the transition state analog saccharo-1,4-lactone to purified rat preputial gland beta-glucuronidase. At pH 4.5 (the pH optimum), the inhibitor induces a difference spectrum indicative of a change in the environment of tryptophyl residues. Based on the magnitude of the induced difference spectrum as a quantitative measure of inhibitor binding, a titration curve for saccharo-1,4-lactone was obtained. A Scatchard plot of the titration data indicates that 4 molecules of inhibitor bind to the enzyme tetramer at a K-I of 4 times 10-7 M. The inhibitor also induces a similar difference spectrum at pH 7.5, although the binding is considerably weaker at this pH than at pH 4.5. When the native enzyme at pH 4.5 is compared with the native enzyme at pH 7.5, a difference spectrum, distinct from that of the binding of saccharo-1,4-lactone, is observed, indicating that the enzyme exists in different conformations at these pH values. The indication that tryptophyl residues are perturbed upon binding of saccharo-1,4-lactone was supported by studies carried out with N-bromosuccinimide. At pH 4.3, this reagent was found to oxidize 6 tryptophyl residues in the native enzyme but only three in the saccharo-1,4-lactone-inhibited enzyme. A spectrophotometric titration of the enzyme indicated that of the 33 tyrosyl residues per subunit, only 5 to 6 ionize at the pK expected for free phenolic groups.     

Characterization of point mutations in the collagen COL1A1 and COL1A2 genes causing lethal perinatal osteogenesis imperfecta

Type I collagen mutations in a group of patients with lethal perinatal osteogenesis imperfecta were identified in fibroblast RNA by a new method which can detect, by chemical modification and cleavage, single mismatched bases in heteroduplexes formed between mRNA and normal cDNA probes. Control cDNA probes spanning the area of the pro-alpha 1(I) and pro-alpha 2(I) chains likely to contain the mutations were radioactively labeled and used to form heteroduplexes with total patient RNA. Treatment of these heteroduplexes with hydroxylamine followed by cleavage of the cDNA strand at reactive bases by piperidine identified mismatches in the pro-alpha 1(I) cDNA in four patients. In the fifth patient a mismatch was detected in the pro-alpha 2(I) cDNA. To characterize these mutations the regions containing the mismatches were amplified by the polymerase chain reaction, cloned, and sequenced. All were heterozygous single base mutations which led to the substitution of glycine residues in the helical region of the pro-alpha-chains. The substitutions were pro-alpha 1(I) Gly973 and Gly1006 to Val, Gly928 to Ala, Gly976 to Arg, and pro-alpha 2(I) Gly865 to Ser. These mutations emphasize the importance of the Gly-X-Y repeating amino acid sequence for normal collagen helix formation and function in the extracellular matrix.     

Perinatal lethal osteogenesis imperfecta (OI type II): a biochemically heterogeneous disorder usually due to new mutations in the genes for type I collagen.

To resolve uncertainty concerning the inheritance of the perinatal lethal form of osteogenesis imperfecta (OI type II), we collected family data and radiographs for 71 probands and analyzed the collagens synthesized by dermal fibroblastic cells cultured from 43 of the probands, 19 parental pairs, and single parents of each of four additional probands. In 65 families for which there were complete data on sibship size, there was recurrence of the OI type II phenotype in five families such that six (8.6%) of 70 sibs were affected. In two families with recurrence, the radiographic phenotype was milder than that for the remainder; and one of those families was consanguinous, suggesting autosomal recessive inheritance. In the remaining three families there was no evidence of consanguinity, but in one of them the structure was compatible with gonadal mosaicism in the mother. In studies of collagens synthesized by cells from 43 infants, we identified two probands with separate rearrangements in an allele of one of the genes of type I collagen; but in the rest there were subtle mutations that disrupted the normal triple-helix structure of type I collagen molecules. In two probands we identified de novo mutations; in 16 additional families cells from the parents made only normal collagens, compatible with new mutations in their offsprings. These findings indicate that the OI type II phenotype is biochemically heterogeneous, that the majority result from new dominant mutations in the genes encoding type I collagen, and that some recurrences can be accounted for by gonadal mosaicism in one of the parents.     

Two novelCOL1A1 mutations in patients with osteogenesis imperfecta (OI) affect the stability of the collagen type I triple-helix

Osteogenesis imperfecta (OI) is a bone dysplasia caused by mutations in theCOL1A1 andCOL1A2 genes. Although the condition has been intensely studied for over 25 years and recently over 800 novel mutations have been published, the relation between the location of mutations and clinical manifestation is poorly understood. Here we report missense mutations inCOL1A1 of several OI patients. Two novel mutations were found in the D1 period. One caused a substitution of glycine 200 by valine at the N-terminus of D1 in OI type I/IV, lowering collagen stability by 50% at 34°C. The other one was a substitution of valine 349 by phenylalanine at the C-terminus of D1 in OI type I, lowering collagen stability at 37.5°C. Two other mutations, reported before, changed amino residues in D4. One was a lethal substitution changing glycine 866 to serine in genetically identical twins with OI type II. That mutated amino acid was near the border of D3 and D4. The second mutation changed glycine 1040 to serine located at the border of D4 and D0.4, in a proband manifesting OI type III, and lowered collagen stability at 39°C (2°C lower than normal). Our results confirm the hypothesis on a critical role of the D1 and D4 regions in stabilization of the collagen triple-helix. The defect in D1 seemed to produce a milder clinical type of OI, whereas the defect in the C-terminal end of collagen type caused the more severe or lethal types of OI.     

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.     

Kinetic and structural studies of the allosteric conformational changes induced by binding of cAMP to the cAMP receptor protein from Escherichia coli

The cAMP receptor protein, allosterically activated by cAMP, regulates the expression of more than 100 genes in Escherichia coli. CRP is a homodimer of two-domain subunits. It has been suggested that binding of cAMP to CRP leads to a long-distance signal transduction from the N-terminal cAMP binding domain to the C-terminal domain of the protein responsible for interaction with specific sequences of DNA. In this study, the stopped-flow and time-resolved fluorescence lifetime measurements were used to observe the kinetics of the distance changes between the N-terminal and C-terminal domain of CRP induced by binding of cAMP to high-affinity binding sites. In these measurements, we used the constructed CRP heterodimer, which possesses a single Trp85 residue localized at the N-terminal domain of one CRP subunit, and fluorescently labeled by 1,5-I-AEDANS Cys178 localized at the C-terminal domain of the same subunit or at the opposite one. The F?rster resonance energy transfer method has been used to study the distance changes, induced by binding of cAMP, between Trp85 (fluorescence donor) and Cys178-AEDANS (fluorescence acceptor) in the CRP structure. The obtained results show that the allosteric transitions of CRP at micromolar cAMP concentrations follow the sequential binding model, in which binding of cAMP to high-affinity sites causes a 4 A movement of the C-terminal domain toward N-terminal domains of the protein, with kinetics faster than 2 ms, and CRP adopts the "closed" conformation. This fast process is followed by the slower reorientation of both CRP subunits.     

Probing subtle conformational changes induced by phosphorylation and point mutations in the TIR domains of TLR2 and TLR3

Extensive research performed on Toll‐like receptor (TLR) signaling has identified residues in the Toll/interleukin‐1 receptor (TIR) domains that are essential for its proper functioning. Among these residues, those in BB loop are particularly significant as single amino acid mutations in this region can cause drastic changes in downstream signaling. However, while the effect of these mutations on the function is well studied (like the P681H mutation in TLR2, the A795P mutation in TLR3, and the P714H mutation in TLR4), their influence on the dynamics and inter‐residue networks is not well understood. The effects of local perturbations induced by these mutations could propagate throughout the TIR domain, influencing interactions with other TIR domain‐containing proteins. The identification of these subtle changes in inter‐residue interactions can provide new insights and structural rationale for how single‐point mutations cause drastic changes in TIR–TIR interactions. We employed molecular dynamics simulations and protein structure network (PSN) analyses to investigate the structural transitions with special emphasis on TLR2 and TLR3. Our results reveal that phosphorylation of the Tyr 759 residue in the TIR domain of TLR3 introduces rigidity to its BB loop. Subtle differences in the intra BB loop hydrogen bonding network between TLR3 and TLR2 are also observed. The PSN analyses indicate that the TIR domain is highly connected and pinpoints key differences in the inter‐residue interactions between the wild‐type and mutant TIR domains, suggesting that TIR domain structure is prone to allosteric effects, consistent with the current view of the influence of allostery on TLR signaling.     

Sequence of ionophore conformational changes induced by a simulated membrane/water interface

We demonstrate how the progressive passage through an interface simulated by a linear increase in dielectric constant may mediate a reversible transconformation of Na⁺/X337A (lasalocid A) from a cyclic to an extended structure, During the passage through the interface, the complex adopts progressively a more extended conformation favorable to the Na⁺ complexation or decomplexation at the interface.     

Defects in the processing of procollagen to collagen are demonstrable in cultured fibroblasts from patients with the Ehlers-Danlos and osteogenesis imperfecta syndromes

This is a study of the processing of procollagen to collagen in cultures of skin and tendon fibroblasts. Processing was markedly increased by growing cells for 2-4 days postconfluence and then adding ascorbate to the medium for 2 days prior to labeling with [3H] proline. With this system, more than two-thirds of the pro-alpha chains of type I procollagen in the culture medium, and more than 90% of those in the cell layer, were rapidly processed to pC-alpha, pN-alpha, or alpha chains. Purified, exogenous procollagen was also rapidly processed in cell-free culture medium. The results showed for the first time that exogenous procollagen can be processed in conditioned cell-free medium. The system was then used to compare the processing of procollagen in the medium of normal fibroblasts, cells from one bovine and four human variants of osteogenesis imperfecta, and those from eight human variants of the Ehlers-Danlos syndrome. The cells could be divided into three groups, based on their ability to process type I procollagen: normal, consistently slow, and very slow. The cause of the decreased processing was shown to be associated with either a mutation causing a shortening of an alpha chain or decreased activity of procollagen N-proteinase in cell-free culture medium. Decreased processing of procollagen to collagen occurred with cultured fibroblasts from patients with different forms of both osteogenesis imperfecta and Ehlers-Danlos syndrome. Both of these disease syndromes are associated with abnormalities in the structure or metabolism of procollagen in fibrous connective tissues, bones, and teeth. The results show that defects in the structure, synthesis, or processing of procollagen are readily demonstrated with cultured fibroblasts.     

Kinetics of conformational changes induced by the binding of various metal ions to bovine alpha-lactalbumin.

The kinetic effects of the binding of various metal ions (Ca(2+), Cd(2+), Co(2+), Mg(2+), Mn(2+), Sr(2+) and Zn(2+)) to apo bovine alpha-lactalbumin has been monitored by means of stopped-flow fluorescence spectroscopy. Our results show that the measured rate constant for the binding of metal ions to the Ca(2+)-site increases with increasing binding constant. This is, however, not the case for metal ions binding to the Zn(2+)-site. The binding experiments performed at different temperatures allowed us to calculate the activation energy for the transition from the metal-free to the metal-loaded state of the protein. These values do not depend on the nature of the metal ion but are correlated with the type of binding site. As a result, we were able to demonstrate that Mg(2+), a metal ion which was thought to bind to the Ca(2+)-site, shows the same binding characteristics as Co(2+) and Zn(2+) and therefore most likely interacts with the residues belonging to the Zn(2+)-binding site.     

Fluorescence quenching and kinetic studies of conformational changes induced by DNA and cAMP binding to cAMP receptor protein from Escherichia coli

Cyclic AMP receptor protein (CRP) regulates the expression of more then 100 genes in Escherichia coli. It is known that the allosteric activation of CRP by cAMP involves a long-distance signal transmission from the N-terminal cAMP-binding domain to the C-terminal domain of CRP responsible for the interactions with specific sequences of DNA. In this report we have used a CRP mutant containing a single Trp13 located in the N-terminal domain of the protein. We applied the iodide and acrylamide fluorescence quenching method in order to study how different DNA sequences and cAMP binding induce the conformational changes in the CRP molecule. The results presented provide evidence for the occurrence of a long-distance conformational signal transduction within the protein from the C-terminal DNA-binding domain to the N-terminal domain of CRP. This conformational signal transmission depends on the promoter sequence. We also used the stopped-flow and Forster resonance energy transfer between labeled Cys178 of CRP and fluorescently labeled DNA sequences to study the kinetics of DNA-CRP interactions. The results thus obtained lead to the conclusion that CRP can exist in several conformational states and that their distribution is affected by binding of both the cAMP and of specific DNA sequences.