A single carboxy-terminal arginine determines the amino-terminal helix conformation of an alanine-based peptide

Sequence rules for gold-binding peptides

Metal-binding peptides play a central role in bionanotechnology, wherein they are responsible for directing growth and influencing the resulting properties of inorganic nanomaterials. One of the key advantages of using peptides to create nanomaterials is their versatility, wherein subtle changes in the sequence can have a dramatic effect on the structure and properties of the nanomaterial. However, precisely knowing which position and which amino acid should be modified within a given sequence to enhance a specific property can be a daunting challenge owing to combinatorial complexity. In this study, classification based on association rules was performed using 860 gold-binding peptides. Using a minimum support threshold of 0.035 and confidence of 0.9, 30 rules with confidence and lift values greater than 0.9 and 1, respectively, were extracted that can differentiate high-binding from low-binding peptides. The test performance of these rules for categorizing the peptides was found to be satisfactory, as characterized by accuracy = 0.942, F₁ = 0.941, MCC = 0.884. What stands out from the extracted rules are the importance of tryptophan and arginine residues in differentiating peptides with high binding affinity from those with low affinity. In addition, the association rules revealed that positions 2 and 4 within a decapeptide are frequently involved in the rules, thus suggesting their importance in influencing peptide binding affinity to AuNPs. Collectively, this study identified sequence rules that may be used to design peptides with high binding affinity.

Growth Reaction of Gold Nanorods in the Presence of Mutated Peptides and Amine-Modified Single-Stranded Nucleic Acids

Conformation of biomolecules like DNA, peptides and amino acids play vital role during nanoparticle growth. Herein, we have experimentally explored the effect of different noncovalent interaction between a 5'-amine modified DNA sequence (NH₂ -C₆ H₁₂ -5'-ACATCAGT-3', PMR) and arginine during the seed-mediated growth reaction of gold nanorods (GNRs). Amino acid-mediated growth reaction of GNRs results in a snowflake-like gold nanoarchitecture. However, in case of Arg, prior incubation of GNRs with PMR selectively produces sea urchin-like gold suprastructures, via strong hydrogen bonding and cation-π interaction between PMR and Arg. This distinctive structure formation strategy has been extended to study the structural modulation caused by two structurally close α-helical RRR (Ac-(AAAAR)₃ A-NH₂ ) peptide and the lysine mutated KKR (Ac-AAAAKAAAAKAAAARA-NH₂ ) peptide with partial helix at the amino terminus. Simulation studies confirm that a greater number of hydrogen bonding and cation-π interaction between the Arg residues and PMR resulted in the gold sea urchin structure for RRR peptide against KKR peptide.

Heterozygosity for bisphosphoglycerate mutase deficiency expressing clinically as congenital erythrocytosis: A case series and literature review

Erythrocytosis is associated with increased red blood cell mass and can be either congenital or acquired. Congenital secondary causes are rare and include germline variants increasing haemoglobin (Hb)-oxygen affinity (e.g., Hb or bisphosphoglycerate mutase (BPGM) variants) or affecting oxygen-sensing pathway proteins. Here, we describe five adults from three kindreds with erythrocytosis associated with heterozygosity for BPGM variants, including one novel. Functional analyses showed partial BPGM deficiency, reduced 2,3-bisphosphoglycerate levels and/or increased Hb-oxygen affinity. We also review currently known BPGM variants. This study contributes to raising awareness of BPGM variants, and in particular that heterozygosity for BPGM deficiency may already manifest clinically.

Cationic Side Chain Identity Directs the Hydrophobically Driven Self-Assembly of Amphiphilic β-Peptides in Aqueous Solution

