Acylation stabilizes a protease-resistant conformation of protoporphyrinogen oxidase, the molecular target of diphenyl ether-type herbicides

Protein acylation is an important way in which a number of proteins with a variety of functions are modified. The physiological role of the acylation of cellular proteins is still poorly understood. Covalent binding of fatty acids to nonintegral membrane proteins is thought to produce transient or permanent enhancement of the association of the polypeptide chains with biological membranes. In this paper, we investigate the functional role for the palmitoylation of an atypical membrane-bound protein, yeast protoporphyrinogen oxidase, which is the molecular target of diphenyl ether-type herbicides. Palmitoylation stabilizes an active heat- and protease-resistant conformation of the protein. Palmitoylation of protoporphyrinogen oxidase has been demonstrated to occur in vivo both in yeast cells and in a heterologous bacterial expression system, where it may be inhibited by cerulenin leading to the accumulation of degradation products of the protein. The thiol ester linking palmitoleic acid to the polypeptide chain was shown to be sensitive to hydrolysis by hydroxylamine and also by the widely used serine-protease inhibitor phenylmethylsulfonyl fluoride.

A peptide that inhibits hydroxyapatite growth is in an extended conformation on the crystal surface

Proteins play an important role in the biological mechanisms controlling hard tissue development, but the details of molecular recognition at inorganic crystal interfaces remain poorly characterized. We have applied a recently developed homonuclear dipolar recoupling solid-state NMR technique, dipolar recoupling with a windowless sequence (DRAWS), to directly probe the conformation of an acidic peptide adsorbed to hydroxyapatite (HAP) crystals. The phosphorylated hexapeptide, DpSpSEEK (N6, where pS denotes phosphorylated serine), was derived from the N terminus of the salivary protein statherin. Constant-composition kinetic characterization demonstrated that, like the native statherin, this peptide inhibits the growth of HAP seed crystals when preadsorbed to the crystal surface. The DRAWS technique was used to measure the internuclear distance between two 13C labels at the carbonyl positions of the adjacent phosphoserine residues. Dipolar dephasing measured at short mixing times yielded a mean separation distance of 3.2 ± 0.1 Å. Data obtained by using longer mixing times suggest a broad distribution of conformations about this average distance. Using a more complex model with discrete α-helical and extended conformations did not yield a better fit to the data and was not consistent with chemical shift analysis. These results suggest that the peptide is predominantly in an extended conformation rather than an α-helical state on the HAP surface. Solid-state NMR approaches can thus be used to determine directly the conformation of biologically relevant peptides on HAP surfaces. A better understanding of peptide and protein conformation on biomineral surfaces may provide design principles useful for the modification of orthopedic and dental implants with coatings and biological growth factors that are designed to enhance biocompatibility with surrounding tissue.

Molecular cloning of a peroxisomal Ca2+-dependent member of the mitochondrial carrier superfamily

A cDNA from a novel Ca2+-dependent member of the mitochondrial solute carrier superfamily was isolated from a rabbit small intestinal cDNA library. The full-length cDNA clone was 3,298 nt long and coded for a protein of 475 amino acids, with four elongation factor-hand motifs located in the N-terminal half of the molecule. The 25-kDa N-terminal polypeptide was expressed in Escherichia coli, and it was demonstrated that it bound Ca2+, undergoing a reversible and specific conformational change as a result. The conformation of the polypeptide was sensitive to Ca2+ which was bound with high affinity (Kd ≈ 0.37 μM), the apparent Hill coefficient for Ca2+-induced changes being about 2.0. The deduced amino acid sequence of the C-terminal half of the molecule revealed 78% homology to Grave disease carrier protein and 67% homology to human ADP/ATP translocase; this sequence homology identified the protein as a new member of the mitochondrial transporter superfamily. Northern blot analysis revealed the presence of a single transcript of about 3,500 bases, and low expression of the transporter could be detected in the kidney but none in the liver. The main site of expression was the colon with smaller amounts found in the small intestine proximal to the ileum. Immunoelectron microscopy localized the transporter in the peroxisome, although a minor fraction was found in the mitochondria. The Ca2+ binding N-terminal half of the transporter faces the cytosol.

Molecular switch in signal transduction: Reaction paths of the conformational changes in ras p21

