The role of plastic beta-hairpin and weak hydrophobic core in the stability and unfolding of a full sequence design protein

In this study, the thermal stability of a designed alpha/beta protein FSD (full sequence design) was studied by explicit solvent simulations at three moderate temperatures, 273 K, 300 K, and 330 K. The average properties of the ten trajectories at each temperature were analyzed. The thermal unfolding, as judged by backbone root-mean-square deviation and percentage of native contacts, was displayed with increased sampling outside of the native basin as the temperature was raised. The positional fluctuation of the hairpin residues was significantly higher than that of the helix residues at all three temperatures. The hairpin segment displayed certain plasticity even at 273 K. Apart from the terminal residues, the highest fluctuation was shown in the turn residues 7-9. Secondary structure analysis manifested the structural heterogeneity of the hairpin segment. It was also revealed by the simulation that the hydrophobic core was vulnerable to thermal denaturation. Consistent with the experiment, the I7Y mutation in the double mutant FSD-EY (FSD with mutations Q1E and I7Y) dramatically increased the protein stability in the simulation, suggesting that the plasticity of the hairpin can be partially compensated by a stronger hydrophobic core. As for the unfolding pathway, the breathing of the hydrophobic core and the separation of the two secondary structure elements (alpha helix and beta hairpin) was the initiation step of the unfolding. The loss of global contacts from the separation further destabilized the hairpin structure and also led to the unwinding of the helix.     

Interplay of secondary structures and side-chain contacts in the denatured state of BBA1

The denatured state of a miniprotein BBA1 is studied under the native condition with the AMBER/Poisson-Boltzmann energy model and with the self-guided enhanced sampling technique. Forty independent trajectories are collected to sample the highly diversified denatured structures. Our simulation data show that the denatured BBA1 contains high percentage of native helix and native turn, but low percentage of native hairpin. Conditional population analysis indicates that the native helix formation and the native hairpin formation are not cooperative in the denatured state. Side-chain analysis shows that the native hydrophobic contacts are more preferred than the non-native hydrophobic contacts in the denatured BBA1. In contrast, the salt-bridge contacts are more or less nonspecific even if their populations are higher than those of hydrophobic contacts. Analysis of the trajectories shows that the native helix mostly initiates near the N terminus and propagates to the C terminus, and mostly forms from 3(10)-helix/turn to alpha helix. The same analysis shows that the native turn is important but not necessary in its formation in the denatured BBA1. In addition, the formations of the two strands in the native hairpin are rather asymmetric, demonstrating the likely influence of the protein environment. Energetic analysis shows that the native helix formation is largely driven by electrostatic interactions in denatured BBA1. Further, the native helix formation is associated with the breakup of non-native salt-bridge contacts and the accumulation of native salt-bridge contacts. However, the native hydrophobic contacts only show a small increase upon the native helix formation while the non-native hydrophobic contacts stay essentially the same, different from the evolution of hydrophobic contacts observed in an isolated helix folding.     

Tertiary DNA structure in the single-stranded hTERT promoter fragment unfolds and refolds by parallel pathways via cooperative or sequential events

The discovery of G-quadruplexes and other DNA secondary elements has increased the structural diversity of DNA well beyond the ubiquitous double helix. However, it remains to be determined whether tertiary interactions can take place in a DNA complex that contains more than one secondary structure. Using a new data analysis strategy that exploits the hysteresis region between the mechanical unfolding and refolding traces obtained by a laser-tweezers instrument, we now provide the first convincing kinetic and thermodynamic evidence that a higher order interaction takes place between a hairpin and a G-quadruplex in a single-stranded DNA fragment that is found in the promoter region of human telomerase. During the hierarchical unfolding or refolding of the DNA complex, a 15-nucleotide hairpin serves as a common species among three intermediates. Moreover, either a mutant that prevents this hairpin formation or the addition of a DNA fragment complementary to the hairpin destroys the cooperative kinetic events by removing the tertiary interaction mediated by the hairpin. The coexistence of the sequential and the cooperative refolding events provides direct evidence for a unifying kinetic partition mechanism previously observed only in large proteins and complex RNA structures. Not only does this result rationalize the current controversial observations for the long-range interaction in complex single-stranded DNA structures, but also this unexpected complexity in a promoter element provides additional justification for the biological function of these structures in cells.     

Proteins that convert from alpha helix to beta sheet: implications for folding and disease

