Detailed folding structures of M-lycotoxin-Hc1a and its mutageneses using 2D HP model

Although the hydrophobic-polar (HP) model was proposed a decade ago, it applies almost to no real-case study because of its intense computation. In this study, a 2D HP model was applied to study the folding structures of M-lycotoxin-Hc1a, an antimicrobial peptide, in order to get full pictures of its numerous folding structures. The normalised hydrophobicity index was used to convert M-lycotoxin-Hc1a and its six mutageneses into HP sequences, and then the 2D HP model was used to compute all the possible folding structures (3²⁴ = 282,429,536,481), and finally the normalised hydrophobicity index was used to distinguish the native state. The results showed that M-lycotoxin-Hc1a had 6 and 138 folding structures at their native state with the minimal energy of ? 13 at pH 2 and pH 7 when glycine served as hydrophobic amino acid. When glycine serves as polar amino acid, M-lycotoxin-Hc1a had 12 and 54 folding structures at their native state with the minimal energy of ? 12 and ? 13 at pH 2 and pH 7, respectively. This study advanced the knowledge on how to apply the HP model to real-life study, and how the mutageneses influenced the folding structures of M-lycotoxin-Hc1a, their native states and minimal energy at different pH levels.     

Detailed folding structures of kappa-conotoxin RIIIJ and its mutageneses obtained from 2-dimensional HP model

Kappa-conotoxin RIIIJ is a conopeptide to inhibit voltage-gated potassium channels, however, its detailed folding structures have yet to be studied. With the advance in computing power, it is possible to use the HP model to analyze all its possible folding structures. In this study, the amino acid sequences of kappa-conotoxin RIIIJ and its four mutageneses were converted into ten HP sequences according to the normalized hydrophobicity index. All 282 429 536 481 possible folding structures in each HP sequence were found using the 2-dimensional HP model, and the detailed folding structures at native state were studied. The results showed that kappa-conotoxin RIIIJ had 180 and 90 folding structures at their native state with minimal energy of -9 and -10 at pH 2 and pH 7; its mutagenesis (6-8) TPP - > SLN increased the numbers of the folding structures to 456 and 564 at pH 2 and pH 7; whereas its mutageneses (6-11) TPPKKH - > SLNLRL, (9- 11) KKH - > LRL, and (10-11) KH - > RL decreased the numbers of the folding structures to 60, 30 and 90 at both pH levels, respectively. Thereafter, the normalized hydrophobicity index was employed to distinguish those native states, and attempts were made to explain the effect of mutageneses on potassium channels in terms of the number of folding structures and numerical native states.     

Tuning a polar molecule for selective cytoplasmic delivery by a pH (Low) insertion peptide

Drug molecules are typically hydrophobic and small in order to traverse membranes to reach cytoplasmic targets, but we have discovered that more polar molecules can be delivered across membranes using water-soluble, moderately hydrophobic membrane peptides of the pHLIP (pH low insertion peptide) family. Delivery of polar cargo molecules could expand the chemical landscape for pharmacological agents that have useful activity but are too polar by normal drug criteria. The spontaneous insertion and folding of the pHLIP peptide across a lipid bilayer seeks a free energy minimum, and insertion is accompanied by a release of energy that can be used to translocate cell-impermeable cargo molecules. In this study, we report our first attempt to tune the hydrophobicity of a polar cargo, phallacidin, in a systematic manner. We present the design, synthesis, and characterization of three phallacidin cargoes, where the hydrophobicity of the cargo was tuned by the attachment of diamines of various lengths of hydrophobic chains. The phallacidin cargoes were conjugated to pHLIP and shown to selectively inhibit the proliferation of cancer cells in a concentration-dependent manner at low pH.     

Modeling the effects of mutations on the denatured states of proteins.

