Structural genomics: computational methods for structure analysis

The success of structural genomics initiatives requires the development and application of tools for structure analysis, prediction, and annotation. In this paper we review recent developments in these areas; specifically structure alignment, the detection of remote homologs and analogs, homology modeling and the use of structures to predict function. We also discuss various rationales for structural genomics initiatives. These include the structure-based clustering of sequence space and genome-wide function assignment. It is also argued that structural genomics can be integrated into more traditional biological research if specific biological questions are included in target selection strategies.     

Therapeutics development for triplet repeat expansion diseases

The underlying genetic mutations for many inherited neurodegenerative disorders have been identified in recent years. One frequent type of mutation is trinucleotide repeat expansion. Depending on the location of the repeat expansion, the mutation might result in a loss of function of the disease gene, a toxic gain of function or both. Disease gene identification has led to the development of model systems for investigating disease mechanisms and evaluating treatments. Examination of experimental findings reveals similarities in disease mechanisms as well as possibilities for treatment.     

Evolution of a triplet repeat in a conifer.

The opportunity to trace the evolution of a triplet repeat is rare, especially for seed-plant lineages with a well-defined fossil record. Microsatellite PtTX2133 sequences from 18 species in 2 conifer genera were used to calibrate the birth of a CAGn repeat, from its protomicrosatellite origins to its repeat expansion. Birth occurred in the hard-pine genome ~ 136 million years ago, or 14 million generations ago, then expanded as a polymorphic triplet repeat 136-100 million years before a major North American vicariance event. Calibration of the triplet-repeat birth and expansion is supported by the shared allelic lineages among Old and New World hard pines and the shared alleles solely among North American diploxylon or hard pines. Five CAGn repeat units appeared to be the expansion threshold for Old and New World diploxylon pines. Haploxylon pine species worldwide did not undergo birth and repeat expansion, remaining monomorphic, with a single imperfect 198-bp allele. A sister genus, Picea, had only a region of cryptic simplicity, preceding a proto-microsatellite region. The polymorphic triplet repeat in hard pines is older than some long-lived microsatellites reported for reptiles, yet younger than those reported for insects. Some cautionary points are raised about phylogenetic applications for this long-lived microsatellite.     

Energetic coupling between clustered lesions modulated by intervening triplet repeat bulge loops: Allosteric implications for DNA repair and triplet repeat expansion

Clusters of closely spaced oxidative DNA lesions present challenges to the cellular repair machinery. When located in opposing strands, base excision repair (BER) of such lesions can lead to double strand DNA breaks (DSB). Activation of BER and DSB repair pathways has been implicated in inducing enhanced expansion of triplet repeat sequences. We show here that energy coupling between distal lesions (8oxodG and/or abasic sites) in opposing DNA strands can be modulated by a triplet repeat bulge loop located between the lesion sites. We find this modulation to be dependent on the identity of the lesions (8oxodG vs. abasic site) and the positions of the lesions (upstream vs. downstream) relative to the intervening bulge loop domain. We discuss how such bulge loop‐mediated lesion crosstalk might influence repair processes, while favoring DNA expansion, the genotype of triplet repeat diseases. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 355–369, 2010. This article was originally published online as an acceptedpreprint. The “Published Online” date corresponds to the preprint version. You can reqest a copy of the preprint byemailing the Biopolymers editorial office at biopolymers@wiley.com     

Structural basis of VDR-DNA interactions on direct repeat response elements

The vitamin D receptor (VDR) forms homo- or heterodimers on response elements composed of two hexameric half-sites separated by 3 bp of spacer DNA. We describe here the crystal structures at 2.7-2.8 A resolution of the VDR DNA-binding region (DBD) in complex with response elements from three different promoters: osteopontin (SPP), canonical DR3 and osteocalcin (OC). These structures reveal the chemical basis for the increased affinity of VDR for the SPP response element, and for the poor stability of the VDR-OC complex, relative to the canonical DR3 response element. The homodimeric protein-protein interface is stabilized by van der Waals interactions and is predominantly non-polar. An extensive alpha-helix at the C-terminal end of the VDR DBD resembles that found in the thyroid hormone receptor (TR), and suggests a mechanism by which VDR and TR discriminate among response elements. Selective structure-based mutations in the asymmetric homodimeric interface result in a VDR DBD protein that is defective in homodimerization but now forms heterodimers with the 9-cis retinoic acid receptor (RXR) DBD.     

Structural and transcriptional analysis of a human subtelomeric repeat

