Stabilization of the long central helix of troponin C by intrahelical salt bridges between charged amino acid side chains

Functional Analysis of Novicidin Peptide: Coordinated Delivery System for Zinc via Schiff Base Ligand

Novicidin (NVC), is a membrane-penetrating peptide, which forms a stable complex with Zn-Schiff base with interesting antitumor selectivity. We studied NVC derivatives to determine functional roles of key amino acids in toxicity, helicity, and binding of the Zn-Schiff base complex. Trimmed derivatives highlighted the role of peptide length and helicity in toxicity and membrane penetration. The removal of Lys from position 1 and 2 strongly increases the ability to disrupt the membranes. The trimming of the N-terminal residues significantly increases the stability of peptide helicity enhancing penetrating properties. Gly residue derivatives undermined a role of peptide bending in membrane penetration and toxicity. After the substitution of the central Gly derivatives with Ile or Lys, the peptides retained toxicity. These results illustrate the minor role of central helix bending in NVC toxicity. Binding-site-peptide derivatives identified His residue as the sole Zn-Schiff base binding site and eliminated the role of other aromatic residues.

Molecular and functional consequences of mutations in the central helix of cardiac troponin C

The objective of this work was to investigate the role of acidic residues within the exposed middle segment of the central helix of cTnC in (1) cTnC-cTnI interactions, (2) Ca(2+) binding and exchange with the regulatory N-domain of cTnC in increasingly complex biochemical systems, and (3) ability of the cTn complex to regulate actomyosin ATPase. In order to achieve this objective, we introduced the D87A/D88A and E94A/E95A/E96A mutations into the central helix of cTnC. The D87A/D88A and E94A/E95A/E96A mutations decreased affinity of cTnC for the regulatory region of cTnI. The Ca(2+) sensitivity of the regulatory N-domain of isolated cTnC was decreased by the D87A/D88A, but not E94A/E95A/E96A mutation. However, both the D87A/D88A and E94A/E95A/E96A mutations desensitized the cTn complex and reconstituted thin filaments to Ca(2+). Decreases in the Ca(2+) sensitivity of the cTn complex and reconstituted thin filaments were, at least in part, due to faster rates of Ca(2+) dissociation. In addition, the D87A/D88A and E94A/E95A/E96A mutations desensitized actomyosin ATPase to Ca(2+), and decreased maximal actomyosin ATPase activity. Thus, our results indicate that conserved acidic residues within the exposed middle segment of the central helix of cTnC are important for the proper regulatory function of the cTn complex.

Dynamics of the rhomboid-like protein RHBDD2 expression in mouse retina and involvement of its human ortholog in retinitis pigmentosa

The novel rhomboid-like protein RHBDD2 is distantly related to rhomboid proteins, a group of highly specialized membrane-bound proteases that catalyze regulated intramembrane proteolysis. In retina, RHBDD2 is expressed from embryonic stages to adulthood, and its levels show age-dependent changes. RHBDD2 is distinctly abundant in the perinuclear region of cells, and it localizes to their Golgi. A glycine zipper motif present in one of the transmembrane domains of RHBDD2 is important for its packing into the Golgi membranes. Its deletion causes dislodgment of RHBDD2 from the Golgi. A specific antibody against RHBDD2 recognizes two forms of the protein, one with low (39 kDa; RHBDD2(L)) and the other with high (117 kDa; RHBDD2H) molecular masses in mouse retinal extracts. RHBDD2(L) seems to be ubiquitously expressed in all retinal cells. In contrast, RHBDD2H seems to be present only in the outer segments of cone photoreceptors and may correspond to a homotrimer of RHBDD2(L). This protein consistently co-localizes with S- and M-types of cone opsins. We identified a homozygous mutation in the human RHBDD2 gene, R85H, that co-segregates with disease in affected members of a family with autosomal recessive retinitis pigmentosa. Our findings suggest that the RHBDD2 protein plays important roles in the development and normal function of the retina.

Unlocking the mystery behind the activation phenomenon of T1 lipase: a molecular dynamics simulations approach

