Time-resolved structural analysis of an RNA-cleaving DNA catalyst.

The 10-23 DNAzyme is one of the most prominent catalytically active DNA sequences1,2. Its ability to cleave a wide range of RNA targets with high selectivity entails a substantial therapeutic and biotechnological potential2. However, the high expectations have not yet been met, a fact that coincides with the lack of high-resolution and time-resolved information about its mode of action3. Here we provide high-resolution NMR characterization of all apparent states of the prototypic 10-23 DNAzyme and present a comprehensive survey of the kinetics and dynamics of its catalytic function. The determined structure and identified metal-ion-binding sites of the precatalytic DNAzyme-RNA complex reveal that the basis of the DNA-mediated catalysis is an interplay among three factors: an unexpected, yet exciting molecular architecture; distinct conformational plasticity; and dynamic modulation by metal ions. We further identify previously hidden rate-limiting transient intermediate states in the DNA-mediated catalytic process via real-time NMR measurements. Using a rationally selected single-atom replacement, we could considerably enhance the performance of the DNAzyme, demonstrating that the acquired knowledge of the molecular structure, its plasticity and the occurrence of long-lived intermediate states constitutes a valuable starting point for the rational design of next-generation DNAzymes.

The cyclophilin A-binding loop of the capsid regulates the human TRIM5α sensitivity of nonpandemic HIV-1.

The rather few cases of humans infected by HIV-1 N, O, or P raise the question of their incomplete adaptation to humans. We hypothesized that early postentry restrictions may be relevant for the impaired spread of these HIVs. One of the best-characterized species-specific restriction factors is TRIM5α. HIV-1 M can escape human (hu) TRIM5α restriction by binding cyclophilin A (CYPA, also known as PPIA, peptidylprolyl isomerase A) to the so-called CYPA-binding loop of its capsid protein. How non-M HIV-1s interact with huTRIM5α is ill-defined. By testing full-length reporter viruses (Δ env) of HIV-1 N, O, P, and SIVgor (simian IV of gorillas), we found that in contrast to HIV-1 M, the nonpandemic HIVs and SIVgor showed restriction by huTRIM5α. Work to identify capsid residues that mediate susceptibility to huTRIM5α revealed that residue 88 in the capsid CYPA-binding loop was important for such differences. There, HIV-1 M uses alanine to resist, while non-M HIV-1s have either valine or methionine, which avail them for huTRIM5α. Capsid residue 88 determines the sensitivity to TRIM5α in an unknown way. Molecular simulations indicated that capsid residue 88 can affect trans-to-cis isomerization patterns on the capsids of the viruses we tested. These differential CYPA usages by pandemic and nonpandemic HIV-1 suggest that the enzymatic activity of CYPA on the viral core might be important for its protective function against huTRIM5α.

Respiratory and C4-photosynthetic NAD-malic enzyme coexist in bundle sheath cell mitochondria and evolved via association of differentially adapted subunits.

In plant mitochondria, nicotinamide adenine dinucleotide-malic enzyme (NAD-ME) has a housekeeping function in malate respiration. In different plant lineages, NAD-ME was independently co-opted in C4 photosynthesis. In the C4 Cleome species, Gynandropsis gynandra and Cleome angustifolia, all NAD-ME genes (NAD-MEα, NAD-MEß1, and NAD-MEß2) were affected by C4 evolution and are expressed at higher levels than their orthologs in the C3 species Tarenaya hassleriana. In T. hassleriana, the NAD-ME housekeeping function is performed by two heteromers, NAD-MEα/ß1 and NAD-MEα/ß2, with similar biochemical properties. In both C4 species, this role is restricted to NAD-MEα/ß2. In the C4 species, NAD-MEα/ß1 is exclusively present in the leaves, where it accounts for most of the enzymatic activity. Gynandropsis gynandra NAD-MEα/ß1 (GgNAD-MEα/ß1) exhibits high catalytic efficiency and is differentially activated by the C4 intermediate aspartate, confirming its role as the C4-decarboxylase. During C4 evolution, NAD-MEß1 lost its catalytic activity; its contribution to the enzymatic activity results from a stabilizing effect on the associated α-subunit and the acquisition of regulatory properties. We conclude that in bundle sheath cell mitochondria of C4 species, the functions of NAD-ME as C4 photosynthetic decarboxylase and as a housekeeping enzyme coexist and are performed by isoforms that combine the same α-subunit with differentially adapted ß-subunits.

A MademoiseLLE domain binding platform links the key RNA transporter to endosomes.