Conformational changes in ras p21 triggered by the hydrolysis of GTP play an essential role in the signal transduction pathway. The path for the conformational change is determined by molecular dynamics simulation with a holonomic constraint directing the system from the known GTP-bound structure (with the γ-phosphate removed) to the GDP-bound structure. The simulation is done with a shell of water molecules surrounding the protein. In the switch I region, the side chain of Tyr-32, which undergoes a large displacement, moves through the space between loop 2 and the rest of the protein, rather than on the outside of the protein. As a result, the charged residues Glu-31 and Asp-33, which interact with Raf in the homologous RafRBD–Raps complex, remain exposed during the transition. In the switch II region, the conformational changes of α2 and loop 4 are strongly coupled. A transient hydrogen bonding complex between Arg-68 and Tyr-71 in the switch II region and Glu-37 in switch I region stabilizes the intermediate conformation of α2 and facilitates the unwinding of a helical turn of α2 (residues 66–69), which in turn permits the larger scale motion of loop 4. Hydrogen bond exchange between the protein and solvent molecules is found to be important in the transition. Possible functional implications of the results are discussed.

Smooth muscle myosin mutants containing a single tryptophan reveal molecular interactions at the actin-binding interface

Elucidation of the molecular details of the cyclic actomyosin interaction requires the ability to examine structural changes at specific sites in the actin-binding interface of myosin. To study these changes dynamically, we have expressed two mutants of a truncated fragment of chicken gizzard smooth muscle myosin, which includes the motor domain and essential light chain (MDE). These mutants were engineered to contain a single tryptophan at (Trp-546) or near (Trp-625) the putative actin-binding interface. Both 546- and 625-MDE exhibited actin-activated ATPase and actin-binding activities similar to wild-type MDE. Fluorescence emission spectra and acrylamide quenching of 546- and 625-MDE suggest that Trp-546 is nearly fully exposed to solvent and Trp-625 is less than 50% exposed in the presence and absence of ATP, in good agreement with the available crystal structure data. The spectrum of 625-MDE bound to actin was quite similar to the unbound spectrum indicating that, although Trp-625 is located near the 50/20-kDa loop and the 50-kDa cleft of myosin, its conformation does not change upon actin binding. However, a 10-nm blue shift in the peak emission wavelength of 546-MDE observed in the presence of actin indicates that Trp-546, located in the A-site of the lower 50-kDa subdomain of myosin, exists in a more buried environment and may directly interact with actin in the rigor acto-S1 complex. This change in the spectrum of Trp-546 constitutes direct evidence for a specific molecular interaction between residues in the A-site of myosin and actin.

Molecular mechanism of use-dependent calcium channel block by phenylalkylamines: Role of inactivation

The role of channel inactivation in the molecular mechanism of calcium (Ca2+) channel block by phenylalkylamines (PAA) was analyzed by designing mutant Ca2+ channels that carry the high affinity determinants of the PAA receptor site [Hockerman, G. H., Johnson, B. D., Scheuer, T., and Catterall, W. A. (1995) J. Biol. Chem. 270, 22119–22122] but inactivate at different rates. Use-dependent block by PAAs was studied after expressing the mutant Ca2+ channels in Xenopus oocytes. Substitution of single putative pore-orientated amino acids in segment IIIS6 by alanine (F-1499-A, F-1500-A, F-1510-A, I-1514-A, and F-1515-A) gradually slowed channel inactivation and simultaneously reduced inhibition of barium currents (IBa) by (−)D600 upon depolarization by 100 ms steps at 0.1 Hz. This apparent reduction in drug sensitivity was only evident if test pulses were applied at a low frequency of 0.1 Hz and almost disappeared at the frequency of 1 Hz. (−)D600 slowed IBa recovery after maintained membrane depolarization (1–3 sec) to a comparable extent in all channel constructs. A drug-induced delay in the onset of IBa recovery from inactivation suggests that PAAs promote the transition to a deep inactivated channel conformation. These findings indicate that apparent PAA sensitivity of Ca2+ channels is not only defined by drug interaction with its receptor site but also crucially dependent on intrinsic gating properties of the channel molecule. A molecular model for PAA-Ca2+ channel interaction that accounts for the relationship between drug induced inactivation and channel block by PAA is proposed.

Structure and function in rhodopsin: Rhodopsin mutants with a neutral amino acid at E134 have a partially activated conformation in the dark state*