Hydrophobic interactions mediated by nonpolar molecular fragments in water are influenced by local chemical and physical contexts in ways that are not yet fully understood. Here, we use globally amphiphilic (GA) β-peptides (GA-Lys and GA-Arg) with stable conformations to explore if replacement of β³-homolysine (βLys) with β³-homoarginine (βArg) influences the hydrophobically driven assembly of these peptides in bulk aqueous solution. The studies were conducted in 10 mM triethanolamine buffer at pH 7, where both βLys (ammonium) and βArg (guanidinium) side chains are substantially protonated. Comparisons of light scattering measurements and cryo-electron micrographs before and after the addition of 60 vol% MeOH indicate very different outcomes of the hydrophobically driven assembly of AcY-GA-Lys versus AcY-GA-Arg (AcY denotes an N-acetylated-β³-homotyrosine (βTyr) at each N-terminus). Nuclear magnetic resonance and analytical ultracentrifugation confirm that AcY-GA-Lys assembles into large aggregates in aqueous buffer, whereas AcY-GA-Arg at comparable concentrations forms only small oligomers. Titration of AcY-GA-Arg into aqueous solutions of AcY-GA-Lys reveals that the driving force for AcY-GA-Lys association is far stronger than that for AcY-GA-Arg association. We discuss these results in the light of past experimental observations involving single molecule force measurements with GA β-peptides and hydrophobically driven dimerization of conventional peptides that form a GA α-helix upon dimerization (but do not display the Lys versus Arg trend predicted by extrapolating from the earlier AFM studies with β-peptides). Overall, our results establish that the identity of proximal cationic groups, ammonium versus guanidinium, profoundly modulates the hydrophobically driven self-assembly of conformationally stable β-peptides in bulk aqueous solution.

Unraveling unidirectional threading of α-cyclodextrin in a [2]rotaxane through spin labeling approach

We present here the results of a CW-ESR investigation of a double spin labeled α-cyclodextrin-based [2]rotaxane that is characterized by the presence of nitroxide labels both at the wheel and at the dumbbell. This was accomplished by synthesizing a spin labeled α-CD (the wheel) that was mechanically blocked on a thread containing a nitroxide unit by a Cu(I) catalyzed azide-alkyne cycloaddition (CuAAC). Both ESI-MS analysis and NMR spectroscopy were used to provide evidence of the threading processes. Because of the unsymmetrical structure of both the wheel and the axle, two different geometrical isomers could be predicted on the basis of the orientation of the CD along the thread. By combining molecular dynamic calculations and information extracted from the CW-ESR spectra, we were able to determine the geometrical nature of the isomer that was isolated as the only species. The ESR spectra showed J-coupling between the two mechanically connected nitroxide units and were analyzed by a model assuming three main molecular states. The intramolecular noncovalent nature of spin exchange was confirmed by reversibly switching the magnetic interaction on-off by changing the pH of the solution in the presence of a competing macrocyclic host.

Water-soluble poly(L-serine)s with elongated and charged side-chains: synthesis, conformations, and cell-penetrating properties

Water-soluble poly(L-serine)s bearing long side-chain with terminal charge groups were synthesized via ring-opening polymerization of O-pentenyl-L-serine N-carboxyanhydride followed by thiol-ene reactions. These side-chain modified poly(L-serine)s adopt β-sheet conformation in aqueous solution with excellent stability against changes in pH and temperature. These water-soluble poly(L-serine) derivatives with charged side-chain functional groups and stable β-sheet conformations showed membrane-penetrating capabilities in different cell lines with low cytotoxicity.

The role of alpha-, 3(10)-, and pi-helix in helix-->coil transitions

