Phytosulfokine-α, a sulfated pentapeptide, stimulates the proliferation of rice cells by means of specific high- and low-affinity binding sites

Peptide growth factors were isolated from conditioned medium derived from rice (Oryza sativa L.) suspension cultures and identified to be a sulfated pentapeptide [H-Tyr(SO3H)-Ile-Tyr(SO3H)-Thr-Gln-OH] and its C-terminal-truncated tetrapeptide [H-Tyr(SO3H)-Ile-Tyr(SO3H)-Thr-OH]. These structures were identical to the phytosulfokines originally found in asparagus (Asparagus officinalis L.) mesophyll cultures. The pentapeptide [phytosulfokine-α (PSK-α)] very strongly stimulated colony formation of rice protoplasts at concentrations above 10−8 M, indicating a similar mode of action in rice of phytosulfokines. Binding assays using 35S-labeled PSK-α demonstrated the existence of both high- and low-affinity specific saturable binding sites on the surface of rice cells in suspension. Analysis of [35S]PSK-α binding in differential centrifugation fractions suggested association of the binding with a plasma membrane-enriched fraction. The apparent Kd values for [35S]PSK-α binding were found to be 1 × 10−9 M for the high-affinity type and 1 × 10−7 M for the low-affinity type, with maximal numbers of binding sites of 1 × 104 sites per cell and 1 × 105 sites per cell, respectively. Competition studies with [35S]PSK-α and several synthetic PSK-α analogs demonstrated that only peptides that possesses mitogenic activity can effectively displace the radioligand. These results suggest that a signal transduction pathway mediated by peptide factors is involved in plant cell proliferation.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Construction of a high-affinity receptor site for dihydropyridine agonists and antagonists by single amino acid substitutions in a non-l-type Ca2+ channel

The activity of l-type Ca2+ channels is increased by dihydropyridine (DHP) agonists and inhibited by DHP antagonists, which are widely used in the therapy of cardiovascular disease. These drugs bind to the pore-forming α1 subunits of l-type Ca2+ channels. To define the minimal requirements for DHP binding and action, we constructed a high-affinity DHP receptor site by substituting a total of nine amino acid residues from DHP-sensitive l-type α1 subunits into the S5 and S6 transmembrane segments of domain III and the S6 transmembrane segment of domain IV of the DHP-insensitive P/Q-type α1A subunit. The resulting chimeric α1A/DHPS subunit bound DHP antagonists with high affinity in radioligand binding assays and was inhibited by DHP antagonists with high affinity in voltage clamp experiments. Substitution of these nine amino acid residues yielded 86% of the binding energy of the l-type α1C subunit and 92% of the binding energy of the l-type α1S subunit for the high-affinity DHP antagonist PN200–110. The activity of chimeric Ca2+ channels containing α1A/DHPS was increased 3.5 ± 0.7-fold by the DHP agonist (−)Bay K8644. The effect of this agonist was stereoselective as in l-type Ca2+ channels since (+) Bay K8644 inhibited the activity of α1A/DHPS. The results show conclusively that DHP agonists and antagonists bind to a single receptor site at which they have opposite effects on Ca2+ channel activity. This site contains essential components from both domains III and IV, consistent with a domain interface model for binding and allosteric modulation of Ca2+ channel activity by DHPs.

Peptide interfacial adsorption is kinetically limited by the thermodynamic stability of self association

We present a study of the adsorption of two peptides at the octane–water interface. The first peptide, Lac21, exists in mixed monomer–tetramer equilibrium in bulk solution with an appreciable monomer concentration. The second peptide, Lac28, exists as a tetramer in solution, with minimal exposed hydrophobic surface. A kinetic limitation to interfacial adsorption exists for Lac28 at moderate to high surface coverage that is not observed for Lac21. We estimate the potential energy barrier for Lac28 adsorption to be 42 kJ/mol and show that this is comparable to the expected free energy barrier for tetramer dissociation. This finding suggests that, at moderate to high surface coverage, adsorption is kinetically limited by the availability of interfacially active monomeric “domains” in the subinterfacial region. We also show how the commonly used empirical equation for protein adsorption dynamics can be used to estimate the potential energy barrier for adsorption. Such an approach is shown to be consistent with a formal description of diffusion–adsorption, provided a large potential energy barrier exists. This work demonstrates that the dynamics of interfacial adsorption depend on protein thermodynamic stability, and hence structure, in a quantifiable way.

