Search

Biomolecular function is realized by recognition, and increasing evidence shows that recognition is determined not only by structure but also by flexibility and dynamics. We explored a biomolecular recognition process that involves a major conformational change – protein folding. In particular, we explore the binding-induced folding of IA3, an intrinsically disordered protein that blocks the active site cleft of the yeast aspartic proteinase saccharopepsin (YPrA) by folding its own N-terminal residues into an amphipathic alpha helix. We developed a multi-scaled approach that explores the underlying mechanism by combining structure-based molecular dynamics simulations at the residue level with a stochastic path method at the atomic level. Both the free energy profile and the associated kinetic paths reveal a common scheme whereby IA3 binds to its target enzyme prior to folding itself into a helix. This theoretical result is consistent with recent time-resolved experiments. Furthermore, exploration of the detailed trajectories reveals the important roles of non-native interactions in the initial binding that occurs prior to IA3 folding. In contrast to the common view that non-native interactions contribute only to the roughness of landscapes and impede binding, the non-native interactions here facilitate binding by reducing significantly the entropic search space in the landscape. The information gained from multi-scaled simulations of the folding of this intrinsically disordered protein in the presence of its binding target may prove useful in the design of novel inhibitors of aspartic proteinases., Author Summary The intrinsically disordered peptide IA3 is the endogenous inhibitor for the enzyme named yeast aspartic proteinase saccharopepsin (YPrA). In the presence of YPrA, IA3 folds itself into an amphipathic helix that blocks the active site cleft of the enzyme. We developed a multi-scaled approach to explore the underlying mechanism of this binding-induced ordering transition. Our approach combines a structure-based molecular dynamics model at the residue level with a stochastic path method at the atomic level. Our simulations suggest that IA3 inhibits YPrA through an induced-fit mechanism where the enzyme (YPrA) induces conformational change of its inhibitor (IA3). This expands the definition of an induced-fit model from its original meaning that the binding of substrate (IA3) drives conformational change in the protein (YPrA). Our result is consistent with recent kinetic experiments and provides a microscopic explanation for the underlying mechanism. We also discuss the important roles of non-native interactions and backtracking. These results enrich our understanding of the enzyme-inhibition mechanism and may have value in the design of drugs.

We applied a novel application of FLIM-FRET to in situ measurement and quantification of protein interactions to explore isoform specific differences in Aβ-ApoE interaction and ApoE tertiary conformation in senile plaques in human Alzheimer brain. ApoE3 interacts more closely with Aβ than ApoE4, but a greater proportion of Aβ molecules within plaques are decorated with ApoE4 than ApoE3, lending strong support to the hypothesis that isoform specific differences in ApoE are linked with Aβ deposition. We found an increased number of ApoE N-terminal fragments in ApoE4 plaques, consistent with the observation that ApoE4 is more easily cleaved than ApoE3. In addition, we measured a small but significant isoform specific difference in ApoE domain interaction. Based on our in situ data, supported by traditional biochemical data, we propose a pathway by which isoform specific conformational differences increase the level of cleavage at the hinge region of ApoE4, leading to a loss of ApoE function to mediate clearance of Aβ and thereby increase the risk of AD for carriers of the APOEε4 allele.

BACKGROUND: The question of how the aggregation of the neuronal protein α-synuclein contributes to neuronal toxicity in Parkinson's disease has been the subject of intensive research over the past decade. Recently, attention has shifted from the amyloid fibrils to soluble oligomeric intermediates in the α-synuclein aggregation process. These oligomers are hypothesized to be cytotoxic and to permeabilize cellular membranes, possibly by forming pore-like complexes in the bilayer. Although the subject of α-synuclein oligomer-membrane interactions has attracted much attention, there is only limited evidence that supports the pore formation by α-synuclein oligomers. In addition the existing data are contradictory.METHODOLOGY/PRINCIPAL FINDINGS: Here we have studied the mechanism of lipid bilayer disruption by a well-characterized α-synuclein oligomer species in detail using a number of in vitro bilayer systems and assays. Dye efflux from vesicles induced by oligomeric α-synuclein was found to be a fast all-or-none process. Individual vesicles swiftly lose their contents but overall vesicle morphology remains unaltered. A newly developed assay based on a dextran-coupled dye showed that non-equilibrium processes dominate the disruption of the vesicles. The membrane is highly permeable to solute influx directly after oligomer addition, after which membrane integrity is partly restored. The permeabilization of the membrane is possibly related to the intrinsic instability of the bilayer. Vesicles composed of negatively charged lipids, which are generally used for measuring α-synuclein-lipid interactions, were unstable to protein adsorption in general.CONCLUSIONS/SIGNIFICANCE: The dye efflux from negatively charged vesicles upon addition of α-synuclein has been hypothesized to occur through the formation of oligomeric membrane pores. However, our results show that the dye efflux characteristics are consistent with bilayer defects caused by membrane instability. These data shed new insights into potential mechanisms of toxicity of oligomeric α-synuclein species.

Prion diseases like Creutzfeldt-Jakob disease in humans, Scrapie in sheep or bovine spongiform encephalopathy are fatal neurodegenerative diseases, which can be of sporadic, genetic, or infectious origin. Prion diseases are transmissible between different species, however, with a variable species barrier. The key event of prion amplification is the conversion of the cellular isoform of the prion protein (PrP(C)) into the pathogenic isoform (PrP(Sc)). We developed a sodiumdodecylsulfate-based PrP conversion system that induces amyloid fibril formation from soluble α-helical structured recombinant PrP (recPrP). This approach was extended applying pre-purified PrP(Sc) as seeds which accelerate fibrillization of recPrP. In the present study we investigated the interspecies coherence of prion disease. Therefore we used PrP(Sc) from different species like Syrian hamster, cattle, mouse and sheep and seeded fibrillization of recPrP from the same or other species to mimic in vitro the natural species barrier. We could show that the in vitro system of seeded fibrillization is in accordance with what is known from the naturally occurring species barriers.

Background Pathogenic bacteria adhere to the host cell surface using a family of outer membrane proteins called Trimeric Autotransporter Adhesins (TAAs). Although TAAs are highly divergent in sequence and domain structure, they are all conceptually comprised of a C-terminal membrane anchoring domain and an N-terminal passenger domain. Passenger domains consist of a secretion sequence, a head region that facilitates binding to the host cell surface, and a stalk region. Methodology/Principal Findings Pathogenic species of Burkholderia contain an overabundance of TAAs, some of which have been shown to elicit an immune response in the host. To understand the structural basis for host cell adhesion, we solved a 1.35 Å resolution crystal structure of a BpaA TAA head domain from Burkholderia pseudomallei, the pathogen that causes melioidosis. The structure reveals a novel fold of an intricately intertwined trimer. The BpaA head is composed of structural elements that have been observed in other TAA head structures as well as several elements of previously unknown structure predicted from low sequence homology between TAAs. These elements are typically up to 40 amino acids long and are not domains, but rather modular structural elements that may be duplicated or omitted through evolution, creating molecular diversity among TAAs. Conclusions/Significance The modular nature of BpaA, as demonstrated by its head domain crystal structure, and of TAAs in general provides insights into evolution of pathogen-host adhesion and may provide an avenue for diagnostics.

The type III secretion system of the Salmonella flagellum consists of 6 integral membrane proteins: FlhA, FlhB, FliO, FliP, FliQ, and FliR. However, in some other type III secretion systems, a homologue of FliO is apparently absent, suggesting it has a specialized role. Deleting the fliO gene from the chromosome of a motile strain of Salmonella resulted in a drastic decrease of motility. Incubation of the ΔfliO mutant strain in motility agar, gave rise to pseudorevertants containing extragenic bypass mutations in FliP at positions R143H or F190L. Using membrane topology prediction programs, and alkaline phosphatase or GFPuv chimeric protein fusions into the FliO protein, we demonstrated that FliO is bitopic with its N-terminus in the periplasm and C-terminus in the cytoplasm. Truncation analysis of FliO demonstrated that overexpression of FliO43–125 or FliO1–95 was able to rescue motility of the ΔfliO mutant. Further, residue leucine 91 in the cytoplasmic domain was identified to be important for function. Based on secondary structure prediction, the cytoplasmic domain, FliO43–125, should contain beta-structure and alpha-helices. FliO43–125-Ala was purified and studied using circular dichroism spectroscopy; however, this domain was disordered, and its structure was a mixture of beta-sheet and random coil. Coexpression of full-length FliO with FliP increased expression levels of FliP, but coexpression with the cytoplasmic domain of FliO did not enhance FliP expression levels. Overexpression of the cytoplasmic domain of FliO further rescued motility of strains deleted for the fliO gene expressing bypass mutations in FliP. These results suggest FliO maintains FliP stability through transmembrane domain interaction. The results also demonstrate that the cytoplasmic domain of FliO has functionality, and it presumably becomes structured while interacting with its binding partners., Author Summary The propeller-like flagella, which some bacteria use to swim, possess a specialized secretion apparatus, which is imbedded in the cell membrane for their formation. The components are highly conserved among flagella systems and also to the Type III secretion apparatus used by some bacteria in conjunction with virulence-associated needle complexes. The ubiquity of these secretion apparatuses and their function as intricate nanomachines has made them fascinating for biologists. The most studied flagellar system is that of Salmonella enterica, which consists of 6 integral membrane proteins: FlhA, FlhB, FliO, FliP, FliQ, and FliR. Among these proteins, FliO shows a sporadic distribution in bacteria, and its function is unknown, suggesting it might have a specialized role to play where it is present. In this study, we show that FliO has an important role in maintaining stability of FliP, which is a highly conserved member of the secretion apparatus. We have characterized the important regions of FliO through mutagenesis. We have shown that it is possible to bypass the effect of not producing the FliO protein, by encoding mutations within FliP or by overexpressing the cytoplasmic domain of FliO only.

