Hard α-keratin IF: A structural model lacking a head-to-tail molecular overlap but having hybrid features characteristic of both epidermal keratin and vimentin IF

In intermediate filaments (IF) both epidermal keratin and vimentin molecules have been shown to have an eight residue head to-tail overlap between the rod domains of similarly directed molecules. In the case of the epidermal keratins this region has also been shown to have particular structural/functional significance since it represents a hot-spot for mutations in the four keratinopathies characterized to date. While there is good evidence that this head-to-tail overlap is present in IF containing Type III, IV, and V chains, as well as in the epidermal keratin IF (Ib/IIb), there are no data currently available for the hard α-keratin IF (Ia/IIa). Using a variety of data derived from X-ray diffraction and crosslinking studies, as well as theoretical modeling, it is now possible to demonstrate that the overlap region is not a feature of hard α-keratin IF. Indeed, it is shown that there is a nine residue gap between consecutive parallel molecules in the IF. An explanation for this observation is presented in terms of compensating disulfide bonds that occur both within the IF, and between the IF and the matrix in which the IF are embedded. © 1995 Wiley-Liss, Inc.     

Packing and hydrophobicity effects on protein folding and stability: effects of beta-branched amino acids, valine and isoleucine, on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers.

The aim of this study was to examine the differences between hydrophobicity and packing effects in specifying the three-dimensional structure and stability of proteins when mutating hydrophobes in the hydrophobic core. In DNA-binding proteins (leucine zippers), Leu residues are conserved at positions "d," and beta-branched amino acids, Ile and Val, often occur at positions "a" in the hydrophobic core. In order to discern what effect this selective distribution of hydrophobes has on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers, three Val or three Ile residues were simultaneously substituted for Leu at either positions "a" (9, 16, and 23) or "d" (12, 19, and 26) in both chains of a model coiled coil. The stability of the resulting coiled coils was monitored by CD in the presence of Gdn.HCl. The results of the mutations of Ile to Val at either positions "a" or "d" in the reduced or oxidized coiled coils showed a significant hydrophobic effect with the additional methylene group in Ile stabilizing the coiled coil (delta delta G values range from 0.45 to 0.88 kcal/mol/mutation). The results of mutations of Leu to Ile or Val at positions "a" in the reduced or oxidized coiled coils showed a significant packing effect in stabilizing the coiled coil (delta delta G values range from 0.59 to 1.03 kcal/mol/mutation). Our results also indicate the subtle control hydrophobic packing can have not only on protein stability but on the conformation adopted by the amphipathic alpha-helices. These structural findings correlate with the observation that in DNA-binding proteins, the conserved Leu residues at positions "d" are generally less tolerant of amino acid substitutions than the hydrophobic residues at positions "a."     

The role of position a in determining the stability and oligomerization state of alpha-helical coiled coils: 20 amino acid stability coefficients in the hydrophobic core of proteins

We describe here a systematic investigation into the role of position a in the hydrophobic core of a model coiled-coil protein in determining coiled-coil stability and oligomerization state. We employed a model coiled coil that allowed the formation of an extended three-stranded trimeric oligomerization state for some of the analogs; however, due to the presence of a Cys-Gly-Gly linker, unfolding occurred from the same two-stranded monomeric oligomerization state for all of the analogs. Denaturation from a two-stranded state allowed us to measure the relative contribution of 20 different amino acid side chains to coiled-coil stability from chemical denaturation profiles. In addition, the relative hydrophobicity of the substituted amino acid side chains was assessed by reversed-phase high-performance liquid chromatography and found to correlate very highly (R = 0.95) with coiled-coil stability. We also determined the effect of position a in specifying the oligomerization state using ultracentrifugation as well as high-performance size-exclusion chromatography. We found that nine of the analogs populated one oligomerization state exclusively at peptide concentrations of 50 microM under benign buffer conditions. The Leu-, Tyr-, Gln-, and His-substituted analogs were found to be exclusively three-stranded trimers, while the Asn-, Lys-, Orn-, Arg-, and Trp-substituted analogs formed exclusively two-stranded monomers. Modeling results for the Leu-substituted analog showed that a three-stranded oligomerization state is preferred due to increased side-chain burial, while a two-stranded oligomerization state was observed for the Trp analog due to unfavorable cavity formation in the three-stranded state.     

Modeling of the metallo-beta-lactamase from B. fragilis: structural and dynamic effects of inhibitor binding

