Co-transcriptional splicing facilitates transcription of gigantic genes.

Although introns are typically tens to thousands of nucleotides, there are notable exceptions. In flies as well as humans, a small number of genes contain introns that are more than 1000 times larger than typical introns, exceeding hundreds of kilobases (kb) to megabases (Mb). It remains unknown why gigantic introns exist and how cells overcome the challenges associated with their transcription and RNA processing. The Drosophila Y chromosome contains some of the largest genes identified to date: multiple genes exceed 4Mb, with introns accounting for over 99% of the gene span. Here we demonstrate that co-transcriptional splicing of these gigantic Y-linked genes is important to ensure successful transcription: perturbation of splicing led to the attenuation of transcription, leading to a failure to produce mature mRNA. Cytologically, defective splicing of the Y-linked gigantic genes resulted in disorganization of transcripts within the nucleus suggestive of entanglement of transcripts, likely resulting from unspliced long RNAs. We propose that co-transcriptional splicing maintains the length of nascent transcripts of gigantic genes under a critical threshold, preventing their entanglement and ensuring proper gene expression. Our study reveals a novel biological significance of co-transcriptional splicing.

Constrained Mutational Sampling of Amino Acids in HIV-1 Protease Evolution.

The evolution of HIV-1 protein sequences should be governed by a combination of factors including nucleotide mutational probabilities, the genetic code, and fitness. The impact of these factors on protein sequence evolution is interdependent, making it challenging to infer the individual contribution of each factor from phylogenetic analyses alone. We investigated the protein sequence evolution of HIV-1 by determining an experimental fitness landscape of all individual amino acid changes in protease. We compared our experimental results to the frequency of protease variants in a publicly available data set of 32,163 sequenced isolates from drug-naïve individuals. The most common amino acids in sequenced isolates supported robust experimental fitness, indicating that the experimental fitness landscape captured key features of selection acting on protease during viral infections of hosts. Amino acid changes requiring multiple mutations from the likely ancestor were slightly less likely to support robust experimental fitness than single mutations, consistent with the genetic code favoring chemically conservative amino acid changes. Amino acids that were common in sequenced isolates were predominantly accessible by single mutations from the likely protease ancestor. Multiple mutations commonly observed in isolates were accessible by mutational walks with highly fit single mutation intermediates. Our results indicate that the prevalence of multiple-base mutations in HIV-1 protease is strongly influenced by mutational sampling.

Molecular Determinants of Epistasis in HIV-1 Protease: Elucidating the Interdependence of L89V and L90M Mutations in Resistance.

Protease inhibitors have the highest potency among antiviral therapies against HIV-1 infections, yet the virus can evolve resistance. Darunavir (DRV), currently the most potent Food and Drug Administration-approved protease inhibitor, retains potency against single-site mutations. However, complex combinations of mutations can confer resistance to DRV. While the interdependence between mutations within HIV-1 protease is key for inhibitor potency, the molecular mechanisms that underlie this control remain largely unknown. In this study, we investigated the interdependence between the L89V and L90M mutations and their effects on DRV binding. These two mutations have been reported to be positively correlated with one another in HIV-1 patient-derived protease isolates, with the presence of one mutation making the probability of the occurrence of the second mutation more likely. The focus of our investigation is a patient-derived isolate, with 24 mutations that we call "KY"; this variant includes the L89V and L90M mutations. Three additional KY variants with back-mutations, KY(V89L), KY(M90L), and the KY(V89L/M90L) double mutation, were used to experimentally assess the individual and combined effects of these mutations on DRV inhibition and substrate processing. The enzymatic assays revealed that the KY(V89L) variant, with methionine at residue 90, is highly resistant, but its catalytic function is compromised. When a leucine to valine mutation at residue 89 is present simultaneously with the L90M mutation, a rescue of catalytic efficiency is observed. Molecular dynamics simulations of these DRV-bound protease variants reveal how the L90M mutation induces structural changes throughout the enzyme that undermine the binding interactions.

The bone-specific Runx2-P1 promoter displays conserved three-dimensional chromatin structure with the syntenic Supt3h promoter.

