Mutational analysis and NMR spectroscopy of quail cysteine and glycine-rich protein CRP2 reveal an intrinsic segmental flexibility of LIM domains

LIM domain-wide comprehensive virtual mutagenesis provides structural rationale for cardiomyopathy mutations in CSRP3

Cardiomyopathies are a severe and chronic cardiovascular burden worldwide, affecting a large cohort in the general population. Cysteine and glycine-rich protein 3 (CSRP3) is one of key proteins implicated in dominant dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM). In this study, we device a rapid in silico screening protocol that creates a mutational landscape map for all possible allowed and disallowed substitutions in the protein of interest. This map provides the structural and functional insights on the stability of LIM domains of CSRP3. Further, the sequence analysis delineates the eukaryotic CSRP3 protein orthologs which complements the mutational map, but provide limited information of amino acid exchanges. Next, we also evaluated the effect of HCM/DCM mutations on these domains. One of highly destabilising mutations—L44P (also disease causing) and a neutral mutation—L44M were further subjected to molecular dynamics (MD) simulations. The results establish that L44P substitution affects the LIM domain structure by altering secondary structure and due to loss of hydrophobic interaction with Phenylananine 35. The present study provides a useful perspective to our understanding of the role of mutations in the CSRP3 LIM domains and their evolution. This study provides a novel computational screening method for quick identification of key mutation sites for specific protein structures that can reduce the burden on experimental research.

Performance comparison of computational methods for modeling alpha-helical structures

Geometry optimization results are reported for secondary structural elements of small proteins and polypeptides. Emphasis is placed on how well molecular mechanics as well as semiempirical, ab initio, and density functional methods describe α-helical and related structures in purely theoretical models (Gly10, Ile10) as well as in realistic models (an α-helical region of calmodulin, and the complete structure of a small protein). Many of the methods examined here were found to provide unsatisfactory descriptions of the hydrogen-bonding interactions within polypeptide-type structures, as the α-helical canonical secondary structure motif was not reproduced accurately. Ab initio and DFT methods provided reasonable results only when solvation models were included, although Hartree-Fock failed even with solvation in one of the test cases; among the semiempirical methods, one of the PM6 implementations performed very well.

Four and a half LIM protein 1 gene mutations cause four distinct human myopathies: a comprehensive review of the clinical, histological and pathological features

Mutations in the four and a half LIM protein 1 (FHL1) gene were recently identified as the cause of four distinct skeletal muscle diseases. Since the initial report outlining the first fhl1 mutation in 2008, over 25 different mutations have been identified in patients with reducing body myopathy, X-linked myopathy characterized by postural muscle atrophy, scapuloperoneal myopathy and Emery-Dreifuss muscular dystrophy. Reducing body myopathy was first described four decades ago, its underlying genetic cause was unknown until the discovery of fhl1 mutations. X-linked myopathy characterized by postural muscle atrophy is a novel disease where fhl1 mutations are the only cause. This review will profile each of the FHL1, with a comprehensive analysis of mutations, a comparison of the clinical and histopathological features and will present several hypotheses for the possible disease mechanism(s).

Phenotypic heterogeneity in British patients with a founder mutation in the FHL1 gene

Mutations in the four-and-a-half LIM domain 1 (FHL1) gene, which encodes a 280-amino-acid protein containing four LIM domains and a single zinc-finger domain in the N-terminal region, have been associated with a broad clinical spectrum of X-linked muscle diseases encompassing a variety of different phenotypes. Patients might present with a scapuloperoneal myopathy, a myopathy with postural muscle atrophy and generalized hypertrophy, an Emery-Dreifuss muscular dystrophy, or an early onset myopathy with reducing bodies. It has been proposed that the phenotypic variability is related to the position of the mutation within the FHL1 gene. Here, we report on three British families with a heterogeneous clinical presentation segregating a single FHL1 gene mutation and haplotype, suggesting that this represents a founder mutation. The underlying FHL1 gene mutation was detected by direct sequencing and the founder effect was verified by haplotype analysis of the FHL1 gene locus. A 3-bp insertion mutation (p.Phe127_Thr128insIle) within the second LIM domain of the FHL1 gene was identified in all available affected family members of the three families. Haplotype analysis of the FHL1 region on Xq26 revealed that the families shared a common haplotype. The p.Phe127_Thr128insIle mutation in the FHL1 gene therefore appears to be a British founder mutation and FHL1 gene screening, in particular of exon 6, should therefore be indicated in British patients with a broad phenotypic spectrum of X-linked muscle diseases.

