Human telomerase model shows the role of the TEN domain in advancing the double helix for the next polymerization step

Mechanistic influence of discreet conformation of human telomerase linker region

Approximately 90% of malignancies have been shown to have human telomerase activity, establishing it as a viable therapeutic target. The crystal structure of telomerase was determined recently. However, the tertiary structure of the non-conserved flexible linker region remains unresolved. This study aims to predict the full-length tertiary structure of the human telomerase reverse transcriptase (hTERT). Two strategies were employed to determine the full-length structure of hTERT (1132 amino acids); iterative threading and a conjoined model generated from machine learning and energy functions. After energy minimization, Ramachandran Plot analysis, and simulation; the conjoined model was considered of better quality and stability. The linker region of the conjoined depicted two helices from approximately 275-284 and 201-211 amino acids respectively in contrast to the iterative threading model which has a single helix. Moreover, the region was observed to undergo major structural changes throughout the simulation. These changes signify its flexibility which might be due to the region having a significant number of glycine and proline and could enhance the clamping movement.Communicated by Ramaswamy H. Sarma.

Conformational plasticity and allosteric communication networks explain Shelterin protein TPP1 binding to human telomerase

The Shelterin complex protein TPP1 interacts with human telomerase (TERT) by means of the TEL-patch region, controlling telomere homeostasis. Aberrations in the TPP1-TERT heterodimer formation might lead to short telomeres and severe diseases like dyskeratosis congenita and Hoyeraal-Hreidarsson syndrome. In the present study, we provide a thorough characterization of the structural properties of the TPP1's OB-domain by combining data coming from microsecond-long molecular dynamics calculations, time-series analyses, and graph-based networks. Our results show that the TEL-patch conformational freedom is influenced by a network of long-range amino acid communications that together determine the proper TPP1-TERT binding. Furthermore, we reveal that in TPP1 pathological variants Glu169Δ, Lys170Δ and Leu95Gln, the TEL-patch plasticity is reduced, affecting the correct binding to TERT and, in turn, telomere processivity, which eventually leads to accelerated aging of affected cells. Our study provides a structural basis for the design of TPP1-targeting ligands with therapeutic potential against cancer and telomeropathies.

Cycloastragenol activation of telomerase improves β-Klotho protein level and attenuates age-related malfunctioning in ovarian tissues

Age-related deterioration in the reproductive capacity of women is directly related to the poor developmental potential of ovarian follicles. Although telomerase plays a key role in female fertility, TERT-targeting therapeutic strategies for age-related female infertility have yet to be investigated. This study elucidated the effect of Telomerase activation on mice ovaries and more specifically on Klb (β-Klotho) gene expression, which is linked to ageing, female hormonal regulation, and cyclicity. The homology-based 3D model of hTERT was used to predict its binding mode of Cycloastragenol (CAG) using molecular docking and molecular dynamics simulations. Based on docking score, simulation behavior, and interaction with hTERT residues it was observed that CAG could bind with the hTERT model. CAG treatment to primary cultured mouse granulosa cells and activation of telomerase was examined via telomerase activity assay (Mouse TE (telomerase) ELISA Kit) and telomere length by quantitative fluorescence in situ hybridization. CAG mediated telomerase also significantly improved β-Klotho protein level in the aged granulosa cells. To demonstrate that β-Klotho is telomerase dependent, the TERT was knocked down via siRNA in granulosa cells and protein level of β-Klotho was examined. Furthermore, CAG-mediated telomerase activation significantly enhanced the level of Klb and recovered ovarian follicles in the D-galactose (D-gal)-induced ovarian ageing mouse model. Moreover, Doxorubicin-induced ovarian damage, which changes ovarian hormones, and inhibit follicular growth was successfully neutralized by CAG activated telomerase and its recovery of β-Klotho level. In conclusion, TERT dependent β-Klotho regulation in ovarian tissues is one of the mechanisms, which can overcome female infertility.

Mechanisms of nucleotide selection by telomerase

Telomerase extends telomere sequences at chromosomal ends to protect genomic DNA. During this process it must select the correct nucleotide from a pool of nucleotides with various sugars and base pairing properties, which is critically important for the proper capping of telomeric sequences by shelterin. Unfortunately, how telomerase selects correct nucleotides is unknown. Here, we determined structures of Tribolium castaneum telomerase reverse transcriptase (TERT) throughout its catalytic cycle and mapped the active site residues responsible for nucleoside selection, metal coordination, triphosphate binding, and RNA template stabilization. We found that TERT inserts a mismatch or ribonucleotide ~1 in 10,000 and ~1 in 14,000 insertion events, respectively. At biological ribonucleotide concentrations, these rates translate to ~40 ribonucleotides inserted per 10 kilobases. Human telomerase assays determined a conserved tyrosine steric gate regulates ribonucleotide insertion into telomeres. Cumulatively, our work provides insight into how telomerase selects the proper nucleotide to maintain telomere integrity.

Targeting telomerase for its advent in cancer therapeutics