Spatiotemporal expression can be achieved by transport and translation of mRNAs at defined subcellular sites. An emerging mechanism mediating mRNA trafficking is microtubule-dependent co-transport on shuttling endosomes. Although progress has been made in identifying various components of the endosomal mRNA transport machinery, a mechanistic understanding of how these RNA-binding proteins are connected to endosomes is still lacking. Here, we demonstrate that a flexible MademoiseLLE (MLLE) domain platform within RNA-binding protein Rrm4 of Ustilago maydis is crucial for endosomal attachment. Our structure/function analysis uncovered three MLLE domains at the C-terminus of Rrm4 with a functionally defined hierarchy. MLLE3 recognises two PAM2-like sequences of the adaptor protein Upa1 and is essential for endosomal shuttling of Rrm4. MLLE1 and MLLE2 are most likely accessory domains exhibiting a variable binding mode for interaction with currently unknown partners. Thus, endosomal attachment of the mRNA transporter is orchestrated by a sophisticated MLLE domain binding platform.

TopEnzyme: a framework and database for structural coverage of the functional enzyme space.

MOTIVATION: TopEnzyme is a database of structural enzyme models created with TopModel and is linked to the SWISS-MODEL repository and AlphaFold Protein Structure Database to provide an overview of structural coverage of the functional enzyme space for over 200 000 enzyme models. It allows the user to quickly obtain representative structural models for 60% of all known enzyme functions. RESULTS: We assessed the models with TopScore and contributed 9039 good-quality and 1297 high-quality structures. Furthermore, we compared these models to AlphaFold2 models with TopScore and found that the TopScore differs only by 0.04 on average in favor of AlphaFold2. We tested TopModel and AlphaFold2 for targets not seen in the respective training databases and found that both methods create qualitatively similar structures. When no experimental structures are available, this database will facilitate quick access to structural models across the currently most extensive structural coverage of the functional enzyme space within Swiss-Prot. AVAILABILITY AND IMPLEMENTATION: We provide a full web interface to the database at https://cpclab.uni-duesseldorf.de/topenzyme/.

Molecular Mechanisms Underlying Single Nucleotide Polymorphism-Induced Reactivity Decrease in CYP2D6.

Cytochrome P450 2D6 (CYP2D6) is one of the most important enzymes involved in drug metabolism. Genetic polymorphism can influence drug metabolism by CYP2D6 such that a therapy is seriously affected by under- or overdosing of drugs. However, a general explanation at the atomistic level for poor activity is missing so far. Here we show for the 20 most common single nucleotide polymorphisms (SNPs) of CYP2D6 that poor metabolism is driven by four mechanisms. We found in extensive all-atom molecular dynamics simulations that the rigidity of the I-helix (central helix), distance between central phenylalanines (stabilizing bound substrate), availability of basic residues on the surface of CYP2D6 (binding of cytochrome P450 reductase), and position of arginine 132 (electron transfer to heme) are essential for an extensive function of the enzyme. These results were applied to SNPs with unknown effects, and potential SNPs that may lead to poor drug metabolism were identified. The revealed molecular mechanisms might be important for other drug-metabolizing cytochrome P450 enzymes.

Thermodynamic profile of mutual subunit control in a heteromeric receptor.

Cyclic nucleotide-gated (CNG) ion channels of olfactory neurons are tetrameric membrane receptors that are composed of two A2 subunits, one A4 subunit, and one B1b subunit. Each subunit carries a cyclic nucleotide-binding domain in the carboxyl terminus, and the channels are activated by the binding of cyclic nucleotides. The mechanism of cooperative channel activation is still elusive. Using a complete set of engineered concatenated olfactory CNG channels, with all combinations of disabled binding sites and fit analyses with systems of allosteric models, the thermodynamics of microscopic cooperativity for ligand binding was subunit- and state-specifically quantified. We show, for the closed channel, that preoccupation of each of the single subunits increases the affinity of each other subunit with a Gibbs free energy (ΔΔG) of ∼-3.5 to ∼-5.5 kJ ⋅ mol-1, depending on the subunit type, with the only exception that a preoccupied opposite A2 subunit has no effect on the other A2 subunit. Preoccupation of two neighbor subunits of a given subunit causes the maximum affinity increase with ΔΔG of ∼-9.6 to ∼-9.9 kJ ⋅ mol-1 Surprisingly, triple preoccupation leads to fewer negative ΔΔG values for a given subunit as compared to double preoccupation. Channel opening increases the affinity of all subunits. The equilibrium constants of closed-open isomerizations systematically increase with progressive liganding. This work demonstrates, on the example of the heterotetrameric olfactory CNG channel, a strategy to derive detailed insights into the specific mutual control of the individual subunits in a multisubunit membrane receptor.

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin.

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

