Solution structure by NMR and molecular dynamics of a duplex containing a guanine opposite a N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide lesion

One- and two-dimensional NMR spectroscopy has been used combined with molecular dynamics to determine the fine structure of the DNA duplex 5'-d(AGGAGCCACG).d(CGTGGFTCCT) where F is the N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide residue which is a ring fragmentation product of thymine. The formamide deoxyribose exists as two isomers with respect to the orientation about the peptide bond. The two isomers (trans and cis) are observed in a ratio 3:2 in solution. For both species, the oligonucleotide adopts a globally B form structure although conformational changes are observed around the mismatch site. The formamide residue, whatever the isomer, is intrahelical and can pair with the guanine on the opposite strand with one hydrogen bond. For the cis isomer, the residue adopts a syn orientation and is able to form a second hydrogen bond with the guanine on the 5' side on the same strand. Off-resonance ROESY experiments have been used to investigate the chemical exchange observed at low temperature of the duplex. Conformational exchange has only been found for the oligonucleotide with the formamide residue in the trans conformation.     

Structure of an oligonucleotide containing a N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide residue facing a guanine

Formamide residue is a major oxidative DNA damage product from ionizing radiation on thymine residues in DNA. We report NMR and molecular modeling studies on a DNA duplex structure which contains guanine opposite formamide residue. Formamide residue exists as either the cis and trans isomer. For the trans and the cis isomers, we find that guanine and formamide are stacked inside the helix and are hydrogen bonded. The oligonucleotide adopts globally a B form structure for the two isomers. Conformational changes are observed between the two isomers.     

Thymine glycols and urea residues in M13 DNA constitute replicative blocks in vitro.

Thymine glycols were produced in M13 DNA in a concentration dependent manner by treating the DNA with osmium tetroxide (OsO4). For the formation of urea-containing M13 DNA, OsO4-oxidized DNA was hydrolyzed in alkali (pH 12) to convert the thymine glycols to urea residues. With both thymine glycol- and urea-containing M13 DNA, DNA synthesis catalyzed by Escherichia coli DNA polymerase I Klenow fragment was decreased in proportion to the number of damages present in the template DNA. Sequencing gel analysis of the products synthesized by E. coli DNA polymerase I and T4 DNA polymerase showed that DNA synthesis terminated opposite the putative thymine glycol site and at one nucleotide before the putative urea site. Substitution of manganese for magnesium in the reaction mix resulted in increased processivity of DNA synthesis so that a base was incorporated opposite urea. With thymine glycol-containing DNA, processivity in the presence of manganese was strongly dependent on the presence of a pyrimidine 5' to the thymine glycol in the template.     

Solution structure of a DNA duplex containing a formamide-adenine base pair

The N-(2-deoxy-beta3-D-erythro-pentofuranosyl) formamide residue results from a ring fragmentation product of thymine or cytosine. The presence of a formamide-adenine base pair in the sequence 5'd(AGGAACCACG).d(CGTGGFTCCT) has been studied by 1H and 31P nuclear magnetic resonance (NMR) and molecular dynamics. There are two possible isomers for the formamide side chain, either cis or trans. For each isomer, we observed an equilibrium in solution between two forms. First, a species where the formamide is intrahelical and paired with the facing adenine. For the cis isomer, the formamide is in a syn conformation and two hydrogen bonds with adenine are formed. The trans isomer is in an anti conformation and a single hydrogen bond is observed. In the second form, whatever the isomer, the formamide is rejected outside the helix, whereas the adenine remains inside.     

Initial denaturing conditions influence the slow folding phase of acylphosphatase associated with proline isomerization

The folding kinetics of human common-type acylphosphatase (cAcP) from its urea- and TFE-denatured states have been determined by stopped-flow fluorescence techniques. The refolding reaction from the highly unfolded state formed in urea is characterized by double exponential behavior that includes a slow phase associated with isomerism of the Gly53-Pro54 peptide bond. However, this slow phase is absent when refolding is initiated by dilution of the highly a-helical denatured state formed in the presence of 40% trifluoroethanol (TFE). NMR studies of a peptide fragment corresponding to residues Gly53-Gly69 of cAcP indicate that only the native-like trans isomer of the Gly-Pro peptide bond is significantly populated in the presence of TFE, whereas both the cis and trans isomers are found in an approximately 1:9 ratio for the peptide bond in aqueous solution. Molecular modeling studies in conjunction with NMR experiments suggest that the trans isomer of the Gly53-Pro54 peptide bond is stabilized in TFE by the formation of a nonnative-like hydrogen bond between the CO group of Gly53 and the NH group of Lys57. These results therefore reveal that a specific nonnative interaction in the denatured state can increase significantly the overall efficiency of refolding.     

