Mutations in WDR62, encoding a centrosomal and nuclear protein, in Indian primary microcephaly families with cortical malformations

Primary microcephaly is an autosomal recessive disorder characterized by smaller than normal brain size and mental retardation. It is genetically heterogeneous with seven loci: MCPH1-MCPH7. We have previously reported genetic analysis of 35 families, including the identification of the MCPH7 gene STIL. Of the 35 families, three families showed linkage to the MCPH2 locus. Recent whole-exome sequencing studies have shown that the WDR62 gene, located in the MCPH2 candidate region, is mutated in patients with severe brain malformations. We therefore sequenced the WDR62 gene in our MCPH2 families and identified two novel homozygous protein truncating mutations in two families. Affected individuals in the two families had pachygyria, microlissencephaly, band heterotopias, gyral thickening, and dysplastic cortex. Using immunofluorescence study, we showed that, as with other MCPH proteins, WDR62 localizes to centrosomes in A549, HepG2, and HaCaT cells. In addition, WDR62 was also localized to nucleoli. Bioinformatics analysis predicted two overlapping nuclear localization signals and multiple WD-40 repeats in WDR62. Two other groups have also recently identified WDR62 mutations in MCPH2 families. Our results therefore add further evidence that WDR62 is the MCPH2 gene. The present findings will be helpful in genetic diagnosis of patients linked to the MCPH2 locus.

Advances in the Mathematical Modelling of Hepatitis C Virus Dynamics

Mathematical models have provided key insights into the pathogenesis of hepatitis C virus (HCV) in vivo, suggested predominant mechanism(s) of drug action, explained confounding patterns of viral load changes in HCV infected patients undergoing therapy, and presented a framework for therapy optimization. In this article, I present an overview of the major advances in the mathematical modeling of HCV dynamics.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

Molecular Docking Studies of Novel Palmitoyl-ligands for Cyclooxygenase-2

An in silico approach was adopted to identify potential cyclooxygenase-2 inhibitors through molecular docking studies. The in vivo studies indicated that synthetic palmitoyl derivatives of salicylic acid, para amino phenol, para amino benzoic acid, and anthranilic acid possessed significant pharmacological activities like anti-inflammatory, analgesic, and antipyretic activities. None of the tested substances produced any significant gastric lesions in experimental animals. In an attempt to understand the ligandprotein interactions in terms of the binding affinity, the above synthetic molecules were subjected to docking analysis using AutoDock. The palmitoyl derivatives palmitoyl anthranilic acid, palmitoyl para amino benzoic acid, palmitoyl para amino phenol, and palmitoyl salicylic acid showed better binding energy than the known inhibitor diclofenac bound to 1PXX. All the palmitoyl derivatives made similar interactions with the binding site residues of cyclooxygenase-2 as compared to that of the known inhibitor. Thus, structure-based drug discovery approach was successfully employed to identify some promising pro-drugs for the treatment of pain and inflammation.

Use of Vertex Index in Structure-Activity Analysis and Design of Molecules

The last few decades have witnessed application of graph theory and topological indices derived from molecular graph in structure-activity analysis. Such applications are based on regression and various multivariate analyses. Most of the topological indices are computed for the whole molecule and used as descriptors for explaining properties/activities of chemical compounds. However, some substructural descriptors in the form of topological distance based vertex indices have been found to be useful in identifying activity related substructures and in predicting pharmacological and toxicological activities of bioactive compounds. Another important aspect of drug discovery e. g. designing novel pharmaceutical candidates could also be done from the distance distribution associated with such vertex indices. In this article, we will review the development and applications of this approach both in activity prediction as well as in designing novel compounds.

Sequence-Structure Similarity: Do Sequentially Identical Peptide Fragments have Similar Three-Dimensional Structures?

