Analyses of the folding properties of ferredoxin‐like fold proteins by means of a coarse‐grained Gō model: Relationship between the free energy profiles and folding cores

The folding mechanisms of proteins with multi‐state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on‐going research issues. In this study, we investigate ferredoxin‐like fold proteins which have a simple topology and multi‐state transitions. We analyze the folding processes by means of a coarse‐grained Gō model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high‐free‐energy intermediate state, and ADA2h and S6, which fold into the native structure through two‐state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54‐55, are reproduced and compared with experimental observations. We show that the ferredoxin‐like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two‐state mechanism. Proteins 2014; 82:954–965. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulations of a winter flounder “antifreeze” polypeptide in aqueous solution

A winter flounder antifreeze polypeptide (HPLC-6) has been studied in vacuo and in aqueous solution using molecular dynamics computer simulation techniques. The helical conformation of this polypeptide was found to be stable both in vacuum and in solution. The major stabilizing interactions were found to be the main-chain hydrogen bonds, a salt-bridge interaction, and solute–solvent hydrogen bonds. A significant bending in the middle of the polypeptide chain was observed both in vacuo and in solvent at 300 K. Possible causes of the bending are discussed. From simulations of mutant polypeptide molecules in vacuo, it is concluded that the bend in the native polypeptide was caused by side chain to backbone hydrogen bond competition involving the Thr 24 side chain and facilitated by strains on the helix resulting from the Lys 18-Glu 22 salt bridge. © 1993 John Wiley & Sons, Inc.     

A rare polyglycine type II‐like helix motif in naturally occurring proteins

Common structural elements in proteins such as α‐helices or β‐sheets are characterized by uniformly repeating, energetically favorable main chain conformations which additionally exhibit a completely saturated hydrogen‐bonding network of the main chain NH and CO groups. Although polyproline or polyglycine type II helices (PP_II or PG_II) are frequently found in proteins, they are not considered as equivalent secondary structure elements because they do not form a similar self‐contained hydrogen‐bonding network of the main chain atoms. In this context our finding of an unusual motif of glycine‐rich PG_II‐like helices in the structure of the acetophenone carboxylase core complex is of relevance. These PG_II‐like helices form hexagonal bundles which appear to fulfill the criterion of a (largely) saturated hydrogen‐bonding network of the main‐chain groups and therefore may be regarded in this sense as a new secondary structure element. It consists of a central PG_II‐like helix surrounded by six nearly parallel PG_II‐like helices in a hexagonal array, plus an additional PG_II‐like helix extending the array outwards. Very related structural elements have previously been found in synthetic polyglycine fibers. In both cases, all main chain NH and CO groups of the central PG_II‐helix are saturated by either intra‐ or intermolecular hydrogen‐bonds, resulting in a self‐contained hydrogen‐bonding network. Similar, but incomplete PG_II‐helix patterns were also previously identified in a GTP‐binding protein and an antifreeze protein.     

Helix‐sheet packing in proteins

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Although much is known about the packing geometry observed between α‐helices and between β‐sheets, there has been little progress on characterizing helix–sheet interactions. We present an analysis of the conformation of αβ₂ motifs in proteins, corresponding to all occurrences of helices in contact with two strands that are hydrogen bonded. The geometry of the αβ₂ motif is characterized by the azimuthal angle θ between the helix axis and an average vector representing the two strands, the elevation angle ψ between the helix axis and the plane containing the two strands, and the distance D between the helix and the strands. We observe that the helix tends to align to the two strands, with a preference for an antiparallel orientation if the two strands are parallel; this preference is diminished for other topologies of the β‐sheet. Side‐chain packing at the interface between the helix and the strands is mostly hydrophobic, with a preference for aliphatic amino acids in the strand and aromatic amino acids in the helix. From the knowledge of the geometry and amino acid propensities of αβ₂ motifs in proteins, we have derived different statistical potentials that are shown to be efficient in picking native‐like conformations among a set of non‐native conformations in well‐known decoy datasets. The information on the geometry of αβ₂ motifs as well as the related statistical potentials have applications in the field of protein structure prediction. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Exhaustive Metropolis Monte Carlo sampling and analysis of polyalanine conformations adopted under the influence of hydrogen bonds

