Prevalence Estimates of Predicted Pathogenic COL4A3-COL4A5 Variants in a Population Sequencing Database and Their Implications for Alport Syndrome.

BACKGROUND: The reported prevalence of Alport syndrome varies from one in 5000 to one in 53,000 individuals. This study estimated the frequencies of predicted pathogenic COL4A3-COL4A5 variants in sequencing databases of populations without known kidney disease. METHODS: Predicted pathogenic variants were identified using filtering steps based on the ACMG/AMP criteria, which considered collagen IV α3-α5 position 1 Gly to be critical domains. The population frequencies of predicted pathogenic COL4A3-COL4A5 variants were then determined per mean number of sequenced alleles. Population frequencies for compound heterozygous and digenic combinations were calculated from the results for heterozygous variants. RESULTS: COL4A3-COL4A5 variants resulting in position 1 Gly substitutions were confirmed to be associated with hematuria (for each, P<0.001). Predicted pathogenic COL4A5 variants were found in at least one in 2320 individuals. p.(Gly624Asp) represented nearly half (16 of 33, 48%) of the variants in Europeans. Most COL4A5 variants (54 of 59, 92%) had a biochemical feature that potentially mitigated the clinical effect. The predicted pathogenic heterozygous COL4A3 and COL4A4 variants affected one in 106 of the population, consistent with the finding of thin basement membrane nephropathy in normal donor kidney biopsy specimens. Predicted pathogenic compound heterozygous variants occurred in one in 88,866 individuals, and digenic variants in at least one in 44,793. CONCLUSIONS: The population frequencies for Alport syndrome are suggested by the frequencies of predicted pathogenic COL4A3-COL4A5 variants, but must be adjusted for the disease penetrance of individual variants and for the likelihood of already diagnosed disease and non-Gly substitutions. Disease penetrance may depend on other genetic and environmental factors.

A systematic survey of the Cys2His2 zinc finger DNA-binding landscape.

Cys2His2 zinc fingers (C2H2-ZFs) comprise the largest class of metazoan DNA-binding domains. Despite this domain's well-defined DNA-recognition interface, and its successful use in the design of chimeric proteins capable of targeting genomic regions of interest, much remains unknown about its DNA-binding landscape. To help bridge this gap in fundamental knowledge and to provide a resource for design-oriented applications, we screened large synthetic protein libraries to select binding C2H2-ZF domains for each possible three base pair target. The resulting data consist of >160 000 unique domain-DNA interactions and comprise the most comprehensive investigation of C2H2-ZF DNA-binding interactions to date. An integrated analysis of these independent screens yielded DNA-binding profiles for tens of thousands of domains and led to the successful design and prediction of C2H2-ZF DNA-binding specificities. Computational analyses uncovered important aspects of C2H2-ZF domain-DNA interactions, including the roles of within-finger context and domain position on base recognition. We observed the existence of numerous distinct binding strategies for each possible three base pair target and an apparent balance between affinity and specificity of binding. In sum, our comprehensive data help elucidate the complex binding landscape of C2H2-ZF domains and provide a foundation for efforts to determine, predict and engineer their DNA-binding specificities.

Predicting DNA recognition by Cys2His2 zinc finger proteins.

MOTIVATION: Cys(2)His(2) zinc finger (ZF) proteins represent the largest class of eukaryotic transcription factors. Their modular structure and well-conserved protein-DNA interface allow the development of computational approaches for predicting their DNA-binding preferences even when no binding sites are known for a particular protein. The 'canonical model' for ZF protein-DNA interaction consists of only four amino acid nucleotide contacts per zinc finger domain. RESULTS: We present an approach for predicting ZF binding based on support vector machines (SVMs). While most previous computational approaches have been based solely on examples of known ZF protein-DNA interactions, ours additionally incorporates information about protein-DNA pairs known to bind weakly or not at all. Moreover, SVMs with a linear kernel can naturally incorporate constraints about the relative binding affinities of protein-DNA pairs; this type of information has not been used previously in predicting ZF protein-DNA binding. Here, we build a high-quality literature-derived experimental database of ZF-DNA binding examples and utilize it to test both linear and polynomial kernels for predicting ZF protein-DNA binding on the basis of the canonical binding model. The polynomial SVM outperforms previously published prediction procedures as well as the linear SVM. This may indicate the presence of dependencies between contacts in the canonical binding model and suggests that modification of the underlying structural model may result in further improved performance in predicting ZF protein-DNA binding. Overall, this work demonstrates that methods incorporating information about non-binding and relative binding of protein-DNA pairs have great potential for effective prediction of protein-DNA interactions. AVAILABILITY: An online tool for predicting ZF DNA binding is available at http://compbio.cs.princeton.edu/zf/.