Telomerase has emerged as an important primary target in anticancer therapy. It is a distinctive reverse transcriptase enzyme, which extends the length of telomere at the 3' chromosomal end, and uses telomerase reverse transcriptase (TERT) and telomerase RNA template-containing domains. Telomerase has a vital role and is a contributing factor in human health, mainly affecting cell aging and cell proliferation. Due to its unique feature, it ensures unrestricted cell proliferation in malignancy and plays a major role in cancer disease. The development of telomerase inhibitors with increased specificity and better pharmacokinetics is being considered to design and develop newer potent anticancer agents. Use of natural and synthetic compounds for the inhibition of telomerase activity can lead to an opening of new vistas in cancer treatment. This review details about the telomerase biochemistry, use of natural and synthetic compounds; vaccines and oncolytic virus in therapy that suppress the telomerase activity. We have discussed structure-activity relationships of various natural and synthetic telomerase inhibitors to help medicinal chemists and chemical biology researchers with a ready reference and updated status of their clinical trials. Suppression of human TERT (hTERT) activity through inhibition of hTERT promoter is an important approach for telomerase inhibition.

In silico design of telomerase inhibitors

Telomerase is a reverse transcriptase enzyme involved in DNA synthesis at the end of linear chromosomes. Unlike in most other cells, telomerase is reactivated most cancerous cells and, therefore, has become a promising new anticancer target. Despite extensive research, direct telomerase inhibitors have yet not been introduced to the clinics because of the complexity of this enzyme. Structures of this protein from simple organisms and human homology models are currently available and have been used in structure-based drug design efforts to find potential inhibitors. Different is silico strategies have been applied and different chemical groups have been explored. Here, we provide an overview of recent discoveries.

Investigation on Quantitative Structure-Activity Relationships of 1,3,4-Oxadiazole Derivatives as Potential Telomerase Inhibitors

BACKGROUND

Telomerase, a reverse transcriptase, maintains telomere and chromosomes integrity of dividing cells, while it is inactivated in most somatic cells. In tumor cells, telomerase is highly activated, and works in order to maintain the length of telomeres causing immortality, hence it could be considered as a potential marker to tumorigenesis.A series of 1,3,4-oxadiazole derivatives showed significant broad-spectrum anticancer activity against different cell lines, and demonstrated telomerase inhibition.

METHODS

This series of 24 N-benzylidene-2-((5-(pyridine-4-yl)-1,3,4-oxadiazol-2yl)thio)acetohydrazide derivatives as telomerase inhibitors has been considered to carry out QSAR studies. The endpoint to build QSAR models is determined by the IC50 values for telomerase inhibition, i.e., the concentration (μM) of inhibitor that produces 50% inhibition. These values were converted to pIC50 (- log IC50) values. We used the most common and transparent method, where models are described by clearly expressed mathematical equations: Multiple Linear Regression (MLR) by Ordinary Least Squares (OLS).

RESULTS

Validated models with high correlation coefficients were developed. The Multiple Linear Regression (MLR) models, by Ordinary Least Squares (OLS), showed good robustness and predictive capability, according to the Multi-Criteria Decision Making (MCDM = 0.8352), a technique that simultaneously enhances the performances of a certain number of criteria. The descriptors selected for the models, such as electrotopological state (E-state) descriptors, and extended topochemical atom (ETA) descriptors, showed the relevant chemical information contributing to the activity of these compounds.

CONCLUSION

The results obtained in this study make sure about the identification of potential hits as prospective telomerase inhibitors.

Single-Molecule Analysis of Reverse Transcriptase Enzymes

The original discovery of enzymes that synthesize DNA using an RNA template appeared to contradict the central dogma of biology, in which information is transferred, in a unidirectional way, from DNA genes into RNA molecules. The paradigm-shifting discovery of RNA-dependent DNA polymerases, also called reverse transcriptases (RTs), reshaped existing views for how cells function; however, the scope of the impact RTs impose on biology had yet to be realized. In the decades of research since the early 1970s, the biomedical and biotechnological significance of retroviral RTs, as well as the evolutionarily related telomerase enzyme, has become exceedingly clear. One common theme that has emerged in the course of RT-related research is the central role of nucleic acid binding and dynamics during enzyme function. However, directly interrogating these dynamic properties is challenging because of the stochastic properties of biological macromolecules. In this review, we describe how the development of single-molecule biophysical techniques has opened new windows through which to observe the dynamic behavior of this remarkable class of enzymes. Specifically, we focus on how the powerful single-molecule Förster resonance energy transfer (FRET) method has been exploited to study the structure and function of the human immunodeficiency virus (HIV) RT and telomerase ribonucleoprotein (RNP) enzymes. These exciting studies have refined our understanding of RT catalysis, have revealed unforeseen structural rearrangements between RTs and their nucleic acid substrates, and have helped to characterize the mode of action of RT-inhibiting drugs. We conclude with a discussion of how the ongoing development of single-molecule technologies will continue to empower researchers to probe RT mechanisms in new and exciting ways.

Structure and function of the N-terminal domain of the yeast telomerase reverse transcriptase