We propose a novel Metropolis Monte Carlo procedure for protein modeling and analyze the influence of hydrogen bonding on the distribution of polyalanine conformations. We use an atomistic model of the polyalanine chain with rigid and planar polypeptide bonds, and elastic alpha carbon valence geometry. We adopt a simplified energy function in which only hard-sphere repulsion and hydrogen bonding interactions between the atoms are considered. Our Metropolis Monte Carlo procedure utilizes local crankshaft moves and is combined with parallel tempering to exhaustively sample the conformations of 16-mer polyalanine. We confirm that Flory's isolated-pair hypothesis (the steric independence between the dihedral angles of individual amino acids) does not hold true in long polypeptide chains. In addition to 3(10)- and alpha-helices, we identify a kink stabilized by 2 hydrogen bonds with a shared acceptor as a common structural motif. Varying the strength of hydrogen bonds, we induce the helix-coil transition in the model polypeptide chain. We compare the propensities for various hydrogen bonding patterns and determine the degree of cooperativity of hydrogen bond formation in terms of the Hill coefficient. The observed helix-coil transition is also quantified according to Zimm-Bragg theory.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

The effect of simple shear flow on the helix–coil transition of poly‐L‐lysine

Studies of the helix‐to‐coil transition in dilute solutions of poly‐L ‐lysine, dissolved in mixtures of water and methanol (MeOH), have been carried under shear flow using flow birefringence and modulated polarimetry. The fraction of helical conformations in a given solution remains independent of shear rate for MeOH concentrations above and below the critical value for the helix‐coil transition (i.e., 87.5% MeOH). For the 87.5% MeOH solutions, a shear‐induced helix‐to‐“stretched” coil transition occurs above a critical shear rate. Induction times for the transition show a temperature and shear rate dependence that can be described in terms of an activated jump process. Measurements of circular birefringence on cessation of flow also show that the transition is reversible, with the stretched coil reverting to the helical state on a time scale of several seconds. The activation energy for the jump process is found to be 16.2 kJ/mole. © 1999 John Wiley & Sons, Inc. Biopoly 50: 589–594, 1999     

Adsorption of polypeptides on the solid surface. II. Effect of secondary chain structure and helix–coil transition

The theory of adsorption of semistiff chains on a planar surface developed by the authors previously has been used to consider the helix–coil transition in single-stranded macromolecule interacting with an adsorbent plane. The cases of nonselective interaction when the adsorption energy is independent of the unit conformation (a) and selective interaction with only helical (b) or coiled (c) sequences active in adsorption were investigated. In case (b) the existence of secondary structure favors chain bonding to the surface. This leads to the increase in the stability of the helical state and complete polypeptide chain spiralization. The profile of the conformational helix–coil transition acquires an asymmetrical shape inherent to the second-order phase transition. In case (c) the bonding of a partially helical chain to the surface is similar to the adsorption of Gaussian coils and is accompanied by the destruction of secondary structure, this destruction being appreciable even if the helical state in space was favorable.     

Comparison of native and non‐native ubiquitin oligomers reveals analogous structures and reactivities

Ubiquitin (Ub) chains regulate a wide range of biological processes, and Ub chain connectivity is a critical determinant of the many regulatory roles that this post‐translational modification plays in cells. To understand how distinct Ub chains orchestrate different biochemical events, we and other investigators have developed enzymatic and non‐enzymatic methods to synthesize Ub chains of well‐defined length and connectivity. A number of chemical approaches have been used to generate Ub oligomers connected by non‐native linkages; however, few studies have examined the extent to which non‐native linkages recapitulate the structural and functional properties associated with native isopeptide bonds. Here, we compare the structure and function of Ub dimers bearing native and non‐native linkages. Using small‐angle X‐ray scattering (SAXS) analysis, we show that scattering profiles for the two types of dimers are similar. Moreover, using an experimental structural library and atomistic simulations to fit the experimental SAXS profiles, we find that the two types of Ub dimers can be matched to analogous structures. An important application of non‐native Ub oligomers is to probe the activity and selectivity of deubiquitinases. Through steady‐state kinetic analyses, we demonstrate that different families of deubiquitinases hydrolyze native and non‐native isopeptide linkages with comparable efficiency and selectivity. Considering the significant challenges associated with building topologically diverse native Ub chains, our results illustrate that chains harboring non‐native linkages can serve as surrogate substrates for explorations of Ub function.     