Y-position cysteine substitution in type I collagen (alpha1(I) R888C/p.R1066C) is associated with osteogenesis imperfecta/Ehlers-Danlos syndrome phenotype.

The most common mutations in type I collagen causing types II-IV osteogenesis imperfecta (OI) result in substitution for glycine in a Gly-Xaa-Yaa triplet by another amino acid. We delineated a Y-position substitution in a small pedigree with a combined OI/Ehlers-Danlos Syndrome (EDS) phenotype, characterized by moderately decreased DEXA z-score (-1.3 to -2.6), long bone fractures, and large-joint hyperextensibility. Affected individuals have an alpha1(I)R888C (p.R1066C) substitution in one COL1A1 allele. Polyacrylamide gel electrophoresis (PAGE) of [(3)H]-proline labeled steady-state collagen reveals slight overmodification of the alpha1(I) monomer band, much less than expected for a substitution of a neighboring glycine residue, and a faint alpha1(I) dimer. Dimers form in about 10% of proband type I collagen. Dimer formation is inefficient compared to a possible 25%, probably because the SH-side chains have less proximity in this Y-position than when substituting for a glycine. Theoretical stability calculations, differential scanning calorimetry (DSC) thermograms, and thermal denaturation curves showed only weak local destabilization from the Y-position substitution in one or two chains of a collagen helix, but greater destabilization is seen in collagen containing dimers. Y-position collagen dimers cause kinking of the helix, resulting in a register shift that is propagated the full length of the helix and causes resistance to procollagen processing by N-proteinase. Collagen containing the Y-position substitution is incorporated into matrix deposited in culture, including immaturely and maturely cross-linked fractions. In vivo, proband dermal fibrils have decreased density and increased diameter compared to controls, with occasional aggregate formation. This report on Y-position substitutions in type I collagen extends the range of phenotypes caused by nonglycine substitutions and shows that, similar to X- and Y-position substitutions in types II and III collagen, the phenotypes resulting from nonglycine substitutions in type I collagen are distinct from those caused by glycine substitutions.

Self-association of collagen triple helic peptides into higher order structures.

Interest in self-association of peptides and proteins is motivated by an interest in the mechanism of physiologically higher order assembly of proteins such as collagen as well as the mechanism of pathological aggregation such as beta-amyloid formation. The triple helical form of (Pro-Hyp-Gly)(10), a peptide that has proved a useful model for molecular features of collagen, was found to self-associate, and its association properties are reported here. Turbidity experiments indicate that the triple helical peptide self-assembles at neutral pH via a nucleation-growth mechanism, with a critical concentration near 1 mM. The associated form is more stable than individual molecules by about 25 degrees C, and the association is reversible. The rate of self-association increases with temperature, supporting an entropically favored process. After self-association, (Pro-Hyp-Gly)(10) forms branched filamentous structures, in contrast with the highly ordered axially periodic structure of collagen fibrils. Yet a number of characteristics of triple helix assembly for the peptide resemble those of collagen fibril formation. These include promotion of fibril formation by neutral pH and increasing temperature; inhibition by sugars; and a requirement for hydroxyproline. It is suggested that these similar features for peptide and collagen self-association are based on common lateral underlying interactions between triple helical molecules mediated by hydrogen-bonded hydration networks involving hydroxyproline.

Prediction of collagen stability from amino acid sequence.