Structure and dynamics of the human muscle LIM protein

The family of cysteine rich proteins (CRP) comprises three closely homologous members that have been reported to interact with alpha-actinin. Muscular LIM protein (MLP/CRP3), the skeletal muscle variant, was originally discovered as a positive regulator of myogenesis and is suggested to be part of the stretch sensor of the myofibril through its interaction with telethonin (T-Cap). We determined the structure of both LIM domains of human MLP by nuclear magnetic resonance spectroscopy. We confirm by (15)N relaxation measurements that both LIM domains act as independent units and that the adjacent linker regions are fully flexible. With the published structures of CRP1 and CRP2, the complete family has now been structurally characterized.

The structural basis of integrin-linked kinase-PINCH interactions

The heterotrimeric complex between integrin-linked kinase (ILK), PINCH, and parvin is an essential signaling platform, serving as a convergence point for integrin and growth-factor signaling and regulating cell adhesion, spreading, and migration. We report a 1.6-A crystal structure of the ILK ankyrin repeat domain bound to the PINCH1 LIM1 domain, revealing the molecular basis of ILK-PINCH interactions and providing a structural description of this region of ILK. This structure identifies 5 ankyrin repeats in ILK, explains previous deletion mutagenesis data, permits identification of ILK and PINCH1 point mutations that disrupt the interaction, shows how zincs are coordinated by PINCH1 LIM1, and suggests that conformational flexibility and twisting between the 2 zinc fingers within the LIM1 domain may be important for ILK binding. These data provide an atomic-resolution description of a key interaction in the ILK-PINCH-parvin scaffolding complex.

A cysteine-rich LIM-only protein mediates regulation of smooth muscle-specific gene expression by cGMP-dependent protein kinase

Vascular smooth muscle cells (VSMCs) undergo phenotypic modulation, changing from a differentiated, contractile to a de-differentiated, synthetic phenotype; the change is associated with decreased expression of smooth muscle (SM)-specific genes and loss of cGMP-dependent protein kinase (PKG), but transfection of PKG into de-differentiated VSMCs restores SM-specific gene expression. We show that small interference RNA-mediated down-regulation or pharmacologic inhibition of PKG reduced SM-specific gene expression in differentiated VSMCs and provide a mechanism for cGMP/PKG regulation of SM-specific genes involving the cysteine-rich LIM-only protein CRP4. PKG associated with CRP4 and phosphorylated the protein in intact cells. CRP4 had no intrinsic transcriptional activity, but exhibited adaptor function, because it acted synergistically with serum response factor (SRF) and GATA6 to activate the SM-alpha-actin promoter. cGMP stimulation of the promoter required PKG and CRP4 co-expression with SRF and GATA6. A phosphorylation-deficient mutant CRP4 and a CRP4 deletion mutant deficient in PKG binding did not support cGMP/PKG stimulation of the SM-alpha-actin promoter. In the presence of wild-type but not mutant CRP4, cGMP/PKG enhanced SRF binding to a probe encoding the distal SM-alpha-actin promoter CArG (CC(AT)(6)GG) element. CRP4 and SRF associated with CArG elements of endogenous SM-specific genes in intact chromatin. Small interference RNA-mediated down-regulation of CRP4 prevented the positive effects of cGMP/PKG on SM-specific gene expression. In the presence of CRP4, cGMP/PKG increased SRF- and GATA6-dependent expression of endogenous SM-specific genes in pluripotent 10T1/2 cells. Thus, CRP4 mediates cGMP/PKG stimulation of SM-specific gene expression, and PKG plays an important role in regulating the phenotype of VSMCs.

Ribosomal cold-adaptation: characterization of the genes encoding the acidic ribosomal P0 and P2 proteins from the Antarctic ciliate Euplotes focardii