Three-dimensional organization of chromatin is fundamental for transcriptional regulation. Tissue-specific transcriptional programs are orchestrated by transcription factors and epigenetic regulators. The RUNX2 transcription factor is required for differentiation of precursor cells into mature osteoblasts. Although organization and control of the bone-specific Runx2-P1 promoter have been studied extensively, long-range regulation has not been explored. In this study, we investigated higher-order organization of the Runx2-P1 promoter during osteoblast differentiation. Mining the ENCODE database revealed interactions between Runx2-P1 and Supt3h promoters in several non-mesenchymal human cell lines. Supt3h is a ubiquitously expressed gene located within the first intron of Runx2. These two genes show shared synteny across species from humans to sponges. Chromosome conformation capture analysis in the murine pre-osteoblastic MC3T3-E1 cell line revealed increased contact frequency between Runx2-P1 and Supt3h promoters during differentiation. This increase was accompanied by enhanced DNaseI hypersensitivity along with RUNX2 and CTCF binding at the Supt3h promoter. Furthermore, interplasmid-3C and luciferase reporter assays showed that the Supt3h promoter can modulate Runx2-P1 activity via direct association. Taken together, our data demonstrate physical proximity between Runx2-P1 and Supt3h promoters, consistent with their syntenic nature. Importantly, we identify the Supt3h promoter as a potential regulator of the bone-specific Runx2-P1 promoter.

Genome-wide co-occupancy of AML1-ETO and N-CoR defines the t(8;21) AML signature in leukemic cells.