The sequence of a protein normally determines which amino acid residues will form alpha helices, and which one beta sheets, to an extent that allows secondary structure prediction to be made with a reasonable reliability. Nevertheless, non-native helical structures are observed during in vitro folding of several model proteins and may even occur during protein biosynthesis within the ribosomal exit tunnel. Moreover, non-native beta sheet structures are common in amyloid fibrils formed by a variety of pathogenic and even non-pathogenic proteins and peptides. In all of these cases, the formation of alpha helix precedes the appearance of beta sheet, which suggests that conversion from the simpler, more local helix structure to the often more convoluted sheet architecture during folding and pathogenic misfolding processes could be a unifying principle of general importance. A better understanding of this switching process, and the ability to design molecular systems which can be induced to switch between these conformations will have a significant impact on fields ranging from fundamental biochemistry through to applied technology and medicine.     

Probing the formation of stable tertiary structure in a model miniprotein at atomic resolution: determinants of stability of a helical hairpin

The minimal model system to study the basic principles of protein folding is the hairpin. The formation of beta-hairpins, which are the basic components of antiparallel beta-sheets, has been studied extensively in the past decade, but much less is known about helical hairpins. Here, we probe hairpin formation between a polyproline type-II helix and an alpha-helix as present in the natural miniprotein peptide YY (PYY). Both turn sequence and interactions of aromatic side chains from the C-terminal alpha-helix with the pockets formed by N-terminal Pro residues are shown by site-directed mutagenesis and solution NMR spectroscopy in different solvent systems to be important determinants of backbone dynamics and hairpin stability, suggesting a close analogy with some beta-hairpin structures. It is shown that multiple relatively weak contacts between the helices are necessary for the formation of the helical hairpin studied here, whereas the type-I beta-turn acts like a hinge, which through certain single amino acid substitutions is destabilized such that hairpin formation is completely abolished. Denaturation and renaturation of tertiary structure by temperature or cosolvents were probed by measuring changes of chemical shifts. Folding of PYY is both reversible and cooperative as inferred from the sigmoidal denaturation curves displayed by residues at the interface of the helical hairpin. Such miniproteins thus feature an important hallmark of globular proteins and should provide a convenient system to study basic aspects of helical hairpin folding that are complementary to those derived from studies of beta-hairpins.     

Design, synthesis and structural investigations of a beta-peptide forming a 314-helix stabilized by electrostatic interactions

Two different strategies have been employed for the synthesis of Fmoc-protected beta(3)-homoarginine; the Arndt-Eistert homologation of alpha-arginine and the guanidinylation of beta(3)-homoornithine. Solid-phase beta-peptide synthesis was used for the preparation of beta-heptapeptide 1, which was designed to form a helix stabilized by electrostatic interactions through positively (beta(3)hArg) and negatively charged (beta(3)hGlu) amino acid residues. CD measurements and corresponding NMR investigations in MeOH and aqueous solutions do indeed show that the beta-peptidic 3(14)-helix can be stabilized by salt-bridge formation.     

Micelle‐Triggered β‐Hairpin to α‐Helix Transition in a 14‐Residue Peptide from a Choline‐Binding Repeat of the Pneumococcal Autolysin LytA

Choline‐binding modules (CBMs) have a ββ‐solenoid structure composed of choline‐binding repeats (CBR), which consist of a β‐hairpin followed by a short linker. To find minimal peptides that are able to maintain the CBR native structure and to evaluate their remaining choline‐binding ability, we have analysed the third β‐hairpin of the CBM from the pneumococcal LytA autolysin. Circular dichroism and NMR data reveal that this peptide forms a highly stable native‐like β‐hairpin both in aqueous solution and in the presence of trifluoroethanol, but, strikingly, the peptide structure is a stable amphipathic α‐helix in both zwitterionic (dodecylphosphocholine) and anionic (sodium dodecylsulfate) detergent micelles, as well as in small unilamellar vesicles. This β‐hairpin to α‐helix conversion is reversible. Given that the β‐hairpin and α‐helix differ greatly in the distribution of hydrophobic and hydrophilic side chains, we propose that the amphipathicity is a requirement for a peptide structure to interact and to be stable in micelles or lipid vesicles. To our knowledge, this “chameleonic” behaviour is the only described case of a micelle‐induced structural transition between two ordered peptide structures.     

Dynamic folding pathway models of the villin headpiece subdomain (HP‐36) structure

