Secondary Structure Conversion Research Articles

Sulfated Polysaccharides as a Fighter with Protein Non-Physiological Aggregation: The Role of Polysaccharide Flexibility and Charge Density.

Proteins can lose native functionality due to non-physiological aggregation. In this work, we have shown the power of sulfated polysaccharides as a natural assistant to restore damaged protein structures. Protein aggregates enriched by cross-β structures are a characteristic of amyloid fibrils related to different health disorders. Our recent studies demonstrated that model fibrils of hen egg white lysozyme (HEWL) can be disaggregated and renatured by some negatively charged polysaccharides. In the current work, using the same model protein system and FTIR spectroscopy, we studied the role of conformation and charge distribution along the polysaccharide chain in the protein secondary structure conversion. The effects of three carrageenans (κ, ι, and λ) possessing from one to three sulfate groups per disaccharide unit were shown to be different. κ-Carrageenan was able to fully eliminate cross-β structures and complete the renaturation process. ι-Carrageenan only initiated the formation of native-like β-structures in HEWL, retaining most of the cross-β structures. In contrast, λ-carrageenan even increased the content of amyloid cross-β structures. Furthermore, κ-carrageenan in rigid helical conformation loses its capability to restore protein native structures, largely increasing the amount of amyloid cross-β structures. Our findings create a platform for the design of novel natural chaperons to counteract protein unfolding.

An in‐situ method to track the quality change of frozen surimi as a whole: Multi‐molecular infrared spectroscopy in combination with LF‐NMR

The rapid overall characterization of quality changes in complex food systems is considered extremely challenging. Herein, a multi-molecular infrared spectroscopy (MM-IR) methodology in combination with low-field nuclear magnetic resonance and chemometrics, to in-situ reveal the quality changes of surimi (protein) products as a whole was developed. Surimi with cryoprotectant (SWC) and without cryoprotectant (SWO) were used as a typical demonstration. Specifically, secondary structures in proteins could transform into each other and the structural variation of SWC was smaller than that of SWO. Cryoprotectant could slow down the water conversion from non-moveable to free, thereby slow down the rate of protein degradation. principal component analysis results of massive samples verified the differences in quality between SWC and SWO at different storage stages. Therefore, the addition of cryoprotectant led to the frozen surimi having less protein secondary structure conversion and water migration. In sum, the developed MM-IR could be applicable for the quality tracking of frozen surimi. Practical applications Low-temperature storage is currently the most mainstream storage method for surimi, and the preservation time of surimi can be extended by adding cryoprotectants. In this study, the overall quality of frozen surimi could be characterized in situ and simultaneous multi-component analysis could be realized, and the changes in the secondary structure of surimi protein and the migration of water could be tracked in situ. The results of this work had the potential to track the quality changes of complex frozen aquatic products as a whole and predicted the frozen storage time.

Toward the Mode of Action of the Clinical Stage All-d-Enantiomeric Peptide RD2 on Aβ42 Aggregation.

The aggregation of amyloid-β (Aβ) into oligomers and fibrillary structures is critical for the pathogenesis of Alzheimer's disease (AD). Recently, research effort has been focused on developing novel agents that can preferentially suppress Aβ oligomer mediated toxicities, for example, by directly targeting these toxic assemblies. The compound RD2 has been developed and optimized for Aβ42 monomer binding and stabilization of the monomer in its native intrinsically disordered conformation. It has been demonstrated to improve and even reverse the cognitive and behavioral deficits in AD mouse models, while the detailed mechanism of action is not fully clarified. Here we focused on exploring the interaction between RD2 and Aβ42 monomers and its consequences for the fibrillation of Aβ42. RD2 binds to Aβ42 monomers with nanomolar affinities, according to microscale thermophoresis and surface plasmon resonance measurements. Complexes between RD2 and Aβ42 monomers are formed at 1:1 and other stoichiometries, as revealed by analytical ultracentrifugation. At substoichiometric levels, RD2 slows down the secondary structure conversion of Aβ42 and significantly delays the fibril formation. Our research provides experimental evidence in supporting that RD2 eliminates toxic Aβ assemblies by stabilizing Aβ monomers in their native intrinsically disordered conformation. The study further supports the promising application of RD2 in counteracting Aβ aggregation related pathologies.