A human subtelomeric repeat (designated as the HST repeat) has been isolated and characterized from a yeast artificial chromosome containing one human telomere. This repeat is located immediately adjacent to the telomeric T2AG3 repeats at the extreme termini of the human chromosomes. The DNA sequence of 3.6 kb of the HST repeat has been determined. The HST repeat spans over 3.6 kb in length, and contains one evolutionarily conserved CpG-rich region. The copy number of the HST repeat varies among telomeres. Genomic hybridization experiments suggest that the HST repeat consists of two distinct segments, and the distal portions of the HST repeat are also distributed elsewhere in the genome. In HeLa cells, the HST repeat sequence appears to be transcribed into a 6 kb polyadenylated RNA and a variety of non-polyadenylated RNA species.     

Structural and computational basis for potent inhibition of glutamate carboxypeptidase II by carbamate-based inhibitors

A series of carbamate-based inhibitors of glutamate carboxypeptidase II (GCPII) were designed and synthesized using ZJ-43, N-[[[(1S)-1-carboxy-3-methylbutyl]amino]carbonyl]-l-glutamic acid, as a molecular template in order to better understand the impact of replacing one of the two nitrogen atoms in the urea-based GCPII inhibitor with an oxygen atom. Compound 7 containing a C-terminal 2-oxypentanedioic acid was more potent than compound 5 containing a C-terminal glutamic acid (2-aminopentanedioic acid) despite GCPII’s preference for peptides containing an N-terminal glutamate as substrates. Subsequent crystallographic analysis revealed that ZJ-43 and its two carbamate analogs 5 and 7 with the same (S,S)-stereochemical configuration adopt a nearly identical binding mode while (R,S)-carbamate analog 8 containing a d-leucine forms a less extensive hydrogen bonding network. QM and QM/MM calculations have identified no specific interactions in the GCPII active site that would distinguish ZJ-43 from compounds 5 and 7 and attributed the higher potency of ZJ-43 and compound 7 to the free energy changes associated with the transfer of the ligand from bulk solvent to the protein active site as a result of the lower ligand strain energy and solvation/desolvation energy. Our findings underscore a broader range of factors that need to be taken into account in predicting ligand-protein binding affinity. These insights should be of particular importance in future efforts to design and develop GCPII inhibitors for optimal inhibitory potency.     

Unstable nucleotide repeat minireview series: a molecular biography of unstable repeat disorders

Expansion of an unstable nucleotide repeat is a mutational mechanism that is apparently unique to humans and is known to cause a variety of neurological disorders. This collection of minireviews examines several of these unstable repeats, focusing on those where there is considerable molecular information on how the mutation alters function.     

scdNet: a computational tool for single-cell differential network analysis

Background

Single-cell RNA sequencing (scRNA-Seq) is an emerging technology that has revolutionized the research of the tumor heterogeneity. However, the highly sparse data matrices generated by the technology have posed an obstacle to the analysis of differential gene regulatory networks.

Results

Addressing the challenges, this study presents, as far as we know, the first bioinformatics tool for scRNA-Seq-based differential network analysis (scdNet). The tool features a sample size adjustment of gene-gene correlation, comparison of inter-state correlations, and construction of differential networks. A simulation analysis demonstrated the power of scdNet in the analyses of sparse scRNA-Seq data matrices, with low requirement on the sample size, high computation efficiency, and tolerance of sequencing noises. Applying the tool to analyze two datasets of single circulating tumor cells (CTCs) of prostate cancer and early mouse embryos, our data demonstrated that differential gene regulation plays crucial roles in anti-androgen resistance and early embryonic development.

Conclusions

Overall, the tool is widely applicable to datasets generated by the emerging technology to bring biological insights into tumor heterogeneity and other studies. MATLAB implementation of scdNet is available at https://github.com/ChenLabGCCRI/scdNet.

     

Structural basis of kinetic variations in Sucrose Phosphate Phosphatase (SPP) of rice and Anabaena unveiled through computational analysis