Thymidine and Thymine Incorporation into Deoxyribonucleic Acid: Inhibition and Repression by Uridine of Thymidine Phosphorylase of Escherichia coli

Thymidine is poorly incorporated into deoxyribonucleic acid (DNA) of Escherichia coli. Its incorporation is greatly increased by uridine, which acts in two ways. Primarily, uridine competitively inhibits thymidine phosphorylase (E.C.2.4.4), and thereby prevents the degradation of thymidine to thymine which is not incorporated into normally growing E. coli. Uridine also inhibits induction of the enzyme by thymidine. It prevents the actual inducer, probably a deoxyribose phosphate, from being formed rather than competing for a site on the repressor. The inhibition of thymidine phosphorylase by uridine also accounts for inhibition by uracil compounds of thymine incorporation into thymine-requiring mutants. Deoxyadenosine also increases the incorporation of thymidine, by competitively inhibiting thymidine phosphorylase. Deoxyadenosine induces the enzyme, in contrast to uridine. But this is offset by a transfer of deoxyribose from deoxyadenosine to thymine. Thus, deoxyadenosine permits incorporation of thymine into DNA, even in cells induced for thymidine phosphorylase. This incorporation of thymine in the presence of deoxyadenosine did not occur in a thymidine phosphorylase-negative mutant; thus, the utilization of thymine seems to proceed by way of thymidine phosphorylase, followed by thymidine kinase. These results are consistent with the data of others in suggesting that wild-type E. coli cells fail to utilize thymine because they lack a pool of deoxyribose phosphates, the latter being necessary for conversion of thymine to thymidine by thymidine phosphorylase.     

A cis-prolyl peptide bond isomerization dominates the folding of the alpha subunit of Trp synthase,a TIM barrel protein

The cis/trans isomerization of prolyl peptide bonds has been suggested to dominate the folding of the alpha subunit of tryptophan synthase from Escherichia coli (alphaTS). To test the role of the unique cis isomer between Asp27 and Pro28, the folding properties of P28A, P28G and G(3)P28G, a three-glycine insertion mutant between Asp27 and Gly28, were investigated using urea as a denaturant. Circular dichroism analysis demonstrated that none of the mutations perturb the secondary structure significantly, although the aromatic side-chain packing is altered for P28A and P28G. All three mutant proteins inherited the three-state thermodynamic behavior observed in wild-type alphaTS, ensuring that the fundamental features of the energy surface are intact. Kinetic studies showed that neither alanine nor glycine substitutions at Pro28 results in the elimination of any slow-refolding phases. By contrast, the G(3)P28G mutant eliminates the fastest of the slow-refolding phases and one of the two unfolding phases. Double-jump experiments on G(3)P28G confirm the assignment of the missing refolding phase to the isomerization of the Asp27-Pro28 peptide bond. These results imply that the local stability conveyed by the tight, overlapping turns containing the cis peptide bond is sufficient to favor the cis isomer for several non-prolyl residues. The free energy required to drive the isomerization reaction is provided by the formation of the stable intermediate, demonstrating that the acquisition of structure and stability is required to induce subsequent rate-limiting steps in the folding of alphaTS.     

Crystal structures of alpha-lytic protease complexes with irreversibly bound phosphonate esters