The initial stage of structural transformation of Aβ42 peptides from the human and mole rat in the presence of Fe2+ and Fe3+: Related to Alzheimer's disease

The early stage of secondary structural conversion of amyloid beta (Aβ) to misfolded aggregations is a key feature of Alzheimer's disease (AD). Under normal physiological conditions, Aβ peptides can protect neurons from the toxicity of highly concentrated metals. However, they become toxic under certain conditions. Under conditions of excess iron, amyloid precursor proteins (APP) become overexpressed. This subsequently increases Aβ production. Experimental studies suggest that Aβ fibrillation (main-pathway) and amorphous (off-pathway) aggregate formations are two competitive pathways driven by factors such as metal binding, pH and temperature. In this study, we performed molecular dynamic (MD) simulations to examine the initial stage of conformational transformations of human Aβ (hAβ) and rat Aβ (rAβ) peptides in the presence of Fe2+ and Fe3+ ions. Our results demonstrated that Fe2+ and Fe3+ play key roles in Aβs folding and aggregation. Fe3+ had a greater effect than Fe2+on Aβs’ folding during intermolecular interactions and subsequently, had a greater effect in decreasing structural diversity. Fe2+ was observed to be more likely than Fe3+ to interact with nitrogen atoms from the residues of imidazole rings of His. rAβ peptides are more energetically favorable than hAβ for intermolecular interactions and amorphous aggregations. We concluded that most hAβ structures were energetically unfavorable. However, hAβs with intermolecular β-sheet formations in the C-terminal were energetically favorable. It is notable that Fe2+ can change the surface charge of hAβ. Furthermore, Fe3+ can promote C-terminal folding by binding to Glu22 and Ala42, and by forming stable β-sheet formations on the C-terminal. Fe3+ can also pause the main-pathway by inducing random aggregations.

Menadione sodium bisulfite inhibits the toxic aggregation of amyloid-β(1–42)

Protein misfolding and aggregation are associated with amyloidosis. The toxic aggregation of amyloid-β 1–42 (Aβ42) may disrupt cell membranes and lead to cell death and is thus regarded as a contributing factor in Alzheimer's disease (AD). 1,4-naphthoquinone (NQ) has been shown to exhibit strong anti-aggregation effects on amyloidogenic proteins such as insulin and α-synuclein; however, its high toxicity and poor solubility limit its clinical application. Menadione sodium bisulfite (MSB, also known as vitamin K3), is used clinically in China to treat hemorrhagic diseases caused by vitamin K deficiency and globally as a vitamin K supplement. We hypothesized that MSB could inhibit amyloid formation since its backbone structure is similar to NQ. To test our hypothesis, we first investigated the effects of MSB on Aβ42 amyloid formation in vitro. We found that MSB inhibited Aβ42 amyloid formation in a dose dependent manner, delayed the secondary structural conversion of Aβ42 from random coil to ordered β-sheet, and attenuated the ability of Aβ42 aggregates to disrupt membranes; moreover, the quinone backbone rather than lipophilicity is esstial for the inhibitory effects of MSB. Next, in cells expressing a pathogenic APP mutation (Osaka mutation) that results in the formation of intraneuronal Aβ oligomers, MSB inhibited the intracellular aggregation of Aβ. Moreover, MSB treatment significantly extended the life span of Caenorhabditis elegans CL2120, a strain that expresses human Aβ42. Together, these results suggest that MSB and its derivatives may be further explored as potential therapeutic agents for the prevention or treatment of AD.

Temperature Dependent Pair-Wise Amino Acid Structural Propensity in ccβ and Implications for α to β Secondary Structural Conversion