BACKGROUND: A number of amyloid diseases involve deposition of extracellular protein aggregates, which are implicated in mechanisms of cell damage and death. However, the mechanisms involved remain poorly understood. METHODOLOGY/PRINCIPAL FINDINGS: Here we use the yeast prion protein Ure2 as a generic model to investigate how amyloid-like protein aggregates can enter mammalian cells and convey cytotoxicity. The effect of three different states of Ure2 protein (native dimer, protofibrils and mature fibrils) was tested on four mammalian cell lines (SH-SY5Y, MES23.5, HEK-293 and HeLa) when added extracellularly to the medium. Immunofluorescence using a polyclonal antibody against Ure2 showed that all three protein states could enter the four cell lines. In each case, protofibrils significantly inhibited the growth of the cells in a dose-dependent manner, fibrils showed less toxicity than protofibrils, while the native state had no effect on cell growth. This suggests that the structural differences between the three protein states lead to their different effects upon cells. Protofibrils of Ure2 increased membrane conductivity, altered calcium homeostasis, and ultimately induced apoptosis. The use of standard inhibitors suggested uptake into mammalian cells might occur via receptor-mediated endocytosis. In order to investigate this further, we used the chicken DT40 B cell line DKOR, which allows conditional expression of clathrin. Uptake into the DKOR cell-line was reduced when clathrin expression was repressed suggesting similarities between the mechanism of PrP uptake and the mechanism observed here for Ure2. CONCLUSIONS/SIGNIFICANCE: The results provide insight into the mechanisms by which amyloid aggregates may cause pathological effects in prion and amyloid diseases.

Genomic correlates of evolutionary adaptation to very low or very high optimal growth temperature (OGT) values have been the subject of many studies. Whereas these provided a protein-structural rationale of the activity and stability of globular proteins/enzymes, the point has been neglected that adaptation to extreme temperatures could also have resulted from an increased use of intrinsically disordered proteins (IDPs), which are resistant to these conditions in vitro. Contrary to these expectations, we found a conspicuously low level of structural disorder in bacteria of very high (and very low) OGT values. This paucity of disorder does not reflect phylogenetic relatedness, i.e. it is a result of genuine adaptation to extreme conditions. Because intrinsic disorder correlates with important regulatory functions, we asked how these bacteria could exist without IDPs by studying transcription factors, known to harbor a lot of function-related intrinsic disorder. Hyperthermophiles have much less transcription factors, which have reduced disorder compared to their mesophilic counterparts. On the other hand, we found by systematic categorization of proteins with long disordered regions that there are certain functions, such as translation and ribosome biogenesis that depend on structural disorder even in hyperthermophiles. In all, our observations suggest that adaptation to extreme conditions is achieved by a significant functional simplification, apparent at both the level of the genome and individual genes/proteins.

To understand how symmetric structures of many proteins are formed from asymmetric sequences, the proteins with two repeated beta-trefoil domains in Plant Cytotoxin B-chain family and all presently known beta-trefoil proteins are analyzed by structure-based multi-sequence alignments. The results show that all these proteins have similar key structural residues that are distributed symmetrically in their structures. These symmetric key structural residues are further analyzed in terms of inter-residues interaction numbers and B-factors. It is found that they can be distinguished from other residues and have significant propensities for structural framework. This indicates that these key structural residues may conduct the formation of symmetric structures although the sequences are asymmetric.

The folding pathway and rate coefficients of the folding of a knotted protein are calculated for a potential energy function with minimal energetic frustration. A kinetic transition network is constructed using the discrete path sampling approach, and the resulting potential energy surface is visualized by constructing disconnectivity graphs. Owing to topological constraints, the low-lying portion of the landscape consists of three distinct regions, corresponding to the native knotted state and to configurations where either the N- or C-terminus is not yet folded into the knot. The fastest folding pathways from denatured states exhibit early formation of the N-terminus portion of the knot and a rate-determining step where the C-terminus is incorporated. The low-lying minima with the N-terminus knotted and the C-terminus free therefore constitute an off-pathway intermediate for this model. The insertion of both the N- and C-termini into the knot occur late in the folding process, creating large energy barriers that are the rate limiting steps in the folding process. When compared to other protein folding proteins of a similar length, this system folds over six orders of magnitude more slowly., 19 pages

As nascent proteins are synthesized by the ribosome, they depart via an exit tunnel running through the center of the large subunit. The exit tunnel likely plays an important part in various aspects of translation. Although water plays a key role in many bio-molecular processes, the nature of water confined to the exit tunnel has remained unknown. Furthermore, solvent in biological cavities has traditionally been characterized as either a continuous dielectric fluid, or a discrete tightly bound molecule. Using atomistic molecular dynamics simulations, we predict that the thermodynamic and kinetic properties of water confined within the ribosome exit tunnel are quite different from this simple two-state model. We find that the tunnel creates a complex microenvironment for the solvent resulting in perturbed rotational dynamics and heterogenous dielectric behavior. This gives rise to a very rugged solvation landscape and significantly retarded solvent diffusion. We discuss how this non-bulk-like solvent is likely to affect important biophysical processes such as sequence dependent stalling, co-translational folding, and antibiotic binding. We conclude with a discussion of the general applicability of these results to other biological cavities., Author Summary The physical properties of water are of critical importance to biomolecular processes such as protein folding and protein-ligand binding. The predominant view of water from a structural biological perspective, however, is that of a two-state model: water is either a tightly bound molecule (as seen in x-ray crystallography) or a continuous bulk fluid of constant dielectric. When water is confined on the nanometer length scale (as it is in a number of critically important biological cavities) the two-state characterization is likely to be inadequate. Using molecular dynamics simulations we report on the thermodynamic and kinetic properties of water confined to the ribosome exit tunnel. We show that this solvent is slow-diffusing and semi-structured, sampling a great variety of states between its characteristic ‘end-points’ of bulk-like and tightly bound. We conclude with a discussion of how the properties of this nano-confined solvent can affect the interpretation of recent experimental results regarding sequence dependant stalling and co-translational folding in the ribosome exit tunnel.

FSD-1, a designed small ultrafast folder with a ββα fold, has been actively studied in the last few years as a model system for studying protein folding mechanisms and for testing of the accuracy of computational models. The suitability of this protein to describe the folding of naturally occurring α/β proteins has recently been challenged based on the observation that the melting transition is very broad, with ill-resolved baselines. Using molecular dynamics simulations with the AMBER protein force field (ff96) coupled with the implicit solvent model (IGB = 5), we shed new light into the nature of this transition and resolve the experimental controversies. We show that the melting transition corresponds to the melting of the protein as a whole, and not solely to the helix-coil transition. The breadth of the folding transition arises from the spread in the melting temperatures (from ∼325 K to ∼302 K) of the individual transitions: formation of the hydrophobic core, β-hairpin and tertiary fold, with the helix formed earlier. Our simulations initiated from an extended chain accurately predict the native structure, provide a reasonable estimate of the transition barrier height, and explicitly demonstrate the existence of multiple pathways and multiple transition states for folding. Our exhaustive sampling enables us to assess the quality of the Amber ff96/igb5 combination and reveals that while this force field can predict the correct native fold, it nonetheless overstabilizes the α-helix portion of the protein (Tm = ∼387K) as well as the denatured structures., Author Summary The protein folding process, in which a linear chain of amino acids reaches its biologically active three-dimensional shape, is fundamental to life. Small “ultrafast” folders, proteins that fold in microseconds, have received considerable attention, because these proteins serve as model systems for the folding of larger proteins, and thus permit a testing of the accuracy of computational models as well as an assessment of protein folding theories. FSD-1, a designed small ultrafast folder with a ββα fold, has been actively studied in the last few years as a model system for mixed α/β fold proteins. The suitability of this protein to describe the folding of naturally occurring proteins has however recently been challenged based on the observation that the melting transition is very broad, with ill-resolved baselines. Prior simulations have not been successful in providing an interpretation of this broad melting transition. In the present study, our extensive molecular dynamics simulations using the AMBER protein force field (ff96) coupled with the implicit solvent model (IGB = 5) shed new light on the nature of the folding transition of this protein, as well as reveal the strengths and weaknesses of the force field in predicting the thermodynamics and kinetics of folding.

The HIV-1 viral infectivity factor (Vif) protein recruits an E3 ubiquitin ligase complex, comprising the cellular proteins elongin B and C (EloBC), cullin 5 (Cul5) and RING-box 2 (Rbx2), to the anti-viral proteins APOBEC3G (A3G) and APOBEC3F (A3F) and induces their polyubiquitination and proteasomal degradation. In this study, we used purified proteins and direct in vitro binding assays, isothermal titration calorimetry and NMR spectroscopy to describe the molecular mechanism for assembly of the Vif-EloBC ternary complex. We demonstrate that Vif binds to EloBC in two locations, and that both interactions induce structural changes in the SOCS box of Vif as well as EloBC. In particular, in addition to the previously established binding of Vif's BC box to EloC, we report a novel interaction between the conserved Pro-Pro-Leu-Pro motif of Vif and the C-terminal domain of EloB. Using cell-based assays, we further show that this interaction is necessary for the formation of a functional ligase complex, thus establishing a role of this motif. We conclude that HIV-1 Vif engages EloBC via an induced-folding mechanism that does not require additional co-factors, and speculate that these features distinguish Vif from other EloBC specificity factors such as cellular SOCS proteins, and may enhance the prospects of obtaining therapeutic inhibitors of Vif function., Author Summary HIV-1 is the etiologic agent of AIDS. Current therapies are based on cocktails of anti-viral drugs that inhibit viral enzymes essential for virus replication, but this strategy has several shortcomings, including the development of drug-resistant virus strains. Consequently, pharmacologic strategies that interfere with additional aspects of HIV-1 replication have the potential to enhance HIV-1 treatments. The HIV-1 Vif protein is a promising target for the development of new anti-HIV-1 therapeutics; it functions to counteract the cellular proteins A3G and A3F, two components of a human anti-viral defence mechanism. Vif accomplishes this by hijacking a cellular complex (comprising the proteins EloB, EloC, Cul5 and Rbx2), which then eliminates A3G and A3F from infected cells by degradation, therefore evading their anti-viral effect. Here, we used purified proteins to reconstitute in vitro the recruitment of this complex by HIV-1 Vif. Using structural and biochemical methods, we dissected the different events involved in Vif's interaction with the EloBC complex. Our results reveal fundamental differences with cellular proteins known to recruit this complex, suggesting that Vif possesses unique features that could be targeted by pharmacologic intervention, without disturbing normal cell functions. The assays reported here could be utilized for the discovery of such inhibitors.