The structure and dynamics of an inhibitor-bound complex of the metallo-beta-lactamase from Bacteroides fragilis are studied by using molecular dynamics. A search of the conformational space was performed to obtain three distinct models of the complex, which were then subjected to solvated molecular dynamics. A solvated molecular dynamics study of the apo protein was performed to serve as a baseline for comparison with the bound simulations. We find loop conformation changes due to binding as well as a decrease in flexibility of the protein as a whole and especially in the major loop of the beta-lactamase. We report the structural and dynamical features of the inhibitor-bound and apo models, as well as experimentally measurable quantities, which should be capable of distinguishing the two binding modes we have determined.     

Fishing new proteins in the twilight zone of genomes: the test case of outer membrane proteins in Escherichia coli K12, Escherichia coli O157:H7, and other Gram-negative bacteria

We address the problem of clustering the whole protein content of genomes into three different categories-globular, all-alpha, and all-beta membrane proteins-with the aim of fishing new membrane proteins in the pool of nonannotated proteins (twilight zone). The focus is then mainly on outer membrane proteins. This is performed by using an integrated suite of programs (Hunter) specifically developed for predicting the occurrence of signal peptides in proteins of Gram-negative bacteria and the topography of all-alpha and all-beta membrane proteins. Hunter is tested on the well and partially annotated proteins (2160 and 760, respectively) of Escherichia coli K 12 scoring as high as 95.6% in the correct assignment of each chain to the category. Of the remaining 1253 nonannotated sequences, 1099 are predicted globular, 136 are all-alpha, and 18 are all-beta membrane proteins. In Escherichia coli 0157:H7 we filtered 1901 nonannotated proteins. Our analysis classifies 1564 globular chains, 327 inner membrane proteins, and 10 outer membrane proteins. With Hunter, new membrane proteins are added to the list of putative membrane proteins of Gram-negative bacteria. The content of outer membrane proteins per genome (nine are analyzed) ranges from 1.5% to 2.4%, and it is one order of magnitude lower than that of inner membrane proteins. The finding is particularly relevant when it is considered that this is the first large-scale analysis based on validated tools that can predict the content of outer membrane proteins in a genome and can allow cross-comparison of the same protein type between different species.     

Computational analysis of molecular basis of 1:1 interactions of NRG-1beta wild-type and variants with ErbB3 and ErbB4

The neuregulin/ErbB system is a growth factor/receptor cascade that has been proven to be essential in the development of the heart and the sympathetic nervous system. However, the basis of the specificity of ligand-receptor recognition remains to be elucidated. In this study, the structures of NRG-1beta/ErbB3 and NRG-1beta/ErbB4 complexes were modeled based on the available structures of the homologous proteins. The binding free energies of NRG-1beta to ErbB3 and ErbB4 were calculated using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) computational method. In addition, computational alanine-scanning mutagenesis was performed in the binding site of NRG-1beta and the difference in the binding free energies between NRG-1beta mutants and the receptors was calculated. The results specify the contribution of each residue at the interaction interfaces to the binding affinity of NRG-1beta with ErbB3 and ErbB4, identifying several important interaction residue pairs that are in agreement with previously acquired experimental data. This indicates that the presented structural models of NRG-1beta/ErbB3 and NRG-1beta/ErbB4 complexes are reliable and could be used to guide future studies, such as performing desirable mutations on NRG-1beta to increase the binding affinity and selectivity to the receptor and discovering new therapeutic agents for the treatment of heart failure.     

Synthesis, folding, and structure of the beta-turn mimic modified B1 domain of streptococcal protein G

The mechanism of beta-sheet formation remains a fundamental issue in our understanding of the protein folding process, but is hampered by the often encountered kinetic competition between folding and aggregation. The role of local versus nonlocal interactions has been probed traditionally by mutagenesis of both turn and strand residues. Recently, rigid organic molecules that impose a correct chain reversal have been introduced in several small peptides to isolate the importance of the long-range interactions. Here, we present the incorporation of a well-studied beta-turn mimic, designated as the dibenzofuran-based (DBF) amino acid, in the B1 domain of streptococcal protein G (B1G), and compare our results with those obtained upon insertion of the same mimic into the N-terminal beta-hairpin of B1G (O Melnyk et al., 1998, Lett Pept Sci 5:147-150). The DBF-B1G domain conserves the structure and the functional and thermodynamical properties of the native protein, whereas the modified peptide does not adopt a native-like conformation. The nature of the DBF flanking residues in the modified B1G domain prevents the beta-turn mimic from acting as a strong beta-sheet nucleator, which reinforces the idea that the native beta-hairpin formation is not driven by the beta-turn formation, but by tertiary interactions.     