Numerous proteins undergo an alpha-helix to beta-sheet conversion, which can then lead to the formation of toxic amyloid fibrils. The resulting fibrous protein aggregation is associated with the pathology of many debilitating illnesses and neurodegenerative disorders, such as type II diabetes, Alzheimer's disease, and Parkinson's disease. Many factors determine the preference that an amino acid has to assume an α-helix versus β-strand secondary structure. The large size of these proteins presents difficulties in investigating the molecular details of the dynamics of the structural transformation. We investigate the small, engineered protein ccβ that undergoes similar structural transformations. We use Replica Exchange Molecular Dynamics computer simulations to study the amino acids that make up the system. Based upon the importance of amino acid pairs, we examine the propensities for the various pairs of amino acids in ccβ to form specific secondary structures as a function of temperature. We find that there is an interplay between the inherent secondary structure preferences of amino acids with various molecular forces necessary to stabilize the helical trimer. Three specific pairs of amino acids in the ccβ primary sequence are especially important in the temperature-dependent structural conversion, and these pairs are located at positions along the sequence that facilitate the formation of the coiled-coil trimer. We show these effects vary with temperature, leading to the temperature dependent conversion of ccβ from helical to beta conformation.

Role of Hairpin-Quadruplex DNA Secondary Structural Conversion in the Promoter of hnRNP K in Gene Transcriptional Regulation.

The promoter of hnRNP K oncogene was found to contain a G/C-rich sequence on the same DNA strand, which can form interconvertible G-quadruplex, i-motif, and hairpin structures. Protein CNBP could bind and stabilize the G-quadruplex, inducing transformation of the hairpin into the G-quadruplex, resulting in down-regulation of hnRNP K transcription. In contrast, Corticosterone could bind and stabilize the hairpin, inducing transformation of the G-quadruplex into the hairpin, resulting in up-regulation of hnRNP K gene transcription.

Phenolic compounds prevent the oligomerization of α-synuclein and reduce synaptic toxicity.

Lewy bodies, mainly composed of α-synuclein (αS), are pathological hallmarks of Parkinson's disease and dementia with Lewy bodies. Epidemiological studies showed that green tea consumption or habitual intake of phenolic compounds reduced Parkinson's disease risk. We previously reported that phenolic compounds inhibited αS fibrillation and destabilized preformed αS fibrils. Cumulative evidence suggests that low-order αS oligomers are neurotoxic and critical species in the pathogenesis of α-synucleinopathies. To develop disease modifying therapies for α-synucleinopathies, we examined effects of phenolic compounds (myricetin (Myr), curcumin, rosmarinic acid (RA), nordihydroguaiaretic acid, and ferulic acid) on αS oligomerization. Using methods such as photo-induced cross-linking of unmodified proteins, circular dichroism spectroscopy, the electron microscope, and the atomic force microscope, we showed that Myr and RA inhibited αS oligomerization and secondary structure conversion. The nuclear magnetic resonance analysis revealed that Myr directly bound to the N-terminal region of αS, whereas direct binding of RA to monomeric αS was not detected. Electrophysiological assays for long-term potentiation in mouse hippocampal slices revealed that Myr and RA ameliorated αS synaptic toxicity by inhibition of αS oligomerization. These results suggest that Myr and RA prevent the αS aggregation process, reducing the neurotoxicity of αS oligomers. To develop disease modifying therapies for α-synucleinopathies, we examined effects of phenolic compounds on α-synuclein (αS) oligomerization. Phenolic compounds, especially Myricetin (Myr) and Rosmarinic acid (RA), inhibited αS oligomerization and secondary structure conversion. Myr and RA ameliorated αS synaptic toxicity on the experiment of long-term potentiation. Our results suggest that Myr and RA prevent αS aggregation process and reduce the neurotoxicity of αS oligomers. Phenolic compounds are good candidates of disease modifying drugs for α-synucleinopathies.

Reversible and rapid pH-regulated self-assembly of a poly(ethylene glycol)-peptide bioconjugate.

