Molecular genetics of tooth agenesis

Congenital missing of teeth, tooth agenesis or hypodontia, is one of the most common developmental anomalies in man. The common forms in which one or a few teeth are absent, may cause occlusal or cosmetic harm, while severe forms which are relatively rare always require clinical attention to support and maintain the dental function. Observation of tooth agenesis is also important for diagnosis of malformation syndromes. Some external factors may cause developmental defects and agenesis in dentition. However, the role of inheritance in the etiology of tooth agenesis is well established by twin and family studies. Studies on familial tooth agenesis as well as mouse null mutants have also identified several genetic factors. However, these explain syndromic or rare dominant forms of tooth agenesis, whereas the genes and defects responsible for the majority of cases of tooth agenesis, especially the common and less severe forms, are largely unknown. In this study it was shown, that a dominant nonsense mutation in PAX9 was responsible for severe tooth agenesis (oligodontia) in a Finnish family. In a study of tooth agenesis associated with Wolf-Hirschhorn syndrome, it was shown that severe tooth agenesis was present if the causative deletion in 4p spanned the MSX1 locus. It was concluded that severe tooth agenesis was caused by haploinsufficiency of these transcription factors. A summary of the phenotypes associated with known defects in MSX1 and PAX9 showed that, despite similarities, they were significantly different, suggesting that the genes, in addition to known interactions, also have independent roles during the development of human dentition. The original aim of this work was to identify gene defects that underlie the common incisor and premolar hypodontia. After excluding several candidate genes, a genome-wide search was conducted in seven Finnish families in which this phenotype was inherited in an autosomal dominant manner. A promising locus for second premolar agenesis was identified in chromosome 18 in one family and this finding was supported by results from other families. The results also implied the existence of other loci both for second premolar agenesis and for incisor agenesis. On the other hand the results did not lend support for comprehensive involvement of the most obvious candidate genes in the etiology of incisor and premolar hypodontia. Rather, they suggest remarkable genetic heterogeneity of tooth agenesis. The available evidence suggests that quantitative defects during tooth development predispose to a failure to overcome a developmental threshold and to agenesis. The results of the study increase the understanding of the etiology and heredity of tooth agenesis. Further studies may lead to identification of novel genes that affect the development of teeth.

New structural insights and molecular-modelling studies of 4-methyl-5-beta-hydroxyethylthiazole kinase from Pyrococcus horikoshii OT3 (PhThiK)

4-Methyl-5-beta-hydroxyethylthiazole kinase (ThiK) catalyses the phosphorylation of the hydroxyl group of 4-methyl-5-beta-hydroxyethylthiazole. This work reports the first crystal structure of an archaeal ThiK: that from Pyrococcus horikoshii OT3 (PhThiK) at 1.85 angstrom resolution with a phosphate ion occupying the position of the beta-phosphate of the nucleotide. The topology of this enzyme shows the typical ribokinase fold of an alpha/beta protein. The overall structure of PhThiK is similar to those of Bacillus subtilis ThiK (BsThiK) and Enterococcus faecalis V583 ThiK (EfThiK). Sequence analysis of ThiK enzymes from various sources indicated that three-quarters of the residues involved in interfacial regions are conserved. It also revealed that the amino-acid residues in the nucleotide-binding, magnesium ion-binding and substrate-binding sites are conserved. Binding of the nucleotide and substrate to the ThiK enzyme do not influence the quaternary association (trimer) as revealed by the crystal structure of PhThiK.

CADASIL: molecular studies on the most common hereditary vascular dementing disorder