The structures of the complexes with alpha-lytic protease of both phosphorus stereoisomers of N-[(2S)-2-[[[(1R)-1-[N-[(tert-butyloxycarbonyl)-L-alanyl-L-alanyl- L-prolyl]amino]-2-methylpropyl]-phenoxyphosphinyl]oxy]propanoyl]- L-alanine methyl ester, an analogue of the peptide Boc-Ala-Ala-Pro-Val-Ala-Ala where Val is replaced with an analogous phosphonate phenyl ester and the subsequent Ala is replaced with lactate, have been determined to high resolution (1.9 A) by X-ray crystallography. Both stereoisomers inactivate the enzyme but differ by a factor of 2 in the second-order rate constant for inactivation [Sampson, N. S., & Bartlett, P. A. (1991) Biochemistry (preceding paper in this issue)]. One isomer (B) forms a tetrahedral adduct in which the phosphonate phenyl ester is displaced by the active site serine (S195) and interacts with the enzyme across seven substrate recognition sites that span both sides of the scissile bond. Seven hydrogen bonds are formed with the enzyme, and 510 A2 of hydrophobic surface area is buried when the inhibitor interacts with the enzyme. Although two hydrogen bonds are gained by incorporation of two residues on the C-terminal side of the scissile bond into the inhibitor, there is very little adjustment in the structure of the enzyme in this region. Surprisingly, the active site histidine (H57) does not interact with the phosphonate, apparently because the phosphonate lacks negative charge in or near the oxyanion hole, and instead, the side chain rotates out of the active site cleft and hydrogen bonds with solvent. The other isomer (A) forms a mixture of two different tetrahedral adducts in the active site, both covalently bonded to Ser 195. One adduct, at approximately 58% occupancy, is exactly the same in structure as the complex formed with isomer B, and the other adduct, at 42% occupancy, has lost the two residues C-terminal to the scissile bond by hydrolysis. In the lower occupancy structure, His 57 does not rotate out of the active site and forms a hydrogen bond with the phosphonate oxygen instead. The structures of both complexes were insensitive to pH. As very little change in structure accompanies the histidine rotation, the complex with isomer B provides an excellent mimic for the structure of the transition state (or high-energy reaction intermediate) that spans both sides of the scissile bond.     

Role of base stacking and sequence context in the inhibition of yeast DNA polymerase eta by pyrene nucleotide

The Y family DNA polymerase yeast pol eta inserts pyrene deoxyribose monophosphate (dPMP) in preference to A opposite an abasic site, the 3'-T of a thymine dimer, and a normal T with almost equal efficiency. In contrast, pol A family polymerases such as Klenow fragment and T7 DNA polymerase only insert dPMP efficiently opposite an abasic site and the 3'-T of a thymine dimer but not opposite undamaged DNA. Pyrene nucleotide is also an efficient chain-terminating inhibitor of DNA synthesis by pol eta but not by Klenow fragment or T7 DNA polymerase. To better understand the origin of the efficiency and sequence specificity of dPMP insertion by pol eta, the kinetics of dPMP insertion opposite various templates have been determined. In one sequence context, the efficiency of dPMP insertion increases 4.6-fold opposite G < A < T < C, suggesting that the templating nucleotide modulates dPMP insertion efficiency by having to destack prior to dPTP binding. The efficiency of insertion of dPMP opposite T in the same sequence context increases 7-fold for primers terminating in G < A < C < T and is similar to that observed for nontemplated blunt-end extension, suggesting that stacking interactions between the pyrene and the primer terminus are also important. On heterogeneous templates, the average selectivity for dPMP insertion relative to the complementary dNMP decreases in the order of dAMP > dGMP > dTMP > dCMP, from a high of 5.8 when dAMP is to be inserted following a T to a low of 0.5 when dCMP is to be inserted following a C. The relative preference for dPMP insertion at a given site can be largely explained by the energetic cost of destacking the templating base and stacking of pyrene nucleotide relative to that of stacking and base pairing the complementary nucleotide. Thus, pyrene nucleotide represents a novel class of nucleotide-based chain-terminating DNA synthesis inhibitors whose base portion consists of a hydrophobic, non-hydrogen bonding, base-pair mimic.     

Sequence dependence for bypass of thymine glycols in DNA by DNA polymerase I.

Single-stranded phage DNAs containing thymine glycols were prepared by oxidation with osmium tetroxide (OsO4) and were used as templates for DNA synthesis by E. coli DNA polymerase I. The induction of thymine glycol lesions in DNA, as measured by immunoassay, quantitatively accounted for an inhibition of in vitro DNA synthesis on modified templates. Analysis of termination sites for synthesis by DNA polymerase I (Klenow fragment) showed that DNA synthesis terminated at most template thymine sites in OsO4-treated DNA, indicating that incorporation occurred opposite putative thymine glycols in DNA. Nucleotides 5' and 3' to putative thymine glycol sites affect the reaction, however, since termination was not observed at thymines in the sequence 5'-CTPur-3'. Conversion of thymine glycols to urea residues in DNA by alkali treatment caused termination of DNA synthesis one nucleotide 3' to template thymine sites, including thymines in the 5'-CTPur-3' sequence, showing that the effect of surrounding sequence is on the elongation reaction by DNA polymerase rather than differential damage induction by OsO4.     

The cysteine residues in the carboxy terminal domain of tropoelastin form an intrachain disulfide bond that stabilizes a loop structure and positively charged pocket.