The use of external triggers to manipulate the secondary structure of self-assembling peptides conjugated to flexible synthetic polymers is a challenging problem, particularly in terms of reversibility. Here, we demonstrate, for the first time, sustained rapid and reversible, pH-regulated self-assembly of the peptide ELELELELELF (EL-5F) and its conjugates with 2 and 5 kDa poly(ethylene glycol) (EL-5F-PEG-2K and EL-5F-PEG-5K). Circular dichroism indicated the formation of β-sheet structures at pH < 5.9, 5.8, and 5.4 and disassembly to random coils above those pH values for EL-5F, EL-5F-PEG-2K, and EL-5F-PEG-5K, respectively. β-sheets were confirmed by the thioflavin T assay, while transmission electron microscopy revealed the existence of extended fibrillar structures below the above pH values. pH-induced secondary structure conversion was reproducible for over 15 cycles, even at salt concentrations of up to 200 mM NaCl, and was quantitatively related to the pH. Self-supporting hydrogelation after self-assembly was observed at concentrations as low as 0.2 wt %, which is 15-fold lower than previously reported concentrations. This simple approach to mediate reversible self-assembly of EL-5F-PEG bioconjugates is expected to offer novel functionality relevant to drug delivery and bioseparation systems.

Hot spot of structural ambivalence in prion protein revealed by secondary structure principal component analysis.

The conformational conversion of proteins into an aggregation-prone form is a common feature of various neurodegenerative disorders including Alzheimer's, Huntington's, Parkinson's, and prion diseases. In the early stage of prion diseases, secondary structure conversion in prion protein (PrP) causing β-sheet expansion facilitates the formation of a pathogenic isoform with a high content of β-sheets and strong aggregation tendency to form amyloid fibrils. Herein, we propose a straightforward method to extract essential information regarding the secondary structure conversion of proteins from molecular simulations, named secondary structure principal component analysis (SSPCA). The definite existence of a PrP isoform with an increased β-sheet structure was confirmed in a free-energy landscape constructed by mapping protein structural data into a reduced space according to the principal components determined by the SSPCA. We suggest a "spot" of structural ambivalence in PrP-the C-terminal part of helix 2-that lacks a strong intrinsic secondary structure, thus promoting a partial α-helix-to-β-sheet conversion. This result is important to understand how the pathogenic conformational conversion of PrP is initiated in prion diseases. The SSPCA has great potential to solve various challenges in studying highly flexible molecular systems, such as intrinsically disordered proteins, structurally ambivalent peptides, and chameleon sequences.

Structures and Free Energy Landscapes of the A53T Mutant-Type α-Synuclein Protein and Impact of A53T Mutation on the Structures of the Wild-Type α-Synuclein Protein with Dynamics

The A53T genetic missense mutation of the wild-type α-synuclein (αS) protein was initially identified in Greek and Italian families with familial Parkinson's disease. Detailed understanding of the structures and the changes induced in the wild-type αS structure by the A53T mutation, as well as establishing the direct relationships between the rapid conformational changes and free energy landscapes of these intrinsically disordered fibrillogenic proteins, helps to enhance our fundamental knowledge and to gain insights into the pathogenic mechanism of Parkinson's disease. We employed extensive parallel tempering molecular dynamics simulations along with thermodynamic calculations to determine the secondary and tertiary structural properties as well as the conformational free energy surfaces of the wild-type and A53T mutant-type αS proteins in an aqueous solution medium using both implicit and explicit water models. The confined aqueous volume effect in the simulations of disordered proteins using an explicit model for water is addressed for a model disordered protein. We also assessed the stabilities of the residual secondary structure component interconversions in αS based on free energy calculations at the atomic level with dynamics using our recently developed theoretical strategy. To the best of our knowledge, this study presents the first detailed comparison of the structural properties linked directly to the conformational free energy landscapes of the monomeric wild-type and A53T mutant-type α-synuclein proteins in an aqueous solution environment. Results demonstrate that the β-sheet structure is significantly more altered than the helical structure upon A53T mutation of the monomeric wild-type αS protein in aqueous solution. The β-sheet content close to the mutation site in the N-terminal region is more abundant while the non-amyloid-β component (NAC) and C-terminal regions show a decrease in β-sheet abundance upon A53T mutation. Obtained results utilizing our new theoretical strategy show that the residual secondary structure conversion stabilities resulting in α-helix formation are not significantly affected by the mutation. Interestingly, the residual secondary structure conversion stabilities show that secondary structure conversions resulting in β-sheet formation are influenced by the A53T mutation and the most stable residual transition yielding β-sheet occurs directly from the coil structure. Long-range interactions detected between the NAC region and the N- or C-terminal regions of the wild-type αS disappear upon A53T mutation. The A53T mutant-type αS structures are thermodynamically more stable than those of the wild-type αS protein structures in aqueous solution. Overall, the higher propensity of the A53T mutant-type αS protein to aggregate in comparison to the wild-type αS protein is related to the increased β-sheet formation and lack of strong intramolecular long-range interactions in the N-terminal region in comparison to its wild-type form. The specific residual secondary structure component stabilities reported herein provide information helpful for designing and synthesizing small organic molecules that can block the β-sheet forming residues, which are reactive toward aggregation.