One of the applications of Molecular Dynamics (MD) simulations is to explore the energetic barriers to mechanical unfolding of proteins such as occurs in response to the mechanical pulling of single molecules in Atomic Force Microscopy (AFM) experiments. Although Steered Molecular Dynamics simulations have provided microscopic details of the unfolding process during the pulling, the simulated forces required for unfolding are typically far in excess of the measured values. To rectify this, we have developed the Pulsed Unconstrained Fluctuating Forces (PUFF) method, which induces constant-momentum motions by applying forces directly to the instantaneous velocity of selected atoms in a protein system. The driving forces are applied in pulses, which allows the system to relax between pulses, resulting in more accurate unfolding force estimations than in previous methods. In the cases of titin, ubiquitin and e2lip3, the PUFF trajectories produce force fluctuations that agree quantitatively with AFM experiments. Another useful property of PUFF is that simulations get trapped if the target momentum is too low, simplifying the discovery and analysis of unfolding intermediates.

International audience; BACKGROUND: Protein fold recognition usually relies on a statistical model of each fold; each model is constructed from an ensemble of natural sequences belonging to that fold. A complementary strategy may be to employ sequence ensembles produced by computational protein design. Designed sequences can be more diverse than natural sequences, possibly avoiding some limitations of experimental databases. METHODOLOGY/PRINCIPAL FINDINGS: WE EXPLORE THIS STRATEGY FOR FOUR SCOP FAMILIES: Small Kunitz-type inhibitors (SKIs), Interleukin-8 chemokines, PDZ domains, and large Caspase catalytic subunits, represented by 43 structures. An automated procedure is used to redesign the 43 proteins. We use the experimental backbones as fixed templates in the folded state and a molecular mechanics model to compute the interaction energies between sidechain and backbone groups. Calculations are done with the Proteins@Home volunteer computing platform. A heuristic algorithm is used to scan the sequence and conformational space, yielding 200,000-300,000 sequences per backbone template. The results confirm and generalize our earlier study of SH2 and SH3 domains. The designed sequences ressemble moderately-distant, natural homologues of the initial templates; e.g., the SUPERFAMILY, profile Hidden-Markov Model library recognizes 85% of the low-energy sequences as native-like. Conversely, Position Specific Scoring Matrices derived from the sequences can be used to detect natural homologues within the SwissProt database: 60% of known PDZ domains are detected and around 90% of known SKIs and chemokines. Energy components and inter-residue correlations are analyzed and ways to improve the method are discussed. CONCLUSIONS/SIGNIFICANCE: For some families, designed sequences can be a useful complement to experimental ones for homologue searching. However, improved tools are needed to extract more information from the designed profiles before the method can be of general use.

To obtain insight into the functions of proteins and their specific roles, it is important to establish efficient procedures for exploring the states that encapsulate their conformational space. Global Protein folding State mapping by multivariate NMR (GPS NMR) is a powerful high-throughput method that provides such an overview. GPS NMR exploits the unique ability of NMR to simultaneously record signals from individual hydrogen atoms in complex macromolecular systems and of multivariate analysis to describe spectral variations from these by a few variables for establishment of, and positioning in, protein-folding state maps. The method is fast, sensitive, and robust, and it works without isotope-labelling. The unique capabilities of GPS NMR to identify different folding states and to compare different unfolding processes are demonstrated by mapping of the equilibrium folding space of bovine alpha-lactalbumin in the presence of the anionic surfactant sodium dodecyl sulfate, SDS, and compare these with other surfactants, acid, denaturants and heat.

Amino acid replacements at dozens of positions in the dimeric protein human, Cu,Zn superoxide dismutase (SOD1) can cause amyotrophic lateral sclerosis (ALS). Although it has long been hypothesized that these mutations might enhance the populations of marginally-stable aggregation-prone species responsible for cellular toxicity, there has been little quantitative evidence to support this notion. Perturbations of the folding free energy landscapes of metal-free versions of five ALS-inducing variants, A4V, L38V, G93A, L106V and S134N SOD1, were determined with a global analysis of kinetic and thermodynamic folding data for dimeric and stable monomeric versions of these variants. Utilizing this global analysis approach, the perturbations on the global stability in response to mutation can be partitioned between the monomer folding and association steps, and the effects of mutation on the populations of the folded and unfolded monomeric states can be determined. The 2- to 10-fold increase in the population of the folded monomeric state for A4V, L38V and L106V and the 80- to 480-fold increase in the population of the unfolded monomeric states for all but S134N would dramatically increase their propensity for aggregation through high-order nucleation reactions. The wild-type-like populations of these states for the metal-binding region S134N variant suggest that even wild-type SOD1 may also be prone to aggregation in the absence of metals.

Multiple variants of the AMBER all-atom force field were quantitatively evaluated with respect to their ability to accurately characterize helix-coil equilibria in explicit solvent simulations. Using a global distributed computing network, absolute conformational convergence was achieved for large ensembles of the capped A(21) and F(s) helical peptides. Further assessment of these AMBER variants was conducted via simulations of a flexible 164-residue five-helix-bundle protein, apolipophorin-III, on the 100 ns timescale. Of the contemporary potentials that had not been assessed previously, the AMBER-99SB force field showed significant helix-destabilizing tendencies, with beta bridge formation occurring in helical peptides, and unfolding of apolipophorin-III occurring on the tens of nanoseconds timescale. The AMBER-03 force field, while showing adequate helical propensities for both peptides and stabilizing apolipophorin-III, (i) predicts an unexpected decrease in helicity with ALA-->ARG(+) substitution, (ii) lacks experimentally observed 3(10) helical content, and (iii) deviates strongly from average apolipophorin-III NMR structural properties. As is observed for AMBER-99SB, AMBER-03 significantly overweighs the contribution of extended and polyproline backbone configurations to the conformational equilibrium. In contrast, the AMBER-99phi force field, which was previously shown to best reproduce experimental measurements of the helix-coil transition in model helical peptides, adequately stabilizes apolipophorin-III and yields both an average gyration radius and polar solvent exposed surface area that are in excellent agreement with the NMR ensemble.

The presence of amyloid deposits consisting primarily of Amyloid-β (Aβ) fibril in the brain is a hallmark of Alzheimer's disease (AD). The morphologies of these fibrils are exquisitely sensitive to environmental conditions. Using molecular dynamics simulations combined with data from previously published solid-state NMR experiments, we propose the first atomically detailed structures of two asymmetric polymorphs of the Aβ9-40 peptide fibril. The first corresponds to synthetic fibrils grown under quiescent conditions and the second to fibrils derived from AD patients' brain-extracts. Our core structure in both fibril structures consists of a layered structure in which three cross-β subunits are arranged in six tightly stacked β-sheet layers with an antiparallel hydrophobic-hydrophobic and an antiparallel polar-polar interface. The synthetic and brain-derived structures differ primarily in the side-chain orientation of one β-strand. The presence of a large and continually exposed hydrophobic surface (buried in the symmetric agitated Aβ fibrils) may account for the higher toxicity of the asymmetric fibrils. Our model explains the effects of external perturbations on the fibril lateral architecture as well as the fibrillogenesis inhibiting action of amphiphilic molecules., Author Summary Amyloid diseases are characterized by the presence of amyloid fibrils on organs and tissue in the body. Alzheimer's disease, Parkinson's diseases and Type II Diabetes are all examples of amyloid diseases. Determining the structure of amyloid fibrils is critical for understanding the mechanism of fibril formation as well as for the design of inhibitor molecules that can prevent aggregation. In the case of the Alzheimer Amyloid-β (Aβ) peptide, the structure of fibrils grown under conditions of mechanical agitation has been elucidated from a combination of simulation and experiments. However, the structures of the asymmetric quiescent Aβ fibrils (grown under conditions akin to physiological conditions) and of Alzheimer's brain–derived fibrils are not known. In this paper, we propose the first atomically detailed structures of these two fibrils, using molecular dynamics simulations combined with data from previously published experiments. In additions, we suggest a unifying lateral growth mechanism that explains the increased toxicity of quiescent Aβ fibrils, the effects of external perturbations on fibril lateral architecture and the inhibition mechanism of the small molecule inhibitors on fibril formation.

Kinetically stable proteins, those whose stability is derived from their slow unfolding kinetics and not thermodynamics, are examples of evolution's best attempts at suppressing unfolding. Especially in highly proteolytic environments, both partially and fully unfolded proteins face potential inactivation through degradation and/or aggregation, hence, slowing unfolding can greatly extend a protein's functional lifetime. The prokaryotic serine protease α-lytic protease (αLP) has done just that, as its unfolding is both very slow (t1/2 ∼1 year) and so cooperative that partial unfolding is negligible, providing a functional advantage over its thermodynamically stable homologs, such as trypsin. Previous studies have identified regions of the domain interface as critical to αLP unfolding, though a complete description of the unfolding pathway is missing. In order to identify the αLP unfolding pathway and the mechanism for its extreme cooperativity, we performed high temperature molecular dynamics unfolding simulations of both αLP and trypsin. The simulated αLP unfolding pathway produces a robust transition state ensemble consistent with prior biochemical experiments and clearly shows that unfolding proceeds through a preferential disruption of the domain interface. Through a novel method of calculating unfolding cooperativity, we show that αLP unfolds extremely cooperatively while trypsin unfolds gradually. Finally, by examining the behavior of both domain interfaces, we propose a model for the differential unfolding cooperativity of αLP and trypsin involving three key regions that differ between the kinetically stable and thermodynamically stable classes of serine proteases., Author Summary Proteins, synthesized as linear polymers of amino acids, fold up into compact native states, burying their hydrophobic amino acids into their interiors. Protein folding minimizes the non-specific interactions that unfolded protein chains can make, which include aggregation with other proteins and degradation by proteases. Unfortunately, even in the native state, proteins can partially unfold, opening up regions of their structure and making these adverse events possible. Some proteins, particularly those in harsh environments full of proteases, have evolved to virtually eliminate partial unfolding, significantly reducing their rate of degradation. This elimination of partial unfolding is termed “cooperative,” because unfolding is an all-or-none process. One class of proteins has diverged into two families, one bacterial and highly cooperative and the other animal and non-cooperative. We have used detailed simulations of unfolding for members of each family, α-lytic protease (bacterial) and trypsin (animal) to understand the unfolding pathways of each and the mechanism for the differential unfolding cooperativity. Our results explain prior biochemical experiments, reproduce the large difference in unfolding cooperativity between the families, and point to the interface between α-lytic protease's two domains as essential to establishing unfolding cooperativity. As seen in an unrelated protein family, generation of a cooperative domain interface may be a common evolutionary response for ensuring the highest protein stability.