BACKGROUND: Many leukemias result from chromosomal rearrangements. The t(8;21) chromosomal translocation produces AML1-ETO, an oncogenic fusion protein that compromises the function of AML1, a transcription factor critical for myeloid cell differentiation. Because of the pressing need for new therapies in the treatment of acute myleoid leukemia, we investigated the genome-wide occupancy of AML1-ETO in leukemic cells to discover novel regulatory mechanisms involving AML-ETO bound genes. RESULTS: We report the co-localization of AML1-ETO with the N-CoR co-repressor to be primarily on genomic regions distal to transcriptional start sites (TSSs). These regions exhibit over-representation of the motif for PU.1, a key hematopoietic regulator and member of the ETS family of transcription factors. A significant discovery of our study is that genes co-occupied by AML1-ETO and N-CoR (e.g., TYROBP and LAPTM5) are associated with the leukemic phenotype, as determined by analyses of gene ontology and by the observation that these genes are predominantly up-regulated upon AML1-ETO depletion. In contrast, the AML1-ETO/p300 gene network is less responsive to AML1-ETO depletion and less associated with the differentiation block characteristic of leukemic cells. Furthermore, a substantial fraction of AML1-ETO/p300 co-localization occurs near TSSs in promoter regions associated with transcriptionally active loci. CONCLUSIONS: Our findings establish a novel and dominant t(8;21) AML leukemia signature characterized by occupancy of AML1-ETO/N-CoR at promoter-distal genomic regions enriched in motifs for myeloid differentiation factors, thus providing mechanistic insight into the leukemic phenotype.

Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors.

Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) has become the dominant technique for mapping transcription factor (TF) binding regions genome-wide. We performed an integrative analysis centered around 457 ChIP-seq data sets on 119 human TFs generated by the ENCODE Consortium. We identified highly enriched sequence motifs in most data sets, revealing new motifs and validating known ones. The motif sites (TF binding sites) are highly conserved evolutionarily and show distinct footprints upon DNase I digestion. We frequently detected secondary motifs in addition to the canonical motifs of the TFs, indicating tethered binding and cobinding between multiple TFs. We observed significant position and orientation preferences between many cobinding TFs. Genes specifically expressed in a cell line are often associated with a greater occurrence of nearby TF binding in that cell line. We observed cell-line-specific secondary motifs that mediate the binding of the histone deacetylase HDAC2 and the enhancer-binding protein EP300. TF binding sites are located in GC-rich, nucleosome-depleted, and DNase I sensitive regions, flanked by well-positioned nucleosomes, and many of these features show cell type specificity. The GC-richness may be beneficial for regulating TF binding because, when unoccupied by a TF, these regions are occupied by nucleosomes in vivo. We present the results of our analysis in a TF-centric web repository Factorbook (http://factorbook.org) and will continually update this repository as more ENCODE data are generated.

MECP2 directly interacts with RNA polymerase II to modulate transcription in human neurons.

Mutations in the methyl-DNA-binding protein MECP2 cause the neurodevelopmental disorder Rett syndrome (RTT). How MECP2 contributes to transcriptional regulation in normal and disease states is unresolved; it has been reported to be an activator and a repressor. We describe here the first integrated CUT&Tag, transcriptome, and proteome analyses using human neurons with wild-type (WT) and mutant MECP2 molecules. MECP2 occupies CpG-rich promoter-proximal regions in over four thousand genes in human neurons, including a plethora of autism risk genes, together with RNA polymerase II (RNA Pol II). MECP2 directly interacts with RNA Pol II, and genes occupied by both proteins showed reduced expression in neurons with MECP2 patient mutations. We conclude that MECP2 acts as a positive cofactor for RNA Pol II gene expression at many neuronal genes that harbor CpG islands in promoter-proximal regions and that RTT is due, in part, to the loss of gene activity of these genes in neurons.

Targeted therapies prime oncogene-driven lung cancers for macrophage-mediated destruction.

Macrophage immune checkpoint inhibitors, such as anti-CD47 antibodies, show promise in clinical trials for solid and hematologic malignancies. However, the best strategies to use these therapies remain unknown, and ongoing studies suggest they may be most effective when used in combination with other anticancer agents. Here, we developed an unbiased, high-throughput screening platform to identify drugs that render lung cancer cells more vulnerable to macrophage attack, and we found that therapeutic synergy exists between genotype-directed therapies and anti-CD47 antibodies. In validation studies, we found that the combination of genotype-directed therapies and CD47 blockade elicited robust phagocytosis and eliminated persister cells in vitro and maximized antitumor responses in vivo. Importantly, these findings broadly applied to lung cancers with various RTK/MAPK pathway alterations - including EGFR mutations, ALK fusions, or KRASG12C mutations. We observed downregulation of ß2-microglobulin and CD73 as molecular mechanisms contributing to enhanced sensitivity to macrophage attack. Our findings demonstrate that dual inhibition of the RTK/MAPK pathway and the CD47/SIRPa axis is a promising immunotherapeutic strategy. Our study provides strong rationale for testing this therapeutic combination in patients with lung cancers bearing driver mutations.

SARS-CoV-2 infection of human pluripotent stem cell-derived vascular cells reveals smooth muscle cells as key mediators of vascular pathology during infection.

Although respiratory symptoms are the most prevalent disease manifestation of infection by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), nearly 20% of hospitalized patients are at risk for thromboembolic events 1 . This prothrombotic state is considered a key factor in the increased risk of stroke, which has been observed clinically during both acute infection and long after symptoms have cleared 2 . Here we developed a model of SARS-CoV-2 infection using human-induced pluripotent stem cell-derived endothelial cells, pericytes, and smooth muscle cells to recapitulate the vascular pathology associated with SARS-CoV-2 exposure. Our results demonstrate that perivascular cells, particularly smooth muscle cells (SMCs), are a specifically susceptible vascular target for SARS-CoV-2 infection. Utilizing RNA sequencing, we characterized the transcriptomic changes accompanying SARS-CoV-2 infection of SMCs, and endothelial cells (ECs). We observed that infected human SMCs shift to a pro-inflammatory state and increase the expression of key mediators of the coagulation cascade. Further, we showed human ECs exposed to the secretome of infected SMCs produce hemostatic factors that can contribute to vascular dysfunction, despite not being susceptible to direct infection. The findings here recapitulate observations from patient sera in human COVID-19 patients and provide mechanistic insight into the unique vascular implications of SARS-CoV-2 infection at a cellular level.

Targeted therapies prime oncogene-driven lung cancers for macrophage-mediated destruction.

Macrophage immune checkpoint inhibitors, such as anti-CD47 antibodies, show promise in clinical trials for solid and hematologic malignancies. However, the best strategies to use these therapies remain unknown and ongoing studies suggest they may be most effective when used in combination with other anticancer agents. Here, we developed a novel screening platform to identify drugs that render lung cancer cells more vulnerable to macrophage attack, and we identified therapeutic synergy exists between genotype-directed therapies and anti-CD47 antibodies. In validation studies, we found the combination of genotype-directed therapies and CD47 blockade elicited robust phagocytosis and eliminated persister cells in vitro and maximized anti-tumor responses in vivo. Importantly, these findings broadly applied to lung cancers with various RTK/MAPK pathway alterations-including EGFR mutations, ALK fusions, or KRASG12C mutations. We observed downregulation of ß2-microglobulin and CD73 as molecular mechanisms contributing to enhanced sensitivity to macrophage attack. Our findings demonstrate that dual inhibition of the RTK/MAPK pathway and the CD47/SIRPa axis is a promising immunotherapeutic strategy. Our study provides strong rationale for testing this therapeutic combination in patients with lung cancers bearing driver mutations.