Inhibition of protein-protein interactions: the discovery of druglike beta-catenin inhibitors by combining virtual and biophysical screening

The interaction between beta-catenin and Tcf family members is crucial for the Wnt signal transduction pathway, which is commonly mutated in cancer. This interaction extends over a very large surface area (4800 A(2)), and inhibiting such interactions using low molecular weight inhibitors is a challenge. However, protein surfaces frequently contain "hot spots," small patches that are the main mediators of binding affinity. By making tight interactions with a hot spot, a small molecule can compete with a protein. The Tcf3/Tcf4-binding surface on beta-catenin contains a well-defined hot spot around residues K435 and R469. A 17,700 compounds subset of the Pharmacia corporate collection was docked to this hot spot with the QXP program; 22 of the best scoring compounds were put into a biophysical (NMR and ITC) screening funnel, where specific binding to beta-catenin, competition with Tcf4 and finally binding constants were determined. This process led to the discovery of three druglike, low molecular weight Tcf4-competitive compounds with the tightest binder having a K(D) of 450 nM. Our approach can be used in several situations (e.g., when selecting compounds from external collections, when no biochemical functional assay is available, or when no HTS is envisioned), and it may be generally applicable to the identification of inhibitors of protein-protein interactions.     

Molecular dynamics studies of a DNA-binding protein: 1. A comparison of the trp repressor and trp aporepressor aqueous simulations.

The results of two 30-ps molecular dynamics simulations of the trp repressor and trp aporepressor proteins are presented in this paper. The simulations were obtained using the AMBER molecular mechanical force field and in both simulations a 6-A shell of TIP3P waters surrounded the proteins. The trp repressor protein is a DNA-binding regulatory protein and it utilizes a helix-turn-helix (D helix-turn-E helix) motif to interact with DNA. The trp aporepressor, lacking two molecules of the L-tryptophan corepressor, cannot bind specifically to DNA. Our simulations show that the N- and C-termini and the residues in and near the helix-turn-helix motifs are the most mobile regions of the proteins, in agreement with the X-ray crystallographic studies. Our simulations also find increased mobility of the residues in the turn-D helix-turn regions of the proteins. We find the average distance separating the DNA-binding motifs to be larger in the repressor as compared to the aporepressor. In addition to examining the protein residue fluctuations and deviations with respect to X-ray structures, we have also focused on backbone dihedral angles and corepressor hydrogen-bonding patterns in this paper.     

MAPKAP kinase 2 is activated by heat shock and TNF-α: In vivo phosphorylation of small heat shock protein results from stimulation of the MAP kinase cascade

The activation of MAPKAP kinase 2 was investigated under heat-shock conditions in mouse Ehrlich ascites tumor cells and after treatment of human MO7 cells with tumor necrosis factor-α (TNF-α). MAPKAP kinase 2 activity was determined using the small heat-shock proteins (sHsps) Hsp25 and Hsp27 as substrates. In both cell types, about a threefold increase in MAPKAP kinase 2 activity could be detected in a time interval of about 10–15 min after stimulation either by heat shock or TNF-α. Phosphorylation of MAPKAP kinase 2, but not the level of MAPKAP kinase 2 mRNA, was increased after heat shock in EAT cells. It is further shown that activation of MAPKAP kinase 2 in MO7 cells is accompanied by increased MAP kinase activity. These data strongly suggest that increased phosphorylation of the sHsps after heat shock or TNF-α treatment results from phosphorylation by MAPKAP kinase 2, which itself is activated by phosphorylation through MAP kinases. Hence, we demonstrate that MAPKAP kinase 2 is responsible not only for phosphorylation of sHsps in vitro but also in vivo. The findings link sHsp phosphorylation to the MAP kinase cascade, explaining the early phosphorylation of sHsp that is stimulated by a variety of inducers such as mitogens, phorbol esters, thrombin, calcium ionophores, and heat shock.     

Host‐mediated induction of α‐amylases by larvae of the Mexican bean weevil Zabrotes subfasciatus (Coleoptera: Chrysomelidae: Bruchinae) is irreversible and observed from the initiation of the feeding period