The conformational equilibrium between 3(10)- and alpha-helical structure has been studied via high-resolution NMR spectroscopy by Millhauser and coworkers using the MW peptide Ac-AMAAKAWAAKA AAARA-NH2. Their 750-MHz nuclear Overhauser effect spectroscopy (NOESY) spectra were interpreted to reflect appreciable populations of 3(10)-helix throughout the peptide, with the greatest contribution at the N and C termini. The presence of simultaneous alphaN(i,i + 2) and alphaN(i,i + 4) NOE cross-peaks was proposed to represent conformational averaging between 3(10)- and alpha-helical structures. In this study, we describe 25-nsec molecular dynamics simulations of the MW peptide at 298 K, using both an 8 A and a 10 A force-shifted nonbonded cutoff. The ensemble averages of both simulations are in reasonable agreement with the experimental helical content from circular dichroism (CD), the (3)J(HNalpha) coupling constants, and the 57 observed NOEs. Analysis of the structures from both simulations revealed very little formation of contiguous i --> i + 3 hydrogen bonds (3(10)-helix); however, there was a large population of bifurcated i --> i + 3 and i --> i + 4 alpha-helical hydrogen bonds. In addition, both simulations contained considerable populations of pi-helix (i --> i + 5 hydrogen bonds). Individual turns formed over residues 1-9, which we predict contribute to the intensities of the experimentally observed alphaN(i,i + 2) NOEs. Here we show how sampling of both folded and unfolded structures can provide a structural framework for deconvolution of the conformational contributions to experimental ensemble averages.

Cation permeability and cation-anion interactions in a mutant GABA-gated chloride channel from Drosophila

To investigate the structural basis of anion selectivity of Drosophila GABA-gated Cl(-) channels, the permeation properties of wild-type and mutant channels were studied in Xenopus oocytes. This work focused on asparagine 319, which by homology is one amino acid away from a putative extracellular ring of charge that regulates cation permeation in nicotinic receptors. Mutation of this residue to aspartate reduced channel conductance, and mutation to lysine or arginine increased channel conductance. These results are consistent with an electrostatic interaction between this site and permeating anions. The lysine mutant, but not the arginine mutant, formed a channel that is permeable to cations, and this cannot be explained in terms of electrostatics. The lysine mutant had a 25-mV reversal potential in solutions with symmetrical Cl(-) and asymmetrical cations. The permeability ratio of K(+) to Cl(-) was determined as 0. 33 from reversal potential measurements in KCl gradients. Experiments with large organic cations and anions showed that cation permeation can only be seen in the presence of Cl(-), but Cl(-) permeation can be seen in the absence of permeant cations. Measurements of permeability ratios of organic anions indicated that the lysine mutant has an increased pore size. The cation permeability of the lysine-containing mutant channel cannot be accounted for by a simple electrostatic interaction with permeating ions. It is likely that lysine substitution causes a structural change that extends beyond this one residue to influence the positions of other channel-forming residues. Thus protein conformation plays an important role in enabling ion channels to distinguish between anions and cations.

Nitroxide side-chain dynamics in a spin-labeled helix-forming peptide revealed by high-frequency (139.5-GHz) EPR spectroscopy

High-frequency electron paramagnetic resonance (EPR) spectroscopy has been performed on a nitroxide spin-labeled peptide in fluid aqueous solution. The peptide, which follows the single letter sequence, was reacted with the methanethiosulfonate spin label at the cysteine sulfur. The spin sensitivity of high-frequency EPR is excellent with less than 20 pmol of sample required to obtain spectra with good signal-to-noise ratios. Simulation of the temperature-dependent spectral lineshapes reveals the existence of local anisotropic motion about the nitroxide N-O bond with a motional anisotropy tau( perpendicular)/tau( parallel) ( identical with N) approaching 2.6 at 306 K. Comparison with previous work on rigidly labeled peptides suggests that the spin label is reorienting about its side-chain tether. This study demonstrates the feasibility of performing 140-GHz EPR on biological samples in fluid aqueous solution.

The effect of mutations on peptide models of the DNA binding helix of p53: evidence for a correlation between structure and tumorigenesis