The activation of lipases has been postulated to proceed by interfacial activation, temperature switch activation, or aqueous activation. Recently, based on molecular dynamics (MD) simulation experiments, the T1 lipase activation mechanism was proposed to involve aqueous activation in addition to a double-flap mechanism. Because the open conformation structure is still unavailable, it is difficult to validate the proposed theory unambiguously to understand the behavior of the enzyme. In this study, we try to validate the previous reports and uncover the mystery behind the activation process using structural analysis and MD simulations. To investigate the effects of temperature and environmental conditions on the activation process, MD simulations in different solvent environments (water and water-octane interface) and temperatures (20, 50, 70, 80, and 100°C) were performed. Based on the structural analysis of the lipases in the same family of T1 lipase (I.5 lipase family), we proposed that the lid domain comprises α6 and α7 helices connected by a loop, thus forming a helix-loop-helix motif involved in interfacial activation. Throughout the MD simulations experiments, lid displacements were only observed in the water-octane interface, not in the aqueous environment with respect to the temperature effect, suggesting that the activation process is governed by interfacial activation coupled with temperature switch activation. Examining the activation process in detail revealed that the large structural rearrangement of the lid domain was caused by the interaction between the hydrophobic residues of the lid with octane, a nonpolar solvent, and this conformation was found to be thermodynamically favorable.

Importance of the spatial display of charged residues in heparin-peptide interactions

Many studies have examined consensus sequences required for protein-glycosaminoglycan interactions. Through the synthesis of helical heparin binding peptides, this study probes the relationship between spatial arrangement of positive charge and heparin binding affinity. Peptides with a linear distribution of positive charge along one face of the alpha-helix had the highest affinity for heparin. Moving the basic residues away from a single face resulted in drastic changes in heparin binding affinity of up to three orders of magnitude. These findings demonstrate that amino acid sequences, different from the known heparin binding consensus sequences, will form high affinity protein-heparin binding interactions when the charged residues are aligned linearly.

Detecting coevolution without phylogenetic trees? Tree-ignorant metrics of coevolution perform as well as tree-aware metrics

BACKGROUND

Identifying coevolving positions in protein sequences has myriad applications, ranging from understanding and predicting the structure of single molecules to generating proteome-wide predictions of interactions. Algorithms for detecting coevolving positions can be classified into two categories: tree-aware, which incorporate knowledge of phylogeny, and tree-ignorant, which do not. Tree-ignorant methods are frequently orders of magnitude faster, but are widely held to be insufficiently accurate because of a confounding of shared ancestry with coevolution. We conjectured that by using a null distribution that appropriately controls for the shared-ancestry signal, tree-ignorant methods would exhibit equivalent statistical power to tree-aware methods. Using a novel t-test transformation of coevolution metrics, we systematically compared four tree-aware and five tree-ignorant coevolution algorithms, applying them to myoglobin and myosin. We further considered the influence of sequence recoding using reduced-state amino acid alphabets, a common tactic employed in coevolutionary analyses to improve both statistical and computational performance.

RESULTS

Consistent with our conjecture, the transformed tree-ignorant metrics (particularly Mutual Information) often outperformed the tree-aware metrics. Our examination of the effect of recoding suggested that charge-based alphabets were generally superior for identifying the stabilizing interactions in alpha helices. Performance was not always improved by recoding however, indicating that the choice of alphabet is critical.

CONCLUSION

The results suggest that t-test transformation of tree-ignorant metrics can be sufficient to control for patterns arising from shared ancestry.

Centrins in retinal photoreceptor cells: regulators in the connecting cilium

Changes in the intracellular Ca2+ concentration regulate the visual signal transduction cascade directly or more often indirectly through Ca2+-binding proteins. Here we focus on centrins, which are members of a highly conserved subgroup of the EF-hand superfamily of Ca2+-binding proteins in photoreceptor cells of the vertebrate retina. Centrins are commonly associated with centrosome-related structures. In mammalian retinal photoreceptor cells, four centrin isoforms are expressed as prominent components in the connecting cilium linking the light-sensitive outer segment compartment with the metabolically active inner segment compartment. Our data indicate that Ca2+-activated centrin isoforms assemble into protein complexes with the visual heterotrimeric G-protein transducin. This interaction of centrins with transducin is mediated by binding to the betagamma-dimer of the heterotrimeric G-protein. More recent findings show that these interactions of centrins with transducin are reciprocally regulated via site-specific phosphorylations mediated by the protein kinase CK2. The assembly of centrin/G-protein complexes is a novel aspect of translocation regulation of signalling proteins in sensory cells, and represents a potential link between molecular trafficking and signal transduction in general.

Sequence comparisons of intermediate filament chains: evidence of a unique functional/structural role for coiled-coil segment 1A and linker L1