Analysis of purified bovine tropoelastin with Ellman's reagent and [14C]iodoacetamide demonstrated that the only two cysteine residues in the molecule form an intrachain disulfide bond. Molecular modeling suggests that the cysteine residues are juxtaposed as the result of a tight turn that produces an antiparallel beta structure. Protruding from the C-terminal end of the turn is the sequence Arg-Lys-Arg-Lys which forms the floor of a positively charged pocket created by the extension of the arginine and lysine side chains on opposite sides of the peptide chain perpendicular to the plane of the turn. The side chain of a conserved lysine residue in the disulfide-bonded loop forms the top of the pocket. This positively charged pocket may define a binding site for acidic microfibrillar proteins that mediate elastic fiber assembly.     

Deletion of a single amino acid changes the folding of an apamin hybrid sequence peptide to that of endothelin

The solution conformations of a hybrid sequence peptide related to the bee venom peptide apamin have been determined using two-dimensional ¹H-nmr. Apamin is an 18 amino acid peptide containing a C-terminal helix that is stabilized by two disulfide bonds. The deletion of one residue (K4) of the N-terminal “scaffold” region of the apamin sequence results in a helical peptide, but with a change in the pairing of cysteines to form the disulfide cross links. The new disulfide arrangement is analogous to that of the vasoconstrictor peptide endothelin. Two sets of nmr resonances were observed for the apamin-deletion (AD) peptide, due to cis-trans isomerism at the A4-P5 peptide bond. The cis isomer of the AD peptide contains a tight turn in residues 3–6, which is required for formation of the α-helix in residues 7–15. Nuclear Overhauser effects observed for the trans AD peptide are not consistent with any single unique fold, indicating the presence of conformational averaging when the peptide adopts the trans form. Distance geometry calculations on the cis AD peptide reveal an α-helical structure that appears to be more like that of apamin than the crystal structure of human endothelin, despite the reversal of the disulfide pattern in the AD peptide from that of apamin to that of endothelin.© 1997 John Wiley & Sons, Inc. Biopoly 41 : 451–460, 1997     

Interaction of a bacterially expressed peptide from the receptor binding domain of Pseudomonas aeruginosa pili strain PAK with a cross-reactive antibody: conformation of the bound peptide

The C-terminal receptor binding region of Pseudomonas aeruginosa pilin protein strain PAK (residues 128-144) has been the target for the design of a vaccine effective against P. aeruginosa infections. We have recently cloned and expressed a (15)N-labeled PAK pilin peptide spanning residues 128-144 of the PAK pilin protein. The peptide exists as a major (trans) and minor (cis) species in solution, arising from isomerization around a central Ile(138)-Pro(139) peptide bond. The trans isomer adopts two well-defined turns in solution, a type I beta-turn spanning Asp(134)-Glu-Gln-Phe(137) and a type II beta-turn spanning Pro(139)-Lys-Gly-Cys(142). The cis isomer adopts only one well-defined type II beta-turn spanning Pro(139)-Lys-Gly-Cys(142) but displays evidence of a less ordered turn spanning Asp(132)-Gln-Asp-Glu(135). These turns have been implicated in cross-reactive antibody recognition. (15)N-edited NMR spectroscopy was used to study the binding of the (15)N-labeled PAK pilin peptide to an Fab fragment of a cross-reactive monoclonal antibody, PAK-13, raised against the intact PAK pilus. The results of these studies are as follows: the trans and cis isomers bind with similar affinity to the Fab, despite their different topologies; both isomers maintain the conformational integrity of their beta-turns when bound; binding leads to the preferential stabilization of the first turn over the second turn in each isomer; and binding leads to the perturbation of resonances within regions of the trans and cis backbone that undergo microsecond to millisecond motions. These slow motions may play a role in induced fit binding of the first turn to Fab PAK-13, which would allow the same antibody combining site to accommodate either trans or cis topology. More importantly for vaccine design, these motions may also play a role in the development of a broad-spectrum vaccine capable of generating an antibody therapeutic effective against the multiple strains of P. aeruginosa.     

Escherichia coli thioredoxin folds into two compact forms of different stability to urea denaturation

