Scalable approaches for generating, validating and incorporating data from high-throughput functional assays to improve clinical variant classification.

As the adoption and scope of genetic testing continue to expand, interpreting the clinical significance of DNA sequence variants at scale remains a formidable challenge, with a high proportion classified as variants of uncertain significance (VUSs). Genetic testing laboratories have historically relied, in part, on functional data from academic literature to support variant classification. High-throughput functional assays or multiplex assays of variant effect (MAVEs), designed to assess the effects of DNA variants on protein stability and function, represent an important and increasingly available source of evidence for variant classification, but their potential is just beginning to be realized in clinical lab settings. Here, we describe a framework for generating, validating and incorporating data from MAVEs into a semi-quantitative variant classification method applied to clinical genetic testing. Using single-cell gene expression measurements, cellular evidence models were built to assess the effects of DNA variation in 44 genes of clinical interest. This framework was also applied to models for an additional 22 genes with previously published MAVE datasets. In total, modeling data was incorporated from 24 genes into our variant classification method. These data contributed evidence for classifying 4043 observed variants in over 57,000 individuals. Genetic testing laboratories are uniquely positioned to generate, analyze, validate, and incorporate evidence from high-throughput functional data and ultimately enable the use of these data to provide definitive clinical variant classifications for more patients.

Structural basis of antibody inhibition and chemokine activation of the human CC chemokine receptor 8.

The C-C motif chemokine receptor 8 (CCR8) is a class A G-protein coupled receptor that has emerged as a promising therapeutic target in cancer. Targeting CCR8 with an antibody has appeared to be an attractive therapeutic approach, but the molecular basis for chemokine-mediated activation and antibody-mediated inhibition of CCR8 are not fully elucidated. Here, we obtain an antagonist antibody against human CCR8 and determine structures of CCR8 in complex with either the antibody or the endogenous agonist ligand CCL1. Our studies reveal characteristic antibody features allowing recognition of the CCR8 extracellular loops and CCL1-CCR8 interaction modes that are distinct from other chemokine receptor - ligand pairs. Informed by these structural insights, we demonstrate that CCL1 follows a two-step, two-site binding sequence to CCR8 and that antibody-mediated inhibition of CCL1 signaling can occur by preventing the second binding event. Together, our results provide a detailed structural and mechanistic framework of CCR8 activation and inhibition that expands our molecular understanding of chemokine - receptor interactions and offers insight into the development of therapeutic antibodies targeting chemokine GPCRs.

Transient water wires mediate selective proton transport in designed channel proteins.

Selective proton transport through proteins is essential for forming and using proton gradients in cells. Protons are conducted along hydrogen-bonded 'wires' of water molecules and polar side chains, which, somewhat surprisingly, are often interrupted by dry apolar stretches in the conduction pathways, inferred from static protein structures. Here we hypothesize that protons are conducted through such dry spots by forming transient water wires, often highly correlated with the presence of the excess protons in the water wire. To test this hypothesis, we performed molecular dynamics simulations to design transmembrane channels with stable water pockets interspersed by apolar segments capable of forming flickering water wires. The minimalist designed channels conduct protons at rates similar to viral proton channels, and they are at least 106-fold more selective for H+ over Na+. These studies inform the mechanisms of biological proton conduction and the principles for engineering proton-conductive materials.

Interaction specificity of clustered protocadherins inferred from sequence covariation and structural analysis.

Clustered protocadherins, a large family of paralogous proteins that play important roles in neuronal development, provide an important case study of interaction specificity in a large eukaryotic protein family. A mammalian genome has more than 50 clustered protocadherin isoforms, which have remarkable homophilic specificity for interactions between cellular surfaces. A large antiparallel dimer interface formed by the first 4 extracellular cadherin (EC) domains controls this interaction. To understand how specificity is achieved between the numerous paralogs, we used a combination of structural and computational approaches. Molecular dynamics simulations revealed that individual EC interactions are weak and undergo binding and unbinding events, but together they form a stable complex through polyvalency. Strongly evolutionarily coupled residue pairs interacted more frequently in our simulations, suggesting that sequence coevolution can inform the frequency of interaction and biochemical nature of a residue interaction. With these simulations and sequence coevolution, we generated a statistical model of interaction energy for the clustered protocadherin family that measures the contributions of all amino acid pairs at the interface. Our interaction energy model assesses specificity for all possible pairs of isoforms, recapitulating known pairings and predicting the effects of experimental changes in isoform specificity that are consistent with literature results. Our results show that sequence coevolution can be used to understand specificity determinants in a protein family and prioritize interface amino acid substitutions to reprogram specific protein-protein interactions.