Fragile X Syndrome Patient-Derived Neurons Developing in the Mouse Brain Show FMR1-Dependent Phenotypes.

BACKGROUND: Fragile X syndrome (FXS) is characterized by physical abnormalities, anxiety, intellectual disability, hyperactivity, autistic behaviors, and seizures. Abnormal neuronal development in FXS is poorly understood. Data on patients with FXS remain scarce, and FXS animal models have failed to yield successful therapies. In vitro models do not fully recapitulate the morphology and function of human neurons. METHODS: To mimic human neuron development in vivo, we coinjected neural precursor cells derived from FXS patient-derived induced pluripotent stem cells and neural precursor cells derived from corrected isogenic control induced pluripotent stem cells into the brain of neonatal immune-deprived mice. RESULTS: The transplanted cells populated the brain and a proportion differentiated into neurons and glial cells. Immunofluorescence and single and bulk RNA sequencing analyses showed accelerated maturation of FXS neurons after an initial delay. Additionally, we found increased percentages of Arc- and Egr-1-positive FXS neurons and wider dendritic protrusions of mature FXS striatal medium spiny neurons. CONCLUSIONS: This transplantation approach provides new insights into the alterations of neuronal development in FXS by facilitating physiological development of cells in a 3-dimensional context.

A maternally programmed intergenerational mechanism enables male offspring to make piRNAs from Y-linked precursor RNAs in Drosophila.

In animals, PIWI-interacting RNAs (piRNAs) direct PIWI proteins to silence complementary targets such as transposons. In Drosophila and other species with a maternally specified germline, piRNAs deposited in the egg initiate piRNA biogenesis in the progeny. However, Y chromosome loci cannot participate in such a chain of intergenerational inheritance. How then can the biogenesis of Y-linked piRNAs be initiated? Here, using Suppressor of Stellate (Su(Ste)), a Y-linked Drosophila melanogaster piRNA locus as a model, we show that Su(Ste) piRNAs are made in the early male germline via 5'-to-3' phased piRNA biogenesis initiated by maternally deposited 1360/Hoppel transposon piRNAs. Notably, deposition of Su(Ste) piRNAs from XXY mothers obviates the need for phased piRNA biogenesis in sons. Together, our study uncovers a developmentally programmed, intergenerational mechanism that allows fly mothers to protect their sons using a Y-linked piRNA locus.

A machine learning approach for the prediction of protein surface loop flexibility.

Proteins often undergo conformational changes when binding to each other. A major fraction of backbone conformational changes involves motion on the protein surface, particularly in loops. Accounting for the motion of protein surface loops represents a challenge for protein-protein docking algorithms. A first step in addressing this challenge is to distinguish protein surface loops that are likely to undergo backbone conformational changes upon protein-protein binding (mobile loops) from those that are not (stationary loops). In this study, we developed a machine learning strategy based on support vector machines (SVMs). Our SVM uses three features of loop residues in the unbound protein structures-Ramachandran angles, crystallographic B-factors, and relative accessible surface area-to distinguish mobile loops from stationary ones. This method yields an average prediction accuracy of 75.3% compared with a random prediction accuracy of 50%, and an average of 0.79 area under the receiver operating characteristic (ROC) curve using cross-validation. Testing the method on an independent dataset, we obtained a prediction accuracy of 70.5%. Finally, we applied the method to 11 complexes that involve members from the Ras superfamily and achieved prediction accuracy of 92.8% for the Ras superfamily proteins and 74.4% for their binding partners.

A combined experimental and theoretical study of ion solvation in liquid N-methylacetamide.