Kinetics of peptide secondary structure conversion during amyloid β-protein fibrillogenesis

Amyloid fibrils are a common component in many debilitating human neurological diseases such as Alzheimer's (AD), Parkinson's, and Creutzfeldt-Jakob, and in animal diseases such as BSE. The role of fibrillar Αβ proteins in AD has stimulated interest in the kinetics of Αβ fibril formation. Kinetic models that include reaction pathways and rate parameters for the various stages of the process can be helpful towards understanding the dynamics on a molecular level. Based upon experimental data, we have developed a mathematical model for the reaction pathways and determined rate parameters for peptide secondary structural conversion and aggregation during the entire fibrillogenesis process from random coil to mature fibrils, including the molecular species that accelerate the conversions. The model and the rate parameters include different molecular structural stages in the nucleation and polymerization processes and the numerical solutions yield graphs of concentrations of different molecular species versus time that are in close agreement with experimental results. The model also allows for the calculation of the time-dependent increase in aggregate size. The calculated results agree well with experimental results, and allow differences in experimental conditions to be included in the calculations. The specific steps of the model and the rate constants that are determined by fitting to experimental data provide insight on the molecular species involved in the fibril formation process.

Design of anti- and pro-aggregation variants to assess the effects of methionine oxidation in human prion protein.

Prion disease is characterized by the alpha-->beta structural conversion of the cellular prion protein (PrP(C)) into the misfolded and aggregated "scrapie" (PrP(Sc)) isoform. It has been speculated that methionine (Met) oxidation in PrP(C) may have a special role in this process, but has not been detailed and assigned individually to the 9 Met residues of full-length, recombinant human PrP(C) [rhPrP(C)(23-231)]. To better understand this oxidative event in PrP aggregation, the extent of periodate-induced Met oxidation was monitored by electrospray ionization-MS and correlated with aggregation propensity. Also, the Met residues were replaced with isosteric and chemically stable, nonoxidizable analogs, i.e., with the more hydrophobic norleucine (Nle) and the highly hydrophilic methoxinine (Mox). The Nle-rhPrP(C) variant is an alpha-helix rich protein (like Met-rhPrP(C)) resistant to oxidation that lacks the in vitro aggregation properties of the parent protein. Conversely, the Mox-rhPrP(C) variant is a beta-sheet rich protein that features strong proaggregation behavior. In contrast to the parent Met-rhPrP(C), the Nle/Mox-containing variants are not sensitive to periodate-induced in vitro aggregation. The experimental results fully support a direct correlation of the alpha-->beta secondary structure conversion in rhPrP(C) with the conformational preferences of Met/Nle/Mox residues. Accordingly, sporadic prion and other neurodegenerative diseases, as well as various aging processes, might also be caused by oxidative stress leading to Met oxidation.

Secondary structure conversions of Alzheimer’s Aβ(1–40) peptide induced by membrane‐mimicking detergents