Larvae of Zabrotes subfasciatus secrete α‐amylases that are insensitive to the α‐amylase inhibitor found in seeds of Phaseolus vulgaris. By analyzing amylase activities during larval development on P. vulgaris, we detected activity of the constitutive amylase and the two inducible amylase isoforms at all stages. When larvae were transferred from the non α‐amylase inhibitor containing seeds of Vigna unguiculata to P. vulgaris, the inducible α‐amylases were expressed at the same level as in control larvae fed on P. vulgaris. Interestingly, when larvae were transferred from seeds of P. vulgaris to those of V. unguiculata, inducible α‐amylases continued to be expressed at a level similar to that found in control larvae fed P. vulgaris continuously. When 10‐day‐old larvae were removed from seeds of V. unguiculata and transferred into capsules containing flour of P. vulgaris cotyledons, and thus maintained until completing 17 days (age when the larvae stopped feeding), we could detect higher activity of the inducible α‐amylases. However, when larvae of the same age were transferred from P. vulgaris into capsules containing flour of V. unguiculata, the inducible α‐amylases remained up‐regulated. These results suggest that the larvae of Z. subfasciatus have the ability to induce insensitive amylases early in their development. A short period of feeding on P. vulgaris cotyledon flour was sufficient to irreversibly induce the inducible α‐amylase isoforms. Incubations of brush border membrane vesicles with the α‐amylase inhibitor 1 from P. vulgaris suggest that the inhibitor is recognized by putative receptors found in the midgut microvillar membranes. © 2010 Wiley Periodicals, Inc.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long‐range hydrophobic interactions in hairpin formation

Two peptides, corresponding to the turn region of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51–56 [IG(51–56)] and 50–57 [IG(50–57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the β‐turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283–323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the β‐hairpin structure of peptides corresponding to the C‐terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long‐range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the β‐hairpin structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Independent tyrosyl contributions to the CD of Ff gene 5 protein and the distinctive effects of Y41H and Y41F mutants on protein-protein cooperative interactions

The gene 5 protein (g5p) of the Ff virus contains five Tyr, individual mutants of which have now all been characterized by CD spectroscopy. The protein has a dominant tyrosyl 229-nm L(a) CD band that is shown to be approximately the sum of the five individual Tyr contributions. Tyr41 is particularly important in contributing to the high cooperativity with which the g5p binds to ssDNA, and Y41F and Y41H mutants are known to differ in dimer-dimer packing interactions in crystal structures. We compared the solution structures and binding properties of the Y41F and Y41H mutants using CD spectroscopy. Secondary structures of the mutants were similar by CD analyses and close to those derived from the crystal structures. However, there were significant differences in the binding properties of the two mutant proteins. The Y41H protein had an especially low binding affinity and perturbed the spectrum of poly[d(A)] in 2 mM Na(+) much less than did Y41F and the wild-type gene 5 proteins. Moreover, a change in the Tyr 229 nm band, assigned to the perturbation of Tyr34 at the dimer-dimer interface, was absent in titrations with the Y41H mutant under low salt conditions. In contrast, titrations with the Y41H mutant in 50 mM Na(+) exhibited typical CD changes of both the nucleic acid and the Tyr 229-nm band. Thus, protein-protein and g5p-ssDNA interactions appeared to be mutually influenced by ionic strength, indicative of correlated changes in the ssDNA binding and cooperativity loops of the protein or of indirect structural constraints.     

Topology and sequence in the folding of a TIM barrel protein: global analysis highlights partitioning between transient off-pathway and stable on-pathway folding intermediates in the complex folding mechanism of a (betaalpha)8 barrel of unknown function from B. subtilis

The relative contributions of chain topology and amino acid sequence in directing the folding of a (betaalpha)(8) TIM barrel protein of unknown function encoded by the Bacillus subtilis iolI gene (IOLI) were assessed by reversible urea denaturation and a combination of circular dichroism, fluorescence and time-resolved fluorescence anisotropy spectroscopy. The equilibrium reaction for IOLI involves, in addition to the native and unfolded species, a stable intermediate with significant secondary structure and stability and self-associated forms of both the native and intermediate states. Global kinetic analysis revealed that the unfolded state partitions between an off-pathway refolding intermediate and the on-pathway equilibrium intermediate early in folding. Comparisons with the folding mechanisms of two other TIM barrel proteins, indole-3-glycerol phosphate synthase from the thermophile Sulfolobus solfataricus (sIGPS) and the alpha subunit of Escherichia coli tryptophan synthase (alphaTS), reveal striking similarities that argue for a dominant role of the topology in both early and late events in folding. Sequence-specific effects are apparent in the magnitudes of the relaxation times and relative stabilities, in the presence of additional monomeric folding intermediates for alphaTS and sIGPS and in rate-limiting proline isomerization reactions for alphaTS.