Computational Reconstruction of Multidomain Proteins Using Atomic Force Microscopy Data

SummaryClassical structural biology techniques face a great challenge to determine the structure at the atomic level of large and flexible macromolecules. We present a novel methodology that combines high-resolution AFM topographic images with atomic coordinates of proteins to assemble very large macromolecules or particles. Our method uses a two-step protocol: atomic coordinates of individual domains are docked beneath the molecular surface of the large macromolecule, and then each domain is assembled using a combinatorial search. The protocol was validated on three test cases: a simulated system of antibody structures; and two experimentally based test cases: Tobacco mosaic virus, a rod-shaped virus; and Aquaporin Z, a bacterial membrane protein. We have shown that AFM-intermediate resolution topography and partial surface data are useful constraints for building macromolecular assemblies. The protocol is applicable to multicomponent structures connected in the polypeptide chain or as disjoint molecules. The approach effectively increases the resolution of AFM beyond topographical information down to atomic-detail structures

Identification of rare DNA variants in mitochondrial disorders with improved array-based sequencing.

A common goal in the discovery of rare functional DNA variants via medical resequencing is to incur a relatively lower proportion of false positive base-calls. We developed a novel statistical method for resequencing arrays (SRMA, sequence robust multi-array analysis) to increase the accuracy of detecting rare variants and reduce the costs in subsequent sequence verifications required in medical applications. SRMA includes single and multi-array analysis and accounts for technical variables as well as the possibility of both low- and high-frequency genomic variation. The confidence of each base-call was ranked using two quality measures. In comparison to Sanger capillary sequencing, we achieved a false discovery rate of 2% (false positive rate 1.2 × 10⁻⁵, false negative rate 5%), which is similar to automated second-generation sequencing technologies. Applied to the analysis of 39 nuclear candidate genes in disorders of mitochondrial DNA (mtDNA) maintenance, we confirmed mutations in the DNA polymerase gamma POLG in positive control cases, and identified novel rare variants in previously undiagnosed cases in the mitochondrial topoisomerase TOP1MT, the mismatch repair enzyme MUTYH, and the apurinic-apyrimidinic endonuclease APEX2. Some patients carried rare heterozygous variants in several functionally interacting genes, which could indicate synergistic genetic effects in these clinically similar disorders

Genotyping with a 198 Mutation Arrayed Primer Extension Array for Hereditary Hearing Loss: Assessment of Its Diagnostic Value for Medical Practice

Molecular diagnostic testing of individuals with congenital sensorineural hearing loss typically begins with DNA sequencing of the GJB2 gene. If the cause of the hearing loss is not identified in GJB2, additional testing can be ordered. However, the step-wise analysis of several genes often results in a protracted diagnostic process. The more comprehensive Hereditary Hearing Loss Arrayed Primer Extension microarray enables analysis of 198 mutations across eight genes (GJB2, GJB6, GJB3, GJA1, SLC26A4, SLC26A5, MTRNR1 and MTTS1) in a single test. To evaluate the added diagnostic value of this microarray for our ethnically diverse patient population, we tested 144 individuals with congenital sensorineural hearing loss who were negative for biallelic GJB2 or GJB6 mutations. The array successfully detected all GJB2 changes previously identified in the study group, confirming excellent assay performance. Additional mutations were identified in the SLC26A4, SLC26A5 and MTRNR1 genes of 12/144 individuals (8.3%), four of whom (2.8%) had genotypes consistent with pathogenicity. These results suggest that the current format of this microarray falls short of adding diagnostic value beyond the customary testing of GJB2, perhaps reflecting the array's limitations on the number of mutations included for each gene, but more likely resulting from unknown genetic contributors to this phenotype. We conclude that mutations in other hearing loss associated genes should be incorporated in the array as knowledge of the etiology of hearing loss evolves. Such future modification of the flexible configuration of the Hereditary Hearing Loss Arrayed Primer Extension microarray would improve its impact as a diagnostic tool

Comparative functional characterization of novel non-syndromic GJB2 gene variant p.Gly45Arg and lethal syndromic variant p.Gly45Glu