The Glu-134–Arg-135 residues in rhodopsin, located near the cytoplasmic end of the C helix, are involved in G protein binding, or activation, or both. Furthermore, the charge-neutralizing mutation Glu-134 to Gln-134 produces hyperactivity in the activated state and produces constitutive activity in opsin. The Glu/Asp-Arg charge pair is highly conserved in equivalent positions in other G protein-coupled receptors. To investigate the structural consequences of charge-neutralizing mutations at Glu-134 and Arg-135 in rhodopsin, single spin-labeled side chains were introduced at sites in the cytoplasmic domains of helices C (140), E (227), F (250), or G (316) to serve as “molecular sensors” of the local helix bundle conformation. In each of the spin-labeled rhodopsins, a Gln substitution was introduced at either Glu-134 or Arg-135, and the electron paramagnetic resonance spectrum of the spin label was used to monitor the structural response of the helix bundle. The results indicate that a Gln substitution at Glu-134 induces a photoactivated conformation around helices C and G even in the dark state, an observation of potential relevance to the hyperactivity and constitutive activity of the mutant. In contrast, little change is induced in helix F, which has been shown to undergo a dominant motion upon photoactivation. This result implies that the multiple helix motions accompanying photoactivation are not strongly coupled and can be induced to take place independently. Gln substitution at Arg-135 produces only minor structural changes in the dark- or light-activated conformation, suggesting that this residue is not a determinant of structure in the regions investigated, although it may be functionally important.

Specific Molecular Chaperone Interactions and an ATP-dependent Conformational Change Are Required during Posttranslational Protein Translocation into the Yeast ER

The posttranslational translocation of proteins across the endoplasmic reticulum (ER) membrane in yeast requires ATP hydrolysis and the action of hsc70s (DnaK homologues) and DnaJ homologues in both the cytosol and ER lumen. Although the cytosolic hsc70 (Ssa1p) and the ER lumenal hsc70 (BiP) are homologous, they cannot substitute for one another, possibly because they interact with specific DnaJ homologues on each side of the ER membrane. To investigate this possibility, we purified Ssa1p, BiP, Ydj1p (a cytosolic DnaJ homologue), and a GST–63Jp fusion protein containing the lumenal DnaJ region of Sec63p. We observed that BiP, but not Ssa1p, is able to associate with GST–63Jp and that Ydj1p stimulates the ATPase activity of Ssa1p up to 10-fold but increases the ATPase activity of BiP by <2-fold. In addition, Ydj1p and ATP trigger the release of an unfolded polypeptide from Ssa1p but not from BiP. To understand further how BiP drives protein translocation, we purified four dominant lethal mutants of BiP. We discovered that each mutant is defective for ATP hydrolysis, fails to undergo an ATP-dependent conformational change, and cannot interact with GST–63Jp. Measurements of protein translocation into reconstituted proteoliposomes indicate that the mutants inhibit translocation even in the presence of wild-type BiP. We conclude that a conformation- and ATP-dependent interaction of BiP with the J domain of Sec63p is essential for protein translocation and that the specificity of hsc70 action is dictated by their DnaJ partners.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

Molecular mechanisms of translation initiation in eukaryotes

Translation initiation is a complex process in which initiator tRNA, 40S, and 60S ribosomal subunits are assembled by eukaryotic initiation factors (eIFs) into an 80S ribosome at the initiation codon of mRNA. The cap-binding complex eIF4F and the factors eIF4A and eIF4B are required for binding of 43S complexes (comprising a 40S subunit, eIF2/GTP/Met-tRNAi and eIF3) to the 5′ end of capped mRNA but are not sufficient to promote ribosomal scanning to the initiation codon. eIF1A enhances the ability of eIF1 to dissociate aberrantly assembled complexes from mRNA, and these factors synergistically mediate 48S complex assembly at the initiation codon. Joining of 48S complexes to 60S subunits to form 80S ribosomes requires eIF5B, which has an essential ribosome-dependent GTPase activity and hydrolysis of eIF2-bound GTP induced by eIF5. Initiation on a few mRNAs is cap-independent and occurs instead by internal ribosomal entry. Encephalomyocarditis virus (EMCV) and hepatitis C virus epitomize distinct mechanisms of internal ribosomal entry site (IRES)-mediated initiation. The eIF4A and eIF4G subunits of eIF4F bind immediately upstream of the EMCV initiation codon and promote binding of 43S complexes. EMCV initiation does not involve scanning and does not require eIF1, eIF1A, and the eIF4E subunit of eIF4F. Initiation on some EMCV-like IRESs requires additional noncanonical initiation factors, which alter IRES conformation and promote binding of eIF4A/4G. Initiation on the hepatitis C virus IRES is even simpler: 43S complexes containing only eIF2 and eIF3 bind directly to the initiation codon as a result of specific interaction of the IRES and the 40S subunit.