Molecular adaptation at low temperature requires specificities represented mainly by modifications in the gene sequence and consequently in the protein primary structure. To characterize the molecular mechanisms responsible for ribosome cold-adaptation, we compared the ribosomal P0 and P2 genes from the Antarctic ciliate Euplotes focardii with homologous genes from mesophilic organisms, including the ciliates Tetrahymena thermophila and non cold-adapted Euplotes species. This analysis revealed the presence of non synonymous mutations unique to E. focardii. In the P0 protein the mutations produced amino acid substitutions that increased the molecular flexibility that may facilitate a conformational adjustment associated with the interaction with the GTPase center of the large subunit rRNA, and increased the hydrophobicity of the region involved in the interaction with P1/P2 heterodimer, probably to keep associated the ribosomal stalk in the cold. In the P2 protein the mutations produced amino acid substitutions that increased the N-terminus flexibility, which may facilitate interactions with P1 protein in the formation of the heterodimer, and reduced the mobility of the C-terminus, to stabilize the stalk during ribosomal activity. Finally, P proteins appeared to be valid markers for investigating the phylogenetic origin of early eukaryotes.

Structural basis for the recognition of ldb1 by the N-terminal LIM domains of LMO2 and LMO4

LMO2 and LMO4 are members of a small family of nuclear transcriptional regulators that are important for both normal development and disease processes. LMO2 is essential for hemopoiesis and angiogenesis, and inappropriate overexpression of this protein leads to T-cell leukemias. LMO4 is developmentally regulated in the mammary gland and has been implicated in breast oncogenesis. Both proteins comprise two tandemly repeated LIM domains. LMO2 and LMO4 interact with the ubiquitous nuclear adaptor protein ldb1/NLI/CLIM2, which associates with the LIM domains of LMO and LIM homeodomain proteins via its LIM interaction domain (ldb1-LID). We report the solution structures of two LMO:ldb1 complexes (PDB: 1M3V and 1J2O) and show that ldb1-LID binds to the N-terminal LIM domain (LIM1) of LMO2 and LMO4 in an extended conformation, contributing a third strand to a beta-hairpin in LIM1 domains. These findings constitute the first molecular definition of LIM-mediated protein-protein interactions and suggest a mechanism by which ldb1 can bind a variety of LIM domains that share low sequence homology.

Structural characterization of En-1, a cold-adapted protein pheromone isolated from the Antarctic ciliate Euplotes nobilii

The second of two diffusible cell signal proteins (pheromones) purified from a wild-type strain of the Antarctic ciliate, Euplotes nobilii, has been determined by automated Edman degradation of the whole molecule and peptides generated by its chymotryptic digestion. The proposed sequence of 52 amino acids of this new pheromone, designated En-1, is: NPEDWFTPDT(10)CAYGDSNTAW(20)TTCTTPGQTC(30)YTCCSSCFDV(40)VGEQACQMSA(50)QC. In common with the previously determined 60-amino-acid sequence of the other pheromone, En-2, it bears eight cysteines in conserved positions (presumably linked into four conserved intrachain disulfide bonds), and physicochemical features of potential significance for cold adaptation, such as a reduced hydrophobicity, an increased solvent accessibility, and an improved local backbone flexibility. However, En-1 diverges from En-2 for having evolved a threonine cluster in the place of a glycine cluster to apparently make more flexible a region that is likely functionally important.

Deuterium spin probes of side-chain dynamics in proteins. 2. Spectral density mapping and identification of nanosecond time-scale side-chain motions