Cerebral Autosomal Dominant Arteriopathy with Subcortical Infarcts and Leukoencephalopathy (CADASIL) is the most common hereditary vascular dementia. CADASIL is a systemic disease of small and medium-sized arteries although the symptoms are almost exclusively neurological, including migraineous headache, recurrent ischemic episodes, cognitive impairment and, finally, subcortical dementia. CADASIL is caused by over 170 different mutations in the NOTCH3 gene, which encodes a receptor expressed in adults predominantly in the vascular smooth muscle cells. The function of NOTCH3 is not crucial for embryonic development but is needed after birth. NOTCH3 directs postnatal arterial maturation and helps to maintain arterial integrity. It is involved in regulation of vascular tone and in the wound healing of a vascular injury. In addition, NOTCH3 promotes cell survival by inducing expression of anti-apoptotic proteins. NOTCH3 is a membrane-spanning protein with a large extracellular domain (N3ECD) containing 34 epidermal growth factor-like (EGF) repeats and a smaller intracellular domain with six ankyrin repeats. All CADASIL mutations are located in the EGF repeats and the majority of the mutations cause gain or loss of one cysteine residue in one of these repeats leading to an odd number of cysteine residues, which in turn leads to misfolding of N3ECD. This misfolding most likely alters the maturation, targetting, degradation and/or function of the NOTCH3 receptor. CADASIL mutations do not seem to affect the canonical NOTCH3 signalling pathway. The main pathological findings are the accumulation of the NOTCH3 extracellular domain on degenerating vascular smooth muscle cells (VSMCs), accumulation of granular osmiophilic material (GOM) in the close vicinity of VSMCs as well as fibrosis and thickening of arterial walls. Narrowing of the arterial lumen and local thrombosis cause insufficient blood flow, mainly in small arteries of the cerebral white matter, resulting in tissue damage and lacunar infarcts. CADASIL is suspected in patients with a suggestive family history and clinical picture as well as characteristic white matter alterations in magnetic resonance imaging. A definitive verification of the diagnosis can be achieved by identifying a pathogenic mutation in the NOTCH3 gene or through the detection of GOM by electron microscopy. To understand the pathology underlying CADASIL, we have generated a unique set of cultured vascular smooth muscle cell (VSMC) lines from umbilical cord, placental, systemic and cerebral arteries of CADASIL patients and controls. Analyses of these VSMCs suggest that mutated NOTCH3 is misfolded, thus causing endoplasmic reticulum stress, activation of the unfolded protein response and increased production of reactive oxygen species. In addition, mutation in NOTCH3 causes alterations in actin cytoskeletal structures and protein expression, increased branching and abnormal node formation. These changes correlate with NOTCH3 expression levels within different VSMCs lines, suggesting that the phenotypic differences of SMCs may affect the vulnerability of the VSMCs and, therefore, the pathogenic impact of mutated NOTCH3 appears to vary in the arteries of different locations. Furthermore, we identified PDGFR- as an immediate downstream target gene of NOTCH3 signalling. Activation of NOTCH induces up-regulation of the PDGFR- expression in control VSMCs, whereas this up-regulation is impaired in CADASIL VSMCs and might thus serve as an alternative molecular mechanism that contributes to CADASIL pathology. In addition, we have established the congruence between NOTCH3 mutations and electron microscopic detection of GOM with a view to constructing a strategy for CADASIL diagnostics. In cases where the genetic analysis is not available or the mutation is difficult to identify, a skin biopsy is an easy-to-perform and highly reliable diagnostic method. Importantly, it is invaluable in setting guidelines concerning how far one should proceed with the genetic analyses.

The Crystal and Molecular Structure of l-(Diphenylmethyl)azetidin-3-ol

1-(Diphenylmethyl)azetidin-3-ol is triclinic, space group P1, with a=8.479(2), b=17.294(4),c = 10.606 (3) A, a = 118.59 (2),/~ = 100.30 (2), y = 89.63 (2) °, Z = 4. The structure was solved by multisolution methods and refined to an R of 0.044 for 2755 reflexions. The four-membered rings in the two independent molecules are puckered with dihedral angles of 156 and 153 ° . The two molecules differ in conformation with respect to rotation of the phenyl rings about the C-C bonds. The structure is stabilized by a network of O-H. • • N intermolecular hydrogen bonds.

Oxazolidine-2-thiones: A molecular modeling study

Two oxazolidine-2-thiones, thio-analogs of linezolid, were synthesized and their antibacterial properties evaluated. Unlike oxazolidinones, the thio-analogs did not inhibit the growth of Gram positive bacteria. A molecular modeling study has been carried out to aid understanding of this unexpected finding.

Phylogeography and Adaptive Divergence of Three-spined Stickleback Populations

Genetic studies on phylogeography and adaptive divergence in Northern Hemisphere fish species such as three-spined stickleback (Gasterosteus aculeatus) provide an excellent opportunity to investigate genetic mechanisms underlying population differentiation. According to the theory, the process of population differentiation results from a complex interplay between random and deterministic processes as well historical factors. The main scope in this thesis was to study how historical factors like the Pleistocene ice ages have shaped the patterns molecular diversity in three-spined stickleback populations in Europe and how this information could be utilized in the conservation genetic context. Furthermore, identifying footprints of natural selection at the DNA level might be used in identifying genes involved in evolutionary change. Overall, the results from phylogeographic studies indicate that the three-spined stickleback has colonized the Atlantic basin relatively recently but constitutes three major evolutionary lineages in Europe. In addition, the colonization of freshwater appears to result from multiple and independent invasions by the marine conspecifics. Molecular data together with morphology suggest that the most divergent freshwater populations are located in the Balkan Peninsula and these populations deserve a special conservation genetic status without warranting further taxonomical classification. In order to investigate the adaptive divergence in Fennoscandian three-spined stickleback populations several approaches were used. First, sequence variability in the Eda-gene, coding for the number of lateral plates, was concordant with the previously observed global pattern. Full plated allele is in high frequencies among marine populations whereas low plated allele dominates in the freshwater populations. Second, a microsatellite based genome scan identified both indications of balancing and directional selection in the three-spined stickleback genome, i.e. loci with unusually similar or unusually different allele frequencies over populations. The directionally selected loci were mainly associated with the adaptation to freshwater. A follow up study conducting a more detailed analysis in a chromosome region containing a putatively selected gene locus identified a fairly large genomic region affected by natural selection. However, this region contained several gene predictions, all of which might be the actual target of natural selection. All in all, the phylogeographic and adaptive divergence studies indicate that most of the genetic divergence has occurred in the freshwater populations whereas the marine populations have remained relatively uniform.