The elongation of single-stranded DNA repeats at the 3′-ends of chromosomes by telomerase is a key process in maintaining genome integrity in eukaryotes. Abnormal activation of telomerase leads to uncontrolled cell division, whereas its down-regulation is attributed to ageing and several pathologies related to early cell death. Telomerase function is based on the dynamic interactions of its catalytic subunit (TERT) with nucleic acids—telomerase RNA, telomeric DNA and the DNA/RNA heteroduplex. Here, we present the crystallographic and NMR structures of the N-terminal (TEN) domain of TERT from the thermotolerant yeast Hansenula polymorpha and demonstrate the structural conservation of the core motif in evolutionarily divergent organisms. We identify the TEN residues that are involved in interactions with the telomerase RNA and in the recognition of the ‘fork’ at the distal end of the DNA product/RNA template heteroduplex. We propose that the TEN domain assists telomerase biological function and is involved in restricting the size of the heteroduplex during telomere repeat synthesis.

Structural, functional, and stability change predictions in human telomerase upon specific point mutations

Overexpression of telomerase is one of the hallmarks of human cancer. Telomerase is important for maintaining the integrity of the ends of chromosomes, which are called telomeres. A growing number of human disease syndromes are associated with organ failure caused by mutations in telomerase (hTERT or hTR). Mutations in telomerase lead to telomere shortening by decreasing the stability of the telomerase complex, reducing its accumulation, or directly affecting its enzymatic activity. In this work, potential human telomerase mutations were identified by a systematic computational approach. Moreover, molecular docking methods were used to predict the effects of these mutations on the affinity of certain ligands (C_9i, C_9k, 16A, and NSC749234). The C_9k inhibitor had the best binding affinity for wild-type (WT) telomerase. Moreover, C_9i and C_9k had improved interactions with human telomerase in most of the mutant models. The R631 and Y717 residues of WT telomerase formed interactions with all studied ligands and these interactions were also commonly found in most of the mutant models. Residues forming stable interactions with ligands in molecular dynamics (MD) were traced, and the MD simulations showed that the C_9k ligand formed different conformations with WT telomerase than the C_9i ligand.

Telomerase Reverse Transcriptase Contains a BH3-Like Motif and Interacts with BCL-2 Family Members

Upregulation of human telomerase reverse transcriptase (hTERT) expression is an important factor in the cellular survival and cancer. Although growing evidence suggests that hTERT inhibits cellular apoptosis by telomere-independent functions, the mechanisms involved are not fully understood. Here, we show that hTERT contains a BH3-like motif, a short peptide sequence found in BCL-2 family proteins, and interacts with anti-apoptotic BCL-2 family proteins MCL-1 and BCL-xL, suggesting a functional link between hTERT and the mitochondrial pathway of apoptosis. Additionally, we propose that hTERT can be categorized into the atypical BH3-only proteins that promote cellular survival, possibly due to the non-canonical interaction between hTERT and antiapoptotic proteins. Although the detailed mechanisms underlying the hTERT BH3-like motif functions and interactions between hTERT and BCL-2 family proteins have not been elucidated, this work proposes a possible connection between hTERT and BCL-2 family members and reconsiders the role of the BH3-like motif as an interaction motif.

Telomerase Inhibition by a New Synthetic Derivative of the Aporphine Alkaloid Boldine

Telomerase, the enzyme responsible for cell immortality, is an important target in anti-cancer drug discovery. Boldine, an abundant aporphine alkaloid of Peumus boldus, is known to inhibit telomerase at non-toxic concentrations. Cytotoxicity of N-benzylsecoboldine hydrochloride (BSB), a synthetic derivative of boldine, was determined using the MTT method in MCF7 and MDA-MB231 cells. Aliquots of cell lysates were incubated with various concentrations of BSB in qTRAP (quantitative telomere repeat amplification protocol)-ligand experiments before substrate elongation by telomerase or amplification by hot-start Taq polymerase. The crystal structure of TERT, the catalytic subunit of telomerase from Tribolium castaneum, was used for docking and molecular dynamics analysis. The qTRAP-ligand data gave an IC₅₀ value of about 0.17 ± 0.1 µM for BSB, roughly 400 times stronger than boldine, while the LD₅₀ in the cytotoxicity assays were 12.5 and 21.88 µM, respectively, in cells treated for 48 h. Although both compounds interacted well with the active site, MD analysis suggests a second binding site with which BSB interacts via two hydrogen bonds, much more strongly than boldine. Theoretical analyses also evaluated the IC₅₀ for BSB as submicromolar. BSB, with greater hydrophobicity and flexibility than boldine, represents a promising structure to inhibit telomerase at non-toxic concentrations.

Direct observation of nucleic acid binding dynamics by the telomerase essential N-terminal domain