The tumor suppresser protein p53 has been called the "guardian of the genome." DNA damage induces p53 to either halt the cell cycle, allowing for repair, or initiate apoptosis. P53 is mutated in over 50% of human tumors and it has been proposed that many tumorigenic mutations are deleterious to p53 because they induce local unfolding. To explore this hypothesis, peptide models have been developed to study tumorigenic mutations in the H2 helix of the p53 core domain. This helix is rich with charged residues and is a key component of the DNA binding region. A 16-residue peptide corresponding to the H2 wild-type sequence extended with an Ala-rich C-terminus was synthesized and studied by 1H-nmr (500 MHz) and CD. The nmr studies demonstrate that this peptide adopts helical structure in solution. Six additional peptides corresponding to subtle tumorigenic mutations were synthesized and CD was used to assess the relative stability of these "mutant analogues." All six mutations studied are destabilizing relative to the wild type, with delta delta G values in the range of 0.26 to 1.35 kcal mol-1. Surprisingly, substitution of Asp 281 with Ala resulted in a peptide with the greatest destabilization even though Ala possesses the largest helix propensity of the common 20 amino acids. Because this helix appears to be stabilized mainly by local electrostatics, we conclude that its structure is susceptible to even the most conservative mutations. These results provide support for the hypothesis that tumorigenic mutations induce local unfolding of p53.

Addition of side-chain interactions to 3(10)-helix/coil and alpha-helix/3(10)-helix/coil theory

An increasing number of experimental and theoretical studies have demonstrated the importance of the 3(10)-helix/ alpha-helix/coil equilibrium for the structure and folding of peptides and proteins. One way to perturb this equilibrium is to introduce side-chain interactions that stabilize or destabilize one helix. For example, an attractive i, i + 4 interaction, present only in the alpha-helix, will favor the alpha-helix over 3(10), while an i, i + 4 repulsion will favor the 3(10)-helix over alpha. To quantify the 3(10)/alpha/coil equilibrium, it is essential to use a helix/coil theory that considers the stability of every possible conformation of a peptide. We have previously developed models for the 3(10)-helix/coil and 3(10)-helix/alpha-helix/ coil equilibria. Here we extend this work by adding i, i + 3 and i, i + 4 side-chain interaction energies to the models. The theory is based on classifying residues into alpha-helical, 3(10)-helical, or nonhelical (coil) conformations. Statistical weights are assigned to residues in a helical conformation with an associated helical hydrogen bond, a helical conformation with no hydrogen bond, an N-cap position, a C-cap position, or the reference coil conformation plus i, i + 3 and i, i + 4 side-chain interactions. This work may provide a framework for quantitatively rationalizing experimental work on isolated 3(10)-helices and mixed 3(10)-/alpha-helices and for predicting the locations and stabilities of these structures in peptides and proteins. We conclude that strong i, i + 4 side-chain interactions favor alpha-helix formation, while the 3(10)-helix population is maximized when weaker i, i + 4 side-chain interactions are present.

An NMR investigation of the conformational effect of nitroxide spin labels on Ala-rich helical peptides

Nitroxide spin labels, in conjunction with electron spin resonance (ESR) experiments, are extensively employed to probe the structure and dynamics of biomolecules. One of the most ubiquitous spin labeling reagents is the methanethiosulfonate spin label which attaches a spin label selectively to Cys residues via a disulfide bond (Cys-SL). However, the actual effect of the nitroxide spin label upon the conformation of the peptide or protein cannot be unambiguously determined by ESR. In this study, a series of 16-residue Ala-rich helical peptides was characterized by nuclear magnetic resonance techniques. The C alpha H chemical shift analysis, NOEs, and 3JNH alpha coupling constants for peptides with no Cys, free Cys, and Cys-SL (with the N-O group reduced) were compared. These results indicate that while replacement of an Ala with a Cys residue causes a loss of overall helical structure, the Cys-SL residue is helix supporting, as would be expected for a non-beta-branched aliphatic amino acid. Thus, the Cys-SL residue does not perturb helical structure and, instead, exhibits helix-stabilizing characteristics similar to that found for Ala, Met, and Leu.

Estimating the relative populations of 3(10)-helix and alpha-helix in Ala-rich peptides: a hydrogen exchange and high field NMR study