A didactic model is presented to illustrate how the effect of macromolecular crowding on protein folding and association is modeled using current analytical theory and discrete molecular dynamics. While analytical treatments of crowding may consider the effect as a potential of average force acting to compress a polypeptide chain into a compact state, the use of simulations enables the presence of crowding reagents to be treated explicitly. Using an analytically solvable toy model for protein folding, an approximate statistical thermodynamic method is directly compared to simulation in order to gauge the effectiveness of current analytical crowding descriptions. Both methodologies are in quantitative agreement under most conditions, indication that both current theory and simulation methods are capable of recapitulating aspects of protein folding even by utilizing a simplistic protein model.

Mutation bias in prokaryotes varies from extreme adenine and thymine (AT) in obligatory endosymbiotic or parasitic bacteria to extreme guanine and cytosine (GC), for instance in actinobacteria. GC mutation bias deeply influences the folding stability of proteins, making proteins on the average less hydrophobic and therefore less stable with respect to unfolding but also less susceptible to misfolding and aggregation. We study a model where proteins evolve subject to selection for folding stability under given mutation bias, population size, and neutrality. We find a non-neutral regime where, for any given population size, there is an optimal mutation bias that maximizes fitness. Interestingly, this optimal GC usage is small for small populations, large for intermediate populations and around 50% for large populations. This result is robust with respect to the definition of the fitness function and to the protein structures studied. Our model suggests that small populations evolving with small GC usage eventually accumulate a significant selective advantage over populations evolving without this bias. This provides a possible explanation to the observation that most species adopting obligatory intracellular lifestyles with a consequent reduction of effective population size shifted their mutation spectrum towards AT. The model also predicts that large GC usage is optimal for intermediate population size. To test these predictions we estimated the effective population sizes of bacterial species using the optimal codon usage coefficients computed by dos Reis et al. and the synonymous to non-synonymous substitution ratio computed by Daubin and Moran. We found that the population sizes estimated in these ways are significantly smaller for species with small and large GC usage compared to species with no bias, which supports our prediction., UB acknowledges financial support from the Spanish Science and Innovation Ministry through the Ramón y Cajal program and through the projects BIO2008-04384 and CSD2006-00023, and a stay at the Aspen Center for Physics where a first version of this work was written. Our collaboration was facilitated through the program ‘‘Acciones Integradas Espana-Alemania’’ of the Spanish Science and Innovation Ministry, project HA2006-0044, and of the Deutscher Akademischer Austauschdienst project D/06/12848. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Polyglutamine (polyQ) expansion in exon1 (XN1) of the huntingtin protein is linked to Huntington's disease. When the number of glutamines exceeds a threshold of approximately 36–40 repeats, XN1 can readily form amyloid aggregates similar to those associated with disease. Many experiments suggest that misfolding of monomeric XN1 plays an important role in the length-dependent aggregation. Elucidating the misfolding of a XN1 monomer can help determine the molecular mechanism of XN1 aggregation and potentially help develop strategies to inhibit XN1 aggregation. The flanking sequences surrounding the polyQ region can play a critical role in determining the structural rearrangement and aggregation mechanism of XN1. Few experiments have studied XN1 in its entirety, with all flanking regions. To obtain structural insights into the misfolding of XN1 toward amyloid aggregation, we perform molecular dynamics simulations on monomeric XN1 with full flanking regions, a variant missing the polyproline regions, which are hypothesized to prevent aggregation, and an isolated polyQ peptide (Qn). For each of these three constructs, we study glutamine repeat lengths of 23, 36, 40 and 47. We find that polyQ peptides have a positive correlation between their probability to form a β-rich misfolded state and their expansion length. We also find that the flanking regions of XN1 affect its probability to^x_page_count=28 form a β-rich state compared to the isolated polyQ. Particularly, the polyproline regions form polyproline type II helices and decrease the probability of the polyQ region to form a β-rich state. Additionally, by lengthening polyQ, the first N-terminal 17 residues are more likely to adopt a β-sheet conformation rather than an α-helix conformation. Therefore, our molecular dynamics study provides a structural insight of XN1 misfolding and elucidates the possible role of the flanking sequences in XN1 aggregation., Author Summary Huntington's Disease is a neurodegenerative disorder associated with protein aggregation in neurons. The aggregates formed are thought to lead to neurotoxicity and cell death. Understanding the molecular structure of these aggregates may lead to strategies to inhibit aggregation. Exon 1 (XN1) of the huntingtin protein is critical for aggregate formation. This polypeptide has a naturally occurring polyglutamine sequence (polyQ), which is elongated in patients afflicted with the disease. The polyQ region in XN1 has several flanking sequences with distinct physicochemical properties, including the N-terminal 17 residues, two polyproline regions, and C-terminal sequences, that may affect its overall structure and aggregation. What is the overall structure of XN1, and what structural effects do the neighboring sequences have on each other and polyQ? We address these questions by studying computational models of various polypeptides, including XN1 and three mutant forms associated with Huntington's Disease. Certain neighboring sequences are found to inhibit aggregation, while others may be recruited by polyQ to form aggregates. Our results suggest the role that the flanking sequences may play in XN1 aggregation and may subsequently guide future structural models of XN1 aggregation.

BACKGROUND: Alpha synuclein (alpha-Syn) is the main component of Lewy bodies which are associated with several neurodegenerative diseases such as Parkinson's disease. While the glycation with D-glucose that results in alpha-Syn misfold and aggregation has been studied, the effects of glycation with D-ribose on alpha-Syn have not been investigated. METHODOLOGY/PRINCIPAL FINDINGS: Here, we show that ribosylation induces alpha-Syn misfolding and generates advanced glycation end products (AGEs) which form protein molten globules with high cytotoxcity. Results from native- and SDS-PAGE showed that D-ribose reacted rapidly with alpha-Syn, leading to dimerization and polymerization. Trypsin digestion and sequencing analysis revealed that during ribosylation the lysinyl residues (K(58), K(60), K(80), K(96), K(97) and K(102)) in the C-terminal region reacted more quickly with D-ribose than those of the N-terminal region. Using Western blotting, AGEs resulting from the glycation of alpha-Syn were observed within 24 h in the presence of D-ribose, but were not observed in the presence of D-glucose. Changes in fluorescence at 410 nm demonstrated again that AGEs were formed during early ribosylation. Changes in the secondary structure of ribosylated alpha-Syn were not clearly detected by CD spectrometry in studies on protein conformation. However, intrinsic fluorescence at 310 nm decreased markedly in the presence of D-ribose. Observations with atomic force microscopy showed that the surface morphology of glycated alpha-Syn looked like globular aggregates. thioflavin T (ThT) fluorescence increased during alpha-Syn incubation regardless of ribosylation. As incubation time increased, ribosylation of alpha-Syn resulted in a blue-shift (approximately 100 nm) in the fluorescence of ANS. The light scattering intensity of ribosylated alpha-Syn was not markedly different from native alpha-Syn, suggesting that ribosylated alpha-Syn is present as molten protein globules. Ribosylated products had a high cytotoxicity to SH-SY5Y cells, leading to LDH release and increase in the levels of reactive oxygen species (ROS). CONCLUSIONS/SIGNIFICANCE: alpha-Syn is rapidly glycated in the presence of D-ribose generating molten globule-like aggregations which cause cell oxidative stress and result in high cytotoxicity.

Stability of membrane protein is crucial during protein purification and crystallization as well as in the fabrication of protein-based devices. Several recent studies have examined how various surfactants can stabilize membrane proteins out of their native membrane environment. However, there is still no single surfactant that can be universally employed for all membrane proteins. Because of the lack of knowledge on the interaction between surfactants and membrane proteins, the choice of a surfactant for a specific membrane protein remains purely empirical. Here we report that a group of short amphiphilic peptides improve the thermal stability of the multi-domain protein complex photosystem-I (PS-I) in aqueous solution and that the peptide surfactants have obvious advantages over other commonly used alkyl chain based surfactants. Of all the short peptides studied, Ac-I(5)K(2)-CONH(2) (I(5)K(2)) showed the best stabilizing effect by enhancing the melting temperature of PS-I from 48.0 degrees C to 53.0 degrees C at concentration of 0.65 mM and extending the half life of isolated PS-I significantly. AFM experiments showed that PS-I/I(5)K(2)/Triton X-100 formed large and stable vesicles and thus provide interfacial environment mimicking that of native membranes, which may partly explain why I(5)K(2) enhanced the thermal stability of PS-I. Hydrophobic and hydrophilic group length of I(x)K(y) had an important influence on the stabilization of PS-I. Our results showed that longer hydrophobic group was more effective in stabilizing PS-I. These simple short peptides therefore exhibit significant potential for applications in membrane protein studies.

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease., PLoS ONE, 5 (10), ISSN:1932-6203

BACKGROUND: The conformational conversion of the host-derived cellular prion protein (PrP(C)) into the disease-associated scrapie isoform (PrP(Sc)) is responsible for the pathogenesis of transmissible spongiform encephalopathies (TSEs). Various single-point mutations in PrP(C)s could cause structural changes and thereby distinctly influence the conformational conversion. Elucidation of the differences between the wild-type rabbit PrP(C) (RaPrP(C)) and various mutants would be of great help to understand the ability of RaPrP(C) to be resistant to TSE agents. METHODOLOGY/PRINCIPAL FINDINGS: We determined the solution structure of the I214V mutant of RaPrP(C)(91-228) and detected the backbone dynamics of its structured C-terminal domain (121-228). The I214V mutant displays a visible shift of surface charge distribution that may have a potential effect on the binding specificity and affinity with other chaperones. The number of hydrogen bonds declines dramatically. Urea-induced transition experiments reveal an obvious decrease in the conformational stability. Furthermore, the NMR dynamics analysis discloses a significant increase in the backbone flexibility on the pico- to nanosecond time scale, indicative of lower energy barrier for structural rearrangement. CONCLUSIONS/SIGNIFICANCE: Our results suggest that both the surface charge distribution and the intrinsic backbone flexibility greatly contribute to species barriers for the transmission of TSEs, and thereby provide valuable hints for understanding the inability of the conformational conversion for RaPrP(C).