Telomerase is a specialized enzyme that maintains telomere length by adding DNA repeats to chromosome ends. The catalytic protein subunit of telomerase utilizes the integral telomerase RNA to direct telomere DNA synthesis. The telomerase essential N-terminal (TEN) domain is required for enzyme function; however, the precise mechanism of the TEN domain during catalysis is not known. We report a single-molecule study of dynamic TEN-induced conformational changes in its nucleic acid substrates. The TEN domain from the yeast Candida parapsilosis (Cp) exhibits a strong binding preference for double-stranded nucleic acids, with particularly high affinity for an RNA-DNA hybrid mimicking the template-product complex. Surprisingly, the telomere DNA repeat sequence from C. parapsilosis forms a DNA hairpin that also binds CpTEN with high affinity. Mutations to several residues in a putative nucleic acid-binding patch of CpTEN significantly reduced its affinity to the RNA-DNA hybrid and telomere DNA hairpin. Substitution of comparable residues in the related Candida albicans TEN domain caused telomere maintenance defects in vivo and decreased primer extension activity in vitro. Collectively, our results support a working model in which dynamic interactions with telomere DNA and the template-product hybrid underlie the functional requirement for the TEN domain during the telomerase catalytic cycle.

A protocol for CABS-dock protein-peptide docking driven by side-chain contact information

Background

The characterization of protein–peptide interactions is a challenge for computational molecular docking. Protein–peptide docking tools face at least two major difficulties: (1) efficient sampling of large-scale conformational changes induced by binding and (2) selection of the best models from a large set of predicted structures. In this paper, we merge an efficient sampling technique with external information about side-chain contacts to sample and select the best possible models.

Methods

In this paper we test a new protocol that uses information about side-chain contacts in CABS-dock protein–peptide docking. As shown in our recent studies, CABS-dock enables efficient modeling of large-scale conformational changes without knowledge about the binding site. However, the resulting set of binding sites and poses is in many cases highly diverse and difficult to score.

Results

As we demonstrate here, information about a single side-chain contact can significantly improve the prediction accuracy. Importantly, the imposed constraints for side-chain contacts are quite soft. Therefore, the developed protocol does not require precise contact information and ensures large-scale peptide flexibility in the broad contact area.

Conclusions

The demonstrated protocol provides the extension of the CABS-dock method that can be practically used in the structure prediction of protein–peptide complexes guided by the knowledge of the binding interface.

Structural Analysis Reveals the Deleterious Effects of Telomerase Mutations in Bone Marrow Failure Syndromes

Naturally occurring mutations in the ribonucleoprotein reverse transcriptase, telomerase, are associated with the bone marrow failure syndromes dyskeratosis congenita, aplastic anemia, and idiopathic pulmonary fibrosis. However, the mechanism by which these mutations impact telomerase function remains unknown. Here we present the structure of the human telomerase C-terminal extension (or thumb domain) determined by the method of single-wavelength anomalous diffraction to 2.31 Å resolution. We also used direct telomerase activity and nucleic acid binding assays to explain how naturally occurring mutations within this portion of telomerase contribute to human disease. The single mutations localize within three highly conserved regions of the telomerase thumb domain referred to as motifs E-I (thumb loop and helix), E-II, and E-III (the FVYL pocket, comprising the hydrophobic residues Phe-1012, Val-1025, Tyr-1089, and Leu-1092). Biochemical data show that the mutations associated with dyskeratosis congenita, aplastic anemia, and idiopathic pulmonary fibrosis disrupt the binding between the protein subunit reverse transcriptase of the telomerase and its nucleic acid substrates leading to loss of telomerase activity and processivity. Collectively our data show that although these mutations do not alter the overall stability or expression of telomerase reverse transcriptase, these rare genetic disorders are associated with an impaired telomerase holoenzyme that is unable to correctly assemble with its nucleic acid substrates, leading to incomplete telomere extension and telomere attrition, which are hallmarks of these diseases.

Single-molecule FRET-Rosetta reveals RNA structural rearrangements during human telomerase catalysis

Maintenance of telomeres by telomerase permits continuous proliferation of rapidly dividing cells, including the majority of human cancers. Despite its direct biomedical significance, the architecture of the human telomerase complex remains unknown. Generating homogeneous telomerase samples has presented a significant barrier to developing improved structural models. Here we pair single-molecule Förster resonance energy transfer (smFRET) measurements with Rosetta modeling to map the conformations of the essential telomerase RNA core domain within the active ribonucleoprotein. FRET-guided modeling places the essential pseudoknot fold distal to the active site on a protein surface comprising the C-terminal element, a domain that shares structural homology with canonical polymerase thumb domains. An independently solved medium-resolution structure of Tetrahymena telomerase provides a blind test of our modeling methodology and sheds light on the structural homology of this domain across diverse organisms. Our smFRET-Rosetta models reveal nanometer-scale rearrangements within the RNA core domain during catalysis. Taken together, our FRET data and pseudoatomic molecular models permit us to propose a possible mechanism for how RNA core domain rearrangement is coupled to template hybrid elongation.

A better prediction of conformational changes of proteins using minimally connected network models