Met/Val129 polymorphism of the full-length human prion protein dictates distinct pathways of amyloid formation.

Methionine/valine polymorphism at position 129 of the human prion protein, huPrP, is tightly associated with the pathogenic phenotype, disease progress, and age of onset of neurodegenerative diseases such as Creutzfeldt-Jakob disease or Fatal Familial Insomnia. This raises the question of whether and how the amino acid type at position 129 influences the structural properties of huPrP, affecting its folding, stability, and amyloid formation behavior. Here, our detailed biophysical characterization of the 129M and 129V variants of recombinant full-length huPrP(23-230) by amyloid formation kinetics, CD spectroscopy, molecular dynamics simulations, and sedimentation velocity analysis reveals differences in their aggregation propensity and oligomer content, leading to deviating pathways for the conversion into amyloid at acidic pH. We determined that the 129M variant exhibits less secondary structure content before amyloid formation and higher resistance to thermal denaturation compared to the 129V variant, whereas the amyloid conformation of both variants shows similar thermal stability. Additionally, our molecular dynamics simulations and rigidity analyses at the atomistic level identify intramolecular interactions responsible for the enhanced monomer stability of the 129M variant, involving more frequent minimum distances between E196 and R156, forming a salt bridge. Removal of the N-terminal half of the 129M full-length variant diminishes its differences compared to the 129V full-length variant and highlights the relevance of the flexible N terminus in huPrP. Taken together, our findings provide insight into structural properties of huPrP and the effects of the amino acid identity at position 129 on amyloid formation behavior.

Atomic Resolution Insights into pH Shift Induced Deprotonation Events in LS-Shaped Aß(1-42) Amyloid Fibrils.

Alzheimer's disease is a neurodegenerative disorder associated with the deposition of misfolded aggregates of the amyloid-ß protein (Aß). Aß(1-42) is one of the most aggregation-prone components in senile plaques of AD patients. We demonstrated that relatively homogeneous Aß(1-42) fibrils with one predominant fold visible in solid-state NMR spectra can be obtained at acidic pH. The structure of these fibrils differs remarkably from some other polymorphs obtained at neutral pH. In particular, the entire N-terminal region is part of the rigid fibril core. Here, we investigate the effects of a pH shift on the stability and the fold of these fibrils at higher pH values. Fibril bundling at neutral pH values renders cryo-EM studies impractical, but solid-state NMR spectroscopy, molecular dynamics simulations, and biophysical methods provide residue-specific structural information under these conditions. The LS-fold of the Aß(1-42) fibrils does not change over the complete pH range from pH 2 to pH 7; in particular, the N-terminus remains part of the fibril core. We observe changes in the protonation state of charged residues starting from pH 5 on a residue-specific level. The deprotonation of the C-terminal carboxyl group of A42 in the intermolecular salt bridge with D1 and K28 is slow on the NMR time scale, with a local pKa of 5.4, and local conformations of the involved residues are affected by deprotonation of A42. Thus, we demonstrate that this fibril form is stable at physiological pH values.

Expanding the clinical and immunological phenotypes of PAX1-deficient SCID and CID patients.

Paired box 1 (PAX1) deficiency has been reported in a small number of patients diagnosed with otofaciocervical syndrome type 2 (OFCS2). We described six new patients who demonstrated variable clinical penetrance. Reduced transcriptional activity of pathogenic variants confirmed partial or complete PAX1 deficiency. Thymic aplasia and hypoplasia were associated with impaired T cell immunity. Corrective treatment was required in 4/6 patients. Hematopoietic stem cell transplantation resulted in poor immune reconstitution with absent naïve T cells, contrasting with the superior recovery of T cell immunity after thymus transplantation. Normal ex vivo differentiation of PAX1-deficient CD34+ cells into mature T cells demonstrated the absence of a hematopoietic cell-intrinsic defect. New overlapping features with DiGeorge syndrome included primary hypoparathyroidism (n = 5) and congenital heart defects (n = 2), in line with PAX1 expression during early embryogenesis. Our results highlight new features of PAX1 deficiency, which are relevant to improving early diagnosis and identifying patients requiring corrective treatment.

Functional characterization of novel or yet uncharacterized ATP7B missense variants detected in patients with clinical Wilson's disease.