We have investigated the folding pathway of the 36‐residue villin headpiece subdomain (HP‐36) by action‐derived molecular dynamics simulations. The folding is initiated by hydrophobic collapse, after which the concurrent formation of full tertiary structure and α‐helical secondary structure is observed. The collapse is observed to be associated with a couple of specific native contacts contrary to the conventional nonspecific hydrophobic collapse model. Stable secondary structure formation after the collapse suggests that the folding of HP‐36 follows neither the framework model nor the diffusion‐collision model. The C‐terminal helix forms first, followed by the N‐terminal helix positioned in its native orientation. The short middle helix is shown to form last. Signs for multiple folding pathways are also observed. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations

Reaching the native states of small proteins, a necessary step towards a comprehensive understanding of the folding mechanisms, has remained a tremendous challenge to ab initio protein folding simulations despite the extensive effort. In this work, the folding process of the B domain of protein A (BdpA) has been simulated by both conventional and replica exchange molecular dynamics using AMBER FF03 all-atom force field. Started from an extended chain, a total of 40 conventional (each to 1.0 micros) and two sets of replica exchange (each to 200.0 ns per replica) molecular dynamics simulations were performed with different generalized-Born solvation models and temperature control schemes. The improvements in both the force field and solvent model allowed successful simulations of the folding process to the native state as demonstrated by the 0.80 A C(alpha) root mean square deviation (RMSD) of the best folded structure. The most populated conformation was the native folded structure with a high population. This was a significant improvement over the 2.8 A C(alpha) RMSD of the best nativelike structures from previous ab initio folding studies on BdpA. To the best of our knowledge, our results demonstrate, for the first time, that ab initio simulations can reach the native state of BdpA. Consistent with experimental observations, including Phi-value analyses, formation of helix II/III hairpin was a crucial step that provides a template upon which helix I could form and the folding process could complete. Early formation of helix III was observed which is consistent with the experimental results of higher residual helical content of isolated helix III among the three helices. The calculated temperature-dependent profile and the melting temperature were in close agreement with the experimental results. The simulations further revealed that phenylalanine 31 may play critical to achieve the correct packing of the three helices which is consistent with the experimental observation. In addition to the mechanistic studies, an ab initio structure prediction was also conducted based on both the physical energy and a statistical potential. Based on the lowest physical energy, the predicted structure was 2.0 A C(alpha) RMSD away from the experimentally determined structure.     

Simulation studies of amide I IR absorption and two-dimensional IR spectra of beta hairpins in liquid water

Amide I IR absorption and two-dimensional (2D) IR photon echo spectra of a model beta hairpin in aqueous solution are theoretically studied and simulated by combining semiempirical quantum chemistry calculations and molecular dynamics simulation methods. The instantaneous normal-mode analysis of the beta hairpin in solution is performed to obtain the density of states and the inverse participation ratios of the one-exciton states. The motional and exchange narrowing processes are taken into account by employing the time-correlation function theory for the linear and nonlinear response functions. Numerically simulated IR absorption and 2D spectra are then found to be determined largely by the amide I normal modes delocalized on the peptides in the two strands. The site-specific isotope-labeling effects on the IR and 2D IR spectra are discussed. The simulation results for the ideal (A17) beta hairpin are directly compared with those of the realistic 16-residue (GB1) beta hairpin from an immunoglobulin G-binding protein. It was found that the characteristic features in IR and 2D spectra of both the ideal (A17) beta hairpin and the GB1 beta hairpin are the same. The simulated IR spectrum of the GB1 beta hairpin is found to be in good agreement with experiment, which demonstrates that the present computational method is quantitatively reliable.     

Folding mechanisms of individual beta-hairpins in a Go model of Pin1 WW domain by all-atom molecular dynamics simulations

This paper examines the folding mechanism of an individual beta-hairpin in the presence of other hairpins by using an off-lattice model of a small triple-stranded antiparallel beta-sheet protein, Pin1 WW domain. The turn zipper model and the hydrophobic collapse model originally developed for a single beta-hairpin in literature is confirmed to be useful in describing beta-hairpins in model Pin1 WW domain. We find that the mechanism for folding a specific hairpin is independent of whether it folds first or second, but the formation process are significantly dependent on temperature. More specifically, beta1-beta2 hairpin folds via the turn zipper model at a low temperature and the hydrophobic collapse model at a high temperature, while the folding of beta2-beta3 hairpin follows the turn zipper model at both temperatures. The change in folding mechanisms is interpreted by the interplay between contact stability (enthalpy) and loop lengths (entropy), the effect of which is temperature dependent.     

The outlined polymorph selection of iPB-1 revealing by incorporation of a dentritic polymer