The amyloid beta peptide (Abeta) with 39-42 residues is the major component of amyloid plaques found in brains of Alzheimer's disease patients, and soluble oligomeric peptide aggregates mediate toxic effects on neurons. The Abeta aggregation involves a conformational change of the peptide structure to beta-sheet. In the present study, we report on the effect of detergents on the structure transitions of Abeta, to mimic the effects that biomembranes may have. In vitro, monomeric Abeta(1-40) in a dilute aqueous solution is weakly structured. By gradually adding small amounts of sodium dodecyl sulfate (SDS) or lithium dodecyl sulfate to a dilute aqueous solution, Abeta(1-40) is converted to beta-sheet, as observed by CD at 3 degrees C and 20 degrees C. The transition is mainly a two-state process, as revealed by approximately isodichroic points in the titrations. Abeta(1-40) loses almost all NMR signals at dodecyl sulfate concentrations giving rise to the optimal beta-sheet content (approximate detergent/peptide ratio = 20). Under these conditions, thioflavin T fluorescence measurements indicate a maximum of aggregated amyloid-like structures. The loss of NMR signals suggests that these are also involved in intermediate chemical exchange. Transverse relaxation optimized spectroscopy NMR spectra indicate that the C-terminal residues are more dynamic than the others. By further addition of SDS or lithium dodecyl sulfate reaching concentrations close to the critical micellar concentration, CD, NMR and FTIR spectra show that the peptide rearranges to form a micelle-bound structure with alpha-helical segments, similar to the secondary structures formed when a high concentration of detergent is added directly to the peptide solution.

Secondary structure transitions and aggregation induced in dynorphin neuropeptides by the detergent sodium dodecyl sulfate

Dynorphins, endogeneous opioid neuropeptides, function as ligands to the opioid kappa receptors and also induce non-opioid effects in neurons, probably related to direct membrane interactions. We have characterized the structure transitions of dynorphins (big dynorphin, dynorphin A and dynorphin B) induced by the detergent sodium dodecyl sulfate (SDS). In SDS titrations monitored by circular dichroism, we observed secondary structure conversions of the peptides from random coil to α-helix with a highly aggregated intermediate. As determined by Fourier transform infrared spectroscopy, this intermediate exhibited β-sheet structure for dynorphin B and big dynorphin. In contrast, aggregated dynorphin A was α-helical without considerable β-sheet content. Hydrophobicity analysis indicates that the YGGFLRR motif present in all dynorphins is prone to be inserted in the membrane. Comparing big dynorphin with dynorphin A and dynorphin B, we suggest that the potent neurotoxicity of big dynorphin could be related to the combination of amino acid sequences and secondary structure propensities of dynorphin A and dynorphin B, which may generate a synergistic effect for big dynorphin membrane perturbing properties. The induced aggregated α-helix of dynorphin A is also correlated with membrane perturbations, whereas the β-sheet of dynorphin B does not correlate with membrane perturbations.

Secondary structure conversions of Mycobacterium tuberculosis ribonucleotide reductase protein R2 under varying pH and temperature conditions

The structural properties of Mycobacterium tuberculosis (Mtb) ribonucleotide reductase R2 protein were studied under varying pH and temperature conditions by circular dichroism (CD) spectroscopy as well as dynamic light scattering (DLS). Under physiological conditions this protein has a high α-helical content, similar to the corresponding protein from other species, e.g. mouse. Decreasing the pH induced significant structure conversions. When pH was below 6.5 an aggregated structure was observed and reached a maximum at pH 4. The aggregated state of this protein was verified by DLS and was found to be rich in β-structure. This amyloid-like structure transformed into a molten globule state with high temperature stability (between 25 and 80 °C) at pH below 3. The corresponding mouse protein R2 under similar conditions showed no evidence of an aggregated state around pH 4.

Learning probabilistic models of tree edit distance

Nowadays, there is a growing interest in machine learning and pattern recognition for tree-structured data. Trees actually provide a suitable structural representation to deal with complex tasks such as web information extraction, RNA secondary structure prediction, computer music, or conversion of semi-structured data (e.g. XML documents). Many applications in these domains require the calculation of similarities over pairs of trees. In this context, the tree edit distance (ED) has been subject of investigations for many years in order to improve its computational efficiency. However, used in its classical form, the tree ED needs a priori fixed edit costs which are often difficult to tune, that leaves little room for tackling complex problems. In this paper, to overcome this drawback, we focus on the automatic learning of a non-parametric stochastic tree ED. More precisely, we are interested in two kinds of probabilistic approaches. The first one builds a generative model of the tree ED from a joint distribution over the edit operations, while the second works from a conditional distribution providing then a discriminative model. To tackle these tasks, we present an adaptation of the expectation–maximization algorithm for learning these distributions over the primitive edit costs. Two experiments are conducted. The first is achieved on artificial data and confirms the interest to learn a tree ED rather than a priori imposing edit costs; The second is applied to a pattern recognition task aiming to classify handwritten digits.