The urea-induced denaturation of Escherichia coli thioredoxin and thioredoxin variants has been examined by electrophoresis on urea gradient slab gels by the method of Creighton [Creighton, T. (1986) Methods Enzymol. 131, 156-172]. Thioredoxin has only two cysteine residues, and these form a redox-active disulfide at the active site. Oxidized thioredoxin-S2 and reduced thioredoxin-(SH)2 each show two folded isomers with a large difference in stability to urea denaturation. The difference in stability is greater for the isomers of oxidized than for the isomers of reduced thioredoxin. At 2 degrees C, the urea concentrations at the denaturation midpoint are approximately 8 and 4.3 M for the oxidized isomers and 4.8 and 3.7 M for the reduced isomers. The difference between the gel patterns of samples applied in native versus denaturing buffer, and at 2 and 25 degrees C, is characteristic for the involvement of a cis-proline-trans-proline isomerization. The data very strongly suggest that the two folded forms of different stabilities correspond to the cis and trans isomers of the highly conserved Pro 76 peptide bond, which is cis in the crystal structure of oxidized thioredoxin. Urea gel experiments with the mutant thioredoxin P76A, with alanine substituted for proline at position 76, corroborate this interpretation. The electrophoretic banding pattern diagnostic for an involvement of proline isomerization in urea denaturation is not observed for oxidized P76A. In broad estimates of delta G degree for the native-denatured transition, the difference in delta G degree (no urea) between the putative cis and trans isomers of the Ile 75-Pro 76 peptide bond is approximately 3 kcal/mol for oxidized thioredoxin and approximately 1.5 kcal/mol for reduced thioredoxin. Since cis oxidized thioredoxin is much more stable than trans, folded oxidized thioredoxin is essentially all cis. In folded reduced thioredoxin, cis and trans interconvert slowly, on the minute time scale at 2 and 25 degrees C. In the absence of urea, the folded reduced thioredoxin is less than a few percent trans. Three additional mutants with additions or substitutions at the active site also show electrophoresis banding patterns consistent with a difference in stability between cis and trans isomers.     

The dominant interaction between peptide and urea is electrostatic in nature: a molecular dynamics simulation study

The conformational equilibrium of a blocked valine peptide in water and aqueous urea solution is studied using molecular dynamics simulations. Pair correlation functions indicate enhanced concentration of urea near the peptide. Stronger hydrogen bonding of urea-peptide compared to water-peptide is observed with preference for helical conformation. The potential of mean force, computed using umbrella sampling, shows only small differences between urea and water solvation that are difficult to quantify. The changes in solvent structure around the peptide are explained by favorable electrostatic interactions (hydrogen bonds) of urea with the peptide backbone. There is no evidence for significant changes in hydrophobic interactions in the two conformations of the peptide in urea solution. Our simulations suggest that urea denatures proteins by preferentially forming hydrogen bonds to the peptide backbone, reducing the barrier for exposing protein residues to the solvent, and reaching the unfolded state.     

Substrate specificity at the P1' site of Escherichia coli OmpT under denaturing conditions

Though OmpT has been reported to mainly cleave the peptide bond between consecutive basic amino acids, we identified more precise substrate specificity by using a series of modified substrates, termed PRX fusion proteins, consisting of 184 residues. The cleavage site of the substrate PRR was Arg140-Arg141 and the modified substrates PRX substituted all 19 natural amino acids at the P1' site instead of Arg141. OmpT under denaturing conditions (in the presence of 4 M urea) cleaved not only between two consecutive basic amino acids but also at the carboxyl side of Arg140 except for the Arg140-Asp141, -Glu141, and -Pro141 pairs. In addition to Arg140 at the P1 site, similar results were obtained when Lys140 was substituted into the P1 site. In the absence of urea, an aspartic acid residue at the P1' site was unfavorable for OmpT cleavage of synthetic decapeptides, the enzyme showed a preference for a dibasic site.     

Use of N-chlorosuccinimide/urea for the selective cleavage of tryptophanyl peptide bonds in proteins. Cytochrome c.

The conditions and utility of the N-chlorosuccinimide/urea (NCS/urea) reagent for the selective cleavage of tryptophanyl peptide bonds in proteins is demonstrated with cytochrome c. At low concentrations of NCS/urea the oxidation of thioether side chains in cytochrome c is the predominant reaction. Methionyl residues are oxidized to sulfoxide and the heme-thioether bridge is partially cleaved. At 10-fold excess of NCS/urea reagent, cleavage of the tryptophanyl peptide bond is optimal at approximately 50% yield in several species of cytochrome c studied. Analytical data on isolated horse cytochrome c peptide fragments demonstrate lack of modification and cleavage at tyrosyl and histidyl residues. However, at high concentrations of NCS/urea reagent (30-fold) unexpected conversions of methionine to sulfone and cysteine to cysteic acid in intact proteins are observed. This is in contradistinction to the absence of sulfone in NCS/urea-reacted amino acid mixtures. The mechanisms of halogenation and cleavage by N-bromosuccinimide, N-iodosuccinimide, and N-chlorosuccinimide are discussed. It is porposed that the selectivity with respect to halogenation by N-chlorosuccinimide is due to the insignificant participation of molecular chlorine in the NCS/urea reaction. A mechanism of halogenation and cleavage by NCS at tryptophan is also offered.     