Platelet endothelial cell adhesion molecule 1 (PECAM-1) and its interactions with glycosaminoglycans: 1. Molecular modeling studies

Platelet endothelial cell adhesion molecule 1 (PECAM-1) has many functions, including its roles in leukocyte extravasation as part of the inflammatory response and in the maintenance of vascular integrity through its contribution to endothelial cell−cell adhesion. PECAM-1 has been shown to mediate cell−cell adhesion through homophilic binding events that involve interactions between domain 1 of PECAM-1 molecules on adjacent cells. However, various heterophilic ligands of PECAM-1 have also been proposed. The possible interaction of PECAM-1 with glycosaminoglycans (GAGs) is the focus of this study. The three-dimensional structure of the extracellular immunoglobulin (Ig) domains of PECAM-1 were constructed using homology modeling and threading methods. Potential heparin/heparan sulfate-binding sites were predicted on the basis of their amino acid consensus sequences and a comparison with known structures of sulfate-binding proteins. Heparin and other GAG fragments have been docked to investigate the structural determinants of their protein-binding specificity and selectivity. The modeling has predicted two regions in PECAM-1 that appear to bind heparin oligosaccharides. A high-affinity binding site was located in Ig domains 2 and 3, and evidence for a low-affinity site in Ig domains 5 and 6 was obtained. These GAG-binding regions were distinct from regions involved in PECAM-1 homophilic interactions.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

Diaryl dihydropyrazole-3-carboxamides with significant in vivo antiobesity activity related to CB1 receptor antagonism: Synthesis, biological evaluation, and molecular modeling in the homology model

A number of analogues of diaryl dihydropyrazole-3-carboxamides have been synthesized. Their activities were evaluated for appetite suppression and body weight reduction in animal models. Depending on the chemical modification of the selected dihydropyrazole scaffold, the lead compoundsthe bisulfate salt of (±)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 26 and the bisulfate salt of (−)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 30showed significant body weight reduction in vivo, which is attributed to their CB1 antagonistic activity and exhibited a favorable pharmacokinetic profile. The molecular modeling studies also showed interactions of two isomers of (±)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 9 with CB1 receptor in the homology model similar to those of N-piperidino-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4-methyl-3-pyrazole-carboxamide (rimonabant) 1 and 4S-(−)-3-(4-chlorophenyl)-N-methyl-N‘-[(4-chlorophenyl)-sulfonyl]-4-phenyl-4,5-dihydro-1H-pyrazole-1-carboxamidine (SLV-319) 2.

Molecular modeling of Bt Cry1Ac (DI–DII)–ASAL (Allium sativum lectin)–fusion protein and its interaction with aminopeptidase N (APN) receptor of Manduca sexta

Genetic engineering of Bacillus thuringiensis (Bt) Cry proteins has resulted in the synthesis of various novel toxin proteins with enhanced insecticidal activity and specificity towards different insect pests. In this study, a fusion protein consisting of the DI–DII domains of Cry1Ac and garlic lectin (ASAL) has been designed in silico by replacing the DIII domain of Cry1Ac with ASAL. The binding interface between the DI–DII domains of Cry1Ac and lectin has been identified using protein–protein docking studies. Free energy of binding calculations and interaction profiles between the Cry1Ac and lectin domains confirmed the stability of fusion protein. A total of 18 hydrogen bonds was observed in the DI–DII–lectin fusion protein compared to 11 hydrogen bonds in the Cry1Ac (DI–DII–DIII) protein. Molecular mechanics/Poisson–Boltzmann (generalized-Born) surface area [MM/PB (GB) SA] methods were used for predicting free energy of interactions of the fusion proteins. Protein–protein docking studies based on the number of hydrogen bonds, hydrophobic interactions, aromatic–aromatic, aromatic–sulphur, cation–pi interactions and binding energy of Cry1Ac/fusion proteins with the aminopeptidase N (APN) of Manduca sexta rationalised the higher binding affinity of the fusion protein with the APN receptor compared to that of the Cry1Ac–APN complex, as predicted by ZDOCK, Rosetta and ClusPro analysis. The molecular binding interface between the fusion protein and the APN receptor is well packed, analogously to that of the Cry1Ac–APN complex. These findings offer scope for the design and development of customized fusion molecules for improved pest management in crop plants.