Structures in multiple conformations reveal distinct transition metal and proton pathways in an Nramp transporter.

Nramp family transporters-expressed in organisms from bacteria to humans-enable uptake of essential divalent transition metals via an alternating-access mechanism that also involves proton transport. We present high-resolution structures of Deinococcus radiodurans (Dra)Nramp in multiple conformations to provide a thorough description of the Nramp transport cycle by identifying the key intramolecular rearrangements and changes to the metal coordination sphere. Strikingly, while metal transport requires cycling from outward- to inward-open states, efficient proton transport still occurs in outward-locked (but not inward-locked) DraNramp. We propose a model in which metal and proton enter the transporter via the same external pathway to the binding site, but follow separate routes to the cytoplasm, which could facilitate the co-transport of two cationic species. Our results illustrate the flexibility of the LeuT fold to support a broad range of substrate transport and conformational change mechanisms.

Applications of sequence coevolution in membrane protein biochemistry.

Recently, protein sequence coevolution analysis has matured into a predictive powerhouse for protein structure and function. Direct methods, which use global statistical models of sequence coevolution, have enabled the prediction of membrane and disordered protein structures, protein complex architectures, and the functional effects of mutations in proteins. The field of membrane protein biochemistry and structural biology has embraced these computational techniques, which provide functional and structural information in an otherwise experimentally-challenging field. Here we review recent applications of protein sequence coevolution analysis to membrane protein structure and function and highlight the promising directions and future obstacles in these fields. We provide insights and guidelines for membrane protein biochemists who wish to apply sequence coevolution analysis to a given experimental system.

Antiparallel protocadherin homodimers use distinct affinity- and specificity-mediating regions in cadherin repeats 1-4.

Protocadherins (Pcdhs) are cell adhesion and signaling proteins used by neurons to develop and maintain neuronal networks, relying on trans homophilic interactions between their extracellular cadherin (EC) repeat domains. We present the structure of the antiparallel EC1-4 homodimer of human PcdhγB3, a member of the γ subfamily of clustered Pcdhs. Structure and sequence comparisons of α, ß, and γ clustered Pcdh isoforms illustrate that subfamilies encode specificity in distinct ways through diversification of loop region structure and composition in EC2 and EC3, which contains isoform-specific conservation of primarily polar residues. In contrast, the EC1/EC4 interface comprises hydrophobic interactions that provide non-selective dimerization affinity. Using sequence coevolution analysis, we found evidence for a similar antiparallel EC1-4 interaction in non-clustered Pcdh families. We thus deduce that the EC1-4 antiparallel homodimer is a general interaction strategy that evolved before the divergence of these distinct protocadherin families.

Investigation of the interactions between Pt(II) and Pd(II) derivatives of 5,10,15,20-tetrakis (N-methyl-4-pyridyl) porphyrin and G-quadruplex DNA.

G-quadruplexes are non-canonical DNA structures formed by guanine-rich DNA sequences that are implicated in cancer and aging. Understanding how small molecule ligands interact with quadruplexes is essential both to the development of novel anticancer therapeutics and to the design of new quadruplex-selective probes needed for elucidation of quadruplex biological functions. In this work, UV-visible, fluorescence, and circular dichroism spectroscopies, fluorescence resonance energy transfer (FRET) melting assays, and resonance light scattering were used to investigate how the Pt(II) and Pd(II) derivatives of the well-studied 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin (TMPyP4) interact with quadruplexes formed by the human telomeric DNA, Tel22, and by the G-rich sequences from oncogene promoters. Our results suggest that Pt- and PdTMPyP4 interact with Tel22 via efficient π-π stacking with a binding affinity of 10(6)-10(7) M(-1). Under porphyrin excess, PtTMPyP4 aggregates using Tel22 as a template; the aggregates reach maximum size at [PtTMPyP4]/[Tel22] ~8 and dissolve at [PtTMPyP4]/[Tel22] ≤ 2. FRET assays reveal that both porphyrins are excellent stabilizers of human telomeric DNA, with stabilization temperature of 30.7 ± 0.6 °C for PtTMPyP4 and 30.9 ± 0.4 °C for PdTMPyP4 at [PtTMPyP4]/[Tel22] = 2 in K(+) buffer, values significantly higher as compared to those for TMPyP4. The porphyrins display modest selectivity for quadruplex vs. duplex DNA, with selectivity ratios of 150 and 330 for Pt- and PdTMPyP4, respectively. This selectivity was confirmed by observed 'light switch' effect: fluorescence of PtTMPyP4 increases significantly in the presence of a variety of DNA secondary structures, yet the strongest effect is produced by quadruplex DNA.