Sucrose is an important storage form of assimilated carbon in many plant species. Unlike other sucrose biosynthetic enzymes, Sucrose Phosphate Phosphatase (SPP), the terminal enzyme in sucrose biosynthetic pathway, is the least understood. SPPs from different organisms have different kinetic properties. The current study focuses on the structural differences among SPP homologues and unveils the probable structural basis of kinetic variations. We have employed computational methods of molecular modeling and structure comparisons and identified structural variations in some of the substrate binding residues, amino acid substitutions in regions that are lining the active site and minute structural differences that can enhance the nucleophilicity of a catalytic nucleophile (Asp 9 ). We report a structurally and hence functionally important amino acid substitution (Asp 159 by Alanine) in one of the rice SPP isoforms, which can result in the disruption of a H-bond that helps in binding of sucrose at the active site of the enzyme. In this paper we discuss the structural basis of enhanced catalytic efficiency of rice SPP in comparison with a cyanobacterium (Anabaena variabilis). The natural mutations identified in our analysis of the SPP catalytic domain would be useful in re-designing the enzyme for enhanced catalytic efficiency and higher sucrose production.     

Structural basis for catalytic activation of a serine recombinase

Sin resolvase is a site-specific serine recombinase that is normally controlled by a complex regulatory mechanism. A single mutation, Q115R, allows the enzyme to bypass the entire regulatory apparatus, such that no accessory proteins or DNA sites are required. Here, we present a 1.86 ? crystal structure of the Sin Q115R catalytic domain, in a tetrameric arrangement stabilized by an interaction between Arg115 residues on neighboring subunits. The subunits have undergone significant conformational changes from the inactive dimeric state previously reported. The structure provides a new high-resolution view of a serine recombinase active site that is apparently fully assembled, suggesting roles for the conserved active site residues. The structure also suggests how the dimer-tetramer transition is coupled to assembly of the active site. The tetramer is captured in a different rotational substate than that seen in previous hyperactive serine recombinase structures, and unbroken crossover site DNA can be readily modeled into its active sites.     

Structural basis for calmodulin as a dynamic calcium sensor

Calmodulin is a prototypical and versatile Ca(2+) sensor with EF hands as its high-affinity Ca(2+) binding domains. Calmodulin is present in all eukaryotic cells, mediating Ca(2+)-dependent signaling. Upon binding Ca(2+), calmodulin changes its conformation to form complexes with a diverse array of target proteins. Despite a wealth of knowledge on calmodulin, little is known on how target proteins regulate calmodulin's ability to bind Ca(2+). Here, we take advantage of two splice variants of SK2 channels, which are activated by Ca(2+)-bound calmodulin but show different sensitivity to Ca(2+) for their activation. Protein crystal structures and other experiments show that, depending on which SK2 splice variant it binds to, calmodulin adopts drastically different conformations with different affinities for Ca(2+) at its C-lobe. Such target protein-induced conformational changes make calmodulin a dynamic Ca(2+) sensor capable of responding to different Ca(2+) concentrations in cellular Ca(2+) signaling.     

Structural basis for paramyxovirus-mediated membrane fusion

Paramyxoviruses are responsible for significant human mortality and disease worldwide, but the molecular mechanisms underlying their entry into host cells remain poorly understood. We have solved the crystal structure of a fragment of the simian parainfluenza virus 5 fusion protein (SV5 F), revealing a 96 A long coiled coil surrounded by three antiparallel helices. This structure places the fusion and transmembrane anchor of SV5 F in close proximity with a large intervening domain at the opposite end of the coiled coil. Six amino acids, potentially part of the fusion peptide, form a segment of the central coiled coil, suggesting that this structure extends into the membrane. Deletion mutants of SV5 F indicate that putative flexible tethers between the coiled coil and the viral membrane are dispensable for fusion. The lack of flexible tethers may couple a final conformational change in the F protein directly to the fusion of two bilayers.     

Structural basis for sensory rhodopsin function

The crystal structure of sensory rhodopsin II from Natronobacterium pharaonis was recently solved at 2.1 A resolution from lipidic cubic phase-grown crystals. A critical analysis of previous structure-function studies is possible within the framework of the high-resolution structure of this photoreceptor. Based on the structure, a molecular understanding emerges of the efficiency and selectivity of the photoisomerization reaction, of the interaction of the sensory receptor and its cognate transducer protein HtrII, and of the mechanism of spectral tuning in photoreceptors. The architecture of the retinal binding pocket is compact, representing a major determinant for the selective binding of the chromophore, all-trans retinal to the apoprotein, opsin. Several chromophore-protein interactions revealed by the structure were not predicted by previous mutagenesis and spectroscopic analyses. The structure suggests likely mechanisms by which photoisomerization triggers the activation of sensory rhodopsin II, and highlights the possibility of a unified mechanism of signaling mediated by sensory receptors, including visual rhodopsins. Future investigations using time-resolved crystallography, structural dynamics, and computational studies will provide the basis to unveil the molecular mechanisms of sensory receptors-mediated transmembrane signaling.