Helix stability of oligoglycine,oligoalanine, and oligo‐β‐alanine dodecamers reflected by hydrogen‐bond persistence

Helices are important structural/recognition elements in proteins and peptides. Stability and conformational differences between helices composed of α‐ and β‐amino acids as scaffolds for mimicry of helix recognition has become a theme in medicinal chemistry. Furthermore, helices formed by β‐amino acids are experimentally more stable than those formed by α‐amino acids. This is paradoxical because the larger sizes of the hydrogen‐bonding rings required by the extra methylene groups should lead to entropic destabilization. In this study, molecular dynamics simulations using the second‐generation force field, AMOEBA (Ponder, J.W., et al., Current status of the AMOEBA polarizable force field. J Phys Chem B, 2010. 114 (8): p. 2549–64.) explored the stability and hydrogen‐bonding patterns of capped oligo‐β‐alanine, oligoalanine, and oligoglycine dodecamers in water. The MD simulations showed that oligo‐β‐alanine has strong acceptor+2 hydrogen bonds, but surprisingly did not contain a large content of 3₁₂‐helical structures, possibly due to the sparse distribution of the 3₁₂‐helical structure and other structures with acceptor+2 hydrogen bonds. On the other hand, despite its backbone flexibility, the β‐alanine dodecamer had more stable and persistent <3.0 Å hydrogen bonds. Its structure was dominated more by multicentered hydrogen bonds than either oligoglycine or oligoalanine helices. The 3₁ (PII) helical structure, prevalent in oligoglycine and oligoalanine, does not appear to be stable in oligo‐β‐alanine indicating its competition with other structures (stacking structure as indicated by MD analyses). These differences are among the factors that shape helical structural preferences and the relative stabilities of these three oligopeptides. Proteins 2014; 82:3043–3061. © 2014 Wiley Periodicals, Inc.     

Carbon-13 and proton NMR studies of helix-coil transition of poly(gamma-benzyl-L-glutamate).

The helix–coil conformational transition undergone by poly(γ-benzyl-L -glutamate) in solutions of trifluoroacetic acid and deuterated chloroform was studied by proton and carbon-13 nmr. The results indicate that in the case of the solvent-induced helix–coil transition, the side chain assumes a helical conformation before the backbone. In the thermally induced helix–coil transition, the results indicate the existence of an intermediate state, which is between the α-helix and random coil and is free from intramolecular hydrogen bonding.     

Predicting novel trophic interactions in a non‐native world

Humans are altering the global distributional ranges of plants, while their co‐evolved herbivores are frequently left behind. Native herbivores often colonise non‐native plants, potentially reducing invasion success or causing economic loss to introduced agricultural crops. We developed a predictive model to forecast novel interactions and verified it with a data set containing hundreds of observed novel plant–insect interactions. Using a food network of 900 native European butterfly and moth species and 1944 native plants, we built an herbivore host‐use model. By extrapolating host use from the native herbivore–plant food network, we accurately forecasted the observed novel use of 459 non‐native plant species by native herbivores. Patterns that governed herbivore host breadth on co‐evolved native plants were equally important in determining non‐native hosts. Our results make the forecasting of novel herbivore communities feasible in order to better understand the fate and impact of introduced plants.     

Free‐energy function based on an all‐atom model for proteins

We have developed a free‐energy function based on an all‐atom model for proteins. It comprises two components, the hydration entropy (HE) and the total dehydration penalty (TDP). Upon a transition to a more compact structure, the number of accessible configurations arising from the translational displacement of water molecules in the system increases, leading to a water‐entropy gain. To fully account for this effect, the HE is calculated using a statistical‐mechanical theory applied to a molecular model for water. The TDP corresponds to the sum of the hydration energy and the protein intramolecular energy when a fully extended structure, which possesses the maximum number of hydrogen bonds with water molecules and no intramolecular hydrogen bonds, is chosen as the standard one. When a donor and an acceptor (e.g., N and O, respectively) are buried in the interior after the break of hydrogen bonds with water molecules, if they form an intramolecular hydrogen bond, no penalty is imposed. When a donor or an acceptor is buried with no intramolecular hydrogen bond formed, an energetic penalty is imposed. We examine all the donors and acceptors for backbone‐backbone, backbone‐side chain, and side chain‐side chain intramolecular hydrogen bonds and calculate the TDP. Our free‐energy function has been tested for three different decoy sets. It is better than any other physics‐based or knowledge‐based potential function in terms of the accuracy in discriminating the native fold from misfolded decoys and the achievement of high Z‐scores. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Non‐native α‐helix formation is not necessary for folding of lipocalin: Comparison of burst‐phase folding between tear lipocalin and β‐lactoglobulin