A comprehensive analysis of the sequences of all types of intermediate filament chains has been undertaken with a particular emphasis on those of segment 1A and linker L1. This has been done to assess whether structural characteristics can be recognized in the sequences that would be consistent with the role of each region in the recently proposed "swinging head" hypothesis. The analyses show that linker L1 is the most flexible rod domain region, that it is the most elongated structure (on a per residue basis), and that it is the most variable region as regards sequence and length. Segment 1A has one of the two most highly conserved regions of sequence in the rod domain (the other being at the end of segment 2B), with seven particular residues conserved across all chain types. It also contains one of the very few potential interchain ionic interactions that could be conserved across all chain types. However, the aggregation of chains in segment 1A is specified less precisely overall by interchain ionic interactions than are the other coiled-coil segments. The apolar residue contents in positions a and d of the heptad substructure are the highest of any coiled-coil segment in the intermediate filament family. Segment 1A also displays an amino acid composition atypical of not only coiled-coil segments 1B and 2B, but indeed of two-stranded coiled coils in general. Nonetheless, molecular modeling based on the crystal structure of the monomeric 1A fragment from human vimentin shows that coiled-coil formation is plausible. The most extensive regions of apolar/aromatic residues lie at the C-terminal end of segment 2B in the helix termination motif and in segment 1A in and close to the helix initiation motif. The predicted stability of the individual alpha-helices in segment 1A is greater than in those comprising segments 1B and 2B, though potential intrachain ionic interactions are either lacking or are minimal in number. Analysis of the 1A sequence and those regions immediately N- and C-terminal to it has shown that the capping residues are near optimal close to the previously predicted ends, thus adding to the likely stability of the alpha-helical structure. However, a second terminating sequence is predicted in 1A (about 10 residues back from the C-terminus). This allows the possibility of some unwinding of the alpha-helical structure of 1A immediately adjacent to linker L1 when the head domains no longer stabilize the coiled-coil structure. All of these data are consistent with the concept of a flexible hinge at L1 and with the ability of the two alpha-helical coiled-coil strands to separate under appropriate conditions and partly unwind at their C-terminal ends to allow the head domains a greater degree of mobility, thus facilitating function.

Crystal structure of methionine aminopeptidase from hyperthermophile, Pyrococcus furiosus

The structure of methionine aminopeptidase from hyperthermophile Pyrococcus furiosus (PfMAP) with an optimal growth temperature of 100 degreesC was determined by the multiple isomorphous replacement method and refined in three different crystal forms, one monoclinic and two hexagonal, at resolutions of 2.8, 2.9, and 3.5 A. The resolution of the monoclinic crystal form was extended to 1.75 A by water-mediated transformation to a low-humidity form, and the obtained diffraction data used for high-resolution structure refinement. This is the first description of a eukaryotic type methionine aminopeptidase structure. The PfMAP molecule is composed of two domains, a catalytic domain and an insertion domain, connected via two antiparallel beta-strands. The catalytic domain, which possesses an internal 2-fold symmetry and contains two cobalt ions in the active site, resembles the structure of a prokaryotic type MAP from Escherichia coli (EcMAP), while the structure of the insertion domain containing three helices has a novel fold and accounts for a major difference between the eukaryotic and prokaryotic types of methionine aminopeptidase. Analysis of the PfMAP structure in comparison with EcMAP and other mesophile proteins reveals several factors which may contribute to the hyperthermostability of PfMAP: (1) a significantly high number of hydrogen bonds and ion-pairs between side-chains of oppositely charged residues involved in the stabilization of helices; (2) an increased number of hydrogen bonds between the positively charged side-chain and neutral oxygen; (3) a larger number of buried water molecules involved in crosslinking the backbone atoms of sequentially separate segments; (4) stabilization of two antiparallel beta-strands connecting the two domains of the molecule by proline residues; (5) shortening of N and C-terminal tails and stabilization of the loop c3E by deletion of three residues.

A single amino acid can switch the oligomerization state of the alpha-helical coiled-coil domain of cartilage matrix protein

We have studied the oligomerization of an alpha-helical coiled-coil using as an example a peptide corresponding to the C-terminal domain of cartilage matrix protein. By replacing one arginine residue, which forms an interchain ionic interaction with a glutamic acid residue, with glutamine, we found that this peptide assembles into a homotetramer at neutral pH in contrast to the native molecule which forms homotrimers. At acidic and basic pH, however, we again observed the trimer conformation. Another arginine, which is probably involved in an intrachain salt bridge, has no effect on the assembly. Our data demonstrate that besides the specific distribution of hydrophobic residues, interchain ionic interactions can be crucial in modulating the association behavior of alpha-helical coiled-coil domains.

Structural and sequence characteristics of long alpha helices in globular proteins