By incorporating a hydrophilic dentritic polyester into polymorphic isotactic polybutene-1 (iPB-1), we successfully decreased its viscosity as designed to trigger the melt crystallisation of form III due to the proposed decrease of the interactions between the iPB-1 chain sequences. In addition to the usual form II formation, the form III formation resulting from cooling the iPB-1 melt at suitable cooling rates has been verified by the DSC, WAXD, and in-situ synchrotron WAXS measurements. It was proposed that the form III melt crystallisation occurs after the 4/1 helix conformation formation in the iPB-1 melt and the following alignment of the resulting 4/1 helices. Once the interactions between the iPB-1 chain sequences are sufficiently strong as in the usual iPB-1 melt, the collapse of the 4/1 into the 11/3 helix conformation would occur and finally trigger the generally observed melt crystallisation of form II. Thus, we first outlined the polymorph selection process of iPB-1 based on the helix conformation formation and the following alignment of the helices which finally result in crystallisation of iPB-1 into corresponding crystal forms. The results are helpful for understanding of both the polymorph selection of polymorphic polymers like iPB-1 and the general polymer crystallisation.     

Probing the helical content of growth hormone-releasing factor analogs using electrospray ionization mass spectrometry

A series of growth hormone-releasing factor analogs have been studied by both circular dichroism and electrospray ionization mass spectrometry (ESI/MS). The peptides are 32 residues long and are known to adopt a random-coil structure in aqueous solution but become increasingly helical as the proportion of organic solvent is increased. Deuterium exchange was observed as an increase in mass of the peptide, as measured by ESI/MS. Rates of exchange were measured and half-lives calculated for analogs containing amino acid substitutions designed to promote or discourage helix formation. Exchange was slower in peptides that are helical (as shown by circular dichroism) than in randomly coiled peptides. Solution conditions that favor helix formation also produced slower exchange rates. These studies suggest that ESI/MS can provide date about the extent and stability of helix formation.     

DNA polymorphism as an origin of adenine-thymine tract length-dependent threading intercalation rate

Binuclear ruthenium complexes that bind DNA by threading intercalation have recently been found to exhibit an exceptional kinetic selectivity for long polymeric adenine-thymine (AT) DNA. A series of oligonucleotide hairpin duplexes containing a central tract of 6-44 alternating AT base pairs have here been used to investigate the nature of the recognition mechanism. We find that, above a threshold AT tract length corresponding to one helix turn of B-DNA, a dramatic increase in threading intercalation rate occurs. In contrast, such length dependence is not observed for rates of unthreading. Intercalation by any mechanism that depends on the open end of the hairpin was found not to be important in the series of oligonucleotides used, as verified by including in the study a hairpin duplex cyclized by a copper-catalyzed "click" reaction. Our observations are interpreted in terms of a conformational pre-equilibrium, determined by the length of the AT tract. We finally find that mismatches or loops in the oligonucleotide facilitate the threading process, of interest for the development of mismatch-recognizing probes.     

Thioether‐derived Macrocycle for Peptide Secondary Structure Fixation

Recently, we developed methods to stabilize peptides into various secondary structures, including α‐helix, type III turn and β‐hairpin via proper thioether based macrocyclization. These conformationally constrained peptidomimetics confer enhanced biophysical properties and provide a valuable avenue towards clinically‐relevant therapeutic molecules. In this personal account, thioether‐derived macrocyclization methods developed by our group for stabilization of α‐helix, type‐III β turn and β‐hairpin conformations are discussed.     

Experimental resolution of early steps in protein folding: testing molecular dynamics simulations

Time-resolved Tyr fluorescence spectroscopy coupled with a laser-induced temperature-jump (T-jump) was employed to follow the folding relaxation dynamics of the B-domain of Staphylococcal protein A. The single Tyr is located in helix 1 (H1) and is a sensitive probe of the structure of this helix and the overall helical bundle structure. The results from this study were compared to those from a complementary infrared T-jump study on this protein [Vu, D. M.; Myers, J. K.; Oas, T. G.; Dyer, R. B. Biochemistry 2004, 43, 3582]. Both methods detect a microsecond process that follows the cooperative relaxation of the helical bundle core. However, a fast process (10-7 s) that follows the relaxation of the individual helices was observed only with the infrared probe. Thus, fast formation of H1 is not observed, but rather H1 forms in the microsecond phase, concomitantly with the docking to (and stabilization by) the other two helices to form the helical bundle structure. This observation validates the results of several previous molecular dynamics simulations that predict H1 formation only in the final assembly of the helix bundle.     

A test on peptide stability of AMBER force fields with implicit solvation