Metal-Coupled Folding of Cys2His2 Zinc-Finger

Zinc-fingers, which widely exist in eukaryotic cell and play crucial roles in life processes, depend on the binding of zinc ion for their proper folding. To computationally study the zinc-coupled folding of the zinc-fingers, charge transfer and metal induced protonation/deprotonation effects have to be considered. Here, by attempting to implicitly account for such effects in classical molecular dynamics and performing intensive simulations with explicit solvent for the peptides with and without zinc binding, we investigate the folding of the Cys2His2-type zinc-finger motif and the coupling between the peptide folding and zinc binding. We find that zinc ion not only stabilizes the native structure but also participates in the whole folding process. It binds to the peptide at an early stage of folding and directs or modulates the folding and stabilizations of the component beta-hairpin and alpha-helix. Such a crucial role of zinc binding is mediated by the packing of the conserved hydrophobic residues. We also find that the packing of the hydrophobic residues and the coordination of the native ligands are coupled. Meanwhile, the processes of zinc binding, mis-ligation, ligand exchange, and zinc induced secondary structure conversion as well as the water behavior due to the involvement of zinc ion are characterized. Our results are in good agreement with related experimental observations and provide significant insight into the general mechanisms of the metal cofactor dependent protein folding and other metal-induced conformational changes of biological importance.

CSSP2: An improved method for predicting contact-dependent secondary structure propensity

The calculation of contact-dependent secondary structure propensity (CSSP) has been reported to sensitively detect non-native β-strand propensities in the core sequences of amyloidogenic proteins. Here we describe a noble energy-based CSSP method implemented on dual artificial neural networks that rapidly and accurately estimate the potential for the non-native secondary structure formation in local regions of protein sequences. In this method, we attempted to quantify long-range interaction patterns in diverse secondary structures by potential energy calculations and decomposition on a pairwise per-residue basis. The calculated energy parameters and seven-residue sequence information were used as inputs for artificial neural networks (ANNs) to predict sequence potential for secondary structure conversion. The trained single ANN using the >( i, i ± 4) interaction energy parameter exhibited 74% accuracy in predicting the secondary structure of test sequences in their native energy state, while the dual ANN-based predictor using ( i, i ± 4) and >( i, i ± 4) interaction energies showed 83% prediction accuracy. The present method provides a simple and accurate tool for predicting sequence potential for secondary structure conversions without using 3D structural information.

Analysis of Chameleon Sequences by Energy Decomposition on a Pairwise Per-residue Basis

The conversion from alpha-helix to beta-strand that has been widely observed in so-called chameleon sequences has received considerable attention since such a structural change may induce many amyloidogenic proteins to self-assemble into fibrils thus causing fatal diseases. Here we report a large scale-analysis of the energetics of secondary structural conversions in a collection of chameleon sequences retrieved from the Protein Data Bank. Major energetic contributions to the secondary structural conversion were analyzed by carrying out energy decomposition on a pairwise per-residue basis, i.e., (i,i), (i,i +/- 1), (i,i +/- 2), (i,i +/- 3), (i,i +/- 4) and > (i,i +/- 4) intra-/inter-residual interactions. While the overall potential energy differences were subtle, individual residue-based interacting energy differences were observed to vary significantly depending on the specific type of secondary structural conversion. The average energy difference between alpha-helix and beta-strand, <DeltaE (alpha-->beta)>, in the chameleon sequences varied significantly in (i,i), (i,i +/- 1) and > (i,i +/- 4) interactions. The major energetic factors in secondary structure conversions were electrostatic interactions and the polar term for solvation energy. In addition, residue-based average energy differences in alpha-helix --> beta-strand conversions were well-correlated to those in alpha-helix --> random coil --> beta-strand conversions (R2 = 0.92). Assuming that three secondary structural elements can transform in either direction, this strong correlation indicates that the present energy decomposition method using database structures of chameleon sequences provides a reliable tool for the characterization of secondary structure fluctuations in amino acid sequences.