Elucidation of the detailed structural features and sequence requirements for alpha helices of various lengths could be very important in understanding secondary structure formation in proteins and, hence, in the protein folding mechanism. An algorithm to characterize the geometry of an alpha helix from its C(alpha) coordinates has been developed and used to analyze the structures of long alpha helices (number of residues > or = 25) found in globular proteins, the crystal structure coordinates of which are available from the Brookhaven Protein Data Bank. All long alpha helices can be unambiguously characterized as belonging to one of three classes: linear, curved, or kinked, with a majority being curved. Analysis of the sequences of these helices reveals that the long alpha helices have unique sequence characteristics that distinguish them from the short alpha helices in globular proteins. The distribution and statistical propensities of individual amino acids to occur in long alpha helices are different from those found in short alpha helices, with amino acids having longer side chains and/or having a greater number of functional groups occurring more frequently in these helices. The sequences of the long alpha helices can be correlated with their gross structural features, i.e., whether they are curved, linear, or kinked, and in case of the curved helices, with their curvature.

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In alpha-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of alpha-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 angstrom. The sch:mch interactions involving ionizable side chains that belong either to beta-strands or to the central part of alpha-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.

Conformational properties of the Arg-Leu-Gly tripeptide--DMSO--water clusters with the combined use of molecular dynamics and energy minimization studies

The Arg-Leu-Gly tripeptide is the repeating fragment of sequential arginine-rich polypeptides capable of interacting with DNA. The conformational influence of solvent molecules (DMSO/H2O) were investigated with the combined use of molecular dynamics and energy minimization. It was found that water molecules greatly contribute to the peptide structure by solvating all its hydrophylic sites even in the presence of DMSO excess, whereas one water molecule links the ammonium and carboxylic ends of the Arg-Leu-Gly. The persistence of residual water, which was confirmed by varying the computer simulation parameters, indicates that pretreatment of peptide segments in aqueous solutions should greatly affect their conformational properties in organic media. A satisfactory agreement between experimental data (1H-NMR and IR spectroscopy) and the presented computational study deserves also to be noted.

Dissecting titin into its structural motifs: identification of an alpha-helix motif near the titin N-terminus

Titin, also known as connectin, is a giant modular protein specifically found in vertebrate striated muscle. Since the huge size of titin does not allow a direct structure determination, we have started a long-term project to characterize the protein by cutting it into smaller domains or structural units. The major part of the titin sequence is assembled by modules approximately 100 amino acids long that belong to two major protein superfamilies. Most of these modules are linked together by stretches of variable length with unique sequence. No direct structural characterization has been achieved so far for any of these linkers. We present here a study of a stretch located in the titin N-terminus and part of a linker between two modules. Our attention was drawn toward this region because it shows 100% probability to form a coiled coil when analyzed by a prediction program. A synthetic 38 amino acid peptide spanning such a sequence was studied in aqueous solution by circular dichroism, nuclear magnetic resonance, and analytical ultracentrifugation at various pH, salt, and peptide concentrations. Under all conditions, it shows a strong tendency to form alpha-helical structures. In the presence of salt, this conformation is associated with the formation of helical bundles below pH 5. Above pH 5, any aggregate breaks, and the titin peptide is a monomeric helix in equilibrium with its random coil conformation. We discuss the factors which stabilize the helical conformation and the possible role of this stretch in vivo.

Isolation of suppressors of temperature-sensitive folding mutations

Mutations in the tailspike gene (gene 9) of Salmonella typhimurium phage P22 have been used to identify amino acid interactions during the folding of a polypeptide chain. Since temperature-sensitive folding (tsf) mutations cause folding defects in the P22 tailspike polypeptide chain, it is likely that mutants derived from these and correcting the original tsf defects (second-site intragenic suppressors) identify interactions during the folding pathway. We report the isolation and identification of second-site revertants to tsf mutants.

Carbon-13 NMR studies of the lysine side chains of calmodulin and its proteolytic fragments

The pH-titration and dynamic behaviour of the seven lysine side chains in bovine calmodulin were studied by carbon-13 NMR. The amino groups of the calcium saturated protein and its proteolytic fragments TR1C (1-75) and TR2C (78-148) were dimethylated with carbon-13 labeled formaldehyde; this modification did not alter the protein's structure or its ability to activate the enzyme cyclic nucleotide phosphodiesterase. Tentative assignments for 5 out of the 7 dimethyl lysine resonances could be obtained by comparing spectra of the fully and partially modified protein, with those of the proteolytic fragments. The pKa values measured for calcium saturated calmodulin ranged between 9.5 (Lys 75) and 10.2 (Lys 13); two residues (Lys 94 and Lys 13) showed a biphasic titration curve suggesting their possible involvement in ion-pairs. The dynamic behavior of the lysine side chains was deduced from spin lattice relaxation measurements. All side chains were flexible and this was not influenced by the removal of calcium, or the addition of the calmodulin antagonist trifluoperazine. The latter data suggest that the lysine side chains are not directly involved in calmodulin's target binding sites.