We develop a model for the reversible denaturation of proteins and for the effects of single-site mutations on the denatured states. The model is based on short chains of sequences of H (hydrophobic) and P (other) monomers configured as self-avoiding walks on the two-dimensional square lattice. The N (native) state is defined as the unique conformation of lowest contact energy, whereas the D (denatured) state is defined as the collection of all other conformations. With this model we are able to determine the exact partition function, and thus the exact native-denatured equilibrium for various solvent conditions, using the computer to exhaustively enumerate every possible configuration. Previous studies confirm that this model shows many aspects of protein-like behavior. The present study attempts to model how the denatured state (1) depends on the amino acid sequence, and (2) is changed by single-site mutations. The model accounts for two puzzling experimental results: (1) the replacement of a polar residue by a hydrophobic amino acid on the surface of a protein can destabilize a native protein, and (2) the "denaturant slope," m = partial delta G/partial c (where c is the concentration of denaturant--urea, guanidine hydrochloride), can sometimes change by as much as 30% due to a single mutation. The principal conclusion of the present study is that, under strong folding conditions, the denatured conformations that are in equilibrium with the native state are not open random configurations. Instead, they are an ensemble of highly compact conformations with a distribution that depends on the residue sequence and that can be substantially altered by single mutations. Most importantly, we conclude that mutations can exert their dominant effects on protein stability by changing the entropy of folding.     

Studying the folding process of the acylphosphatase from Sulfolobus solfataricus. A comparative analysis with other proteins from the same superfamily

The folding process of the acylphosphatase from Sulfolobus solfataricus (Sso AcP) has been followed, starting from the fully unfolded state, using a variety of spectroscopic probes, including intrinsic fluorescence, circular dichroism, and ANS binding. The results indicate that an ensemble of partially folded or misfolded species form rapidly on the submillisecond time scale after initiation of folding. This conformational ensemble produces a pronounced downward curvature in the Chevron plot, appears to possess a content of secondary structure similar to that of the native state, as revealed by far-UV circular dichroism, and appears to have surface-exposed hydrophobic clusters, as indicated by the ability of this ensemble to bind to 8-anilino-1-naphthalenesulfonic acid (ANS). Sso AcP folds from this conformational state with a rate constant of ca. 5 s(-1) at pH 5.5 and 37 degrees C. A minor slow exponential phase detected during folding (rate constant of 0.2 s(-1) under these conditions) is accelerated by cyclophilin A and is absent in a mutant of Sso AcP in which alanine replaces the proline residue at position 50. This indicates that for a lower fraction of Sso AcP molecules the folding process is rate-limited by the cis-trans isomerism of the peptide bond preceding Pro50. A comparative analysis with four other homologous proteins from the acylphosphatase superfamily shows that sequence hydrophobicity is an important determinant of the conformational stability of partially folded states that may accumulate during folding of a protein. A low net charge and a high propensity to form alpha-helical structure also emerge as possibly important determinants of the stability of partially folded states. A significant correlation is also observed between folding rate and hydrophobic content of the sequence within this superfamily, lending support to the idea that sequence hydrophobicity, in addition to relative contact order and conformational stability of the native state, is a key determinant of folding rate.     

Positive and negative design in stability and thermal adaptation of natural proteins

The aim of this work is to elucidate how physical principles of protein design are reflected in natural sequences that evolved in response to the thermal conditions of the environment. Using an exactly solvable lattice model, we design sequences with selected thermal properties. Compositional analysis of designed model sequences and natural proteomes reveals a specific trend in amino acid compositions in response to the requirement of stability at elevated environmental temperature: the increase of fractions of hydrophobic and charged amino acid residues at the expense of polar ones. We show that this “from both ends of the hydrophobicity scale” trend is due to positive (to stabilize the native state) and negative (to destabilize misfolded states) components of protein design. Negative design strengthens specific repulsive non-native interactions that appear in misfolded structures. A pressure to preserve specific repulsive interactions in non-native conformations may result in correlated mutations between amino acids that are far apart in the native state but may be in contact in misfolded conformations. Such correlated mutations are indeed found in TIM barrel and other proteins.     

Geometrical model for the native-state folds of proteins