The rapidly growing structure databases enhance the probability of finding identical sequences sharing structural similarity. Structure prediction methods are being used extensively to abridge the gap between known protein sequences and the solved structures which is essential to understand its specific biochemical and cellular functions. In this work, we plan to study the ambiguity between sequence-structure relationships and examine if sequentially identical peptide fragments adopt similar three-dimensional structures. Fragments of varying lengths (five to ten residues) were used to observe the behavior of sequence and its three-dimensional structures. The STAMP program was used to superpose the three-dimensional structures and the two parameters (Sequence Structure Similarity Score (Sc) and Root Mean Square Deviation value) were employed to classify them into three categories: similar, intermediate and dissimilar structures. Furthermore, the same approach was carried out on all the three-dimensional protein structures solved in the two organisms, Mycobacterium tuberculosis and Plasmodium falciparum to validate our results.

Suite of tools for statistical N-gram language modeling for pattern mining in whole genome sequences

Genome sequences contain a number of patterns that have biomedical significance. Repetitive sequences of various kinds are a primary component of most of the genomic sequence patterns. We extended the suffix-array based Biological Language Modeling Toolkit to compute n-gram frequencies as well as n-gram language-model based perplexity in windows over the whole genome sequence to find biologically relevant patterns. We present the suite of tools and their application for analysis on whole human genome sequence.

Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

HSPIR: a manually annotated heat shock protein information resource

Heat shock protein information resource (HSPIR) is a concerted database of six major heat shock proteins (HSPs), namely, Hsp70, Hsp40, Hsp60, Hsp90, Hsp100 and small HSP. The HSPs are essential for the survival of all living organisms, as they protect the conformations of proteins on exposure to various stress conditions. They are a highly conserved group of proteins involved in diverse physiological functions, including de novo folding, disaggregation and protein trafficking. Moreover, their critical role in the control of disease progression made them a prime target of research. Presently, limited information is available on HSPs in reference to their identification and structural classification across genera. To that extent, HSPIR provides manually curated information on sequence, structure, classification, ontology, domain organization, localization and possible biological functions extracted from UniProt, GenBank, Protein Data Bank and the literature. The database offers interactive search with incorporated tools, which enhances the analysis. HSPIR is a reliable resource for researchers exploring structure, function and evolution of HSPs.

Ferroelectric-ferroelectric phase coexistence in Na1/2Bi1/2TiO3

Morphotropic phase boundary (MPB) systems are characterized by the coexistence of two ferroelectric phases and are associated with anomalous piezoelectric properties. In general, such coexistence is brought about by composition induced ferroelectric-ferroelectric instability. Here we demonstrate that a pure ferroelectric compound Na1/2Bi1/2TiO3 (NBT) exhibits the coexistence of two ferroelectric phases, rhombohedral (R3c) and monoclinic (Cc), in its equilibrium state at room temperature. This was unravelled by adopting a unique strategy of comparative structural analysis of electrically poled and thermally annealed specimens using high resolution synchrotron x-ray powder diffraction data. The relative fraction of the coexisting phases was found to be highly sensitive to thermal, mechanical, and electrical stimuli. The coexistence of ferroelectric phases in the ground state of the pure compound will have significant bearing on the way MPB is induced in NBT-based lead-free piezoceramics. DOI: 10.1103/PhysRevB.87.060102

Identification of mycobacterial lectins from genomic data

Sixty-four sequences containing lectin domains with homologs of known three-dimensional structure were identified through a search of mycobacterial genomes. They appear to belong to the -prism II, the C-type, the Microcystis virdis (MV), and the -trefoil lectin folds. The first three always occur in conjunction with the LysM, the PI-PLC, and the -grasp domains, respectively while mycobacterial -trefoil lectins are unaccompanied by any other domain. Thirty heparin binding hemagglutinins (HBHA), already annotated, have also been included in the study although they have no homologs of known three-dimensional structure. The biological role of HBHA has been well characterized. A comparison between the sequences of the lectin from pathogenic and nonpathogenic mycobacteria provides insights into the carbohydrate binding region of the molecule, but the structure of the molecule is yet to be determined. A reasonable picture of the structural features of other mycobacterial proteins containing one or the other of the four lectin domains can be gleaned through the examination of homologs proteins, although the structure of none of them is available. Their biological role is also yet to be elucidated. The work presented here is among the first steps towards exploring the almost unexplored area of the structural biology of mycobacterial lectins. Proteins 2013. (c) 2012 Wiley Periodicals, Inc.