Molecular markers reveal that population structure of the human pathogen Candida albicans exhibits both clonality and recombination.

The life history of Candida albicans presents an enigma: this species is thought to be exclusively asexual, yet strains show extensive phenotypic variation. To address the population genetics of C. albicans, we developed a genetic typing method for codominant single-locus markers by screening randomly amplified DNA for single-strand conformation polymorphisms. DNA fragments amplified by arbitrary primers were initially screened for single-strand conformation polymorphisms and later sequenced using locus-specific primers. A total of 12 single base mutations and insertions were detected from six out of eight PCR fragments. Patterns of sequence-level polymorphism observed for individual strains detected considerable heterozygosity at the DNA sequence level, supporting the view that most C. albicans strains are diploid. Population genetic analyses of 52 natural isolates from Duke University Medical Center provide evidence for both clonality and recombination in C. albicans. Evidence for clonality is supported by the presence of several overrepresented genotypes, as well as by deviation of genotypic frequencies from random (Hardy-Weinberg) expectations. However, tests for nonrandom association of alleles across loci reveal less evidence for linkage disequilibrium than expected for strictly clonal populations. Although C. albicans populations are primarily clonal, evidence for recombination suggests that sexual reproduction or some other form of genetic exchange occurs in this species.

Molecular mechanics calculations of the riboacetal internucleotide linkage in double and triple helices.

Structures of Watson-Crick base paired 15-nucleobase oligomer strands in A-type or B-type conformation in which one strand [a strand of alternating nucleotide and riboacetal thymidine nucleoside (RT) units, RP] is DNA and the other is composed of alternating nucleotides and riboacetal nucleosides have been studied by molecular mechanics. Analogously, oligomer strands of RNA in place of DNA have been modeled. The calculations indicate that the RP strand is more stable when complexed in an A-type duplex relative to a B-type form and that this conformational preference is presumably due to the more uniform nature of the former. Nearly planar ribose rings were more commonly observed in the minimized structures of the B-type DNA.RP duplexes as compared with A-type duplexes, despite the fact that planar ribofuranose rings are known to be energetically unfavorable in oligonucleotides. Computed relative stabilities of all duplexes containing the RP strand suggest that such heteroduplexes are less stable than the corresponding double-stranded DNA and double-stranded RNA species. These findings are in agreement with experimental results which show, when equivalent sequences were compared, that a DNA.RNA control forms a more stable duplex than RP hound to a complementary single-stranded RNA strand. In contrast, molecular mechanics studies of complementary triple-helical (DNA)2.RP, (DNA)2.DNA, and (DNA)2.RNA structures indicate that the binding of RP as a Hoogsteen strand stabilizes the underlying duplex to a greater extent compared with native oligonucleotides. These calculations suggest that puckering of the ribose ring in the riboacetal linkage leads to a more favorable interaction with a complementary nucleic acid target than the proposed planar geometry and that this puckering may account for the enhanced binding of RP to a double-stranded target.

Molecular markers reveal cryptic sex in the human pathogen Coccidioides immitis.

Coccidioides immitis, cause of a recent epidemic of "Valley fever" in California, is typical of many eukaryotic microbes in that mating and meiosis have yet to be reported, but it is not clear whether sex is truly absent or just cryptic. To find out, we have undertaken a population genetic study using PCR amplification, screening for single-strand conformation polymorphisms, and direct DNA sequencing to find molecular markers with nucleotide-level resolution. Both population genetic and phylogenetic analyses indicate that C. immitis is almost completely recombining. To our knowledge, this study is the first to find molecular evidence for recombination in a fungus for which no sexual stage has yet been described. These results motivate a directed search for mating and meiosis and illustrate the utility of single-strand conformation polymorphism and sequencing with arbitrary primer pairs in molecular population genetics.

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.

Growth hormone secretagogues: characterization, efficacy, and minimal bioactive conformation.

Another class of growth hormone (GH) secretagogues has been discovered by altering the backbone structure of a flexible linear GH-releasing peptide (GHRP). In vitro and in vivo characterization confirms these GH secretagogues as the most potent and smallest (M(r) < 500) reported. Anabolic efficacy is demonstrated in rodents with intermittent delivery. A convergent model of the bioactive conformation of GHRPs is developed and is supported by the NMR structure of a highly potent cyclic analog of GHRP-2. The model and functional data provide a logical framework for the further design of low-molecular weight secretagogues and illustrate the utility of an interdisciplinary approach to elucidating potential bound-state conformations of flexible peptide ligands.