Elastic network models have recently been used for studying low-frequency collective motions of proteins. These models simplify the complexity that arises from normal mode analysis by considering a simplified potential involving a few parameters. Two common parameters in most of the elastic network models are cutoff radius and force constant. Although the latter has been studied extensively and even elaborate new models were introduced, for the former usually an ad-hoc cutoff radius is considered. Moreover, the quality of the network models is usually assessed by evaluating their prediction against experimental B-factors. In this work, we consider various common elastic network models with different cutoff radii and assess them by their ability to predict conformational changes of proteins in complexes from unbound to bound state. This prediction is performed by perturbing the unbound structure using a number of low-frequency normal modes of its network model to optimally fit the bound structure. We evaluated a dataset of 30 proteins with distinct unbound and bound structures using this criterion. The results showed that, opposed to the common calibration process based on B-factors, a meaningful relationship exists between the quality of the prediction and model parameters. It was shown that the cutoff radius has a major role in this prediction and minimally connected network models, which use the shortest cutoff radius for which the network is stable, give the best results. It was also shown that by considering the first ten normal modes, the conformational changes can be predicted by about 25 percent. Hence, the evaluation process was extended to the case of considering the contribution of all normal modes in the prediction. The results indicated that minimally connected network models are superior, despite their simplicity, when any number of modes are considered in the prediction.

One-Dimensional Structural Properties of Proteins in the Coarse-Grained CABS Model

Despite the significant increase in computational power, molecular modeling of protein structure using classical all-atom approaches remains inefficient, at least for most of the protein targets in the focus of biomedical research. Perhaps the most successful strategy to overcome the inefficiency problem is multiscale modeling to merge all-atom and coarse-grained models. This chapter describes a well-established CABS coarse-grained protein model. The CABS (C-Alpha, C-Beta, and Side chains) model assumes a 2-4 united-atom representation of amino acids, knowledge-based force field (derived from the statistical regularities seen in known protein sequences and structures) and efficient Monte Carlo sampling schemes (MC dynamics, MC replica-exchange, and combinations). A particular emphasis is given to the unique design of the CABS force-field, which is largely defined using one-dimensional structural properties of proteins, including protein secondary structure. This chapter also presents CABS-based modeling methods, including multiscale tools for de novo structure prediction, modeling of protein dynamics and prediction of protein-peptide complexes. CABS-based tools are freely available at http://biocomp.chem.uw.edu.pl/tools.

Telomerase Regulation from Beginning to the End

The vast body of literature regarding human telomere maintenance is a true testament to the importance of understanding telomere regulation in both normal and diseased states. In this review, our goal was simple: tell the telomerase story from the biogenesis of its parts to its maturity as a complex and function at its site of action, emphasizing new developments and how they contribute to the foundational knowledge of telomerase and telomere biology.

The lighthouse at the end of the chromosome

Fluorescence microscopy can be used to assess the dynamic localization and intensity of single entities in vitro or in living cells. It has been applied with aplomb to many different cellular processes and has significantly enlightened our understanding of the heterogeneity and complexity of biological systems. Recently, high-resolution fluorescence microscopy has been brought to bear on telomeres, leading to new insights into telomere spatial organization and accessibility, and into the mechanistic nuances of telomere elongation. We provide a snapshot of some of these recent advances with a focus on mammalian systems, and show how three-dimensional, time-lapse microscopy and single-molecule fluorescence shine a new light on the end of the chromosome.

Dihydropyrazole derivatives as telomerase inhibitors: Structure-based design, synthesis, SAR and anticancer evaluation in vitro and in vivo

It is of our interest to generate and identify novel compounds with regulation telomerase for cancer therapy. In order to carry out more rational design, based on structure-based drug design, several series of N-substituted-dihydropyrazole derivatives, totally 78 compounds as potential human telomerase inhibitors were designed and synthesized. The results demonstrated that some compounds had potent anticancer activity against four tumor cell lines, and showed good selectivity on tumor cells over somatic cells. By the modified TRAP assay, compound 13i exhibited the most potent inhibitory activity against telomerase with an IC50 value of 0.98 μM. In vivo evaluation results indicated that compound 13i could inhibit growth of S180 and HepG2 tumor-bearing mice, and it also significantly enhanced the survival rate of EAC tumor-bearing mice. The further results in vivo confirmed that it could significantly improve pathological changes of N,N-diethylnitrosamine (DEN)-induced rat hepatic tumor. These data support further studies to assess rational design of more efficient telomerase inhibitors in the future.

Mechanism of Folding and Binding of an Intrinsically Disordered Protein As Revealed by ab Initio Simulations

A complex of the phosphorylated kinase-inducible domain (pKID) with its interacting domain (KIX) is a model system for studies of mechanisms by which intrinsically unfolded proteins perform their functions. These mechanisms are not fully understood. Using an efficient coarse-grained model, ab initio simulations were performed of the coupled folding and binding of the pKID to the KIX. The simulations start from an unbound, randomly positioned and disordered pKID structure. During the simulations the pKID chain and its position remain completely unrestricted, while the KIX backbone is limited to near-native fluctuations. Ab initio simulations of such large-scale conformational transitions, unaffected by any knowledge about the bound pKID structure, remain inaccessible to classical simulations. Our simulations recover an ensemble of transient encounter complexes in good agreement with experimental results. We find that a key folding and binding step is linked to the formation of weak native interactions between a preformed nativelike fragment of a pKID helix and KIX surface. Once that nucleus forms, the pKID chain may condense from a largely disordered encounter ensemble to a natively bound and ordered conformation. The observed mechanism is reminiscent of a nucleation-condensation model, a common scenario for folding of globular proteins.

