Large electrostatic differences in the binding thermodynamics of a cationic peptide to oligomeric and polymeric DNA.

Results presented here demonstrate that the thermodynamics of oligocation binding to polymeric and oligomeric DNA are not equivalent because of long-range electrostatic effects. At physiological cation concentrations (0.1-0.3 M) the binding of an oligolysine octacation KWK6-NH2 (+8 charge) to single-stranded poly(dT) is much stronger per site and significantly more salt concentration dependent than the binding of the same ligand to an oligonucleotide, dT(pdT)10 (-10 charge). These large differences are consistent with Poisson-Boltzmann calculations for a model that characterizes the charge distributions with key preaveraged structural parameters. Therefore, both the experimental and the theoretical results presented here show that the polyelectrolyte character of a polymeric nucleic acid makes a large contribution to both the magnitude and the salt concentration dependence of its binding interactions with simple oligocationic ligands.

Interresidue hydrogen bonding in a peptide nucleic acid.RNA heteroduplex.

Previous molecular mechanics calculations suggest that strands of peptide nucleic acids (PNAs) and complementary oligonucleotides form antiparallel duplexes stabilized by interresidue hydrogen bonds. In the computed structures, the amide carbonyl oxygen nearest the nucleobase (O7') forms an interresidue hydrogen bond with the backbone amide proton of the following residue, (n + 1)H1'. Of the 10 published two dimensional 1H NMR structures of a hexameric PNA.RNA heteroduplex. PNA(GAACTC).r(GAGUUC), 9 exhibit two to five potential interresidue hydrogen bonds. In our minimized average structure, created from the coordinates of these 10 NMR structures, three of the five possible interresidue hydrogen bond sites within the PNA backbone display the carbonyl oxygen (O7') and the amide proton (n + 1)H1' distances and N1'-H1'-(n - 1)O7' angles optimal for hydrogen bond formation. The finding of these interresidue hydrogen bonds supports the results of our previous molecular mechanics calculations.

Enhanced conformational sampling in Monte Carlo simulations of proteins: application to a constrained peptide.

A Monte Carlo simulation method for globular proteins, called extended-scaled-collective-variable (ESCV) Monte Carlo, is proposed. This method combines two Monte Carlo algorithms known as entropy-sampling and scaled-collective-variable algorithms. Entropy-sampling Monte Carlo is able to sample a large configurational space even in a disordered system that has a large number of potential barriers. In contrast, scaled-collective-variable Monte Carlo provides an efficient sampling for a system whose dynamics is highly cooperative. Because a globular protein is a disordered system whose dynamics is characterized by collective motions, a combination of these two algorithms could provide an optimal Monte Carlo simulation for a globular protein. As a test case, we have carried out an ESCV Monte Carlo simulation for a cell adhesive Arg-Gly-Asp-containing peptide, Lys-Arg-Cys-Arg-Gly-Asp-Cys-Met-Asp, and determined the conformational distribution at 300 K. The peptide contains a disulfide bridge between the two cysteine residues. This bond mimics the strong geometrical constraints that result from a protein's globular nature and give rise to highly cooperative dynamics. Computation results show that the ESCV Monte Carlo was not trapped at any local minimum and that the canonical distribution was correctly determined.

Three gonadotropin-releasing hormone genes in one organism suggest novel roles for an ancient peptide.

Gonadotropin-releasing hormone (GnRH) is known and named for its essential role in vertebrate reproduction. Release of this decapeptide from neurons in the hypothalamus controls pituitary gonadotropin levels which, in turn, regulate gonadal state. The importance of GnRH is underscored by its widespread expression and conservation across vertebrate taxa: five amino acids are invariant in all nine known forms, whereas two others show only conservative changes. In most eutherian mammals, only one form, expressed in the hypothalamus, is thought to exist, although in a recent report, antibody staining in developing primates suggests an additional form. In contrast, multiple GnRH forms and expression loci have been reported in many non-mammalian vertebrates. However, evidence based on immunological discrimination does not always agree with analysis of gene expression, since GnRH forms encoded by different genes may not be reliably distinguished by antibodies. Here we report the expression of three distinct GnRH genes in a teleost fish brain, including the sequence encoding a novel GnRH preprohormone. Using in situ hybridization, we show that this form is found only in neurons that project to the pituitary and exhibit changes in soma size depending on social and reproductive state. The other two GnRH genes are expressed in other, distinct cell populations. All three genes share the motif of encoding a polypeptide consisting of GnRH and a GnRH-associated peptide. Whereas the GnRH moiety is highly conserved, the GnRH-associated peptides are not, reflecting differential selective pressure on different parts of the gene. GnRH forms expressed in nonhypothalamic regions may serve to coordinate reproductive activities of the animal.

Inhibitors of human heart chymase based on a peptide library.

We have synthesized two sets of noncleavable peptide-inhibitor libraries to map the S and S' subsites of human heart chymase. Human heart chymase is a chymotrypsin-like enzyme that converts angiotensin I to angiotensin II. The first library consists of peptides with 3-fluorobenzylpyruvamides in the P1 position. (Amino acid residues of substrates numbered P1, P2, etc., are toward the N-terminal direction, and P'1, P'2, etc., are toward the C-terminal direction from the scissile bond.) The P'1 and P'2 positions were varied to contain each one of the 20 naturally occurring amino acids and P'3 was kept constant as an arginine. The second library consists of peptides with phenylalanine keto-amides at P1, glycine in P'1, and benzyloxycarbonyl (Z)-isoleucine in P4. The P2 and P3 positions were varied to contain each of the naturally occurring amino acids, except for cysteine and methionine. The peptides of both libraries are attached to a solid support (pins). The peptides are evaluated by immersing the pins in a solution of the target enzyme and evaluating the amount of enzyme absorbed. The pins with the best inhibitors will absorb most enzyme. The libraries select the best and worst inhibitors within each group of peptides and provide an approximate ranking of the remaining peptides according to Ki. Through this library, we determined that Z-Ile-Glu-Pro-Phe-CO2Me and (F)-Phe-CO-Glu-Asp-ArgOMe should be the best inhibitors of chymase in this collection of peptide inhibitors. We synthesized the peptides and found Ki values were 1 nM and 1 microM, respectively. The corresponding Ki values for chymotrypsin were 10 nM and 100 microM. The use of libraries of inhibitors has advantages over the classical method of synthesis of potential inhibitors in solution: the libraries are reusable, the same libraries can be used with a variety of different serine proteases, and the method allows the screening of hundreds of compounds in short periods of time.

Mapping sites of interaction of p47-phox and flavocytochrome b with random-sequence peptide phage display libraries.

During assembly of the phagocyte NADPH oxidase, cytosolic p47-phox translocates to the plasma membrane and binds to flavocytochrome b, and binding domains for p47-phox have been identified on the C-terminal tails of both flavocytochrome b subunits. In the present report, we further examine the interaction of these two oxidase components by using random-sequence peptide phage display library analysis. Screening p47-phox with the peptide libraries identified five potential sites of interaction with flavocytochrome b, including three previously reported regions of interaction and two additional regions of interaction of p47-phox with gp91-phox and p22-phox. The additional sites were mapped to a domain on the first predicted cytosolic loop of gp91-phox encompassing residues S86TRVRRQL93 and to a domain near the cytosolic C-terminal tail of gp91-phox encompassing residues F450EWFADLL457. The mapping also confirmed a previously reported binding domain on gp91-phox (E554SGPRGVHFIF564) and putative Src homology 3 domain binding sites on p22-phox (P156PRPP160 and G177GPPGGP183). To demonstrate that the additional regions identified were biologically significant, peptides mimicking the gp91-phox sequences F77LRGSSACCSTRVRRQL93 and E451WFADLLQLLESQ463 were synthesized and assayed for their ability to inhibit NADPH oxidase activity. These peptides had EC50 values of 1 microM and 230 microM, respectively, and inhibited activation when added prior to assembly but did not affect activity of the preassembled oxidase. Our data demonstrate the usefulness of phage display library analysis for the identification of biologically relevant sites of protein-protein interaction and show that the binding of p47-phox to flavocytochrome b involves multiple binding sites along the C-terminal tails of both gp91- and p22-phox and other regions of gp91-phox nearer to the N terminus.

Kinetics of hydrogen bond breakage in the process of unfolding of ribonuclease A measured by pulsed hydrogen exchange.

A sensitive test for kinetic unfolding intermediates in ribonuclease A (EC 3.1.27.5) is performed under conditions where the enzyme unfolds slowly (10 degrees C, pH 8.0, 4.5 M guanidinium chloride). Exchange of peptide NH protons (2H-1H) is used to monitor structural opening of individual hydrogen bonds during unfolding, and kinetic models are developed for hydrogen exchange during the process of protein unfolding. The analysis indicates that the kinetic process of unfolding can be monitored by EX1 exchange (limited by the rate of opening) for ribonuclease A in these conditions. Of the 49 protons whose unfolding/exchange kinetics was measured, 47 have known hydrogen bond acceptor groups. To test whether exchange during unfolding follows the EX2 (base-catalyzed) or the EX1 (uncatalyzed) mechanism, unfolding/exchange was measured both at pH 8.0 and at pH 9.0. A few faster-exchanging protons were found that undergo exchange by both EX1 and EX2 processes, but the 43 slower-exchanging protons at pH 8 undergo exchange only by the EX1 mechanism, and they have closely similar rates. Thus, it is likely that all 49 protons undergo EX1 exchange at the same rate. The results indicate that a single rate-limiting step in unfolding breaks the entire network of peptide hydrogen bonds and causes the overall unfolding of ribonuclease A. The additional exchange observed for some protons that follows the EX2 mechanism probably results from equilibrium unfolding intermediates and will be discussed elsewhere.