Molecular dynamics simulations of a winter flounder "antifreeze" polypeptide in aqueous solution

A winter flounder antifreeze polypeptide (HPLC-6) has been studied in vacuo and in aqueous solution using molecular dynamics computer simulation techniques. The helical conformation of this polypeptide was found to be stable both in vacuum and in solution. The major stabilizing interactions were found to be the main-chain hydrogen bonds, a salt-bridge interaction, and solute-solvent hydrogen bonds. A significant bending in the middle of the polypeptide chain was observed both in vacuo and in solvent at 300 K. Possible causes of the bending are discussed. From simulations of mutant polypeptide molecules in vacuo, it is concluded that the bend in the native polypeptide was caused by side chain to backbone hydrogen bond competition involving the Thr 24 side chain and facilitated by strains on the helix resulting from the Lys 18-Glu 22 salt bridge.

Design and synthesis of a helix heparin-binding peptide

Elaboration of heparin-protein-binding interactions is necessary to understand how heparin modulates protein function. The heparin-binding domain of some proteins is postulated to be a helix structure which presents a surface of high positive charge density. Thus, a synthetic 19-residue peptide designed to be alpha-helical in character was synthesized, and its interaction with heparin was studied. The peptide was shown to be 75% helix by circular dichroism (CD) spectrometry in neutral pH buffer (at 2 degrees C); helicity increased to nearly 85% under high ionic strength conditions or to nearly 100% in 75% ethanol. Increasing the temperature of the solution caused a change in the spectral envelope consistent with a coil transition of the peptide. The midpoint of the transition (i.e., the temperature at which the helix content was determined to be 50%) was 25 degrees C, and the determined van't Hoff enthalpy change (delta HvH) was 3.2 kcal/mol of peptide. By CD, heparin increases the helix content of the peptide to 100% and increases the apparent thermal stability of the peptide by about 1 kcal/mol. The melting point for the helix/coil transition of the heparin-peptide complex was 50 degrees C. The thermal coefficient of the transition (approximately 300 deg.cm2.dmol-1.degree C-1) was essentially the same for the peptide alone or the peptide-heparin complex. Dissociation of the complex under high ionic strength conditions was also observed in the CD experiment. Biological assays showed less heparin-binding activity than expected (micromolar KD values), but this was attributed to the absence of critical lysyl residues in the peptide.(ABSTRACT TRUNCATED AT 250 WORDS)

The human 64-kDa polyadenylylation factor contains a ribonucleoprotein-type RNA binding domain and unusual auxiliary motifs

Cleavage stimulation factor is one of the multiple factors required for 3'-end cleavage of mammalian pre-mRNAs. We have shown previously that this factor is composed of three subunits with estimated molecular masses of 77, 64, and 50 kDa and that the 64-kDa subunit can be UV-crosslinked to RNA in a polyadenylylation signal (AAUAAA)-dependent manner. We have now isolated cDNAs encoding the 64-kDa subunit of human cleavage stimulation factor. The 64-kDa subunit contains a ribonucleoprotein-type RNA binding domain in the N-terminal region and a repeat structure in the C-terminal region in which a pentapeptide sequence (consensus MEARA/G) is repeated 12 times and the formation of a long alpha-helix stabilized by salt bridges is predicted. An approximately 270-amino acid segment surrounding this repeat structure is highly enriched in proline and glycine residues (approximately 20% for each). When cloned 64-kDa subunit was expressed in Escherichia coli, an N-terminal fragment containing the RNA binding domain bound to RNAs in a polyadenylylation-signal-independent manner, suggesting that the RNA binding domain is directly involved in the binding of the 64-kDa subunit to pre-mRNAs.

The molecular anatomy of caldesmon

Images

A long helix from the central region of smooth muscle caldesmon

The central region of smooth muscle caldesmon is predicted to form alpha-helices on the basis of its primary structure. We have isolated a fragment (CT54) that contains this region. The hydrodynamic properties and the electron microscopic images suggest that CT54 is an elongated (35 nm), monomeric molecule. The circular dichroic spectrum yields an overall alpha-helical content of 55-58%. These results are consistent with the model that the middle portion of CT54 forms a long stretch of single-stranded alpha-helix. Such a structure, if it in fact exists, is thought to be stabilized by numerous salt bridges between charged residues at positions i and i + 4. The structural characteristics of this fragment not only represent an unusual protein configuration but also provide information about the functional role of caldesmon in smooth muscle contraction.