We characterized a novel GJB2 missense variant, c.133G>A, p.Gly45Arg, and compared it with the only other variant at the same amino acid position of the connexin 26 protein (Cx26) reported to date: c.134G>A, p.Gly45Glu. Whereas both variants are associated with hearing loss and are dominantly inherited, p.Gly45Glu has been implicated in the rare fatal keratitis-ichthyosis-deafness (KID) syndrome, which results in cutaneous infections and septicemia with premature demise in the first year of life. In contrast, p.Gly45Arg appears to be non-syndromic. Subcellular localization experiments in transiently co-transfected HeLa cells demonstrated that Cx26-WT (wild-type) and p.Gly45Arg form gap junctions, whereas Cx26-WT with p.Gly45Glu protein does not. The substitution of a nonpolar amino acid glycine in wildtype Cx26 at position 45 with a negatively charged glutamic acid (acidic) has previously been shown to interfere with Ca2+ regulation of hemichannel gating and to inhibit the formation of gap junctions, resulting in cell death. The novel variant p.Gly45Arg, however, changes this glycine to a positively charged arginine (basic), resulting in the formation of dysfunctional gap junctions that selectively affect the permeation of negatively charged inositol 1,4,5-trisphosphate (IP3) and contribute to hearing loss. Cx26 p.Gly45Arg transfected cells, unlike cells transfected with p.Gly45Glu, thrived at physiologic Ca2+ concentrations, suggesting that Ca2+ regulation of hemichannel gating is unaffected in Cx26 p.Gly45Arg transfected cells. Thus, the two oppositely charged amino acids that replace the highly conserved uncharged glycine in p.Gly45Glu and p.Gly45Arg, respectively, produce strikingly different effects on the structure and function of the Cx26 protein

Computational Reconstruction of Multidomain Proteins Using Atomic Force Microscopy Data

International audienceClassical structural biology techniques face a great challenge to determine the structure at the atomic level of large and flexible macromolecules. We present a novel methodology that combines high-resolution AFM topographic images with atomic coordinates of proteins to assemble very large macromolecules or particles. Our method uses a two-step protocol: atomic coordinates of individual domains are docked beneath the molecular surface of the large macromolecule, and then each domain is assembled using a combinatorial search. The protocol was validated on three test cases: a simulated system of antibody structures; and two experimentally based test cases: Tobacco mosaic virus, a rod-shaped virus; and Aquaporin Z, a bacterial membrane protein. We have shown that AFM-intermediate resolution topography and partial surface data are useful constraints for building macromolecular assemblies. The protocol is applicable to multicomponent structures connected in the polypeptide chain or as disjoint molecules. The approach effectively increases the resolution of AFM beyond topographical information down to atomic-detail structures

Combining H/D exchange mass spectroscopy and computational docking reveals extended DNA-binding surface on uracil-DNA glycosylase

X-ray crystallography provides excellent structural data on protein–DNA interfaces, but crystallographic complexes typically contain only small fragments of large DNA molecules. We present a new approach that can use longer DNA substrates and reveal new protein–DNA interactions even in extensively studied systems. Our approach combines rigid-body computational docking with hydrogen/deuterium exchange mass spectrometry (DXMS). DXMS identifies solvent-exposed protein surfaces; docking is used to create a 3-dimensional model of the protein–DNA interaction. We investigated the enzyme uracil-DNA glycosylase (UNG), which detects and cleaves uracil from DNA. UNG was incubated with a 30 bp DNA fragment containing a single uracil, giving the complex with the abasic DNA product. Compared with free UNG, the UNG–DNA complex showed increased solvent protection at the UNG active site and at two regions outside the active site: residues 210–220 and 251–264. Computational docking also identified these two DNA-binding surfaces, but neither shows DNA contact in UNG–DNA crystallographic structures. Our results can be explained by separation of the two DNA strands on one side of the active site. These non-sequence-specific DNA-binding surfaces may aid local uracil search, contribute to binding the abasic DNA product and help present the DNA product to APE-1, the next enzyme on the DNA-repair pathway