BACKGROUND: To examine the structural organisation of heterodimeric soluble guanylyl cyclase (sGC) Förster resonance energy transfer (FRET) was measured between fluorescent proteins fused to the amino- and carboxy-terminal ends of the sGC beta1 and alpha subunits. METHODOLOGY/PRINCIPAL FINDINGS: Cyan fluorescent protein (CFP) was used as FRET donor and yellow fluorescent protein (YFP) as FRET acceptor. After generation of recombinant baculovirus, fluorescent-tagged sGC subunits were co-expressed in Sf9 cells. Fluorescent variants of sGC were analyzed in vitro in cytosolic fractions by sensitized emission FRET. Co-expression of the amino-terminally tagged alpha subunits with the carboxy-terminally tagged beta1 subunit resulted in an enzyme complex that showed a FRET efficiency of 10% similar to fluorescent proteins separated by a helix of only 48 amino acids. Because these findings indicated that the amino-terminus of the alpha subunits is close to the carboxy-terminus of the beta1 subunit we constructed fusion proteins where both subunits are connected by a fluorescent protein. The resulting constructs were not only fluorescent, they also showed preserved enzyme activity and regulation by NO. CONCLUSIONS/SIGNIFICANCE: Based on the ability of an amino-terminal fragment of the beta1 subunit to inhibit activity of an heterodimer consisting only of the catalytic domains (alphacatbetacat), Winger and Marletta (Biochemistry 2005, 44:4083-90) have proposed a direct interaction of the amino-terminal region of beta1 with the catalytic domains. In support of such a concept of "trans" regulation of sGC activity by the H-NOX domains our results indicate that the domains within sGC are organized in a way that allows for direct interaction of the amino-terminal regulatory domains with the carboxy-terminal catalytic region. In addition, we constructed "fluorescent-conjoined" sGC's by fusion of the alpha amino-terminus to the beta1 carboxy-terminus leading to a monomeric, fluorescent and functional enzyme complex. To our knowledge this represents the first example where a fluorescent protein links two different subunits of a higher ordered complex to yield a stoichometrically fixed functionally active monomer.

Characterization of the length dependence of end-to-end loop-closure kinetics in unfolded polypeptide chains provides an understanding of early steps in protein folding. Here, loop-closure in poly-glycine-serine peptides is investigated by combining single-molecule fluorescence spectroscopy with molecular dynamics simulation. For chains containing more than 10 peptide bonds loop-closing rate constants on the 20–100 nanosecond time range exhibit a power-law length dependence. However, this scaling breaks down for shorter peptides, which exhibit slower kinetics arising from a perturbation induced by the dye reporter system used in the experimental setup. The loop-closure kinetics in the longer peptides is found to be determined by the formation of intra-peptide hydrogen bonds and transient β-sheet structure, that accelerate the search for contacts among residues distant in sequence relative to the case of a polypeptide chain in which hydrogen bonds cannot form. Hydrogen-bond-driven polypeptide-chain collapse in unfolded peptides under physiological conditions found here is not only consistent with hierarchical models of protein folding, that highlights the importance of secondary structure formation early in the folding process, but is also shown to speed up the search for productive folding events., Author Summary In studies of protein folding evidence exists for early compaction in the unfolded state, although it is unclear whether these compact conformations contain specific secondary structures (through hydrophilic interactions) or whether compaction is a non-specific hydrophobic-driven effect. Here we combine single-molecule fluorescence spectroscopy and molecular dynamics simulation to demonstrate peptide hydrogen-bond-driven polypeptide-chain collapse involving secondary structure formation as the key process in the early stage of folding. Partial structuring in unfolded polypeptide chains is shown to lead to faster contact formation kinetics than would be expected if the unfolded state were populated by featureless random-coils.

BACKGROUND: The aggregation of the baker's yeast prion Ure2p is at the origin of the [URE3] trait. The Q- and N-rich N-terminal part of the protein is believed to drive Ure2p assembly into fibrils of amyloid nature and the fibrillar forms of full-length Ure2p and its N-terminal part generated in vitro have been shown to induce [URE3] occurrence when introduced into yeast cells. This has led to the view that the fibrillar form of the N-terminal part of the protein is sufficient for the recruitment of constitutive Ure2p and that it imprints its amyloid structure to full-length Ure2p. RESULTS: Here we generate a set of Ure2p N-terminal fragments, document their assembly and structural properties and compare them to that of full-length Ure2p. We identify the minimal region critical for the assembly of Ure2p N-terminal part into amyloids and show that such fibrils are unable to seed the assembly of full length Ure2p unlike fibrils made of intact Ure2p. CONCLUSION: Our results clearly indicate that fibrillar Ure2p shares no structural similarities with the amyloid fibrils made of Ure2p N-terminal part. Our results further suggest that the induction of [URE3] by fibrils made of full-length Ure2p is likely the consequence of fibrils growth by depletion of cytosolic Ure2p while it is the consequence of de novo formation of prion particles following, for example, titration within the cells of a specific set of molecular chaperones when fibrils made of Ure2p N-terminal domain are introduced within the cytoplasm.