Kinetic equilibrium of forces and molecular events in muscle contraction

In biomolecular systems, the mechanical transfer of free energy occurs with both high efficiency and high speed. It is shown here that such a transfer can be achieved only if the participating free-energy-storing elements exhibit opposing relationships between their content of free energy and the force they exert in the transfer direction. A kinetic equilibrium of forces (KEF) results, in which the transfer of free energy is mediated essentially by thermal molecular motion. On the basis of present evidence, KEF is used as a guiding principle in developing a mechanical model of the crossbridge cycle in muscle contraction. The model allows the basic features of molecular events to be visualized in terms of plausible structures. Real understanding of the process will require identification of the elements that perform the functions described here. Besides chemomechanical energy transduction, KEF may have a role in other biomolecular processes in which free energy is transferred mechanically over large distances.

Copper binding to the prion protein: Structural implications of four identical cooperative binding sites

Evidence is growing to support a functional role for the prion protein (PrP) in copper metabolism. Copper ions appear to bind to the protein in a highly conserved octapeptide repeat region (sequence PHGGGWGQ) near the N terminus. To delineate the site and mode of binding of Cu(II) to the PrP, the copper-binding properties of peptides of varying lengths corresponding to 2-, 3-, and 4-octarepeat sequences have been probed by using various spectroscopic techniques. A two-octarepeat peptide binds a single Cu(II) ion with Kd ≈ 6 μM whereas a four-octarepeat peptide cooperatively binds four Cu(II) ions. Circular dichroism spectra indicate a distinctive structuring of the octarepeat region on Cu(II) binding. Visible absorption, visible circular dichroism, and electron spin resonance spectra suggest that the coordination sphere of the copper is identical for 2, 3, or 4 octarepeats, consisting of a square-planar geometry with three nitrogen ligands and one oxygen ligand. Consistent with the pH dependence of Cu(II) binding, proton NMR spectroscopy indicates that the histidine residues in each octarepeat are coordinated to the Cu(II) ion. Our working model for the structure of the complex shows the histidine residues in successive octarepeats bridged between two copper ions, with both the Nɛ2 and Nδ1 imidazole nitrogen of each histidine residue coordinated and the remaining coordination sites occupied by a backbone amide nitrogen and a water molecule. This arrangement accounts for the cooperative nature of complex formation and for the apparent evolutionary requirement for four octarepeats in the PrP.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

Stretching single-domain proteins: Phase diagram and kinetics of force-induced unfolding

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise “quantized” manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.

Role of cysteines in Plasmodium falciparum circumsporozoite protein: Interactions with heparin can rejuvenate inactive protein mutants

Various pathogenic bacteria, viruses, and protozoan bind to glycosaminoglycan-based receptors on host cells and initiate an infection. Sporozoites of Plasmodium predominantly express circumsporozoite (CS) protein on their surface, which binds to heparan sulfate proteoglycans on liver cell surface that subsequently leads to malaria. Here we show that the interaction of free heparin with this parasite ligand has the potential to be a critical component of invasion. CS protein of P. falciparum contains four cysteines at positions 361, 365, 396, and 401. In this study, all four cysteine residues were mutagenized to alanine both individually and in different combinations. Conversion of cysteine 396 to alanine (protein CS3) led to a 10-fold increase in the binding activity of the protein to HepG2 cells. Replacement of cysteines at positions 361, 365, and 401 either alone or in different combinations led to a near total loss of binding. Surprisingly, activity in these inactive mutants could be effectively restored in the presence of submolar concentrations of heparin. Heparin also up-regulated binding of CS3 at submolar concentrations with respect to the protein but down-regulated binding when present in excess. Given the significantly different concentrations of heparin in different organs of the host and the in vitro results described here one can consider in vivo ramifications of this phenomenon for pathogen targeting of specific organs and for the functional effects of antigenic variation on receptor ligand interaction.

Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