Potent Human Telomerase Inhibitors: Molecular Dynamic Simulations, Multiple Pharmacophore-Based Virtual Screening, and Biochemical Assays

Telomere maintenance is a universal cancer hallmark, and small molecules that disrupt telomere maintenance generally have anticancer properties. Since the vast majority of cancer cells utilize telomerase activity for telomere maintenance, the enzyme has been considered as an anticancer drug target. Recently, rational design of telomerase inhibitors was made possible by the determination of high resolution structures of the catalytic telomerase subunit from a beetle and subsequent molecular modeling of the human telomerase complex. A hybrid strategy including docking, pharmacophore-based virtual screening, and molecular dynamics simulations (MDS) were used to identify new human telomerase inhibitors. Docking methodology was applied to investigate the ssDNA telomeric sequence and two well-known human telomerase inhibitors' (BIBR1532 and MST-312) modes of interactions with hTERT TEN domain. Subsequently molecular dynamic simulations were performed to monitor and compare hTERT TEN domain, TEN-ssDNA, TEN-BIBR1532, TEN-MST-312, and TEN-ssDNA-BIBR1532 behavior in a dynamic environment. Pharmacophore models were generated considering the inhibitors manner in the TEN domain anchor site. These exploratory studies identified several new potent inhibitors whose IC50 values were generated experimentally in a low micromolar range with the aid of biochemical assays, including both the direct telomerase and the telomeric repeat amplification protocol (TRAP) assays. The results suggest that the current models of human telomerase are useful templates for rational inhibitor design.

Structural Basis of Telomerase Inhibition by the Highly Specific BIBR1532

BIBR1532 is a highly specific telomerase inhibitor, although the molecular basis for inhibition is unknown. Here we present the crystal structure of BIBR1532 bound to Tribolium castaneum catalytic subunit of telomerase (tcTERT). BIBR1532 binds to a conserved hydrophobic pocket (FVYL motif) on the outer surface of the thumb domain. The FVYL motif is near TRBD residues that bind the activation domain (CR4/5) of hTER. RNA binding assays show that the human TERT (hTERT) thumb domain binds the P6.1 stem loop of CR4/5 in vitro. hTERT mutations of the FVYL pocket alter wild-type CR4/5 binding and cause telomere attrition in cells. Furthermore, the hTERT FVYL mutations V1025F, N1028H, and V1090M are implicated in dyskeratosis congenita and aplastic anemia, further supporting the biological and clinical relevance of this novel motif. We propose that CR4/5 contacts with the telomerase thumb domain contribute to telomerase ribonucleoprotein assembly and promote enzymatic activity.

Modeling of protein-peptide interactions using the CABS-dock web server for binding site search and flexible docking

Protein-peptide interactions play essential functional roles in living organisms and their structural characterization is a hot subject of current experimental and theoretical research. Computational modeling of the structure of protein-peptide interactions is usually divided into two stages: prediction of the binding site at a protein receptor surface, and then docking (and modeling) the peptide structure into the known binding site. This paper presents a comprehensive CABS-dock method for the simultaneous search of binding sites and flexible protein-peptide docking, available as a user's friendly web server. We present example CABS-dock results obtained in the default CABS-dock mode and using its advanced options that enable the user to increase the range of flexibility for chosen receptor fragments or to exclude user-selected binding modes from docking search. Furthermore, we demonstrate a strategy to improve CABS-dock performance by assessing the quality of models with classical molecular dynamics. Finally, we discuss the promising extensions and applications of the CABS-dock method and provide a tutorial appendix for the convenient analysis and visualization of CABS-dock results. The CABS-dock web server is freely available at http://biocomp.chem.uw.edu.pl/CABSdock/.

Molecular modeling of the binding modes of the iron-sulfur protein to the Jac1 co-chaperone from Saccharomyces cerevisiae by all-atom and coarse-grained approaches

The iron-sulfur protein 1 (Isu1) and the J-type co-chaperone Jac1 from yeast are part of a huge ATP-dependent system, and both interact with Hsp70 chaperones. Interaction of Isu1 and Jac1 is a part of the iron-sulfur cluster biogenesis system in mitochondria. In this study, the structure and dynamics of the yeast Isu1-Jac1 complex has been modeled. First, the complete structure of Isu1 was obtained by homology modeling using the I-TASSER server and YASARA software and thereafter tested for stability in the all-atom force field AMBER. Then, the known experimental structure of Jac1 was adopted to obtain initial models of the Isu1-Jac1 complex by using the ZDOCK server for global and local docking and the AutoDock software for local docking. Three most probable models were subsequently subjected to the coarse-grained molecular dynamics simulations with the UNRES force field to obtain the final structures of the complex. In the most probable model, Isu1 binds to the left face of the Γ-shaped Jac1 molecule by the β-sheet section of Isu1. Residues L105 , L109 , and Y163 of Jac1 have been assessed by mutation studies to be essential for binding (Ciesielski et al., J Mol Biol 2012; 417:1-12). These residues were also found, by UNRES/molecular dynamics simulations, to be involved in strong interactions between Isu1 and Jac1 in the complex. Moreover, N(95), T(98), P(102), H(112), V(159), L(167), and A(170) of Jac1, not yet tested experimentally, were also found to be important in binding.