Wilson's disease (WD, MIM#277900) is an autosomal recessive disorder resulting in copper excess caused by biallelic variants in the ATP7B gene (MIM#606882) encoding a copper transporting P-type ATPase. ATP7B variants of unknown significance (VUS) are detected frequently, sometimes impeding a clear diagnosis. Functional analyses can help to classify these variants as benign or pathogenic. Additionally, variants already classified as (likely) pathogenic benefit from functional analyses to understand their pathomechanism, thus contribute to the development of personalized treatment approaches in the future. We described clinical features of six WD patients and functionally characterized five ATP7B missense variants (two VUS, three yet uncharacterized likely pathogenic variants), detected in these patients. We determined the protein level, copper export capacity, and cellular localization in an in vitro model and potential structural consequences using an ATP7B protein model based on AlphaFold. Our analyses give insight into the pathomechanism and allowed reclassification for the two VUS to likely pathogenic and for two of the three likely pathogenic variants to pathogenic.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions.

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

LISTER: Semiautomatic Metadata Extraction from Annotated Experiment Documentation in eLabFTW.

The availability of scientific methods, code, and data is key for reproducing an experiment. Research data should be made available following the FAIR principle (findable, accessible, interoperable, and reusable). For that, the annotation of research data with metadata is central. However, existing research data management workflows often require that metadata be created by the corresponding researchers, which takes effort and time. Here, we developed LISTER as a methodological and algorithmic solution to create and extract metadata from annotated, template-based experimental documentation using minimum effort. We focused on tailoring the integration between existing platforms by using eLabFTW as the electronic lab notebook and adopting the ISA (investigation, study, assay) model as the abstract data model framework. LISTER consists of four components: annotation language to support metadata extraction; customized eLabFTW entries using specific hierarchies, templates, and tags to structure reusable scientific documentation; a "container" concept in eLabFTW, making metadata of a particular container content extractable along with its underlying, related experiments via a single click; a Python-based app to enable easy-to-use, semiautomated metadata extraction from eLabFTW entries. LISTER outputs metadata in machine-readable .json and human-readable .xlsx formats, and Material and Methods (MM) descriptions in .docx format that could be used in a thesis or manuscript. The metadata can be used as a basis to create or extend ontologies, which, when applied to the published research data, will significantly enhance its value. DSpace is used as a data cataloging platform for hosting the extracted metadata and research data. We applied LISTER to computational biophysical chemistry, protein biochemistry, and molecular biology, and our concept should be extendable to other life science areas.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Spectroscopic and QM/MM studies of the Cu(I) binding site of the plant ethylene receptor ETR1.

Herein, we present, to our knowledge, the first spectroscopic characterization of the Cu(I) active site of the plant ethylene receptor ETR1. The x-ray absorption (XAS) and extended x-ray absorption fine structure (EXAFS) spectroscopies presented here establish that ETR1 has a low-coordinate Cu(I) site. The EXAFS resolves a mixed first coordination sphere of N/O and S scatterers at distances consistent with potential histidine and cysteine residues. This finding agrees with the coordination of residues C65 and H69 to the Cu(I) site, which are critical for ethylene activity and well conserved. Furthermore, the Cu K-edge XAS and EXAFS of ETR1 exhibit spectroscopic changes upon addition of ethylene that are attributed to modifications in the Cu(I) coordination environment, suggestive of ethylene binding. Results from umbrella sampling simulations of the proposed ethylene binding helix of ETR1 at a mixed quantum mechanics/molecular mechanics level agree with the EXAFS fit distance changes upon ethylene binding, particularly in the increase of the distance between H69 and Cu(I), and yield binding energetics comparable with experimental dissociation constants. The observed changes in the copper coordination environment might be the triggering signal for the transmission of the ethylene response.

Posttranslational Modification of the NADP-Malic Enzyme Involved in C4 Photosynthesis Modulates the Enzymatic Activity during the Day.

Evolution of the C4 photosynthetic pathway involved in some cases recruitment of housekeeping proteins through gene duplication and their further neofunctionalization. NADP-malic enzyme (ME), the most widespread C4 decarboxylase, has increased its catalytic efficiency and acquired regulatory properties that allowed it to participate in the C4 pathway. Here, we show that regulation of maize (Zea mays) C4-NADP-ME activity is much more elaborate than previously thought. Using mass spectrometry, we identified phosphorylation of the Ser419 residue of C4-NADP-ME in protein extracts of maize leaves. The phosphorylation event increases in the light, with a peak at Zeitgeber time 2. Phosphorylation of ZmC4-NADP-ME drastically decreases its activity as shown by the low residual activity of the recombinant phosphomimetic mutant. Analysis of the crystal structure of C4-NADP-ME indicated that Ser419 is involved in the binding of NADP at the active site. Molecular dynamics simulations and effective binding energy computations indicate a less favorable binding of the cofactor NADP in the phosphomimetic and the phosphorylated variants. We propose that phosphorylation of ZmC4-NADP-ME at Ser419 during the first hours in the light is a cellular mechanism that fine tunes the enzymatic activity to coordinate the carbon concentration mechanism with the CO2 fixation rate, probably to avoid CO2 leakiness from bundle sheath cells.