Recent experimental and theoretical work suggests that alanine-rich peptides fold as a mixture of 3(10)-helix (i --> i + 3 hydrogen bonding) and alpha-helix (i --> i + 4 hydrogen bonding). In order to assess the relative proportions of the two conformers, NMR studies were performed on the 16 residue sequences: Ac-AAAAKAAAAKAAAAKA-NH2 (3K) and Ac-AMAAKAWAAKAAAARA-NH2 (MW). Hydrogen/deuterium-exchange kinetics measured for the first three amide protons of the 3K peptide indicate that the NH of Ala3 is partially protected from exchange. This result is consistent with the presence of an i --> i + 3 hydrogen bond between the carbonyl group of the acetyl blocking group and the NH group of Ala3. The MW peptide is a modified version of the 3K peptide, designed to increase alphaH signal dispersion. 1H NMR spectra of the MW peptide at 750 MHz reveal a series of intermediate range (NOEs) consistent with a mixture of 3(10)-helix and alpha-helix. The relative intensities of the alphaN(i,i + 3) and alphabeta(i,i + 3) (nuclear Overhauser enhancements) NOEs suggest that 3(10)-helix is present throughout the peptide, but with the greatest contribution at the termini. A model was developed to determine the relative contributions of 3(10)-helix and alpha-helix. Lower bounds for the population of 3(10)-helix are approximately 50% at the termini and 25% in the middle of the peptide. The greatest alpha-helical content is between the middle of the peptide and the N terminus.

Lysine as helix C-capping residue in a synthetic peptide

The structure of the synthetic peptide CH3CO(Leu-Ser-Leu-Leu-Leu-Ser-Leu)3Lys-NH2 in trifluoroethanol/water 60/40 (volume ratio) was characterized by two-dimensional nmr spectroscopy. The peptide, closely related to the amphiphilic helix models designed by W. F. De-Grado and co-workers to mimic protein ion channels [(1988) Science, Vol. 240, p. 1177-1181], folds into a regular helix spanning residues 1-20. Evidence for a helix C-terminal capping conformation, involving the terminal lysine residue, was observed from Overhauser effects and checked for consistency by restrained molecular dynamics simulations. The side-chain amino group of Lys22 forms a hydrogen bond with the carbonyl of Leu18, and the distorted helical geometry of the terminal dipeptide allows the inclusion of a water bridge between the backbone NH of the Lys22 residue and the carbonyls of Leu19 and Ser20.

Models for the 3(10)-helix/coil, pi-helix/coil, and alpha-helix/3(10)-helix/coil transitions in isolated peptides

Models for the 3(10)-helix/coil and pi-helix/coil equilibria have been derived. The theory is based on classifying residues into helical or nonhelical (coil) conformations. Statistical weights are assigned to residues in a helical conformation with an associated helical hydrogen bond, a helical conformation with no hydrogen bond, an N-cap position, a C-cap position, or the reference coil conformation. The models for alpha-helix formation and 3(10)-helix formation have also been combined to describe a three-state equilibrium in which alpha-helical, 3(10)-helical, and coil conformations are populated. The results are compared with the modified Lifson-Roig theory for the alpha-helix/coil equilibrium. The comparison accounts for the experimental observations that 3(10)-helices tend to be short and pi-helices are not favored for any length. This work may provide a framework for quantitatively rationalizing experimental work on isolated 3(10)-helices and mixed 3(10)-/alpha-helices.

Determination of the structure of the N-terminal splice region of the cyclic AMP-specific phosphodiesterase RD1 (RNPDE4A1) by 1H NMR and identification of the membrane association domain using chimeric constructs