Vibrational modes of hemoglobin in red blood cells

Equine red blood cells were washed in saline heavy water (2H2O) to exchange the hydrogen atoms of the non-hemoglobin components with deuterons. This led to novel neutron scattering measurements of protein vibrations within a cellular system and permitted a comparison with inelastic neutron scattering measurements on purified horse hemoglobin, either dry or wetted with 2H2O. As a function of wavevector transfer Q and the frequency transfer v the neutron response typified by the dynamic structure factor S(Q, v) was found to be similar for extracted and cellular hemoglobin at low and high temperatures. At 77 K, in the cells, a peak in S(Q, v) due to the protein was found near 0.7 THz, approximately half the frequency of a strong peak in the aqueous medium. Measurements at higher temperatures (170 and 230 K) indicated similar small shifts downwards in the peak frequencies of both components. At 260 K the low frequency component became predominantly quasielastic, but a significant inelastic component could still be ascribed to the aqueous scattering. Near 295 K the frequency responses of both components were similar and centered near zero. When scattering due to water is taken into account it appears that the protein neutron response in, or out of, red blood cells is little affected by hydration in the low frequency regime where Van der Waals forces are thought to be effective.

Design and synthesis of the pseudo-EF hand in calbindin D9K: effect of amino acid substitutions in the alpha-helical regions

A series of 37-residue analogues of the pseudo-EF hand in bovine calbindin D9K has been synthesized by the solid phase method. In the presence of calcium an alpha-helical induction of up to 44% was observed for the peptide with the native sequence with a Kd for calcium binding of 0.35 mM. A number of amino acid substitutions have been carried out to study the packing of the two alpha-helices based on the crystal structure of the entire protein. Three strategies were employed: (1) replacement of the Leu residues, which in the crystal structure do not contribute to the hydrophobic interaction between the two helices, by Gln or Ala in order to control the orientation of the helix packing, (2) stabilization of the individual helix by introducing a Glu-...Lys+ salt bridge or by changing the N-terminal charge to compensate for the helix dipole moment, and (3) introduction of a disulfide bond between the two helices to help the packing of the helices. The mutants with the substitution of (Leu-30, Leu-32) to (Gln-30, Gln-32), (Gln-30, Ala-32), and (Ala-30,Ala-32) designed based on the strategy 1 do not show any affinity for calcium and have low alpha-helicity. The Leu-30 to Lys-30 mutant designed to form a salt bridge between the side chains of Glu-26 and Lys-30 has an apparent Kd for calcium of 6.8 mM. Kd of the N-terminal acetylated and succinylated mutants are 0.41 and 0.45 mM, respectively, and no increase in the alpha-helix content relative to that of the natural sequence peptide is observed. The disulfide containing mutants, namely Tyr-13, Leu-31 to Cys-13, Cys-31 and Tyr-13, Leu-31 to Cys-13, hCys-31, show apparent Kd values of 0.93 and 2.1 mM, respectively. The former mutant shows the highest alpha-helix content among the peptides studied in the presence and absence of calcium. While it is difficult to construct an isolated and rigid helix-loop-helix motif with peptides of this size, introduction of a disulfide bond proved to be effective for this purpose.

The Glu 2- ... Arg 10+ side-chain interaction in the C-peptide helix of ribonuclease A

Previous studies have identified Lys 1, Glu 2, and His 12 as the charged residues responsible for the pH-dependent stability of the helix formed by the isolated C-peptide (residues 1-13 of ribonuclease A). Here we examine whether the helix-stabilizing behavior of Glu 2- results from a Glu 2- ... Arg 10+ interaction, which is known to be present in the crystal structure of ribonuclease A. The general approach is to measure the helix content of C-peptide analogs as a function of three variables: pH (titration of ionizing groups), amino acid identity (substitution test), and NaCl concentration (ion screening test). In order to interpret the results of residue replacement, several factors in addition to the putative Glu 2- ... Arg 10+ interaction have been studied: intrinsic helix-forming tendencies of amino acids; interactions of charged residues with the alpha-helix macrodipole; and helix-lengthening effects. The results provide strong evidence that the Glu 2- ... Arg 10+ interaction is linked to helix formation and contributes to the stability of the isolated C-peptide helix. NMR evidence supports these conclusions and suggests that this interaction also acts as the N-terminal helix stop signal. The implications of this work for protein folding and stability are discussed.

Inhibition of mutant troponin C activity by an intra-domain disulphide bond

Triggering of contraction in striated muscles involves a conformational transition in the N-terminal domain of troponin C, the calcium-binding component of thin filaments. We have designed a mutant troponin C in which the key conformational transition and the calcium-regulatory activity are reversibly blocked by the formation of a disulphide bridge. Our results may be applicable to other proteins of the same family of calcium-binding proteins.