Tear lipocalin and β‐lactoglobulin are members of the lipocalin superfamily. They have similar tertiary structures but unusually low overall sequence similarity. Non‐native helical structures are formed during the early stage of β‐lactoglobulin folding. To address whether the non‐native helix formation is found in the folding of other lipocalin superfamily proteins, the folding kinetics of a tear lipocalin variant were investigated by stopped‐flow methods measuring the time‐dependent changes in circular dichroism (CD) spectrum and small‐angle X‐ray scattering (SAXS). CD spectrum showed that extensive secondary structures are not formed during a burst‐phase (within a measurement dead time). The SAXS data showed that the radius of gyration becomes much smaller than in the unfolded state during the burst‐phase, indicating that the molecule is collapsed during an early stage of folding. Therefore, non‐native helix formation is not general for folding of all lipocalin family members. The non‐native helix content in the burst‐phase folding appears to depend on helical propensities of the amino acid sequence. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Community‐level plant–soil feedbacks explain landscape distribution of native and non‐native plants

Plant–soil feedbacks (PSFs) have gained attention for their potential role in explaining plant growth and invasion. While promising, most PSF research has measured plant monoculture growth on different soils in short‐term, greenhouse experiments. Here, five soil types were conditioned by growing one native species, three non‐native species, or a mixed plant community in different plots in a common‐garden experiment. After 4 years, plants were removed and one native and one non‐native plant community were planted into replicate plots of each soil type. After three additional years, the percentage cover of each of the three target species in each community was measured. These data were used to parameterize a plant community growth model. Model predictions were compared to native and non‐native abundance on the landscape. Native community cover was lowest on soil conditioned by the dominant non‐native, Centaurea diffusa, and non‐native community cover was lowest on soil cultivated by the dominant native, Pseudoroegneria spicata. Consistent with plant growth on the landscape, the plant growth model predicted that the positive PSFs observed in the common‐garden experiment would result in two distinct communities on the landscape: a native plant community on native soils and a non‐native plant community on non‐native soils. In contrast, when PSF effects were removed, the model predicted that non‐native plants would dominate all soils, which was not consistent with plant growth on the landscape. Results provide an example where PSF effects were large enough to change the rank‐order abundance of native and non‐native plant communities and to explain plant distributions on the landscape. The positive PSFs that contributed to this effect reflected the ability of the two dominant plant species to suppress each other's growth. Results suggest that plant dominance, at least in this system, reflects the ability of a species to suppress the growth of dominant competitors through soil‐mediated effects.     

In vitro folding of methionine‐arginine human lyspro‐proinsulin S‐sulfonate—Disulfide formation pathways and factors controlling yield