Atypical abasic sites generated by neocarzinostatin at sequence-specific cytidylate residues in oligodeoxynucleotides

Neocarzinostatin chromophore produces alkali-labile, abasic sites at cytidylate residues in AGC sequences in oligonucleotides in their duplex form. Glutathione is the preferred thiol activator of the drug in the formation of these lesions. The phosphodiester linkages on each side of the abasic site are intact, but when treated with alkali, breaks are formed with phosphate moieties at each end. Similar properties are exhibited by the abasic lesions produced at the purine residue to which the C in AGC is base-paired on the complementary strand. The abasic sites at C residues differ from those produced by acid-induced depurination in the much greater lability of the phosphodiester linkages on both sides of the deoxyribose, in the inability of NaBH4 to prevent alkali-induced cleavage, and in the relative resistance to apurinic/apyrimidinic endonucleases. The importance of DNA microstructure in determining attack site specificity in abasic site formation at C residues is shown not only by the requirement for the sequence AGC but also by the findings that substitution of G by I 5' to the C decreases the attack at C, whereas placement of an I opposite the C markedly enhances the reaction. Quantitation of the abstraction of 3H into the drug from C residues in AGC specifically labeled in the deoxyribose at C-5' or C-1',2' suggests that, in contrast to the attack at C-5' in the induction of direct strand breaks at T residues, abasic site formation at C residues may involve attack at C-1'. Each type of lesion may exist on the complementary strands of the same DNA molecule, forming a double-stranded lesion.     

The importance of anchorage in determining a strained protein loop conformation.

We examine the role of the conformational restriction imposed by constrained ends of a protein loop on the determination of a strained loop conformation. The Lys 116-Pro 117 peptide bond of staphylococcal nuclease A exists in equilibrium between the cis and trans isomers. The folded protein favors the strained cis isomer with an occupancy of 90%. This peptide bond is contained in a solvent-exposed, flexible loop of residues 112-117 whose ends are anchored by Val 111 and Asn 118. Asn 118 is constrained by 2 side-chain hydrogen bonds. We investigate the importance of this constraint by replacing Asn 118 with aspartate, alanine, and glycine. We found that removing 1 or more of the hydrogen bonds observed in Asn 118 stabilizes the trans configuration over the cis configuration. By protonating the Asp 118 side chain of N118D through decreased pH, the hydrogen bonding character of Asp 118 approached that of Asn 118 in nuclease A, and the cis configuration was stabilized relative to the trans configuration. These data suggest that the rigid anchoring of the loop end is important in establishing the strained cis conformation. The segment of residues 112-117 in nuclease A provides a promising model system for study of the basic principles that determine polypeptide conformations. Such studies could be useful in the rational design or redesign of protein molecules.     

Severity of osteogenesis imperfecta and structure of a collagen-like peptide modeling a lethal mutation site

We show that there are correlations between the severities of osteogenesis imperfecta (OI) phenotypes and changes in the residues near the mutation site. Our results show the correlations between the severity of various forms of the inherited disease OI and alteration of residues near the site of OI causing mutations. Among our many observed correlations are particularly striking ones between the presence of nearby proline residues and lethal mutations, and the presence of nearby alanines residues and nonlethal mutations. We investigated the possibility that these correlations have a structural basis using molecular dynamics simulations of collagen-like molecules designed to mimic the site of a lethal OI mutation in collagen type I. Our significant finding is that interchain hydrogen bonding is greatly affected by variations in residue type. We found that the strength of hydrogen bond networks between backbone atoms on different chains depends on the local residue sequence and is weaker in proline-rich regions of the molecule. We also found that an alanine at a site near an OI mutation causes less structural disruption than a proline, and that residue side chains also form interchain hydrogen bonds with frequencies that are dependent on residue type. For example, arginine side chains form strong hydrogen bonds with the backbone of the subsequent peptide chain, while lysine and glutamine less frequently form similar hydrogen bonds. This decrease in the observed hydrogen bond frequency correlates with a decrease in the experimentally determined thermal stability. We contrasted general structural properties of model collagen peptides with and without the mutation to examine the effect of the single-point mutation on the surrounding residues.