Novel myricetin derivatives: Design, synthesis and anticancer activity

Telomere and telomerase were closely related to occurrence and development of some cancers. To enhance ability of myricetin moiety for inhibiting telomerase, we designed a series of novel myricetin derivatives based on reasonable analysis. The telomerase inhibition assay showed that compound 6d displayed the most potent inhibitory activity with IC50 value of 0.91 μM. The anticancer activity assay showed that 6d exhibited high activity against human breast cells MDA-MB-231. The docking simulation of compound 6d was performed to get the probable binding model, the results demonstrated that the furan ring inserted into the active site deeply and had hydrophobic interactions with residues of Phe 568, Pro 627, four methoxy groups had hydrophobic interactions with residues of Phe 568, Pro 627, Lys 902, Val 904 and Pro 929. Western blot results showed that expression of p65 and TERT protein was clearly down-regulated by compound 6d. These data support further studies for the rational design of more efficient p65 and TERT modulators.

Design and synthesis of celastrol derivatives as anticancer agents

A series of celastrol derivatives as potential telomerase inhibitors were designed and synthesized. The bioassays demonstrated that title compounds displayed potent anticancer activities against SGC-7901, SMMC-7721, MGC-803 and HepG-2 cell lines, among them, compounds 3c and 3d which containing hydrophilicity moieties exhibited high anti-proliferative activities (IC50 = 0.10-1.22 μM). The preliminary mechanism of antitumor action indicated that title compound 3c could induce significant SMMC-7721 cells apoptosis. A modified TRAP assay showed that compounds 3c and 3d displayed the most potent inhibitory activity with IC50 values at 0.11 and 0.34 μM, respectively. And there was a good correlation between telomerase inhibition and anti-proliferative inhibition of SMMC-7721 cells. Moreover, molecular docking indicated that the active compound 3c was nicely bound into the telomerase hTERT active site, hydrophobic, van der Waals and two hydrogen bond interactions with conserved residues ASP 628 and TYR 949 were found.

Structure prediction of the second extracellular loop in G-protein-coupled receptors

G-protein-coupled receptors (GPCRs) play key roles in living organisms. Therefore, it is important to determine their functional structures. The second extracellular loop (ECL2) is a functionally important region of GPCRs, which poses significant challenge for computational structure prediction methods. In this work, we evaluated CABS, a well-established protein modeling tool for predicting ECL2 structure in 13 GPCRs. The ECL2s (with between 13 and 34 residues) are predicted in an environment of other extracellular loops being fully flexible and the transmembrane domain fixed in its x-ray conformation. The modeling procedure used theoretical predictions of ECL2 secondary structure and experimental constraints on disulfide bridges. Our approach yielded ensembles of low-energy conformers and the most populated conformers that contained models close to the available x-ray structures. The level of similarity between the predicted models and x-ray structures is comparable to that of other state-of-the-art computational methods. Our results extend other studies by including newly crystallized GPCRs.

Progress in structural studies of telomerase

Telomerase is the ribonucleoprotein (RNP) reverse transcriptase responsible for synthesizing the 3' ends of linear chromosomes. It plays critical roles in tumorigenesis, cellular aging, and stem cell renewal. The past two years have seen exciting progress in determining telomerase holoenzyme architecture and the structural basis of telomerase activity. Notably, the first electron microscopy structures of telomerase were reported, of the Tetrahymena thermophila telomerase holoenzyme and a human telomerase dimer. In addition to new structures of TERT and TER domains, the first structures of telomerase protein domains beyond TERT, and their complexes with TER or telomeric single-stranded DNA, were reported. Together these studies provide the first glimpse into the organization of the proteins and RNA in the telomerase RNP.

The architecture of Tetrahymena telomerase holoenzyme

Telomerase adds telomeric repeats to chromosome ends using an internal RNA template and a specialized telomerase reverse transcriptase (TERT), thereby maintaining genome integrity. Little is known about the physical relationships among protein and RNA subunits within a biologically functional holoenzyme. Here we describe the architecture of Tetrahymena thermophila telomerase holoenzyme determined by electron microscopy. Six of the seven proteins and the TERT-binding regions of telomerase RNA (TER) have been localized by affinity labelling. Fitting with high-resolution structures reveals the organization of TERT, TER and p65 in the ribonucleoprotein (RNP) catalytic core. p50 has an unanticipated role as a hub between the RNP catalytic core, p75-p19-p45 subcomplex, and the DNA-binding Teb1. A complete in vitro holoenzyme reconstitution assigns function to these interactions in processive telomeric repeat synthesis. These studies provide the first view of the extensive network of subunit associations necessary for telomerase holoenzyme assembly and physiological function.