Protein engineering for feedback resistance in 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase.

The shikimate pathway delivers aromatic amino acids (AAAs) in prokaryotes, fungi, and plants and is highly utilized in the industrial synthesis of bioactive compounds. Carbon flow into this pathway is controlled by the initial enzyme 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase (DAHPS). AAAs produced further downstream, phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp), regulate DAHPS by feedback inhibition. Corynebacterium glutamicum, the industrial workhorse for amino acid production, has two isoenzymes of DAHPS, AroF (Tyr sensitive) and AroG (Phe and Tyr sensitive). Here, we introduce feedback resistance against Tyr in the class I DAHPS AroF (AroFcg). We pursued a consensus approach by drawing on structural modeling, sequence and structural comparisons, knowledge of feedback-resistant variants in E. coli homologs, and computed folding free energy changes. Two types of variants were predicted: Those where substitutions putatively either destabilize the inhibitor binding site or directly interfere with inhibitor binding. The recombinant variants were purified and assessed in enzyme activity assays in the presence or absence of Tyr. Of eight AroFcg variants, two yielded > 80% (E154N) and > 50% (P155L) residual activity at 5 mM Tyr and showed > 50% specific activity of the wt AroFcg in the absence of Tyr. Evaluation of two and four further variants at positions 154 and 155 yielded E154S, completely resistant to 5 mM Tyr, and P155I, which behaves similarly to P155L. Hence, feedback-resistant variants were found that are unlikely to evolve by point mutations from the parental gene and, thus, would be missed by classical strain engineering. KEY POINTS: • We introduce feedback resistance against Tyr in the class I DAHPS AroF • Variants at position 154 (155) yield > 80% (> 50%) residual activity at 5 mM Tyr • The variants found are unlikely to evolve by point mutations from the parental gene.

Recurrent Germline Variant in RAD21 Predisposes Children to Lymphoblastic Leukemia or Lymphoma.

Somatic loss of function mutations in cohesin genes are frequently associated with various cancer types, while cohesin disruption in the germline causes cohesinopathies such as Cornelia-de-Lange syndrome (CdLS). Here, we present the discovery of a recurrent heterozygous RAD21 germline aberration at amino acid position 298 (p.P298S/A) identified in three children with lymphoblastic leukemia or lymphoma in a total dataset of 482 pediatric cancer patients. While RAD21 p.P298S/A did not disrupt the formation of the cohesin complex, it altered RAD21 gene expression, DNA damage response and primary patient fibroblasts showed increased G2/M arrest after irradiation and Mitomycin-C treatment. Subsequent single-cell RNA-sequencing analysis of healthy human bone marrow confirmed the upregulation of distinct cohesin gene patterns during hematopoiesis, highlighting the importance of RAD21 expression within proliferating B- and T-cells. Our clinical and functional data therefore suggest that RAD21 germline variants can predispose to childhood lymphoblastic leukemia or lymphoma without displaying a CdLS phenotype.

Allosteric signaling in C-linker and cyclic nucleotide-binding domain of HCN2 channels.

Opening of hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels is controlled by membrane hyperpolarization and binding of cyclic nucleotides to the tetrameric cyclic nucleotide-binding domain (CNBD), attached to the C-linker (CL) disk. Confocal patch-clamp fluorometry revealed pronounced cooperativity of ligand binding among protomers. However, by which pathways allosteric signal transmission occurs remained elusive. Here, we investigate how changes in the structural dynamics of the CL-CNBD of mouse HCN2 upon cAMP binding relate to inter- and intrasubunit signal transmission. Applying a rigidity-theory-based approach, we identify two intersubunit and one intrasubunit pathways that differ in allosteric coupling strength between cAMP-binding sites or toward the CL. These predictions agree with results from electrophysiological and patch-clamp fluorometry experiments. Our results map out distinct routes within the CL-CNBD that modulate different cAMP-binding responses in HCN2 channels. They signify that functionally relevant submodules may exist within and across structurally discernable subunits in HCN channels.