Most current biomolecular simulations are based on potential energy functions that treat the electrostatic energy as a sum of pairwise Coulombic interactions between effective fixed atomic charges. This approximation, in which many-body induced polarization effects are included in an average way, is expected to be satisfactory for a wide range of systems, but less accurate for processes involving the transfer and partition of ions among heterogeneous environments. The limitations of these potential energy functions are perhaps most obvious in studies of ion permeation through membrane channels. In many cases, the pore is so narrow that the permeating ion must shed most of its surrounding water molecules and the large energetic loss due to dehydration must be compensated by coordination with protein atoms. Interactions of cations with protein backbone carbonyl oxygens, in particular, play a critical role in several important biological channels. As a first step toward meeting the challenge of developing an accurate explicit accounting for induced polarization effects, the present work combines experiments and computation to characterize the interactions of alkali and halide ions with N-methylacetamide chosen to represent the peptide bond. From solubility measurements, we extract the solvation free energies of KCl and NaCl in liquid N-methylacetamide. Polarizable models based on the Drude oscillator are then developed and compared with available experimental and ab initio data. The good agreement for a range of structural and thermodynamic properties in the gas and condensed phases suggests that the polarizable models provide an accurate representation of ion-amide interactions in biological systems.

Characterizing Protein-Ligand Binding Using Atomistic Simulation and Machine Learning: Application to Drug Resistance in HIV-1 Protease.

Over the past several decades, atomistic simulations of biomolecules, whether carried out using molecular dynamics or Monte Carlo techniques, have provided detailed insights into their function. Comparing the results of such simulations for a few closely related systems has guided our understanding of the mechanisms by which changes such as ligand binding or mutation can alter the function. The general problem of detecting and interpreting such mechanisms from simulations of many related systems, however, remains a challenge. This problem is addressed here by applying supervised and unsupervised machine learning techniques to a variety of thermodynamic observables extracted from molecular dynamics simulations of different systems. As an important test case, these methods are applied to understand the evasion by human immunodeficiency virus type-1 (HIV-1) protease of darunavir, a potent inhibitor to which resistance can develop via the simultaneous mutation of multiple amino acids. Complex mutational patterns have been observed among resistant strains, presenting a challenge to developing a mechanistic picture of resistance in the protease. In order to dissect these patterns and gain mechanistic insight into the role of specific mutations, molecular dynamics simulations were carried out on a collection of HIV-1 protease variants, chosen to include highly resistant strains and susceptible controls, in complex with darunavir. Using a machine learning approach that takes advantage of the hierarchical nature in the relationships among the sequence, structure, and function, an integrative analysis of these trajectories reveals key details of the resistance mechanism, including changes in the protein structure, hydrogen bonding, and protein-ligand contacts.

Synonymous Mutations at the Beginning of the Influenza A Virus Hemagglutinin Gene Impact Experimental Fitness.

The fitness effects of synonymous mutations can provide insights into biological and evolutionary mechanisms. We analyzed the experimental fitness effects of all single-nucleotide mutations, including synonymous substitutions, at the beginning of the influenza A virus hemagglutinin (HA) gene. Many synonymous substitutions were deleterious both in bulk competition and for individually isolated clones. Investigating protein and RNA levels of a subset of individually expressed HA variants revealed that multiple biochemical properties contribute to the observed experimental fitness effects. Our results indicate that a structural element in the HA segment viral RNA may influence fitness. Examination of naturally evolved sequences in human hosts indicates a preference for the unfolded state of this structural element compared to that found in swine hosts. Our overall results reveal that synonymous mutations may have greater fitness consequences than indicated by simple models of sequence conservation, and we discuss the implications of this finding for commonly used evolutionary tests and analyses.

Elucidating the Interdependence of Drug Resistance from Combinations of Mutations.

HIV-1 protease is responsible for the cleavage of 12 nonhomologous sites within the Gag and Gag-Pro-Pol polyproteins in the viral genome. Under the selective pressure of protease inhibition, the virus evolves mutations within (primary) and outside of (secondary) the active site, allowing the protease to process substrates while simultaneously countering inhibition. The primary protease mutations impede inhibitor binding directly, while the secondary mutations are considered accessory mutations that compensate for a loss in fitness. However, the role of secondary mutations in conferring drug resistance remains a largely unresolved topic. We have shown previously that mutations distal to the active site are able to perturb binding of darunavir (DRV) via the protein's internal hydrogen-bonding network. In this study, we show that mutations distal to the active site, regardless of context, can play an interdependent role in drug resistance. Applying eigenvalue decomposition to collections of hydrogen bonding and van der Waals interactions from a series of molecular dynamics simulations of 15 diverse HIV-1 protease variants, we identify sites in the protease where amino acid substitutions lead to perturbations in nonbonded interactions with DRV and/or the hydrogen-bonding network of the protease itself. While primary mutations are known to drive resistance in HIV-1 protease, these findings delineate the significant contributions of accessory mutations to resistance. Identifying the variable positions in the protease that have the greatest impact on drug resistance may aid in future structure-based design of inhibitors.