The formation of insoluble protein deposits in human tissues is linked to the onset of more than 40 different disorders, ranging from dementia to diabetes. In these diseases, the proteins usually self-assemble into ordered β-sheet enriched aggregates known as amyloid fibrils. Here we study the structure of the inclusions formed by maize transglutaminase (TGZ) in the chloroplasts of tobacco transplastomic plants and demonstrate that they have an amyloid-like nature. Together with the evidence of amyloid structures in bacteria and fungi our data argue that amyloid formation is likely a ubiquitous process occurring across the different kingdoms of life. The discovery of amyloid conformations inside inclusions of genetically modified plants might have implications regarding their use for human applications., This work was supported by grants BIO2007-68046, BFU2006-15115-CO2-01/BMC and BFU2009-08575 from Ministerio de Ciencia e Innovacion (Spain), by grant 2009-SGR 760 from AGAUR(Agencia de Gestio d'Ajuts Universitaris i de Recerca-Generalitat de Catalunya). AV was beneficiary of a PIF fellowship from Universidad Autonoma de Barcelona (UAB). RS was beneficiary of an I3 contract (UAB-Generalitat de Catalunya). JG was beneficiary of a FPI fellowship awarded by Generalitat de Catalunya. SV has been granted with an ICREA ACADEMIA award (ICREA-Institucio Catalana de Recerca i Estudis Avancats). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Aptamers are useful molecular recognition tools in research, diagnostics, and therapy. Despite promising results in other fields, aptamer use has remained scarce in amyloid research, including Alzheimer's disease (AD). AD is a progressive neurodegenerative disease believed to be caused by neurotoxic amyloid beta-protein (Abeta) oligomers. Abeta oligomers therefore are an attractive target for development of diagnostic and therapeutic reagents. We used covalently-stabilized oligomers of the 40-residue form of Abeta (Abeta40) for aptamer selection. Despite gradually increasing the stringency of selection conditions, the selected aptamers did not recognize Abeta40 oligomers but reacted with fibrils of Abeta40, Abeta42, and several other amyloidogenic proteins. Aptamer reactivity with amyloid fibrils showed some degree of protein-sequence dependency. Significant fibril binding also was found for the naïve library and could not be eliminated by counter-selection using Abeta40 fibrils, suggesting that aptamer binding to amyloid fibrils was RNA-sequence-independent. Aptamer binding depended on fibrillogenesis and showed a lag phase. Interestingly, aptamers detected fibril formation with > or =15-fold higher sensitivity than thioflavin T (ThT), revealing substantial beta-sheet and fibril formation undetected by ThT. The data suggest that under physiologic conditions, aptamers for oligomeric forms of amyloidogenic proteins cannot be selected due to high, non-specific affinity of oligonucleotides for amyloid fibrils. Nevertheless, the high sensitivity, whereby aptamers detect beta-sheet formation, suggests that they can serve as superior amyloid recognition tools.

A robust marker to describe mass, hydrophobicity and polarizability distribution holds the key to deciphering structural and folding constraints within proteins. Since each of these distributions is inhomogeneous in nature, the construct should be sensitive in describing the patterns therein. We show, for the first time, that the hydrophobicity and polarizability distributions in protein interior follow fractal scaling. It is found that (barring 'all-alpha') all the major structural classes of proteins have an amount of unused hydrophobicity left in them. This amount of untapped hydrophobicity is observed to be greater in thermophilic proteins, than that in their (structurally aligned) mesophilic counterparts. 'All-beta'(thermophilic, mesophilic alike) proteins are found to have maximum amount of unused hydrophobicity, while 'all-alpha' proteins have been found to have minimum polarizability. A non-trivial dependency is observed between dielectric constant and hydrophobicity distributions within (alpha+beta) and 'all-alpha' proteins, whereas absolutely no dependency is found between them in the 'all-beta' class. This study proves that proteins are not as optimally packed as they are supposed to be. It is also proved that origin of alpha-helices are possibly not hydrophobic but electrostatic; whereas beta-sheets are predominantly hydrophobic in nature. Significance of this study lies in protein engineering studies; because it quantifies the extent of packing that ensures protein functionality. It shows that myths regarding protein interior organization might obfuscate our knowledge of actual reality. However, if the later is studied with a robust marker of strong mathematical basis, unknown correlations can still be unearthed; which help us to understand the nature of hydrophobicity, causality behind protein folding, and the importance of anisotropic electrostatics in stabilizing a highly complex structure named 'proteins'.

Chromosomal translocations, which often generate chimeric proteins by fusing segments of two distinct genes, represent the single major genetic aberration leading to cancer. We suggest that the unifying theme of these events is a high level of intrinsic structural disorder, enabling fusion proteins to evade cellular surveillance mechanisms that eliminate misfolded proteins. Predictions in 406 translocation-related human proteins show that they are significantly enriched in disorder (43.3% vs. 20.7% in all human proteins), they have fewer Pfam domains, and their translocation breakpoints tend to avoid domain splitting. The vicinity of the breakpoint is significantly more disordered than the rest of these already highly disordered fusion proteins. In the unlikely event of domain splitting in fusion it usually spares much of the domain or splits at locations where the newly exposed hydrophobic surface area approximates that of an intact domain. The mechanisms of action of fusion proteins suggest that in most cases their structural disorder is also essential to the acquired oncogenic function, enabling the long-range structural communication of remote binding and/or catalytic elements. In this respect, there are three major mechanisms that contribute to generating an oncogenic signal: (i) a phosphorylation site and a tyrosine-kinase domain are fused, and structural disorder of the intervening region enables intramolecular phosphorylation (e.g., BCR-ABL); (ii) a dimerisation domain fuses with a tyrosine kinase domain and disorder enables the two subunits within the homodimer to engage in permanent intermolecular phosphorylations (e.g., TFG-ALK); (iii) the fusion of a DNA-binding element to a transactivator domain results in an aberrant transcription factor that causes severe misregulation of transcription (e.g. EWS-ATF). Our findings also suggest novel strategies of intervention against the ensuing neoplastic transformations., Author Summary Chromosomal translocations generate chimeric proteins by fusing segments of two distinct genes and are frequently associated with cancer. The proteins involved are large and fairly heterogeneous in sequence and typically have only a few dispersed structural domains connected by long uncharacterized regions. It has never been studied from a structural perspective how these chimeras survive losing significant portions of the original proteins and acquire new oncogenic functions. By analyzing a collection of 406 human translocation proteins we show here that the answer to both questions lies to a large extent in the high level of structural disorder in the fusion partner proteins (on average, they are twice as disordered as all human proteins). The translocation breakpoints usually avoid globular domains. In rare cases when a globular domain is truncated by the fusion, it happens at a location in the domain where the hydrophobicity exposed by the split is favorable (i.e., not too high). Disorder on average is significantly higher in the vicinity of the breakpoint than in the rest of the fusion proteins. Disorder also plays a pivotal role in the acquired oncogenic function by bringing distant/disparate fusion segments together that enables novel intra- and/or intermolecular interactions.

Non-local hydrogen bonding interactions between main chain amide hydrogen atoms and polar side chain acceptors that bracket consecutive betaalpha or alphabeta elements of secondary structure in alphaTS from E. coli, a TIM barrel protein, have previously been found to contribute 4-6 kcal mol(-1) to the stability of the native conformation. Experimental analysis of similar betaalpha-hairpin clamps in a homologous pair of TIM barrel proteins of low sequence identity, IGPS from S. solfataricus and E. coli, reveals that this dramatic enhancement of stability is not unique to alphaTS. A survey of 71 TIM barrel proteins demonstrates a 4-fold symmetry for the placement of betaalpha-hairpin clamps, bracing the fundamental betaalphabeta building block and defining its register in the (betaalpha)(8) motif. The preferred sequences and locations of betaalpha-hairpin clamps will enhance structure prediction algorithms and provide a strategy for engineering stability in TIM barrel proteins.

At the core of amyloid fibrils is the cross-β spine, a long tape of β-sheets formed by the constituent proteins. Recent high-resolution x-ray studies show that the unit of this filamentous structure is a β-sheet bilayer with side chains within the bilayer forming a tightly interdigitating “steric zipper” interface. However, for a given peptide, different bilayer patterns are possible, and no quantitative explanation exists regarding which pattern is selected or under what condition there can be more than one pattern observed, exhibiting molecular polymorphism. We address the structural selection mechanism by performing molecular dynamics simulations to calculate the free energy of incorporating a peptide monomer into a β-sheet bilayer. We test filaments formed by several types of peptides including GNNQQNY, NNQQ, VEALYL, KLVFFAE and STVIIE, and find that the patterns with the lowest binding free energy correspond to available atomistic structures with high accuracy. Molecular polymorphism, as exhibited by NNQQ, is likely because there are more than one most stable structures whose binding free energies differ by less than the thermal energy. Detailed analysis of individual energy terms reveals that these short peptides are not strained nor do they lose much conformational entropy upon incorporating into a β-sheet bilayer. The selection of a bilayer pattern is determined mainly by the van der Waals and hydrophobic forces as a quantitative measure of shape complementarity among side chains between the β-sheets. The requirement for self-complementary steric zipper formation supports that amyloid fibrils form more easily among similar or same sequences, and it also makes parallel β-sheets generally preferred over anti-parallel ones. But the presence of charged side chains appears to kinetically drive anti-parallel β-sheets to form at early stages of assembly, after which the bilayer formation is likely driven by energetics., Author Summary Accumulation of amyloid fibrils is a salient feature of various protein misfolding diseases. Recent advances in precision experiments have begun to reveal their atomistic structures. Quantitative elucidation of how the observed structures are selected over other possible filament patterns would provide much insight into the formation and properties of amyloid fibrils. Using computer simulations and structural modeling, we demonstrate that the most stable filament pattern corresponds to the experimentally observed structure, and molecular polymorphism, selection of two or more patterns, is possible when there are more than one most stable structures. Ability to predict the structure allows for more detailed analysis, so that, for example, we can identify the most important residue for stabilizing the structure that could be therapeutically targeted. Our analysis will be useful for comparing different amyloid structures formed by the same protein or when delineating roles of different intermolecular forces in filament formation.

Background: KaiC, a central clock protein in cyanobacteria, undergoes circadian oscillations between hypophosphorylated and hyperphosphorylated forms in vivo and in vitro. Structural analyses of KaiC crystals have identified threonine and serine residues in KaiC at three residues (T426, S431, and T432) as potential sites at which KaiC is phosphorylated; mutation of any of these three sites to alanine abolishes rhythmicity, revealing an essential clock role for each residue separately and for KaiC phosphorylation in general. Mass spectrometry studies confirmed that the S431 and T432 residues are key phosphorylation sites, however, the role of the threonine residue at position 426 was not clear from the mass spectrometry measurements. Methodology and Principal Findings: Mutational approaches and biochemical analyses of KaiC support a key role for T426 in control of the KaiC phosphorylation status in vivo and in vitro and demonstrates that alternative amino acids at residue 426 dramatically affect KaiC’s properties in vivo and in vitro, especially genetic dominance/recessive relationships, KaiC dephosphorylation, and the formation of complexes of KaiC with KaiA and KaiB. These mutations alter key circadian properties, including period, amplitude, robustness, and temperature compensation. Crystallographic analyses indicate that the T426 site is phosphorylatible under some conditions, and in vitro phosphorylation assays of KaiC demonstrate labile phosphorylation of KaiC when the primary S431 and T432 sites are blocked. Conclusions and Significance: T426 is a crucial site that regulates KaiC phosphorylation status in vivo and in vitro and these studies underscore the importance of KaiC phosphorylation status in the essential cyanobacterial circadian functions. The regulatory roles of these phosphorylation sites–including T426–within KaiC enhance our understanding of the molecular mechanism underlying circadian rhythm generation in cyanobacteria.

BACKGROUND:Kinase-inducible domain (KID) as transcriptional activator can stimulate target gene expression in signal transduction by associating with KID interacting domain (KIX). NMR spectra suggest that apo-KID is an unstructured protein. After post-translational modification by phosphorylation, KID undergoes a transition from disordered to well folded protein upon binding to KIX. However, the mechanism of folding coupled to binding is poorly understood. METHODOLOGY:To get an insight into the mechanism, we have performed ten trajectories of explicit-solvent molecular dynamics (MD) for both bound and apo phosphorylated KID (pKID). Ten MD simulations are sufficient to capture the average properties in the protein folding and unfolding. CONCLUSIONS:Room-temperature MD simulations suggest that pKID becomes more rigid and stable upon the KIX-binding. Kinetic analysis of high-temperature MD simulations shows that bound pKID and apo-pKID unfold via a three-state and a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound pKID folds in the order of KIX access, initiation of pKID tertiary folding, folding of helix alpha(B), folding of helix alpha(A), completion of pKID tertiary folding, and finalization of pKID-KIX binding. Our data show that the folding pathways of apo-pKID are different from the bound state: the foldings of helices alpha(A) and alpha(B) are swapped. Here we also show that Asn139, Asp140 and Leu141 with large Phi-values are key residues in the folding of bound pKID. Our results are in good agreement with NMR experimental observations and provide significant insight into the general mechanisms of binding induced protein folding and other conformational adjustment in post-translational modification.

Protein loops, the flexible short segments connecting two stable secondary structural units in proteins, play a critical role in protein structure and function. Constructing chemically sensible conformations of protein loops that seamlessly bridge the gap between the anchor points without introducing any steric collisions remains an open challenge. A variety of algorithms have been developed to tackle the loop closure problem, ranging from inverse kinematics to knowledge-based approaches that utilize pre-existing fragments extracted from known protein structures. However, many of these approaches focus on the generation of conformations that mainly satisfy the fixed end point condition, leaving the steric constraints to be resolved in subsequent post-processing steps. In the present work, we describe a simple solution that simultaneously satisfies not only the end point and steric conditions, but also chirality and planarity constraints. Starting from random initial atomic coordinates, each individual conformation is generated independently by using a simple alternating scheme of pairwise distance adjustments of randomly chosen atoms, followed by fast geometric matching of the conformationally rigid components of the constituent amino acids. The method is conceptually simple, numerically stable and computationally efficient. Very importantly, additional constraints, such as those derived from NMR experiments, hydrogen bonds or salt bridges, can be incorporated into the algorithm in a straightforward and inexpensive way, making the method ideal for solving more complex multi-loop problems. The remarkable performance and robustness of the algorithm are demonstrated on a set of protein loops of length 4, 8, and 12 that have been used in previous studies., Author Summary Protein loops play an important role in protein function, such as ligand binding, recognition, and allosteric regulation. However, due to their flexibility, it is notoriously difficult to determine their 3D structures using traditional experimental techniques. As a result, one can often find protein structures with missing loops in the Protein Data Bank. Their sequence variability also presents a particular challenge for homology modeling methods, which can only yield good overall structures given sufficient sequence identity and good experimental reference structures. Despite extensive research, the construction of protein loop 3D structures remains an open problem, since a sensible conformation should seamlessly bridge the anchor points without introducing steric clashes within the loop itself or between the loop and its surroundings environment. Here, we present a conceptually simple, mathematically straightforward, numerically robust and computationally efficient approach for building protein loop conformations that simultaneously satisfy end-point, steric, planar and chiral constraints. More importantly, additional constraints derived from experimental sources can be incorporated in a straightforward manner, allowing the processing of more complex structures involving multiple interlocking loops.

Knowledge of the Free Energy Landscape topology is the essential key to understanding many biochemical processes. The determination of the conformers of a protein and their basins of attraction takes a central role for studying molecular isomerization reactions. In this work, we present a novel framework to unveil the features of a Free Energy Landscape answering questions such as how many meta-stable conformers there are, what the hierarchical relationship among them is, or what the structure and kinetics of the transition paths are. Exploring the landscape by molecular dynamics simulations, the microscopic data of the trajectory are encoded into a Conformational Markov Network. The structure of this graph reveals the regions of the conformational space corresponding to the basins of attraction. In addition, handling the Conformational Markov Network, relevant kinetic magnitudes as dwell times and rate constants, or hierarchical relationships among basins, completes the global picture of the landscape. We show the power of the analysis studying a toy model of a funnel-like potential and computing efficiently the conformers of a short peptide, dialanine, paving the way to a systematic study of the Free Energy Landscape in large peptides., Author Summary A complete description of complex polymers, such as proteins, includes information about their structure and their dynamics. In particular it is of utmost importance to answer the following questions: What are the structural conformations possible? Is there any relevant hierarchy among these conformers? What are the transition paths between them? These and other questions can be addressed by analyzing in an efficient way the Free Energy Landscape of the system. With this knowledge, several problems about biomolecular reactions (such as enzymatic activity, protein folding, protein deposition diseases, etc.) can be tackled. In this article we show how to efficiently describe the Free Energy Landscape for small and large peptides. By mapping the trajectories of molecular dynamics simulations into a graph (the Conformational Markov Network) and unveiling its structural organization, we obtain a coarse grained description of the protein dynamics across the Free Energy Landscape in terms of the relevant kinetic magnitudes of the system. Therefore, we show the way to bridge the gap between the microscopic dynamics and the macroscopic kinetics by means of a mesoscopic description of the associated Conformational Markov Network. Along this path the compromise between the physical nature of the process and the magnitudes that characterize the network is carefully kept to assure the reliability of the results shown.

Background Here we continue our efforts to use methods developed in the folding mechanism community to both better understand and improve structure prediction. Our previous work demonstrated that Rosetta's coarse-grained potentials may actually impede accurate structure prediction at full-atom resolution. Based on this work we postulated that it may be time to work completely at full-atom resolution but that doing so may require more careful attention to the kinetics of convergence. Methodology/Principal Findings To explore the possibility of working entirely at full-atom resolution, we apply enhanced sampling algorithms and the free energy theory developed in the folding mechanism community to full-atom protein structure prediction with the prominent Rosetta package. We find that Rosetta's full-atom scoring function is indeed able to recognize diverse protein native states and that there is a strong correlation between score and Cα RMSD to the native state. However, we also show that there is a huge entropic barrier to folding under this potential and the kinetics of folding are extremely slow. We then exploit this new understanding to suggest ways to improve structure prediction. Conclusions/Significance Based on this work we hypothesize that structure prediction may be improved by taking a more physical approach, i.e. considering the nature of the model thermodynamics and kinetics which result from structure prediction simulations.

Since the outbreak of severe acute respiratory syndrome (SARS) in 2003, the three-dimensional structures of several of the replicase/transcriptase components of SARS coronavirus (SARS-CoV), the non-structural proteins (Nsps), have been determined. However, within the large Nsp3 (1922 amino-acid residues), the structure and function of the so-called SARS-unique domain (SUD) have remained elusive. SUD occurs only in SARS-CoV and the highly related viruses found in certain bats, but is absent from all other coronaviruses. Therefore, it has been speculated that it may be involved in the extreme pathogenicity of SARS-CoV, compared to other coronaviruses, most of which cause only mild infections in humans. In order to help elucidate the function of the SUD, we have determined crystal structures of fragment 389–652 (“SUDcore”) of Nsp3, which comprises 264 of the 338 residues of the domain. Both the monoclinic and triclinic crystal forms (2.2 and 2.8 Å resolution, respectively) revealed that SUDcore forms a homodimer. Each monomer consists of two subdomains, SUD-N and SUD-M, with a macrodomain fold similar to the SARS-CoV X-domain. However, in contrast to the latter, SUD fails to bind ADP-ribose, as determined by zone-interference gel electrophoresis. Instead, the entire SUDcore as well as its individual subdomains interact with oligonucleotides known to form G-quadruplexes. This includes oligodeoxy- as well as oligoribonucleotides. Mutations of selected lysine residues on the surface of the SUD-N subdomain lead to reduction of G-quadruplex binding, whereas mutations in the SUD-M subdomain abolish it. As there is no evidence for Nsp3 entering the nucleus of the host cell, the SARS-CoV genomic RNA or host-cell mRNA containing long G-stretches may be targets of SUD. The SARS-CoV genome is devoid of G-stretches longer than 5–6 nucleotides, but more extended G-stretches are found in the 3′-nontranslated regions of mRNAs coding for certain host-cell proteins involved in apoptosis or signal transduction, and have been shown to bind to SUD in vitro. Therefore, SUD may be involved in controlling the host cell's response to the viral infection. Possible interference with poly(ADP-ribose) polymerase-like domains is also discussed., Author Summary The genome of the SARS coronavirus codes for 16 non-structural proteins that are involved in replicating this huge RNA (approximately 29 kilobases). The roles of many of these in replication (and/or transcription) are unknown. We attempt to derive conclusions concerning the possible functions of these proteins from their three-dimensional structures, which we determine by X-ray crystallography. Non-structural protein 3 contains at least seven different functional modules within its 1922-amino-acid polypeptide chain. One of these is the so-called SARS-unique domain, a stretch of about 338 residues that is completely absent from any other coronavirus. It may thus be responsible for the extraordinarily high pathogenicity of the SARS coronavirus, compared to other viruses of this family. We describe here the three-dimensional structure of the SARS-unique domain and show that it consists of two modules with a known fold, the so-called macrodomain. Furthermore, we demonstrate that these domains bind unusual nucleic-acid structures formed by consecutive guanosine nucleotides, where four strands of nucleic acid are forming a superhelix (so-called G-quadruplexes). SUD may be involved in binding to viral or host-cell RNA bearing this peculiar structure and thereby regulate viral replication or fight the immune response of the infected host cell.

Horses are unusual in producing protein-rich sweat for thermoregulation, a major component of which is latherin, a highly surface-active, non-glycosylated protein. The amino acid sequence of latherin, determined from cDNA analysis, is highly conserved across four geographically dispersed equid species (horse, zebra, onager, ass), and is similar to a family of proteins only found previously in the oral cavity and associated tissues of mammals. Latherin produces a significant reduction in water surface tension at low concentrations (< or = 1 mg ml(-1)), and therefore probably acts as a wetting agent to facilitate evaporative cooling through a waterproofed pelt. Neutron reflection experiments indicate that this detergent-like activity is associated with the formation of a dense protein layer, about 10 A thick, at the air-water interface. However, biophysical characterization (circular dichroism, differential scanning calorimetry) in solution shows that latherin behaves like a typical globular protein, although with unusual intrinsic fluorescence characteristics, suggesting that significant conformational change or unfolding of the protein is required for assembly of the air-water interfacial layer. RT-PCR screening revealed latherin transcripts in horse skin and salivary gland but in no other tissues. Recombinant latherin produced in bacteria was also found to be the target of IgE antibody from horse-allergic subjects. Equids therefore may have adapted an oral/salivary mucosal protein for two purposes peculiar to their lifestyle, namely their need for rapid and efficient heat dissipation and their specialisation for masticating and processing large quantities of dry food material.

The influence of lipid molecules on the aggregation of a highly amyloidogenic segment of human islet amyloid polypeptide, hIAPP20–29, and the corresponding sequence from rat has been studied by all-atom replica exchange molecular dynamics (REMD) simulations with explicit solvent model. hIAPP20–29 fragments aggregate into partially ordered β-sheet oligomers and then undergo large conformational reorganization and convert into parallel/antiparallel β-sheet oligomers in mixed in-register and out-of-register patterns. The hydrophobic interaction between lipid tails and residues at positions 23–25 is found to stabilize the ordered β-sheet structure, indicating a catalysis role of lipid molecules in hIAPP20–29 self-assembly. The rat IAPP variants with three proline residues maintain unstructured micelle-like oligomers, which is consistent with non-amyloidogenic behavior observed in experimental studies. Our study provides the atomic resolution descriptions of the catalytic function of lipid molecules on the aggregation of IAPP peptides., Author Summary People diagnosed with diabetes have increased from 30 million to 246 million over the last two decades. One hallmark of type 2 diabetes is the formation of amyloid in the pancreatic islet, which is composed of human islet amyloid polypeptide (90%) and lipid molecules (10%). In the long-lasting endeavors against the disease, it is important to understand, at the atomic level, the interaction between peptide aggregation and lipid molecules. In this study, we use molecular dynamics simulations to explore the influence of lipid molecules on the self-assembly process of toxic peptide segments. Moreover, a negative control simulation, employing the non-amyloidogenic rodent sequence, is also performed to evaluate the robustness of the simulation protocol. Our study provides a generic picture of the catalytic role of lipid molecules in the process of amyloidogenesis.

Genetic background exerts a strong modulatory effect on the toxicity of aggregation-prone proteins in conformational diseases. In addition to influencing the misfolding and aggregation behavior of the mutant proteins, polymorphisms in putative modifier genes may affect the molecular processes leading to the disease phenotype. Mutations in SOD1 in a subset of familial amyotrophic lateral sclerosis (ALS) cases confer dominant but clinically variable toxicity, thought to be mediated by misfolding and aggregation of mutant SOD1 protein. While the mechanism of toxicity remains unknown, both the nature of the SOD1 mutation and the genetic background in which it is expressed appear important. To address this, we established a Caenorhabditis elegans model to systematically examine the aggregation behavior and genetic interactions of mutant forms of SOD1. Expression of three structurally distinct SOD1 mutants in C. elegans muscle cells resulted in the appearance of heterogeneous populations of aggregates and was associated with only mild cellular dysfunction. However, introduction of destabilizing temperature-sensitive mutations into the genetic background strongly enhanced the toxicity of SOD1 mutants, resulting in exposure of several deleterious phenotypes at permissive conditions in a manner dependent on the specific SOD1 mutation. The nature of the observed phenotype was dependent on the temperature-sensitive mutation present, while its penetrance reflected the specific combination of temperature-sensitive and SOD1 mutations. Thus, the specific toxic phenotypes of conformational disease may not be simply due to misfolding/aggregation toxicity of the causative mutant proteins, but may be defined by their genetic interactions with cellular pathways harboring mildly destabilizing missense alleles., Author Summary Correct folding and stability are essential for protein function. In cells, a network of molecular chaperones and degradative enzymes facilitate folding, prevent aggregation and ensure degradation of the misfolded proteins, thus maintaining protein homeostasis. In many diseases, including Amyotrophic Lateral Sclerosis (ALS), expression of a single mutant protein that misfolds and aggregates causes cellular toxicity that is strongly dependent on the genetic background. To address the influence of genetic background on the toxicity of aggregation-prone proteins, we established a C. elegans model of misfolding and aggregation of several distinct ALS-related mutants of superoxide dismutase 1 (SOD1). In one wild type genetic background (N2), these proteins exhibited only mild cellular toxicity despite strong, mutant-specific aggregation phenotypes. However, when SOD1 mutants were expressed in the background of mildly destabilized protein polymorphisms, their toxicity was enhanced and a number of distinct phenotypes were exposed. These synthetic phenotypes reflected the loss-of-function of the destabilized polymorphic proteins. Furthermore, the degree to which each of these phenotypes was exposed depended on the nature of the SOD1 mutation. These data suggest that the presence of mildly destabilizing polymorphisms in the genetic background may modulate and direct the specific toxic phenotypes in protein aggregation diseases.

We perform a large-scale study of intrinsically disordered regions in proteins and protein complexes using a non-redundant set of hundreds of different protein complexes. In accordance with the conventional view that folding and binding are coupled, in many of our cases the disorder-to-order transition occurs upon complex formation and can be localized to binding interfaces. Moreover, analysis of disorder in protein complexes depicts a significant fraction of intrinsically disordered regions, with up to one third of all residues being disordered. We find that the disorder in homodimers, especially in symmetrical homodimers, is significantly higher than in heterodimers and offer an explanation for this interesting phenomenon. We argue that the mechanisms of regulation of binding specificity through disordered regions in complexes can be as common as for unbound monomeric proteins. The fascinating diversity of roles of disordered regions in various biological processes and protein oligomeric forms shown in our study may be a subject of future endeavors in this area., Author Summary Traditionally, protein structure is believed to determine function. Recently, it was observed that many proteins contain regions without well-defined structure (intrinsically disordered regions), including a large fraction of eukaryotic proteins. Intrinsic disorder has been associated with particular functions including cell regulation; signaling; and protein, DNA, and ligand binding. Many proteins are intrinsically disordered in native form and fold upon binding, following the conventional paradigm. Accordingly, disorder in a protein may facilitate binding to multiple partners. However, in some cases disorder has also been found in the bound state. To gain clearer insight into the functional importance of disorder regions in protein complexes, we perform a large-scale analysis of disorder using protein structures in complex and in unbound forms. We show that disorder in protein complexes is rather common and pinpoint changes that occur upon protein binding at interaction interfaces. By illustrating a variety of functional roles for disorder in specific proteins, we emphasize the versatility and importance of this phenomenon.

Peptide inhibitors derived from HIV-gp41 envelope protein play a pivotal role in deciphering the molecular mechanism of HIV-cell fusion. According to accepted models, N-heptad repeat (NHR) peptides can bind two targets in an intermediate fusion conformation, thereby inhibiting progression of the fusion process. In both cases the orientation towards the endogenous intermediate conformation should be important. To test this, we anchored NHR to the cell membrane by conjugating fatty acids with increasing lengths to the N- or C-terminus of N36, as well as to two known N36 mutants; one that cannot bind C-heptad repeat (CHR) but can bind NHR (N36 MUTe,g), and the second cannot bind to either NHR or CHR (N36 MUTa,d). Importantly, the IC50 increased up to 100-fold in a lipopeptide-dependent manner. However, no preferred directionality was observed for the wild type derived lipopeptides, suggesting a planar orientation of the peptides as well as the endogenous NHR region on the cell membrane. Furthermore, based on: (i) specialized analysis of the inhibition curves, (ii) the finding that N36 conjugates reside more on the target cells that occupy the receptors, and (iii) the finding that N36 MUTe,g acts as a monomer both in its soluble form and when anchored to the cell membrane, we suggest that anchoring N36 to the cell changes the inhibitory mode from a trimer which can target both the endogenous NHR and CHR regions, to mainly monomeric lipopetides that target primarily the internal NHR. Besides shedding light on the mode of action of HIV-cell fusion, the similarity between functional regions in the envelopes of other viruses suggests a new approach for developing potent HIV-1 inhibitors., Author Summary Acquired immunodeficiency syndrome (AIDS) is a major global health problem, and its causative agent, human immunodeficiency virus (HIV), is extensively studied. To start an infectious cycle HIV must fuse its membrane with that of its host cell. A specific protein on the virus surface facilitates this process by undergoing major conformational changes. Several virus-cell fusion inhibitors target transiently exposed regions during the conformational changes, thereby preventing progression of the fusion process. Here, we focused on a specific fusion inhibitor peptide having two distinct binding sites and modes of inhibitions. By simple chemical modifications we demonstrate a shift between these two modes of inhibition. Most importantly, we reveal novel details regarding the conformational changes during the fusion process. Furthermore, the chemical modifications extremely enhanced the fusion inhibitory potency of the peptide. Lastly, since the fusion process of HIV shares common features with diverse biological processes, our results might contribute to their research and therapeutic efforts as well.

It has long been proposed that much of the information encoding how a protein folds is contained locally in the peptide chain. Here we present a large-scale simulation study designed to examine the extent to which conformations of peptide fragments in water predict native conformations in proteins. We perform replica exchange molecular dynamics (REMD) simulations of 872 8-mer, 12-mer, and 16-mer peptide fragments from 13 proteins using the AMBER 96 force field and the OBC implicit solvent model. To analyze the simulations, we compute various contact-based metrics, such as contact probability, and then apply Bayesian classifier methods to infer which metastable contacts are likely to be native vs. non-native. We find that a simple measure, the observed contact probability, is largely more predictive of a peptide's native structure in the protein than combinations of metrics or multi-body components. Our best classification model is a logistic regression model that can achieve up to 63% correct classifications for 8-mers, 71% for 12-mers, and 76% for 16-mers. We validate these results on fragments of a protein outside our training set. We conclude that local structure provides information to solve some but not all of the conformational search problem. These results help improve our understanding of folding mechanisms, and have implications for improving physics-based conformational sampling and structure prediction using all-atom molecular simulations., Author Summary Proteins must fold to unique native structures in order to perform their functions. To do this, proteins must solve a complicated conformational search problem, the details of which remain difficult to study experimentally. Predicting folding pathways and the mechanisms by which proteins fold is thus central to understanding how proteins work. One longstanding question is the extent to which proteins solve the search problem locally, by folding into sub-structures that are dictated primarily by local sequence. Here, we address this question by conducting a large-scale molecular dynamics simulation study of protein fragments in water. The simulation data was then used to optimize a statistical model that predicted native and non-native contacts. The performance of the resulting model suggests that local structuring provides some but not all of the information to solve the folding problem, and that molecular dynamics simulation of fragments can be useful for protein structure prediction and design.

Amyloids and prion proteins are clinically and biologically important β-structures, whose supersecondary structures are difficult to determine by standard experimental or computational means. In addition, significant conformational heterogeneity is known or suspected to exist in many amyloid fibrils. Recent work has indicated the utility of pairwise probabilistic statistics in β-structure prediction. We develop here a new strategy for β-structure prediction, emphasizing the determination of β-strands and pairs of β-strands as fundamental units of β-structure. Our program, BETASCAN, calculates likelihood scores for potential β-strands and strand-pairs based on correlations observed in parallel β-sheets. The program then determines the strands and pairs with the greatest local likelihood for all of the sequence's potential β-structures. BETASCAN suggests multiple alternate folding patterns and assigns relative a priori probabilities based solely on amino acid sequence, probability tables, and pre-chosen parameters. The algorithm compares favorably with the results of previous algorithms (BETAPRO, PASTA, SALSA, TANGO, and Zyggregator) in β-structure prediction and amyloid propensity prediction. Accurate prediction is demonstrated for experimentally determined amyloid β-structures, for a set of known β-aggregates, and for the parallel β-strands of β-helices, amyloid-like globular proteins. BETASCAN is able both to detect β-strands with higher sensitivity and to detect the edges of β-strands in a richly β-like sequence. For two proteins (Aβ and Het-s), there exist multiple sets of experimental data implying contradictory structures; BETASCAN is able to detect each competing structure as a potential structure variant. The ability to correlate multiple alternate β-structures to experiment opens the possibility of computational investigation of prion strains and structural heterogeneity of amyloid. BETASCAN is publicly accessible on the Web at http://betascan.csail.mit.edu., Author Summary Amyloid is a highly ordered form of protein aggregation that a wide variety of proteins can form. While the earliest discovered amyloids were associated with systemic and neurodegenerative diseases, recent findings indicate amyloids may have myriad roles and functions ranging from learning and memory, to yeast epigenetics, to biofilm and melanin production. In this study, we expand the range and flexibility of our ability to understand how amyloid properties arise from their polypeptide sequence. By taking advantage of the intrinsic properties of a characteristic amyloid structure—parallel β-strands—and data from available protein structures, we construct and test an algorithm to predict the probability that particular portions of a protein will form amyloid. Our method has the advantage of more accurate detection of the edges of such zones, as well as the ability to consider and evaluate the likelihood of multiple folding patterns.

We consider a recently suggested "equation of state" for natively folded proteins, and verify its validity for a set of about 5800 proteins. The equation is based on a fractal viewpoint of proteins, on a generalization of the Landau-Peierls instability, and on a marginal stability criterion. The latter allows for coexistence of stability and flexibility of proteins, which is required for their proper function. The equation of state relates the protein fractal dimension d(f), its spectral dimension d(s), and the number of amino acids N. Using structural data from the protein data bank (PDB) and the Gaussian network model (GNM), we compute d(f) and d(s) for the entire set and demonstrate that the equation of state is well obeyed. Addressing the fractal properties and making use of the equation of state may help to engineer biologically inspired catalysts.