In the previous paper in this issue we have demonstrated that it is possible to measure the five different relaxation rates of a deuteron in (13)CH(2)D methyl groups of (13)C-labeled, fractionally deuterated proteins. The extensive set of data acquired in these experiments provides an opportunity to investigate side-chain dynamics in proteins at a level of detail that heretofore was not possible. The data, acquired on the B1 domain of peptostreptococcal protein L, include 16 (9) relaxation measurements at 4 (2) different magnetic field strengths, 25 degrees C (5 degrees C). These data are shown to be self-consistent and are analyzed using a spectral density mapping procedure which allows extraction of values of the spectral density function at a number of frequencies with no assumptions about the underlying dynamics. Dynamics data from 31 of 35 methyls in the protein for which data could be obtained were well-fitted using the two-parameter Lipari-Szabo model (Lipari, G.; Szabo, A. J. Am. Chem. Soc. 1982, 104, 4546). The data from the remaining 4 methyls can be fitted using a three-parameter version of the Lipari-Szabo model that takes into account, in a simple manner, additional nanosecond time-scale local dynamics. This interpretation is supported by analysis of a molecular dynamics trajectory where spectral density profiles calculated for side-chain methyl sites reflect the influence of slower (nanosecond) time-scale motions involving jumps between rotameric wells. A discussion of the minimum number of relaxation measurements that are necessary to extract the full complement of dynamics information is presented along with an interpretation of the extracted dynamics parameters.

A unified picture of protein hydration: prediction of hydrodynamic properties from known structures

Hydration is essential for the structural and functional integrity of globular proteins. How much hydration water is required for that integrity? A number of techniques such as X-ray diffraction, nuclear magnetic resonance (NMR) spectroscopy, calorimetry, infrared spectroscopy, and molecular dynamics (MD) simulations indicate that the hydration level is 0.3-0.5 g of water per gram of protein for medium sized proteins. Hydrodynamic properties, when accounted for by modeling proteins as ellipsoids, appear to give a wide range of hydration levels. In this paper we describe an alternative numerical technique for hydrodynamic calculations that takes account of the detailed protein structures. This is made possible by relating hydrodynamic properties (translational and rotational diffusion constants and intrinsic viscosity) to electrostatic properties (capacitance and polarizability). We show that the use of detailed protein structures in predicting hydrodynamic properties leads to hydration levels in agreement with other techniques. A unified picture of protein hydration emerges. There are preferred hydration sites around a protein surface. These sites are occupied nearly all the time, but by different water molecules at different times. Thus, though a given water molecule may have a very short residence time (approximately 100-500 ps from NMR spectroscopy and MD simulations) in a particular site, the site appears fully occupied in experiments in which time-averaged properties are measured.

JAC, a direct target of oncogenic transcription factor Jun, is involved in cell transformation and tumorigenesis

Using subtractive hybridization techniques, we have isolated a gene termed JAC that is strongly and specifically activated in avian fibroblasts transformed by the v-jun oncogene of avian sarcoma virus 17 (ASV17), but not in cells transformed by other oncogenic agents. Furthermore, JAC is highly expressed in cell lines derived from jun-induced avian fibrosarcomas. Kinetic analysis using a doxycycline-controlled conditional cell transformation system showed that expression of the 0.8-kb JAC mRNA is induced rapidly upon activation of the oncogenic v-jun allele. Nucleotide sequence analysis and transcriptional mapping revealed that the JAC gene contains two exons, with the longest ORF confined to exon 2. The deduced 68-amino acid chicken JAC protein is rich in cysteine residues and displays 37% sequence identity to mammalian high-sulfur keratin-associated proteins. The promoter region of JAC contains a consensus (5'-TGACTCA-3') and a nonconsensus (5'-TGAGTAA-3') AP-1 binding site in tandem, which are both specifically bound by the Gag-Jun hybrid protein encoded by ASV17. Mutational analysis revealed that the two AP-1 sites confer strong transcriptional activation by Gag-Jun in a synergistic manner. Ectopic expression of JAC in avian fibroblasts leads to anchorage-independent growth, strongly suggesting that deregulation of JAC is an essential event in jun-induced cell transformation and tumorigenesis.

Application of cross-correlated NMR spin relaxation to the zinc-finger protein CRP2(LIM2): evidence for collective motions in LIM domains