Design of crystalline helices of short oligopeptides as a possible model for nucleation of alpha-helix: role of water molecules in stabilizing helices

We have designed, synthesized, crystallized, and performed x-ray analysis of several hydrophobic tripeptides that show an extended near alpha-helical structure in the crystalline state. All of the tripeptides that show this remarkably stable helix crystallize with two or three water molecules; they all have glycine at the N terminus and have increasing hydrophobicity as one moves from the N to C terminus. Even though three residues in the oligomer are not sufficient to complete a turn, one of the water molecules acts as an added residue and links up adjacent tripeptide segments along the helix axis so that in the crystal, the helix appears effectively as one long continuous helix. Two of these tripeptides are stabilized by two water molecules that enable the peptides to complete a turn of the helix and extend the helical structure throughout the crystal by linking translationally related peptides by hydrogen bonds. In two other peptides, these roles are played by three rather than two water molecules. Though these tripeptides have different crystal symmetry, they all show the basic pattern of hydrated helix and packing, indicating the strong conformational preference for a stable structure even for these tripeptides. Such conformationally stable hydrated structures for short specific related sequences illustrate their possible importance in nucleating protein folding and in the role water molecules play in such events.

The major component of the paraflagellar rod of Trypanosoma brucei is a helical protein that is encoded by two identical, tandemly linked genes

The flagellum of the parasitic hemoflagellate Trypanosoma brucei contains two major structures: (a) the microtubule axoneme, and (b) a highly ordered, filamentous array, the paraflagellar rod (PFR). This is a complex, three-dimensional structure, of yet unknown function, that extends along most of the axoneme and is closely linked to it. Its major structural component is a single protein of 600 amino acids. This PFR protein can assume two different conformations, resulting in two distinct bands of apparent molecular masses of 73 and 69 kD in SDS-gel electrophoresis. Secondary structure predictions indicate a very high helix content. Despite its biochemical similarity to the intermediate filament proteins (solubility properties, amino acid composition, and high degree of helicity), the PFR protein does not belong in this class of cytoskeletal proteins. The PFR protein is coded for by two tandemly linked genes of identical nucleotide sequence. Both genes are transcribed into stable mRNAs of very similar length that carry the mini-exon sequence at their 5' termini.

Mapping myosin light chains by immunoelectron microscopy. Use of anti-fluorescyl antibodies as structural probes

The two classes of light chains in vertebrate fast muscle myosin have been selectively labeled with the thiol specific reagent 5-(iodoacetamido) fluorescein to determine their location in the myosin head. The alkali light chains (A1 and A2) were labeled at a single cysteine residue near the COOH terminus, whereas the regulatory light chain (LC2) was reacted at either cysteine 125 or 154. The two cysteines of LC2 appear to be near each other in the tertiary structure as evidenced by the ease of formation of an intramolecular disulfide bond. Besides having favorable spectral properties, fluorescein is a potent haptenic immunogen for raising high affinity antibodies. When anti-fluorescyl antibodies were added to the fluorescein-labeled light chains, the fluorescence was quenched by greater than 90%, thereby providing a simple method for determining an association constant. The interaction with antibody was the same for light chains exchanged into myosin as for free light chains. Complexes of antibody bound to light chain could be visualized in the electron microscope by rotary shadowing with platinum. By this approach we have shown that the COOH-terminal regions of the two classes of light chains are widely separated in myosin: the cysteine residues of LC2 lie close to the head/rod junction, whereas the single cysteine of A1 or A2 is located approximately 90 A distal to the junction. These sites correspond to the positions of the NH2 termini of the light chains mapped in earlier studies (Winkelmann, D. A., and S. Lowey. 1986. J. Mol. Biol. 188:595-612; Tokunaga, M., M. Suzuki, K. Saeki, and T. Wakabayashi. 1987b. J. Mol. Biol. 194:245-255). We conclude that the two classes of light chains do not lie in a simple colinear arrangement, but instead have a more complex organization in distinct regions of the myosin head.

pH dependent conformation of physalaemin

The undecapeptide physalaemin was investigated by n.m.r. spectroscopy in DMSO solution under acidic and neutral conditions. Large changes of the NH chemical shifts and the temperature gradients of the NH protons occurred on going from pH 3.5 to pH 7.0 for residues around the charged amino acids Asp and Lys. At pH 3.5 the data are in accord with a flexible conformation of the peptide. The results at neutral pH are interpreted in terms of a folded structure having two interresidue and one intraresidue hydrogen bond. They include a beta turn with proline in position i + 1 and asparagine in position i + 2.