Structure and Sequence Analyses of Clustered Protocadherins Reveal Antiparallel Interactions that Mediate Homophilic Specificity.

Clustered protocadherin (Pcdh) proteins mediate dendritic self-avoidance in neurons via specific homophilic interactions in their extracellular cadherin (EC) domains. We determined crystal structures of EC1-EC3, containing the homophilic specificity-determining region, of two mouse clustered Pcdh isoforms (PcdhγA1 and PcdhγC3) to investigate the nature of the homophilic interaction. Within the crystal lattices, we observe antiparallel interfaces consistent with a role in trans cell-cell contact. Antiparallel dimerization is supported by evolutionary correlations. Two interfaces, located primarily on EC2-EC3, involve distinctive clustered Pcdh structure and sequence motifs, lack predicted glycosylation sites, and contain residues highly conserved in orthologs but not paralogs, pointing toward their biological significance as homophilic interaction interfaces. These two interfaces are similar yet distinct, reflecting a possible difference in interaction architecture between clustered Pcdh subfamilies. These structures initiate a molecular understanding of clustered Pcdh assemblies that are required to produce functional neuronal networks.

N-methylmesoporphyrin IX fluorescence as a reporter of strand orientation in guanine quadruplexes.

Guanine quadruplexes (GQ) are four-stranded DNA structures formed by guanine-rich DNA sequences. The formation of GQs inhibits cancer cell growth, although the detection of GQs in vivo has proven difficult, in part because of their structural diversity. The development of GQ-selective fluorescent reporters would enhance our ability to quantify the number and location of GQs, ultimately advancing biological studies of quadruplex relevance and function. N-methylmesoporphyrin IX (NMM) interacts selectively with parallel-stranded GQs; in addition, its fluorescence is sensitive to the presence of DNA, making this ligand a possible candidate for a quadruplex probe. In the present study, we investigated the effect of DNA secondary structure on NMM fluorescence. We found that NMM fluorescence increases by about 60-fold in the presence of parallel-stranded GQs and by about 40-fold in the presence of hybrid GQs. Antiparallel GQs lead to lower than 10-fold increases in NMM fluorescence. Single-stranded DNA, duplex, or i-motif, induce no change in NMM fluorescence. We conclude that NMM shows promise as a 'turn-on' fluorescent probe for detecting quadruplex structures, as well as for differentiating them on the basis of strand orientation.

Optimized end-stacking provides specificity of N-methyl mesoporphyrin IX for human telomeric G-quadruplex DNA.

N-methyl mesoporphyrin IX (NMM) is exceptionally selective for G-quadruplexes (GQ) relative to duplex DNA and, as such, has found a wide range of applications in biology and chemistry. In addition, NMM is selective for parallel versus antiparallel GQ folds, as was recently demonstrated in our laboratory. Here, we present the X-ray crystal structure of a complex between NMM and human telomeric DNA dAGGG(TTAGGG)(3), Tel22, determined in two space groups, P2(1)2(1)2 and P6, at 1.65 and 2.15 Å resolution, respectively. The former is the highest resolution structure of the human telomeric GQ DNA reported to date. The biological unit contains a Tel22 dimer of 5'-5' stacked parallel-stranded quadruplexes capped on both ends with NMM, supporting the spectroscopically determined 1:1 stoichiometry. NMM is capable of adjusting its macrocycle geometry to closely match that of the terminal G-tetrad required for efficient π-π stacking. The out-of-plane N-methyl group of NMM fits perfectly into the center of the parallel GQ core where it aligns with potassium ions. In contrast, the interaction of the N-methyl group with duplex DNA or antiparallel GQ would lead to steric clashes that prevent NMM from binding to these structures, thus explaining its unique selectivity. On the basis of the biochemical data, binding of NMM to Tel22 does not rely on relatively nonspecific electrostatic interactions, which characterize most canonical GQ ligands, but rather it is hydrophobic in nature. The structural features observed in the NMM-Tel22 complex described here will serve as guidelines for developing new quadruplex ligands that have excellent affinity and precisely defined selectivity.

Interaction of human telomeric DNA with N-methyl mesoporphyrin IX.