In vitro and in vivo characterization of designed immunogens derived from the CD-helix of the stem of influenza hemagglutinin

The conserved stem domain of influenza virus hemagglutinin (HA) is a target for broadly neutralizing antibodies and a potential vaccine antigen for induction of hetero-subtypic protection. The epitope of 12D1, a previously reported bnAb neutralizing several H3 subtype influenza strains, was putatively mapped to residues 76-106 of the CD-helix, also referred to as long alpha helix (LAH) of the HA stem. A peptide derivative consisting of wt-LAH residues 76-130 conjugated to keyhole limpet hemocyanin was previously shown to confer robust protection in mice against challenge with influenza strains of subtypes H3, H1, and H5 which motivated the present study. We report the design of multiple peptide derivatives of LAH with or without heterologous trimerization sequences and show that several of these are better folded than wt-LAH. However, in contrast to the previous study immunization of mice with wt-LAH resulted in negligible protection against a lethal homologous virus challenge, while some of the newly designed immunogens could confer weak protection. Combined with structural analysis of HA, our data suggest that in addition to LAH, other regions of HA are likely to significantly contribute to the epitope for 12D1 and will be required to elicit robust protection. In addition, a dynamic, flexible conformation of isolated LAH peptide may be required for eliciting a functional anti-viral response. Proteins 2013; 81:1759-1775. (c) 2013 Wiley Periodicals, Inc.

Information content of molecular graph and prediction of gas phase thermal entropy of organic compounds

Entropy is a fundamental thermodynamic property that has attracted a wide attention across domains, including chemistry. Inference of entropy of chemical compounds using various approaches has been a widely studied topic. However, many aspects of entropy in chemical compounds remain unexplained. In the present work, we propose two new information-theoretical molecular descriptors for the prediction of gas phase thermal entropy of organic compounds. The descriptors reflect the bulk and size of the compounds as well as the gross topological symmetry in their structures, all of which are believed to determine entropy. A high correlation () between the entropy values and our information-theoretical indices have been found and the predicted entropy values, obtained from the corresponding statistically significant regression model, have been found to be within acceptable approximation. We provide additional mathematical result in the form of a theorem and proof that might further help in assessing changes in gas phase thermal entropy values with the changes in molecular structures. The proposed information-theoretical molecular descriptors, regression model and the mathematical result are expected to augment predictions of gas phase thermal entropy for a large number of chemical compounds.

G-Quadruplex Structures Formed at the HOX11 Breakpoint Region Contribute to Its Fragility during t(10;14) Translocation in T-Cell Leukemia

The t(10;14) translocation involving the HOX11 gene is found in several T-cell leukemia patients. Previous efforts to determine the causes of HOX11 fragility were not successful. The role of non-B DNA structures is increasingly becoming an important cause of genomic instability. In the present study, bioinformatics analysis revealed two G-quadruplex-forming motifs at the HOX11 breakpoint cluster. Gel shift assays showed formation of both intra- and intermolecular G-quadruplexes, the latter being more predominant. The structure formation was dependent on four stretches of guanines, as revealed by mutagenesis. Circular dichroism analysis identified parallel conformations for both quadruplexes. The non-B DNA structure could block polymerization during replication on a plasmid, resulting in consistent K K+-dependent pause sites, which were abolished upon mutation of G-motifs, thereby demonstrating the role of the stretches of guanines even on double-stranded DNA. Extrachromosomal assays showed that the G-quadruplex motifs could block transcription, leading to reduced expression of green fluorescent protein (GFP) within cells. More importantly, sodium bisulfite modification assay showed the single-stranded character at regions I and II of HOX11 in the genome. Thus, our findings suggest the occurrence of G-quadruplex structures at the HOX11 breakpoint region, which could explain its fragility during the t(10;14) translocation.