Acid-induced dimerization of skeletal troponin C

We have investigated pH-dependent changes of the properties of troponin C from rabbit skeletal muscle. At pH 7.5 this protein is a monomer and at pH 5.2 it is a dimer. In contrast, bovine cardiac troponin C remains essentially monomeric at pH 5.2. Bovine brain calmodulin is not a dimer, but significantly aggregated at the same acidic pH. The dimerization of skeletal troponin C was demonstrated by low-speed (16,000 rpm) sedimentation equilibrium measurements carried out at 20 degrees C and by polyacrylamide gel electrophoresis under nondenaturing conditions. Dimer formation was significantly inhibited in the ultracentrifuge at rotor speeds of 30,000 and 40,000 rpm at 20 degrees C, and was completely prevented at a rotor speed of 40,000 rpm and 4 degrees C. This temperature and pressure dependence of dimerization strongly suggests that hydrophobic bonding is a major factor in promoting skeletal troponin C association at pH 5.2. The intramolecular distance between Met-25 and Cys-98 of rabbit skeletal troponin C deduced from fluorescence resonance energy transfer measurements increased by a factor of two upon lowering the pH from 7.5 to 5.2, indicating a pH-dependent transition in which the protein changed from a relatively compact conformation to an elongated conformation. The proton-induced increase in the energy transfer distance is related to the acid-induced dimerization of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Salt or ion bridges in biological systems: a study employing quantum and molecular mechanics

Equilibrium geometries and binding energies of model "salt" or "ion" bridge systems have been computed by ab initio quantum chemistry techniques (GAUSSIAN82) and by empirical force field techniques (AMBER2.0). Formate and dimethyl phosphate served as anions in the model compounds while interacting with several organic cations, including methyl ammonium, methyl guanidinium, and divalent metal ion (either Mg2+ or Ca2+) without and with an additional chloride; and a divalent metal ion (either Mg2+ or Ca2+), chloride, and four water molecules of hydration about the metal ion. The majority of the quantum chemical computations were performed using a split-valence basis set. For the model compounds studied we find that the ab initio optimized geometries are in remarkably good agreement with the molecular mechanics geometries. Several calculations were also performed using diffuse fractions. The formate anion binds these model cations more strongly than does dimethyl phosphate, while the organic cation methyl ammonium binds model anions more strongly than does methyl guanidinium. Finally, in model compounds including organic anions, Mg2+ or Ca2+ and four molecules of water, and a chloride anion, we find that the equilibrium structure of the magnesium complex involves a solvent separated ion pair (the magnesium ion is six coordinate), whereas the calcium ion complex remains seven coordinate. Molecular mechanics overestimates binding energies, but the estimates may be close enough to actual binding energies to give useful insight into the details of salt bridges in biological systems.

Molecular modeling of calmodulin: a comparison with crystallographic data

Two methods of side-chain placement on a modeled protein have been examined. Two molecular models of calmodulin were constructed that differ in the treatment of side chains prior to optimization of the molecule. A virtual bond analysis program developed by Purisima and Scheraga was used to determine the backbone conformation based on 2.2 angstroms resolution C alpha coordinates for the molecules. In the first model, side chains were initially constructed in an extended conformation. In the second model, a conformational grid search technique was employed. Calcium ions were treated explicitly during energy optimization using CHARMM. The models are compared to a recently published refined crystal structure of calmodulin. The results indicate that the initial choices for side-chains, but also significant effects on the main-chain conformation and supersecondary structure. The conformational differences are discussed. Analysis of these and other methods makes possible the formulation of a methodology for more appropriate side-chain placement in modeled proteins.

The leucine zipper: a hypothetical structure common to a new class of DNA binding proteins

A 30-amino-acid segment of C/EBP, a newly discovered enhancer binding protein, shares notable sequence similarity with a segment of the cellular Myc transforming protein. Display of these respective amino acid sequences on an idealized alpha helix revealed a periodic repetition of leucine residues at every seventh position over a distance covering eight helical turns. The periodic array of at least four leucines was also noted in the sequences of the Fos and Jun transforming proteins, as well as that of the yeast gene regulatory protein, GCN4. The polypeptide segments containing these periodic arrays of leucine residues are proposed to exist in an alpha-helical conformation, and the leucine side chains extending from one alpha helix interdigitate with those displayed from a similar alpha helix of a second polypeptide, facilitating dimerization. This hypothetical structure is referred to as the "leucine zipper," and it may represent a characteristic property of a new category of DNA binding proteins.