The remarkable selectivity of N-methyl mesoporphyrin IX (NMM) for G-quadruplexes (GQs) is long known, however its ability to stabilize and bind GQs has not been investigated in detail. Through the use of circular dichroism, UV-visible spectroscopy and fluorescence resonance energy transfer (FRET) melting assay we have shown that NMM stabilizes human telomeric DNA dAG(3)(TTAG(3))(3) (Tel22) and is selective for its parallel conformation to which it binds in 1:1 stoichiometry with a binding constant of ≈ 1.0 × 10(5)M(-1). NMM does not interact with an antiparallel conformation of Tel22 in sodium buffer and is the second example in the literature, after TOxaPy, of a ligand with an excellent selectivity for a specific GQ structure. NMM's stabilizing ability toward predominantly parallel GQ conformation is universal: it stabilizes a variety of biologically relevant G-rich sequences including telomeres and oncogene promoters. The N-methyl group is integral for selectivity and stabilization, as the unmethylated analogue, mesoporphyrin IX, does not stabilize GQ DNA in FRET melting assays. Finally, NMM induces the isomerization of Tel22 into a structure with increased parallel component in K(+) but not in Na(+) buffer. The ability of NMM to cause structural rearrangement and efficient stabilization of Tel22 may bear biological significance.

Induction of G-quadruplex DNA structure by Zn(II) 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin.

G-quadruplexes (GQ) are formed by the association of guanine-rich stretches of DNA. Certain small molecules can influence kinetics and thermodynamics of this association. Understanding the mechanism of ligand-assisted GQ folding is necessary for the design of more efficient cancer therapeutics. The oligonucleotide d(TAGGG)(2) forms parallel bimolecular GQ in the presence of ≥66 mM K(+); GQs are not formed under Na(+), Li(+) or low K(+) conditions. The thermodynamic parameters for GQ folding at 60 µM oligonucleotide and 100 mM KCl are ΔH = -35 ± 2 kcal mol(-1) and ΔG(310) = -1.4 kcal mol(-1). Quadruplex [d(TAGGG)(2)](2) binds 2-3 K(+) ions with K(d) of 0.5 ± 0.2 mM. Our work addresses the question of whether metal free 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin (TMPyP4) and its Zn(II), Cu(II), and Pt(II) derivatives are capable of facilitating GQ folding of d(TAGGG)(2) from single stranded, or binding to preformed GQ, using UV-vis and circular dichroism (CD) spectroscopies. ZnTMPyP4 is unique among other porphyrins in its ability to induce GQ structure of d(TAGGG)(2), which also requires at least a low amount of potassium. ZnTMPyP4 binds with 2:1 stoichiometry possibly in an end-stacking mode with a ~10(6) M(-1) binding constant, determined through UV-vis and ITC titrations. This process is entropically driven and has ΔG(298) of -8.0 kcal mol(-1). TMPyP4 binds with 3:1 stoichiometry and K(a) of ~10(6) M(-1). ZnTMPyP4 and TMPyP4 are efficient stabilizers of [d(TAGGG)(2)](2) displaying ΔT(1/2) of 13.5 and 13.8 °C, respectively, at 1:2 GQ to porphyrin ratio; CuTMPyP4 shows a much weaker effect (ΔT(1/2) = 4.7 °C) and PtTMPyP4 is weakly destabilizing (ΔT(1/2) = -2.9 °C). The selectivity of ZnTMPyP4 for GQ versus dsDNA is comparable to that of TMPyP4. The ability of ZnTMPyP4 to bind and stabilize GQ, to induce GQ formation, and speed up its folding may suggest an important biological activity for this molecule.

Rudimentary G-quadruplex-based telomere capping in Saccharomyces cerevisiae.

Telomere capping conceals chromosome ends from exonucleases and checkpoints, but the full range of capping mechanisms is not well defined. Telomeres have the potential to form G-quadruplex (G4) DNA, although evidence for telomere G4 DNA function in vivo is limited. In budding yeast, capping requires the Cdc13 protein and is lost at nonpermissive temperatures in cdc13-1 mutants. Here, we use several independent G4 DNA-stabilizing treatments to suppress cdc13-1 capping defects. These include overexpression of three different G4 DNA binding proteins, loss of the G4 DNA unwinding helicase Sgs1, or treatment with small molecule G4 DNA ligands. In vitro, we show that protein-bound G4 DNA at a 3' overhang inhibits 5'→3' resection of a paired strand by exonuclease I. These findings demonstrate that, at least in the absence of full natural capping, G4 DNA can play a positive role at telomeres in vivo.