We recently introduced a physical model [T.X. Hoang, A. Trovato, F. Seno, J.R. Banavar, A. Maritan, Geometry and symmetry pre-sculpt the free energy landscape of proteins. Proc. Natl. Acad. Sci. U. S. A. 101 (2004) 7960-7964, J.R. Banavar, T.X. Hoang, A. Maritan, F. Seno, A. Trovato, A unified perspective on proteins-a physics approach. Phys. Rev., E 70 (2004) 041905] for proteins which incorporates, in an approximate manner, several key features such as the inherent anisotropy of a chain molecule, the geometrical and energetic constraints placed by the hydrogen bonds and sterics, and the role played by hydrophobicity. Within this framework, marginally compact conformations resembling the native state folds of proteins emerge as broad competing minima in the free energy landscape even for a homopolymer. Here we show how the introduction of sequence heterogeneity using a simple scheme of just two types of amino acids, hydrophobic (H) and polar (P), and sequence design allows a selected putative native fold to become the free energy minimum at low temperature. The folding transition exhibits thermodynamic cooperativity, if one neglects the degeneracy between two different low energy conformations sharing the same fold topology.     

The endoplasmic reticulum glucosyltransferase recognizes nearly native glycoprotein folding intermediates

The UDP-Glc:glycoprotein glucosyltransferase (GT), a key player in the endoplasmic reticulum (ER) quality control of glycoprotein folding, only glucosylates glycoproteins displaying non-native conformations. To determine whether GT recognizes folding intermediates or irreparably misfolded species with nearly native structures, we generated and tested as GT substrates neoglycoprotein fragments derived from chymotrypsin inhibitor 2 (GCI2) bearing from 53 to 64 (full-length) amino acids. Fragment conformations mimicked the last stage-folding structures adopted by a glycoprotein entering the ER lumen. GT catalytic efficiency (V(max)/K(m)) remained constant from GCI2-(1-53) to GCI2-(1-58) and then steadily declined to reach a minimal value with GCI2-(1-64). The same parameter showed a direct hyperbolic relationship with solvent accessibility of the single Trp residue but only in fragments exposing hydrophobic amino acid patches. Mutations introduced (GCI2-(1-63)V63S and GCI2-(1-64)V63S) produced slight structural destabilizations but increased GT catalytic efficiency. This parameter presented an inverse exponential relationship with the free energy of unfolding of canonical and mutant fragments. Moreover, the catalytic efficiency showed a linear relationship with the fraction of unfolded species in water. It was concluded that the GT-derived quality control may be operative with nearly native conformers and that no alternative ER-retaining mechanisms are required when glycoproteins approach their proper folding.     

Hydrophobic Core Formation and Dehydration in Protein Folding Studied by Generalized-Ensemble Simulations

Despite its small size, chicken villin headpiece subdomain HP36 folds into the native structure with a stable hydrophobic core within several microseconds. How such a small protein keeps up its conformational stability and fast folding in solution is an important issue for understanding molecular mechanisms of protein folding. In this study, we performed multicanonical replica-exchange simulations of HP36 in explicit water, starting from a fully extended conformation. We observed at least five events of HP36 folding into nativelike conformations. The smallest backbone root mean-square deviation from the crystal structure was 1.1 Å. In the nativelike conformations, the stably formed hydrophobic core was fully dehydrated. Statistical analyses of the simulation trajectories show the following sequential events in folding of HP36: 1), Helix 3 is formed at the earliest stage; 2), the backbone and the side chains near the loop between Helices 2 and 3 take nativelike conformations; and 3), the side-chain packing at the hydrophobic core and the dehydration of the core side chains take place simultaneously at the later stage of folding. This sequence suggests that the initial folding nucleus is not necessarily the same as the hydrophobic core, consistent with a recent experimental ϕ-value analysis.     

Molecular dynamics simulations of beta-hairpin folding

Molecular dynamics simulations of beta-hairpin folding have been carried out with a solvent-referenced potential at 274 K. The model peptide V4DPGV4 formed stable beta-hairpin conformations and the beta-hairpin ratio calculated by the DSSP algorithm was about 56% in the 50-ns simulation. Folding into beta-hairpin conformations is independent of the initial conformations. The simulations provided insights into the folding mechanism. The hydrogen bond often formed in a beta-turn first, and then propagated by forming more hydrogen bonds along the strands. Unfolding and refolding occurred repeatedly during the simulations. Both the hydrogen bonding and the hydrophobic interaction played important roles in forming the ordered structure. Without the hydrophobic effect, stable beta-hairpin conformations did not form in the simulations. With the same energy functions, the alanine-based peptide (AAQAA)3Y folded into helical conformations, in agreement with experiments. Folding into an alpha-helix or a beta-hairpin is amino acid sequence-dependent.