A 25-residue peptide representing the membrane targeting N-terminal splice region of the cyclic AMP phosphodiesterase RD1 (RNPDE4A1) was synthesized, and its structure was determined by 1H NMR. Two independently folding helical regions were identified, separated by a highly mobile "hinge" region. The first helical region was formed by an N-terminal amphipathic alpha-helix, and the second consisted of multiple overlapping turns and contained a distinct compact, hydrophobic, tryptophan-rich domain (residues 14-20). Chimeric molecules, formed between the N-terminal region of RD1 and the soluble bacterial protein chloramphenicol acetyltransferase, were used in an in vitro system to determine the features within the splice region that were required for membrane association. The ability of RD1-chloramphenicol acetyltransferase chimera to become membrane-associated was not affected by deletion of any of the following regions: the apolar section (residues 2-7) of the first helical region, the polar part of this region together with the hinge region (residues 8-13), or the polar end of the C-terminal helical region (residues 21-25). In marked contrast, deletion of the compact, hydrophobic tryptophan-rich domain (residues 14-20) found in the second helical region obliterated membrane association. Replacement of this domain with a hydrophobic cassette of seven alanine residues also abolished membrane association, indicating that membrane-association occurred by virtue of specific hydrophobic interactions with residues within the compact, tryptophan-rich domain. The structure of this domain is well defined in the peptide, and although the region is helical, both the backbone and the distribution of side chains are somewhat distorted as compared with an ideal alpha-helix. Hydrophobic interactions, such as the "stacked" rings of residues Pro14 and Trp15, stabilize this domain with the side chain of residue Leu16 adopting a central position, interacting with the side chains of all three tryptophan residues 15, 19, and 20. These bulky side chains thus form a hydrophobic cluster. In contrast, the side chain of residue Val17 is relatively exposed, pointing out from the opposite "face" of the peptide. Although it appears that this compact, tryptophan-rich domain is responsible for membrane association, at present the target site and hence the specific interactions involved in membrane targeting by the RD1 splice region remain unidentified.

Detailed description of an alpha helix-->pi bulge transition detected by molecular dynamics simulations of the p185c-erbB2 V659G transmembrane domain

Molecular dynamics simulations of a 29-residue peptide including the transmembrane domain of the V659G mutant of the c-erbB2 protein demonstrate important dynamical behavior. Although the alpha helix is the structure commonly assumed for transmembrane hydrophobic segments, we found that hydrogen bond rearrangements can occur, giving rise to a structural deformation termed pi bulge stabilized by successive hydrogen bonds of pi helix type. A series of simulations enables us to give a detailed description, at the atomic level, of the alpha helix->pi bulge transition. The major consequence of this deformation covering one and a half turn of helix results in a noticeable shift around the helix axis of the C-Terminal residues relatively to those of the N-terminus. Such a deformation closely related to structural motifs described in the literature, induces a change in the distribution of the residues along the helix faces which could modulate the protein activity mediated by a dimerization process.

The stability of transmembrane helices: a molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembrane helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of 7 helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge on the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP.

Competing interactions contributing to alpha-helical stability in aqueous solution

The stability of a 15-residue peptide has been investigated using CD spectroscopy and molecular simulation techniques. The sequence of the peptide was designed to include key features that are known to stabilize alpha-helices, including ion pairs, helix dipole capping, peptide bond capping, and aromatic interactions. The degree of helicity has been determined experimentally by CD in three solvents (aqueous buffer, methanol, and trifluoroethanol) and at two temperatures. Simulations of the peptide in the aqueous system have been performed over 500 ps at the same two temperatures using a fully explicit solvent model. Consistent with the CD data, the degree of helicity is decreased at the higher temperature. Our analysis of the simulation results has focused on competition between different side-chain/side-chain and side-chain/main-chain interactions, which can, in principle, stabilize the helix. The unfolding in aqueous solution occurs at the amino terminus because the side-chain interactions are insufficient to stabilize both the helix dipole and the peptide hydrogen bonds. Loss of capping of the peptide backbone leads to water insertion within the first peptide hydrogen bond and hence unfolding. In contrast, the carboxy terminus of the alpha-helix is stable in both simulations because the C-terminal lysine residue stabilizes the helix dipole, but at the expense of an ion pair.