Alteration of the C-terminal ligand specificity of the Erbin PDZ domain by allosteric mutational effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.

Groucho homologue Grg5 interacts with the transcription factor Runx2-Cbfa1 and modulates its activity during postnatal growth in mice

Runx2-Cbfal, a Runt transcription factor, plays important roles during skeletal development. It is required for differentiation and function of osteoblasts. In its absence, chondrocyte hypertrophy is severely impaired and there is no vascularization of cartilage templates during skeletal development. These tissue-specific functions of Runx2 are likely to be dependent on its interaction with other proteins. We have therefore searched for proteins that may modulate the activity of Runx2. The yeast two-hybrid system was used to identify a groucho homologue, Grg5, as a Runx2-interacting protein. Grg5 enhances Runx2 activity in a cell culture-based assay and by analyses of postnatal growth in mice we demonstrate that Grg5 and Runx2 interact genetically. We also show that Runx2 haploinsufficiency in the absence of Grg5 results in a more severe delay in ossification of cranial sutures and fontanels than occurs with Runx2 haploinsufficiency on a wild-type background. Finally, we find shortening of the proliferative and hypertrophic zones, and expansion of the resting zone in the growth plates of Runx2(+/-)Grg5(-/-) mice that are associated with reduced Ihh expression and Indian hedgehog (Ihh) signaling. We therefore conclude that Grg5 enhances Runx2 activity in vivo.

Androgenetic alopecia: Identification of four genetic risk loci and evidence for the contribution of WNT signaling to its etiology

The pathogenesis of androgenetic alopecia (AGA, male-pattern baldness) is driven by androgens, and genetic predisposition is the major prerequisite. Candidate gene and genome-wide association studies have reported that single-nucleotide polymorphisms (SNPs) at eight different genomic loci are associated with AGA development. However, a significant fraction of the overall heritable risk still awaits identification. Furthermore, the understanding of the pathophysiology of AGA is incomplete, and each newly associated locus may provide novel insights into contributing biological pathways. The aim of this study was to identify unknown AGA risk loci by replicating SNPs at the 12 genomic loci that showed suggestive association (5 x 10(-8)signaling in AGA development.

Proinflammatory cytokines regulate cementogenic differentiation of periodontal ligament cells by Wnt/Ca2+ signaling pathway

Periodontal inflammation can inhibit cell differentiation of periodontal ligament cells (PDLCs), resulting in decreased bone/cementum regeneration ability. The Wnt signaling pathway, including canonical Wnt/β-catenin signaling and noncanonical Wnt/Ca2+ signaling, plays essential roles in cell proliferation and differentiation during tooth development. However, little is still known whether noncanonical Wnt/Ca2+ signaling cascade could regulate cementogenic/osteogenic differentiation capability of PDLCs within an inflammatory environment. Therefore, in this study, human PDLCs (hPDLCs) and their cementogenic differentiation potential were investigated in the presence of cytokines. The data demonstrated that both cytokines interleukin-6 (IL-6) and tumor necrosis factor alpha (TNF-α) inhibited cell proliferation, relative alkaline phosphatase activity, bone/cementum-related gene/protein expression, and canonical Wnt pathway-related gene/protein expression in hPDLCs. Interestingly, both cytokines upregulated the noncanonical Wnt/Ca2+ signaling-related gene and protein expression in hPDLCs. When the Wnt/Ca2+ pathway was blocked by Ca2+/calmodulin-dependent protein kinase II inhibitor KN93, even in the presence of IL-6 and TNF-α, cementogenesis could be stimulated in hPDLCs. Our data indicate that the Wnt/Ca2+ pathway plays an inhibitory role on PDLC cementogenic differentiation in inflammatory microenvironments. Therefore, targeting the Wnt/Ca2+ pathway may provide a novel therapeutic approach to improve periodontal regeneration for periodontal diseases.

Higher-resolution structure of the human insulin receptor ectodomain: Multi-modal inclusion of the insert domain

Insulin receptor (IR) signaling is critical to controlling nutrient uptake and metabolism. However, only a low-resolution (3.8 Å) structure currently exists for the IR ectodomain, with some segments ill-defined or unmodeled due to disorder. Here, we revise this structure using new diffraction data to 3.3 Å resolution that allow improved modeling of the N-linked glycans, the first and third fibronectin type III domains, and the insert domain. A novel haptic interactive molecular dynamics strategy was used to aid fitting to low-resolution electron density maps. The resulting model provides a foundation for investigation of structural transitions in IR upon ligand binding.

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.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.