An algorithm was derived to relate the amino acid sequence of a collagen triple helix to its thermal stability. This calculation is based on the triple helical stabilization propensities of individual residues and their intermolecular and intramolecular interactions, as quantitated by melting temperature values of host-guest peptides. Experimental melting temperature values of a number of triple helical peptides of varying length and sequence were successfully predicted by this algorithm. However, predicted T(m) values are significantly higher than experimental values when there are strings of oppositely charged residues or concentrations of like charges near the terminus. Application of the algorithm to collagen sequences highlights regions of unusually high or low stability, and these regions often correlate with biologically significant features. The prediction of stability from sequence indicates an understanding of the major forces maintaining this protein motif. The use of highly favorable KGE and KGD sequences is seen to complement the stabilizing effects of imino acids in modulating stability and may become dominant in the collagenous domains of bacterial proteins that lack hydroxyproline. The effect of single amino acid mutations in the X and Y positions can be evaluated with this algorithm. An interactive collagen stability calculator based on this algorithm is available online.

Equilibrium thermal transitions of collagen model peptides.

The folding of collagen in vitro is very slow and presents difficulties in reaching equilibrium, a feature that may have implications for in vivo collagen function. Peptides serve as good model systems for examining equilibrium thermal transitions in the collagen triple helix. Investigations were carried out to ascertain whether a range of synthetic triple-helical peptides of varying sequences can reach equilibrium, and whether the triple helix to unfolded monomer transition approximates a two-state model. The thermal transitions for all peptides studied are fully reversible given sufficient time. Isothermal experiments were carried out to obtain relaxation times at different temperatures. The slowest relaxation times, on the order of 10-15 h, were observed at the beginning of transitions, and were shown to result from self-association limited by the low concentration of free monomers, rather than cis-trans isomerization. Although the fit of the CD equilibrium transition curves and the concentration dependence of T(m) values support a two-state model, the more rigorous comparison of the calorimetric enthalpy to the van't Hoff enthalpy indicates the two-state approximation is not ideal. Previous reports of melting curves of triple-helical host-guest peptides are shown to be a two-state kinetic transition, rather than an equilibrium transition.

Fluorescence determination of tryptophan side-chain accessibility and dynamics in triple-helical collagen-like peptides.

We report tryptophan fluorescence measurements of emission intensity, iodide quenching, and anisotropy that describe the environment and dynamics at X and Y sites in stable collagen-like peptides of sequence (Gly-X-Y)(n). About 90% of tryptophans at both sites have similar solvent exposed fluorescence properties and a lifetime of 8.5-9 ns. Analysis of anisotropy decays using an associative model indicates that these long lifetime populations undergo rapid depolarizing motion with a 0.5 ns correlation time; however, the extent of fast motion at the Y site is considerably less than the essentially unrestricted motion at the X site. About 10% of tryptophans at both sites have a shorter ( approximately 3 ns) lifetime indicating proximity to a protein quenching group; these minor populations are immobile on the peptide surface, depolarizing only by overall trimer rotation. Iodide quenching indicates that tryptophans at the X site are more accessible to solvent. Side chains at X sites are more solvent accessible and considerably more mobile than residues at Y sites and can more readily fluctuate among alternate intermolecular interactions in collagen fibrils. This fluorescence analysis of collagen-like peptides lays a foundation for studies on the structure, dynamics, and function of collagen and of triple-helical junctions in gelatin gels.

Peptide investigations of pairwise interactions in the collagen triple-helix.

Pairwise interactions have been studied for the major secondary structures in proteins. The present work extends the characterization of interactions between side-chains to the context of a collagen triple-helix. In this study, the most frequent Gly-X-Y tripeptide sequences in collagen are characterized in terms of interchain interactions between non-imino acid X and Y residues, through the use of host-guest peptides and statistical frequency analysis. Stabilities predicted on the basis of additivity show good agreement with experimental values for almost half of the peptides, indicating a lack of interaction. A small number of peptides have a stability lower than predicted, while a larger number are more stable than expected. Of all triplets containing residues of opposite charge, only Gly-Lys-Asp and Gly-Arg-Asp exhibit stabilizing electrostatic interactions, and these pairs are found together preferentially in collagens. Repulsion of like charges is observed in Gly-Arg-Lys, Gly-Lys-Arg, and Gly-Glu-Asp sequences, and a small degree of hydrophobic stabilization was observed for the Gly-Leu-Leu guest triplet. The data reported here help clarify basic principles of triple-helix stability. In addition, the experimentally determined stabilities of the tripeptide units found most frequently in collagens constitute a database useful for predicting triple-helix stability in peptides, collagens and other triple-helix-containing proteins.