Chromatin dynamics regulate mesenchymal stem cell lineage specification and differentiation to osteogenesis.

Multipotent mesenchymal stromal cells (MSCs) are critical for regeneration of multiple tissues. Epigenetic mechanisms are fundamental regulators of lineage specification and cell fate, and as such, we addressed the question of which epigenetic modifications characterize the transition of nascent MSCs to a tissue specific MSC-derived phenotype. By profiling the temporal changes of seven histone marks correlated to gene expression during proliferation, early commitment, matrix deposition, and mineralization stages, we identified distinct epigenetic mechanisms that regulate transcriptional programs necessary for tissue-specific phenotype development. Patterns of stage-specific enrichment of histone modifications revealed distinct modes of repression and activation of gene expression that would not be detected using single endpoint analysis. We discovered that at commitment, H3K27me3 is removed from genes that are upregulated and is not acquired on downregulated genes. Additionally, we found that the absence of H3K4me3 modification at promoters defined a subset of osteoblast-specific upregulated genes, indicating that acquisition of acetyl modifications drive activation of these genes. Significantly, loss or gain of H3K36me3 was the primary predictor of dynamic changes in temporal gene expression. Using unsupervised pattern discovery analysis the signature of osteogenic-related histone modifications identified novel functional cis regulatory modules associated with enhancer regions that control tissue-specific genes. Our work provides a cornerstone to understand the epigenetic regulation of transcriptional programs that are important for MSC lineage commitment and lineage, as well as insights to facilitate MSC-based therapeutic interventions.

Demethylated HSATII DNA and HSATII RNA Foci Sequester PRC1 and MeCP2 into Cancer-Specific Nuclear Bodies.

This study reveals that high-copy satellite II (HSATII) sequences in the human genome can bind and impact distribution of chromatin regulatory proteins and that this goes awry in cancer. In many cancers, master regulatory proteins form two types of cancer-specific nuclear bodies, caused by locus-specific deregulation of HSATII. DNA demethylation at the 1q12 mega-satellite, common in cancer, causes PRC1 aggregation into prominent Cancer-Associated Polycomb (CAP) bodies. These loci remain silent, whereas HSATII loci with reduced PRC1 become derepressed, reflecting imbalanced distribution of UbH2A on these and other PcG-regulated loci. Large nuclear foci of HSATII RNA form and sequester copious MeCP2 into Cancer-Associated Satellite Transcript (CAST) bodies. Hence, HSATII DNA and RNA have an exceptional capacity to act as molecular sponges and sequester chromatin regulatory proteins into abnormal nuclear bodies in cancer. The compartmentalization of regulatory proteins within nuclear structure, triggered by demethylation of "junk" repeats, raises the possibility that this contributes to further compromise of the epigenome and neoplastic progression.

Structural Basis of the Disorder in the Tandem Zinc Finger Domain of the RNA-Binding Protein Tristetraprolin.

Tristetraprolin (TTP) and TIS11d are two human RNA-binding proteins that belong to the CCCH-type tandem zinc finger family. In the RNA-free state, TIS11d coordinates a zinc ion in each of its two fingers, while TTP coordinates a single zinc ion with the N-terminal zinc finger. We have previously identified three residues, located in the C-terminal half of a short α-helix in the second zinc finger, that control how structured the RNA-binding domain is in these two proteins: Y151, L152, and Q153 in TTP and H201, T202, and I203 in TIS11d. Here, we have used molecular dynamics, NMR spectroscopy, and other biochemical methods to investigate the role of these three residues in the stability of the RNA-binding domain. We found that the intrahelical hydrogen bond formed by the T202 hydroxyl group in the C-terminal zinc finger of TIS11d is necessary to allow for π-π stacking between the side chains of a conserved phenylalanine and the zinc-coordinating histidine. We demonstrated that the lack of this hydrogen bond in TTP is responsible for the reduced zinc affinity of the C-terminal zinc finger.