The vibrational energy flow transition in organic molecules: Theory meets experiment

Most large dynamical systems are thought to have ergodic dynamics, whereas small systems may not have free interchange of energy between degrees of freedom. This assumption is made in many areas of chemistry and physics, ranging from nuclei to reacting molecules and on to quantum dots. We examine the transition to facile vibrational energy flow in a large set of organic molecules as molecular size is increased. Both analytical and computational results based on local random matrix models describe the transition to unrestricted vibrational energy flow in these molecules. In particular, the models connect the number of states participating in intramolecular energy flow to simple molecular properties such as the molecular size and the distribution of vibrational frequencies. The transition itself is governed by a local anharmonic coupling strength and a local state density. The theoretical results for the transition characteristics compare well with those implied by experimental measurements using IR fluorescence spectroscopy of dilution factors reported by Stewart and McDonald [Stewart, G. M. & McDonald, J. D. (1983) J. Chem. Phys. 78, 3907–3915].

The nature of the excited state of the reaction center of photosystem II of green plants: A high-resolution fluorescence spectroscopy study

We studied the electronically excited state of the isolated reaction center of photosystem II with high-resolution fluorescence spectroscopy at 5 K and compared the obtained spectral features with those obtained earlier for the primary electron donor. The results show that there is a striking resemblance between the emitting and charge-separating states in the photosystem II reaction center, such as a very similar shape of the phonon wing with characteristic features at 19 and 80 cm−1, almost identical frequencies of a number of vibrational modes, a very similar double-Gaussian shape of the inhomogeneous distribution function, and relatively strong electron-phonon coupling for both states. We suggest that the emission at 5 K originates either from an exciton state delocalized over the inactive branch of the photosystem or from a fraction of the primary electron donor that is long-lived at 5 K. The latter possibility can be explained by a distribution of the free energy difference of the primary charge separation reaction around zero. Both possibilities are in line with the idea that the state that drives primary charge separation in the reaction center of photosystem II is a collective state, with contributions from all chlorophyll molecules in the central part of the complex.

Thermostable and site-specific DNA binding of the gene product ORF56 from the Sulfolobus islandicus plasmid pRN1, a putative archael plasmid copy control protein

There is still a lack of information on the specific characteristics of DNA-binding proteins from hyperthermophiles. Here we report on the product of the gene orf56 from plasmid pRN1 of the acidophilic and thermophilic archaeon Sulfolobus islandicus. orf56 has not been characterised yet but low sequence similarily to several eubacterial plasmid-encoded genes suggests that this 6.5 kDa protein is a sequence-specific DNA-binding protein. The DNA-binding properties of ORF56, expressed in Escherichia coli, have been investigated by EMSA experiments and by fluorescence anisotropy measurements. Recombinant ORF56 binds to double-stranded DNA, specifically to an inverted repeat located within the promoter of orf56. Binding to this site could down-regulate transcription of the orf56 gene and also of the overlapping orf904 gene, encoding the putative initiator protein of plasmid replication. By gel filtration and chemical crosslinking we have shown that ORF56 is a dimeric protein. Stoichiometric fluorescence anisotropy titrations further indicate that ORF56 binds as a tetramer to the inverted repeat of its target binding site. CD spectroscopy points to a significant increase in ordered secondary structure of ORF56 upon binding DNA. ORF56 binds without apparent cooperativity to its target DNA with a dissociation constant in the nanomolar range. Quantitative analysis of binding isotherms performed at various salt concentrations and at different temperatures indicates that approximately seven ions are released upon complex formation and that complex formation is accompanied by a change in heat capacity of –6.2 kJ/mol. Furthermore, recombinant ORF56 proved to be highly thermostable and is able to bind DNA up to 85°C.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.