We used replica exchange molecular dynamics (REMD) simulations to evaluate four different AMBER force fields and three different implicit solvent models. Our aim was to determine if these physics-based models captured the correct secondary structures of two alpha-helical and two beta-peptides: the 14-mer EK helix of Baldwin and co-workers, the C-terminal helix of ribonuclease, the 16-mer C-terminal hairpin of protein G, and the trpzip2 miniprotein. The different models gave different results, but generally we found that AMBER ff96 plus the implicit solvent model of Onufriev, Bashford, and Case gave reasonable structures, and is fairly well-balanced between helix and sheet. We also observed differences in the strength of ion pairing in the solvent models, we but found that the native secondary structures were retained even when salt bridges were prevented in the conformational sampling. Overall, this work indicates that some of these all-atom physics-based force fields may be good starting points for protein folding and protein structure prediction.     

Modulation of the Thermodynamic Signatures of an RNA Thermometer by Osmolytes and Salts

Folding of ribonucleic acids (RNAs) is driven by several factors, such as base pairing and stacking, chain entropy, and ion‐mediated electrostatics, which have been studied in great detail. However, the power of background molecules in the cellular milieu is often neglected. Herein, we study the effect of common osmolytes on the folding equilibrium of a hairpin‐structured RNA and, using pressure perturbation, provide novel thermodynamic and volumetric insights into the modulation mechanism. The presence of TMAO causes an increased thermal stability and a more positive volume change for the helix‐to‐coil transition, whereas urea destabilizes the hairpin and leads to an increased expansibility of the unfolded state. Further, we find a strong interplay between water, salt, and osmolyte in driving the thermodynamics and defining the temperature and pressure stability limit of the RNA. Our results support a universal working mechanism of TMAO and urea to (de)stabilize proteins and the RNA.     

Explanation of Ionic Sequences in Various Phenomena. VI. The Formation and Destruction of Helices in Bovine Plasma Albumin in Neutral and Acidic Solutions

Abstract

The structure of bovine plasma albumin (BPA) was examined by optical rotatory dispersion studies at both low (pH 1.5 and 2.0) and high (pH 9.0) pH values in various aqueous salt solutions. The resulting cationic sequences were compared to those observed by Pedersen for values of the sedimentation constant. At pH 9.0 the destruction of the “helix” produces an acidic sequence. The relative pH position of the “helix” transition, the fact that addition of salt increases the apparent helical content of BPA, and the observed acidic-type of sequence rule out the possibility of (1)ionic bonds between carboxylate and ?-amino groups, (2) hydrophobic bonds, or (3) hydrogen bonds between peptide linkages as major contributing forces in the formation of the helix. The stability of the “helix” in BPA between pH 3.0 and 9.0 must therefore be due to hydrogen bonds between carboxylate ions and hydroxyl groups such as those of serine, threonine, and tyrosine. Repulsive forces between the positively charged groups on BPA strengthen these bonds by preventing the expanded form of BPA from collapsing. At pH 2.0 two types of sequences were observed: The s⁰ _{20, w}, [α]_D and a₀ values gave an acidic-type cationic sequence. The b₀ (helix content), λ_c and [λ]₂₃₃ values gave essentially a nonpolar sequence. The nonpolar or hydrophobic salting-out sequences show that the formation of hydrophobic bonds at pH 2.0 hinders the formation of the helix or folded structure. The acidic sequences show that hydrogen bonds between carboxylic acid groups stabilize both the apparent helix or helices and the intermolecular aggregation of the BPA molecules. From a comparison of the S⁰ ₂₀,w values and the helical content of BPA at pH 9.0 it is also concluded that the formation of these apparent helices or folded structures expands or stiffens the BPA molecule.     

Supramolecular Chemistry of Cyclodextrin-Peptide Hybrids: Azobenzene-Tagged Peptides

AC17, which is composed of 17 amino acids and has an azobenzene moiety but has no cyclodextrin (CD) unit in the side chain, exhibits 54% helix content. However, AC17, which has both trans-azobenzene and -CD, shows 82% helix content. This result suggests that the helix structure is stabilized by host (CD)-guest (azobenzene) bridge in the side chain of the peptide. The helix content changed by trans-cis photoisomerization as shown by 64% helix content for AC17 in its cis form. This result suggests that cis-azobenzene unit is excluded from the -CD cavity, thus resulting in the smaller helix content. The helix contents for AC17, which has both azobenzene and -CD, are 94% in the cis form and 87% in the trans form, suggesting that the cis form is included in the -CD cavity. Azobenzene-tagged CD-peptide hybrids with histidine unit were also prepared and photoregulation of catalytic activity in ester hydrolysis was examined.