The solution structure of quail CRP2(LIM2) was significantly improved by using an increased number of NOE constraints obtained from a 13C,15N-labeled protein sample and by applying a recently developed triple-resonance cross-correlated relaxation experiment for the determination of the backbone dihedral angle psi. Additionally, the relative orientation of the 15N(i)-1HN(i) dipole and the 13CO(i) CSA tensor, which is related to both backbone angles phi and psi, was probed by nitrogen-carbonyl multiple-quantum relaxation and used as an additional constraint for the refinement of the local geometry of the metal-coordination sites in CRP2(LIM2). The backbone dynamics of residues located in the folded part of CRP2(LIM2) have been characterized by proton-detected 13C'(i-1)-15N(i) and 15N(i)-1HN(i) multiple-quantum relaxation, respectively. We show that regions having cross-correlated time modulation of backbone isotropic chemical shifts on the millisecond to microsecond time scale correlate with residues that are structurally altered in the mutant protein CRP2(LIM2)R122A (disruption of the CCHC zinc-finger stabilizing side-chain hydrogen bond) and that these residues are part of an extended hydrogen-bonding network connecting the two zinc-binding sites. This indicates the presence of long-range collective motions in the two zinc-binding subdomains. The conformational plasticity of the LIM domain may be of functional relevance for this important protein recognition motif.

Structure, function, and dynamics of the dimerization and DNA-binding domain of oncogenic transcription factor v-Myc

The protein product (c-Myc) of the protooncogene c-myc is a transcriptional regulator playing a key role in cellular growth, differentiation, and apoptosis. Deregulated myc genes, like the transduced retroviral v-myc allele, are oncogenic and cause cell transformation. The C-terminal bHLHZip domain of v-Myc, encompassing protein dimerization (helix-loop-helix, leucine zipper) and DNA contact (basic region) surfaces, was expressed in bacteria as a highly soluble p15(v-myc )recombinant protein. Dissociation constants (K(d)) for the heterodimer formed with the recombinant bHLHZip domain of the Myc binding partner Max (p14(max)) and for the Myc-Max-DNA complex were estimated using circular dichroism (CD) spectroscopy and quantitative electrophoretic mobility shift assay (EMSA). Multi-dimensional NMR spectroscopy was used to characterize the solution structural and dynamic properties of the v-Myc bHLHZip domain. Significant secondary chemical shifts indicate the presence of two separated alpha-helical regions. The C-terminal leucine zipper region forms a compact alpha-helix, while the N-terminal basic region exhibits conformational averaging with substantial alpha-helical content. Both helices lack stabilizing tertiary side-chain interactions and represent exceptional examples for loosely coupled secondary structural segments in a native protein. These results and CD thermal denaturation data indicate a monomeric state of the v-Myc bHLHZip domain. The (15)N relaxation data revealed backbone mobilities which corroborate the existence of a partially folded state, and suggest a "beads-on-a-string" motional behaviour of the v-Myc bHLHZip domain in solution. The preformation of alpha-helical regions was confirmed by CD thermal denaturation studies, and quantification of the entropy changes caused by the hydrophobic effect and the reduction of conformational entropy upon protein dimerization. The restricted conformational space of the v-Myc bHLHZip domain reduces the entropy penalty associated with heterodimerization and allows rapid and accurate recognition by the authentic Myc binding partner Max.

Solution structure of the focal adhesion adaptor PINCH LIM1 domain and characterization of its interaction with the integrin-linked kinase ankyrin repeat domain

PINCH is a recently identified adaptor protein that comprises an array of five LIM domains. PINCH functions through LIM-mediated protein-protein interactions that are involved in cell adhesion, growth, and differentiation. The LIM1 domain of PINCH interacts with integrin-linked kinase (ILK), thereby mediating focal adhesions via a specific integrin/ILK signaling pathway. We have solved the NMR structure of the PINCH LIM1 domain and characterized its binding to ILK. LIM1 contains two contiguous zinc fingers of the CCHC and CCCH types and adopts a global fold similar to that of functionally distinct LIM domains from cysteine-rich protein and cysteine-rich intestinal protein families with CCHC and CCCC zinc finger types. Gel-filtration and NMR experiments demonstrated a 1:1 complex between PINCH LIM1 and the ankyrin repeat domain of ILK. A chemical shift mapping experiment identified regions in PINCH LIM1 that are important for interaction with ILK. Comparison of surface features between PINCH LIM1 and other functionally different LIM domains indicated that the LIM motif might have a highly variable mode in recognizing various target proteins.