We investigated the in vitro folding of an oxidized proinsulin (methionine‐arginine human lyspro‐proinsulin S‐sulfonate), using cysteine as a reducing agent at 5°C and high pH (10.5–11). Folding intermediates were detected and characterized by means of matrix‐assisted laser desorption ionization mass spectrometry (MALDI‐MS), reversed‐phase chromatography (RPC), size‐exclusion chromatography, and gel electrophoresis. The folding kinetics and yield depended on the protein and cysteine concentrations. RPC coupled with MALDI‐MS analyses indicated a sequential formation of intermediates with one, two, and three disulfide bonds. The MALDI‐MS analysis of Glu‐C digested, purified intermediates indicated that an intra‐A‐chain disulfide bond formed first among A6, A7, and A11. Various non‐native intra‐A (A20 with A6, A7, or A11), intra‐B (between B7 and B19), and inter‐A‐B disulfide bonds were observed in the intermediates with two disulfide bonds. The intermediates with three disulfide bonds had mainly the non‐native intra‐A and intra‐B bonds. At a cysteine‐to‐proinsulin‐SH ratio of 3.5, all intermediates with the non‐native disulfide bonds were converted to properly folded proinsulin via disulfide bond reshuffling, which was the slowest step. Aggregation via the formation of intermolecular disulfide bonds of early intermediates was the major cause of yield loss. At a higher cysteine‐to‐proinsulin‐SH ratio, some intermediates and folded MR‐KPB‐hPI were reduced to proteins with thiolate anions, which caused unfolding and even more yield loss than what resulted from aggregation of the early intermediates. Reducing protein concentration, while keeping an optimal cysteine‐to‐protein ratio, can improve folding yield significantly. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Differential effects of phe19 and phe20 on fibril formation by amyloidogenic peptide Aβ16–22 (Ac‐KLVFFAE‐NH2)

The sequence KLVFFAE (Aβ16–22) in Alzheimer's β‐amyloid is thought to be a core β‐structure that could act as a template for folding other parts of the polypeptide or molecules into fibrillar assemblies rich in β‐sheet. To elucidate the mechanism of the initial folding process, we undertook combined X‐ray fiber/powder diffraction and infrared (IR) spectroscopy to analyze lyophilized Aβ16–22 and solubilized/dried peptide containing nitrile probes at F19 and/or F20. Solubilized/dried wild‐type (WT) Aβ16–22 and the peptide containing cyanophenylalanine at F19 (19CN) or at F20 (20CN) gave fiber patterns consistent with slab‐like β‐crystallites that were cylindrically averaged around the axis parallel to the polypeptide chain direction. The WT and 19CN assemblies showed 30‐Å period arrays arising from the stacking of the slabs along the peptide chain direction, whereas the 20CN assemblies lacked any such stacking. The electron density projection along the peptide chain direction indicated similar side‐chain dispositions for WT and 20CN, but not for 19CN. These X‐ray results and modeling imply that in the assembly of WT Aβ16–22 the F19 side chain is localized within the intersheet space and is involved in hydrophobic contact with amino acids across the intersheet space, whereas the F20 side chain localized near the slab surface is less important for the intersheet interaction, but involved in slab stacking. IR observations for the same peptides in dilute solution showed a greater degree of hydrogen bonding for the nitrile groups in 20CN than in 19CN, supporting this interpretation. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Direct analysis of backbone-backbone hydrogen bond formation in protein folding transition states

Here we investigate the role of backbone-backbone hydrogen bonding interactions in stabilizing the protein folding transition states of two model protein systems, the B1 domain of protein L (ProtL) and the P22 Arc repressor. A backbone modified analogue of ProtL containing an amide-to-ester bond substitution between residues 105 and 106 was prepared by total chemical synthesis, and the thermodynamic and kinetic parameters associated with its folding reaction were evaluated. Ultimately, these parameters were used in a Phi-value analysis to determine if the native backbone-backbone hydrogen bonding interaction perturbed in this analogue (i.e. a hydrogen bond in the first beta-turn of ProtL's beta-beta-alpha-beta-beta fold) was formed in the transition state of ProtL's folding reaction. Also determined were the kinetic parameters associated with the folding reactions of two Arc repressor analogues, each containing an amide-to-ester bond substitution in the backbone of their polypeptide chains. These parameters were used together with previously established thermodynamic parameters for the folding of these analogues in Phi-value analyses to determine if the native backbone-backbone hydrogen bonding interactions perturbed in these analogues (i.e. a hydrogen bond at the end of the intersubunit beta-sheet interface and hydrogen bonds at the beginning of the second alpha-helix in Arc repressor's beta-alpha-alpha structure) were formed in the transition state of Arc repressor's folding reaction. Our results reveal that backbone-backbone hydrogen bonding interactions are formed in the beta-turn and alpha-helical transition state structures of ProtL and Arc repressor, respectively; and they were not formed in the intersubunit beta-sheet interface of Arc repressor, a region of Arc repressor's polypeptide chain previously shown to have other non-native-like conformations in Arc's protein folding transition state.