Coarse grained normal mode analysis vs. refined Gaussian Network Model for protein residue-level structural fluctuations

We investigate several approaches to coarse grained normal mode analysis on protein residual-level structural fluctuations by choosing different ways of representing the residues and the forces among them. Single-atom representations using the backbone atoms C(α), C, N, and C(β) are considered. Combinations of some of these atoms are also tested. The force constants between the representative atoms are extracted from the Hessian matrix of the energy function and served as the force constants between the corresponding residues. The residue mean-square-fluctuations and their correlations with the experimental B-factors are calculated for a large set of proteins. The results are compared with all-atom normal mode analysis and the residue-level Gaussian Network Model. The coarse-grained methods perform more efficiently than all-atom normal mode analysis, while their B-factor correlations are also higher. Their B-factor correlations are comparable with those estimated by the Gaussian Network Model and in many cases better. The extracted force constants are surveyed for different pairs of residues with different numbers of separation residues in sequence. The statistical averages are used to build a refined Gaussian Network Model, which is able to predict residue-level structural fluctuations significantly better than the conventional Gaussian Network Model in many test cases.

Elastic network normal modes provide a basis for protein structure refinement

It is well recognized that thermal motions of atoms in the protein native state, the fluctuations about the minimum of the global free energy, are well reproduced by the simple elastic network models (ENMs) such as the anisotropic network model (ANM). Elastic network models represent protein dynamics as vibrations of a network of nodes (usually represented by positions of the heavy atoms or by the C(α) atoms only for coarse-grained representations) in which the spatially close nodes are connected by harmonic springs. These models provide a reliable representation of the fluctuational dynamics of proteins and RNA, and explain various conformational changes in protein structures including those important for ligand binding. In the present paper, we study the problem of protein structure refinement by analyzing thermal motions of proteins in non-native states. We represent the conformational space close to the native state by a set of decoys generated by the I-TASSER protein structure prediction server utilizing template-free modeling. The protein substates are selected by hierarchical structure clustering. The main finding is that thermal motions for some substates, overlap significantly with the deformations necessary to reach the native state. Additionally, more mobile residues yield higher overlaps with the required deformations than do the less mobile ones. These findings suggest that structural refinement of poorly resolved protein models can be significantly enhanced by reduction of the conformational space to the motions imposed by the dominant normal modes.

Single-molecule analysis of telomerase structure and function

The telomerase ribonucleoprotein is a specialized reverse transcriptase required to maintain protective chromosome end-capping structures called telomeres. In most cells, telomerase is not active and the natural shortening of telomeres with each round of DNA replication ultimately triggers cell growth arrest. In contrast, the presence of telomerase confers a high level of renewal capacity upon rapidly dividing cells. Telomerase is aberrantly activated in 90% of human cancers and thus represents an important target for anticancer therapeutics. However, the naturally low abundance of telomerase has hampered efforts to obtain high-resolution models for telomerase structure and function. To circumvent these challenges, single-molecule techniques have recently been employed to investigate telomerase assembly, structure, and catalysis.

The TPR-containing domain within Est1 homologs exhibits species-specific roles in telomerase interaction and telomere length homeostasis

Background

The first telomerase-associated protein (Est1) was isolated in yeast due to its essential role in telomere maintenance. The human counterparts EST1A, EST1B, and EST1C perform diverse functions in nonsense-mediated mRNA decay (NMD), telomere length homeostasis, and telomere transcription. Although Est1 and EST1A/B interact with the catalytic subunit of yeast and human telomerase (Est2 and TERT, respectively), the molecular determinants of these interactions have not been elaborated fully.

Results

To investigate the functional conservation of the EST1 protein family, we performed protein-protein interaction mapping and structure-function analysis. The domain in hEST1A most conserved between species, containing a TPR (tricotetrapeptide repeat), was sufficient for interaction of hEST1A with multiple fragments of hTERT including the N-terminus. Two mutations within the hTERT N-terminus that perturb in vivo function (NAAIRS₉₂, NAAIRS₁₂₂) did not affect this protein interaction. ScEst1 hybrids containing the TPR of hEST1A, hEST1B, or hEST1C were expressed in yeast strains lacking EST1, yet they failed to complement senescence. Point mutations within and outside the cognate ScEst1 TPR, chosen to disrupt a putative protein interaction surface, resulted in telomere lengthening or shortening without affecting recruitment to telomeres.

Conclusions

These results identify a domain encompassing the TPR of hEST1A as an hTERT interaction module. The TPR of S. cerevisiae Est1 is required for telomerase-mediated telomere length maintenance in a manner that appears separable from telomere recruitment. Discrete residues in or adjacent to the TPR of Est1 also regulate telomere length homeostasis.