Search

Amyloid-beta (Aβ) toxic oligomers are critical early players in the molecular pathology of Alzheimer's disease (AD). We have developed a Soluble Oligomer Binding Assay (SOBA-AD) for detection of these Aβ oligomers that contain α-sheet secondary structure that discriminates plasma samples from patients on the AD continuum from non-AD controls. We tested 265 plasma samples from two independent cohorts to investigate the performance of SOBA-AD. Testing was performed at two different sites, with different personnel, reagents, and instrumentation. Across two cohorts, SOBA-AD discriminated AD patients from cognitively unimpaired (CU) subjects with 100% sensitivity, > 95% specificity, and > 98% area under the curve (AUC) (95% CI 0.95–1.00). A SOBA-AD positive readout, reflecting α-sheet toxic oligomer burden, was found in AD patients, and not in controls, providing separation of the two populations, aside from 5 SOBA-AD positive controls. Based on an earlier SOBA-AD study, the Aβ oligomers detected in these CU subjects may represent preclinical cases of AD. The results presented here support the value of SOBA-AD as a promising blood-based tool for the detection and confirmation of AD. [ABSTRACT FROM AUTHOR]

Colocalization of microbial pathogens and the β-amyloid peptide (Aβ) in the brain of Alzheimer’s disease (AD) patients suggests that microbial infection may play a role in sporadic AD. Aβ exhibits antimicrobial activity against numerous pathogens, supporting a potential role for Aβ in the innate immune response. While mammalian amyloid is associated with disease, many bacteria form amyloid fibrils to fortify the biofilm that protects the cells from the surrounding environment. In the microbial AD hypothesis, Aβ aggregates in response to infection to combat the pathogen. We hypothesize that this occurs through toxic Aβ oligomers that contain α-sheet structure and form prior to fibrillization. De novo designed α-sheet peptides specifically bind to the α-sheet structure present in the oligomers of both bacterial and mammalian amyloidogenic proteins to neutralize toxicity and inhibit aggregation. Here, we measure the effect of E. coli on Aβ, including upregulation, aggregation, and toxicity. Additionally, we determined the effect of Aβ structure on E. coli amyloid fibrils, or curli comprised of the CsgA protein, and biofilm formation. We found that curli formation by E. coli increased Aβ oligomer production, and Aβ oligomers inhibited curli biogenesis and reduced biofilm cell density. Further, curli and biofilm inhibition by Aβ oligomers increased E. coli susceptibility to gentamicin. Toxic oligomers of Aβ and CsgA interact via α-sheet interactions, neutralizing their toxicity. These results suggest that exposure to toxic oligomers formed by microbial pathogens triggers Aβ oligomer upregulation and aggregation to combat infection via selective interactions between α-sheet oligomers to neutralize toxicity of both species with subsequent inhibition of fibrillization. [ABSTRACT FROM AUTHOR]

Type 2 diabetes (T2D) results from insulin secretory dysfunction arising in part from the loss of pancreatic islet β‐cells. Several factors contribute to β‐cell loss, including islet amyloid formation, which is observed in over 90% of individuals with T2D. The amyloid is comprised of human islet amyloid polypeptide (hIAPP). Here we provide evidence that early in aggregation, hIAPP forms toxic oligomers prior to formation of amyloid fibrils. The toxic oligomers contain α‐sheet secondary structure, a nonstandard secondary structure associated with toxic oligomers in other amyloid diseases. De novo, synthetic α‐sheet compounds designed to be nontoxic and complementary to the α‐sheet structure in the toxic oligomers inhibit hIAPP aggregation and neutralize oligomer‐mediated cytotoxicity in cell‐based assays. In vivo administration of an α‐sheet design to mice for 4 weeks revealed no evidence of toxicity nor did it elicit an immune response. Furthermore, the α‐sheet designs reduced endogenous islet amyloid formation and mitigation of amyloid‐associated β‐cell loss in cultured islets isolated from an hIAPP transgenic mouse model of islet amyloidosis. Characterization of the involvement of α‐sheet in early aggregation of hIAPP and oligomer toxicity contributes to elucidation of the molecular mechanisms underlying amyloid‐associated β‐cell loss. [ABSTRACT FROM AUTHOR]

Uropathogenic Escherichia coli account for the largest proportion of nosocomial infections in the United States. Nosocomial infections are a major source of increased costs and treatment complications. Many infections are biofilm associated, rendering antibiotic treatments ineffective or cause additional complications (e.g., microbiome depletion). This work presents a potentially complementary non-antibiotic strategy to fight nosocomial infections by inhibiting the formation of amyloid fibrils, a proteinaceous structural reinforcement known as curli in E. coli biofilms. Despite extensive characterization of the fibrils themselves and their associated secretion system, mechanistic details of curli assembly in vivo remain unclear. We hypothesized that, like other amyloid fibrils, curli polymerization involves a unique secondary structure termed "α-sheet". Biophysical studies herein confirmed the presence of α-sheet structure in prefibrillar species of CsgA, the major component of curli, as it aggregated. Binding of synthetic α-sheet peptides to the soluble α-sheet prefibrillar species inhibited CsgA aggregation in vitro and suppressed amyloid fibril formation in biofilms. Application of synthetic α-sheet peptides also enhanced antibiotic susceptibility and dispersed biofilm-resident bacteria for improved uptake by phagocytic cells. The ability of synthetic α-sheet peptides to reduce biofilm formation, improve antibiotic susceptibility, and enhance clearance by macrophages has broad implications for combating biofilm-associated infections. [ABSTRACT FROM AUTHOR]

The formation of toxic Amyloid β-peptide (Aβ) oligomers is one of the earliest events in the molecular pathology of Alzheimer's Disease (AD). These oligomers lead to a variety of downstream effects, including impaired neuronal signaling, neuroinflammation, tau phosphorylation, and neurodegeneration, and it is estimated that these events begin 10 to 20 y before the presentation of symptoms. Toxic Aβ oligomers contain a nonstandard protein structure, termed α-sheet, and designed α-sheet peptides target this main-chain structure in toxic oligomers independent of sequence. Here we show that a designed α-sheet peptide inhibits the deleterious effects on neuronal signaling and also serves as a capture agent in our soluble oligomer binding assay (SOBA). Pre-incubated synthetic α-sheet-containing Aβ oligomers produce strong SOBA signals, while monomeric and β-sheet protofibrillar Aβ do not. α-sheet containing oligomers were also present in cerebrospinal fluid (CSF) from an AD patient versus a noncognitively impaired control. For the detection of toxic oligomers in plasma, we developed a plate coating to increase the density of the capture peptide. The proof of concept was achieved by testing 379 banked human plasma samples. SOBA detected Aβ oligomers in patients on the AD continuum, including controls who later progressed to mild cognitive impairment. In addition, SOBA discriminated AD from other forms of dementia, yielding sensitivity and specificity of 99% relative to clinical and neuropathological diagnoses. To explore the broader potential of SOBA, we adapted the assay for α-synuclein oligomers and confirmed their presence in CSF from patients with Parkinson's disease and Lewy body dementia. [ABSTRACT FROM AUTHOR]

Alzheimer's Disease (AD) is a neurodegenerative disorder that is characterized clinically by progressive cognitive decline and pathologically by the β-sheet rich fibril plaque deposition of the amyloid-β (Aβ) peptide in the brain. While plaques are a hallmark of AD, plaque burden is not correlated with cognitive impairment. Instead, Aβ oligomers formed during the aggregation process represent the main agents of neurotoxicity, which occurs 10–20 years before patients begin to show symptoms. These oligomers are dynamic in nature and represented by a heterogeneous distribution of aggregates ranging from low- to high-molecular weight, some of which are toxic while others are not. A major difficulty in determining the pathological mechanism(s) of Aβ, developing reliable diagnostic markers for early-stage detection, as well as effective therapeutics for AD are the differentiation and characterization of oligomers formed throughout disease propagation based on their molecular features, effects on biological function, and relevance to disease propagation and pathology. Thus, it is critical to methodically identify the mechanisms of Aβ aggregation and toxicity, as well as describe the roles of different oligomers and aggregates in disease progression and molecular pathology. Here, we describe a variety of biophysical techniques used to isolate and characterize a range of Aβ oligomer populations, as well as discuss proposed mechanisms of toxicity and therapeutic interventions aimed at specific assemblies formed during the aggregation process. The approaches being used to map the misfolding and aggregation of Aβ are like what was done during the fundamental early studies, mapping protein folding pathways using combinations of biophysical techniques in concert with protein engineering. Such information is critical to the design and molecular engineering of future diagnostics and therapeutics for AD. [ABSTRACT FROM AUTHOR]

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long‐range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures. [ABSTRACT FROM AUTHOR]

Key points: Skeletal muscle relaxation has been primarily studied by assessing the kinetics of force decay. Little is known about the resultant dynamics of structural changes in myosin heads during relaxation.The naturally occurring nucleotide 2‐deoxy‐ATP (dATP) is a myosin activator that enhances cross‐bridge binding and kinetics.X‐ray diffraction data indicate that with elevated dATP, myosin heads were extended closer to actin in relaxed muscle and myosin heads return to an ordered, resting state after contraction more quickly.Molecular dynamics simulations of post‐powerstroke myosin suggest that dATP induces structural changes in myosin heads that increase the surface area of the actin‐binding regions promoting myosin interaction with actin, which could explain the observed delays in the onset of relaxation.This study of the dATP‐induced changes in myosin may be instructive for determining the structural changes desired for other potential myosin‐targeted molecular compounds to treat muscle diseases. Here we used time‐resolved small‐angle X‐ray diffraction coupled with force measurements to study the structural changes in FVB mouse skeletal muscle sarcomeres during relaxation after tetanus contraction. To estimate the rate of myosin deactivation, we followed the rate of the intensity recovery of the first‐order myosin layer line (MLL1) and restoration of the resting spacing of the third and sixth order of meridional reflection (SM3 and SM6) following tetanic contraction. A transgenic mouse model with elevated skeletal muscle 2‐deoxy‐ATP (dATP) was used to study how myosin activators may affect soleus muscle relaxation. X‐ray diffraction evidence indicates that with elevated dATP, myosin heads were extended closer to actin in resting muscle. Following contraction, there is a slight but significant delay in the decay of force relative to WT muscle while the return of myosin heads to an ordered resting state was initially slower, then became more rapid than in WT muscle. Molecular dynamics simulations of post‐powerstroke myosin suggest that dATP induces structural changes in myosin that increase the surface area of the actin‐binding regions, promoting myosin interaction with actin. With dATP, myosin heads may remain in an activated state near the thin filaments following relaxation, accounting for the delay in force decay and the initial delay in recovery of resting head configuration, and this could facilitate subsequent contractions. Key points: Skeletal muscle relaxation has been primarily studied by assessing the kinetics of force decay. Little is known about the resultant dynamics of structural changes in myosin heads during relaxation.The naturally occurring nucleotide 2‐deoxy‐ATP (dATP) is a myosin activator that enhances cross‐bridge binding and kinetics.X‐ray diffraction data indicate that with elevated dATP, myosin heads were extended closer to actin in relaxed muscle and myosin heads return to an ordered, resting state after contraction more quickly.Molecular dynamics simulations of post‐powerstroke myosin suggest that dATP induces structural changes in myosin heads that increase the surface area of the actin‐binding regions promoting myosin interaction with actin, which could explain the observed delays in the onset of relaxation.This study of the dATP‐induced changes in myosin may be instructive for determining the structural changes desired for other potential myosin‐targeted molecular compounds to treat muscle diseases. [ABSTRACT FROM AUTHOR]

The p53 protein is a commonly studied cancer target because of its role in tumor suppression. Unfortunately, it is susceptible to mutation‐associated loss of function; approximately 50% of cancers are associated with mutations to p53, the majority of which are located in the central DNA‐binding domain. Here, we report molecular dynamics simulations of wild‐type (WT) p53 and 20 different mutants, including a stabilized pseudo‐WT mutant. Our findings indicate that p53 mutants tend to exacerbate latent structural‐disruption tendencies, or vulnerabilities, already present in the WT protein, suggesting that it may be possible to develop cancer therapies by targeting a relatively small set of structural‐disruption motifs rather than a multitude of effects specific to each mutant. In addition, α‐sheet secondary structure formed in almost all of the proteins. α‐Sheet has been hypothesized and recently demonstrated to play a role in amyloidogenesis, and its presence in the reported p53 simulations coincides with the recent re‐consideration of cancer as an amyloid disease. [ABSTRACT FROM AUTHOR]

The Dynameomics project contains native state and unfolding simulations of 807 protein domains, where each domain is representative of a different metafold; these metafolds encompass ~97% of protein fold space. There is a long-standing question in structural biology as to whether proteins in the same fold family share the same folding/unfolding characteristics. Using molecular dynamics simulations from the Dynameomics project, we conducted a detailed study of protein unfolding/folding pathways for 5 protein domains from the immunoglobulin (Ig)-like β-sandwich metafold (the highest ranked metafold in our database). The domains have sequence similarities ranging from 4 to 15% and are all from different SCOP superfamilies, yet they share the same overall Ig-like topology. Despite having very different amino acid sequences, the dominant unfolding pathway is very similar for the 5 proteins, and the secondary structures that are peripheral to the aligned, shared core domain add variability to the unfolding pathway. Aligned residues in the core domain display consensus structure in the transition state primarily through conservation of hydrophobic positions. Commonalities in the obligate folding nucleus indicate that insights into the major events in the folding/unfolding of other domains from this metafold may be obtainable from unfolding simulations of a few representative proteins. [ABSTRACT FROM AUTHOR]

The naturally occurring nucleotide 2-deoxy-adenosine 5Œ-triphosphate (dATP) can be used by cardiac muscle as an alternative energy substrate for myosin chemomechanical activity. We and others have previously shown that dATP increases contractile force in normal hearts and models of depressed systolic function, but the structural basis of these effects has remained unresolved. In this work, we combine multiple techniques to provide structural and functional information at the angstrom-nanometer and millisecond time scales, demonstrating the ability to make both structural measurements and quantitative kinetic estimates of weak actin. myosin interactions that underpin sarcomere dynamics. Exploiting dATP as a molecular probe, we assess how small changes in myosin structure translate to electrostatic-based changes in sarcomere function to augment contractility in cardiac muscle. Through Brownian dynamics simulation and computational structural analysis, we found that deoxy-hydrolysis products [2-deoxy-adenosine 5Œ-diphosphate (dADP) and inorganic phosphate (Pi)] bound to prepowerstroke myosin induce an allosteric restructuring of the actin-binding surface on myosin to increase the rate of crossbridge formation. We then show experimentally that this predicted effect translates into increased electrostatic interactions between actin and cardiac myosin in vitro. Finally, using small-angle X-ray diffraction analysis of sarcomere structure, we demonstrate that the proposed increased electrostatic affinity of myosin for actin causes a disruption of the resting conformation of myosin motors, resulting in their repositioning toward the thin filament before activation. The dATP-mediated structural alterations in myosin reported here may provide insight into an improved criterion for the design or selection of small molecules to be developed as therapeutic agents to treat systolic dysfunction. [ABSTRACT FROM AUTHOR]

Alzheimer's disease (AD) is characterized by the deposition of β-sheet-rich, insoluble amyloid β-peptide (Aβ) plaques; however, plaque burden is not correlated with cognitive impairment in AD patients; instead, it is correlated with the presence of toxic soluble oligomers. Here, we show, by a variety of different techniques, that these Aβ oligomers adopt a nonstandard secondary structure, termed "α-sheet." These oligomers form in the lag phase of aggregation, when Aβ-associated cytotoxicity peaks, en route to forming nontoxic β-sheet fibrils. De novo-designed α-sheet peptides specifically and tightly bind the toxic oligomers over monomeric and fibrillar forms of Aβ, leading to inhibition of aggregation in vitro and neurotoxicity in neuroblastoma cells. Based on this specific binding, a soluble oligomer-binding assay (SOBA) was developed as an indirect probe of α-sheet content. Combined SOBA and toxicity experiments demonstrate a strong correlation between α-sheet content and toxicity. The designed α-sheet peptides are also active in vivo where they inhibit Aβ-induced paralysis in a transgenic Aβ Caenorhabditis elegans model and specifically target and clear soluble, toxic oligomers in a transgenic APPsw mouse model. The α-sheet hypothesis has profound implications for further understanding the mechanism behind AD pathogenesis. [ABSTRACT FROM AUTHOR]

Abstract: Metamorphic proteins are biomolecules prone to adopting alternative conformations. Because of this feature, they represent ideal systems to investigate the general rules allowing primary structure to dictate protein topology. A comparative molecular dynamics study was performed on the denatured states of two proteins, sharing nearly identical amino‐acid sequences (88 %) but different topologies, namely an all‐α‐helical bundle protein named GA88 and an α+β‐protein named GB88. The analysis allowed successful design of and experimental validation of a site‐directed mutant that promotes, at least in part, the switch in folding from GB88 to GA88. The mutated position, in which a glutamic acid was replaced by a glutamine, does not make any intramolecular interactions in the native state of GA88, such that its stabilization can be explained by considering the effects on the denatured state. The results represent a direct demonstration of the role of the denatured state in sculpting native structure. [ABSTRACT FROM AUTHOR]

Copyright of Angewandte Chemie is the property of Wiley-Blackwell and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Computational resources have contributed to the design and engineering of novel proteins by integrating genomic, structural and dynamic aspects of proteins. Non-canonical amino acids, such as d -amino acids, expand the available sequence space for designing and engineering proteins; however, the rotamer libraries for d -amino acids are usually constructed as the mirror images of l -amino acid rotamer libraries, an assumption that has not been tested. To this end, we have performed molecular dynamics (MD) simulations of model host–guest peptide systems containing d -amino acids. Our simulations systematically address the applicability of the mirror image convention as well as the effects of neighboring residue chirality. Rotamer libraries derived from these systems provide realistic rotamer distributions suitable for use in both rational and computational design workflows. Our simulations also address the impact of chirality on the intrinsic conformational preferences of amino acids, providing fundamental insights into the relationship between chirality and biomolecular dynamics. While d -amino acids are rare in naturally occurring proteins, they are used in designed proteins to stabilize a desired conformation, increase bioavailability or confer favorable biochemical and physical attributes. Here, we present d -amino acid rotamer libraries derived from MD simulations of alanine-based host–guest pentapeptides and show how certain residues can deviate from mirror image symmetry. Our simulations directly model d -amino acids as guest residues within the chiral l -Ala and d -Ala pentapeptide series to explicitly incorporate any contributions resulting from the chiralities of neighboring residues. [ABSTRACT FROM AUTHOR]

Prion diseases are associated with the misfolding of the prion protein (PrP) from its normal cellular form (Pr PC) to its infectious scrapie form (Pr PSc). Post-translational modifications in PrP in vivo can play an important role in modulating the process of misfolding. To gain more insight into the effects of post-translational modifications in PrP structure and dynamics and to test the hypothesis that such modifications can interact with the protein, we have performed molecular dynamics simulations of diglycosylated human Pr PC bound to a lipid bilayer via a glycophosphatidylinositol anchor. Multiple simulations were performed at three different pH ranges to explore pH effects on structure and dynamics. In contrast to simulations of protein-only Pr PC, no large effects were observed upon lowering the pH of the system. The protein tilted toward the membrane surface in all of the simulations and the putative Pr PSc oligomerization sites became inaccessible, thereby offering a possible protective mechanism against Pr PSc-induced misfolding of PrPC. [ABSTRACT FROM AUTHOR]

Myosin activation is a viable approach to treat systolic heart failure. We previously demonstrated that striated muscle myosin is a promiscuous ATPase that can use most nucleoside triphosphates as energy substrates for contraction. When 2-deoxy ATP (dATP) is used, it acts as a myosin activator, enhancing cross-bridge binding and cycling. In vivo, we have demonstrated that elevated dATP levels increase basal cardiac function and rescues function of infarcted rodent and pig hearts. Here we investigate the molecular mechanism underlying this physiological effect. We show with molecular dynamics simulations that the binding of dADP.Pi (dATP hydrolysis products) to myosin alters the structure and dynamics of the nucleotide binding pocket, myosin cleft conformation, and actin binding sites, which collectively yield a myosin conformation that we predict favors weak, electrostatic binding to actin. In vitro motility assays at high ionic strength were conducted to test this prediction and we found that dATP increased motility. These results highlight alterations to myosin that enhance cross-bridge formation and reveal a potential mechanism that may underlie dATP-induced improvements in cardiac function. [ABSTRACT FROM AUTHOR]

ABSTRACT Amyloids have long been associated with protein dysfunction and neurodegenerative diseases, but recent research has demonstrated that some organisms utilize the unique properties of the amyloid fold to create functional structures with important roles in biological processes. Additionally, new engineering approaches have taken advantage of amyloid structures for implementation in a wide variety of materials and devices. In this review, the role of amyloid in human disease is discussed and compared to the functional amyloids, which serve a largely structural purpose. We then consider the use of amyloid constructs in engineering applications, including their utility as building blocks for synthetic biology and molecular engineering. Biotechnol. Bioeng. 2017;114: 7-20. © 2016 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

The p53 tumor suppressor protein performs a critical role in stimulating apoptosis and cell cycle arrest in response to oncogenic stress. The function of p53 can be compromised by mutation, leading to increased risk of cancer; approximately 50% of cancers are associated with mutations in the p53 gene, the majority of which are in the core DNA-binding domain. The Y220C mutation of p53, for example, destabilizes the core domain by 4 kcal/mol, leading to rapid denaturation and aggregation. The associated loss of tumor suppressor functionality is associated with approximately 75 000 new cancer cases every year. Destabilized p53 mutants can be 'rescued' and their function restored; binding of a small molecule into a pocket on the surface of mutant p53 can stabilize its wild-type structure and restore its function. Here, we describe an in silico algorithm for identifying potential rescue pockets, including the algorithm's integration with the Dynameomics molecular dynamics data warehouse and the DIVE visual analytics engine. We discuss the results of the application of the method to the Y220C p53 mutant, entailing finding a putative rescue pocket through MD simulations followed by an in silico search for stabilizing ligands that dock into the putative rescue pocket. The top three compounds from this search were tested experimentally and one of them bound in the pocket, as shown by nuclear magnetic resonance, and weakly stabilized the mutant. [ABSTRACT FROM AUTHOR]

The conformational propensities of amino acids are an amalgamation of sequence effects, environmental effects and underlying intrinsic behavior. Many have attempted to investigate neighboring residue effects to aid in our understanding of protein folding and improve structure prediction efforts, especially with respect to difficult to characterize states, such as disordered or unfolded states. Host-guest peptide series are a useful tool in examining the propensities of the amino acids free from the surrounding protein structure. Here, we compare the distributions of the backbone dihedral angles (ϕ/ψ) of the 20 proteogenic amino acids in two different sequence contexts using the AAXAA and GGXGG host-guest pentapeptide series. We further examine their intrinsic behaviors across three environmental contexts: water at 298 K, water at 498 K, and 8 M urea at 298 K. The GGXGG systems provide the intrinsic amino acid propensities devoid of any conformational context. The alanine residues in the AAXAA series enforce backbone chirality, thereby providing a model of the intrinsic behavior of amino acids in a protein chain. Our results show modest differences in ϕ/ψ distributions due to the steric constraints of the Ala side chains, the magnitudes of which are dependent on the denaturing conditions. One of the strongest factors modulating φ/ψ distributions was the protonation of titratable side chains, and the largest differences observed were in the amino acid propensities for the rarely sampled αL region. [ABSTRACT FROM AUTHOR]

Protein function is intimately linked to protein structure and dynamics yet experimentally determined structures frequently omit regions within a protein due to indeterminate data, which is often due protein dynamics. We propose that atomistic molecular dynamics simulations provide a diverse sampling of biologically relevant structures for these missing segments (and beyond) to improve structural modeling and structure prediction. Here we make use of the Dynameomics data warehouse, which contains simulations of representatives of essentially all known protein folds. We developed novel computational methods to efficiently identify, rank and retrieve small peptide structures, or fragments, from this database. We also created a novel data model to analyze and compare large repositories of structural data, such as contained within the Protein Data Bank and the Dynameomics data warehouse. Our evaluation compares these structural repositories for improving loop predictions and analyzes the utility of our methods and models. Using a standard set of loop structures, containing 510 loops, 30 for each loop length from 4 to 20 residues, we find that the inclusion of Dynameomics structures in fragment-based methods improves the quality of the loop predictions without being dependent on sequence homology. Depending on loop length, ∼25-75% of the best predictions came from the Dynameomics set, resulting in lower main chain root-mean-square deviations for all fragment lengths using the combined fragment library. We also provide specific cases where Dynameomics fragments provide better predictions for NMR loop structures than fragments from crystal structures. Online access to these fragment libraries is available at . [ABSTRACT FROM AUTHOR]

d-amino acids are useful building blocks for de novo peptide design and they play a role in aging-related diseases associated with gradual protein racemization. For amino acids with achiral side chains, one should be able to presume that the conformational propensities of l- and d-amino acids are a reflection of one another due to the straightforward geometric inversion at the Cα atom. However, this presumption does not account for the directionality of the backbone dipole and the inverted propensities have never been definitively confirmed in this context. Furthermore, there is little known of how alternative side chain chirality affects the backbone conformations of isoleucine and threonine. Using a GGXGG host–guest pentapeptide system, we have completed exhaustive sampling of the conformational propensities of the d-amino acids, including d-allo-isoleucine and d-allo-threonine, using atomistic molecular dynamics simulations. Comparison of these simulations with the same systems hosting the cognate l-amino acids verifies that the intrinsic backbone conformational propensities of the d-amino acids are the inverse of their cognate l-enantiomers. Where amino acids have a chiral center in their side chain (Thr, Ile) the β-configuration affects the backbone sampling, which in turn can confer different biological properties. [ABSTRACT FROM PUBLISHER]

Integrative structural biology attempts to model the structures of protein complexes that are challenging or intractable by classical structural methods (due to size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such experimental method is chemical crosslinking mass spectrometry (XL-MS), in which protein complexes are crosslinked and characterized using liquid chromatography-mass spectrometry to pinpoint specific amino acid residues in close structural proximity. The commonly used lysine-reactive N-hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS3) have a linker arm that is 11.4 Å long when fully extended, allowing Cα (alpha carbon of protein backbone) atoms of crosslinked lysine residues to be up to ∼24 Å apart. However, XL-MS studies on proteins of known structure frequently report crosslinks that exceed this distance. Typically, a tolerance of ∼3 Å is added to the theoretical maximum to account for this observation, with limited justification for the chosen value. We used the Dynameomics database, a repository of high-quality molecular dynamics simulations of 807 proteins representative of diverse protein folds, to investigate the relationship between lysine-lysine distances in experimental starting structures and in simulation ensembles. We conclude that for DSS/BS3, a distance constraint of 26-30 Å between Cα atoms is appropriate. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL-MS results to structures or in modeling. We also discuss the comparison of XL-MS results to MD simulations and known structures as a means to test and validate experimental XL-MS methods. [ABSTRACT FROM AUTHOR]

The need for data-centric scientific tools is growing; domains such as biology, chemistry, and physics are increasingly adopting computational approaches. So, scientists must deal with the challenges of big data. To address these challenges, researchers built a visual-analytics platform named DIVE (Data Intensive Visualization Engine). DIVE is a data-agnostic, ontologically expressive software framework that can stream large datasets at interactive speeds. In particular, DIVE makes novel contributions to structured-data-model manipulation and high-throughput streaming of large, structured datasets. [ABSTRACT FROM AUTHOR]

Bovine spongiform encephalopathy (BSE), or mad cow disease, is a fatal neurodegenerative disease that is transmissible to humans and that is currently incurable. BSE is caused by the prion protein (PrP), which adopts two conformers; PrPC is the native innocuous form, which is a-helix rich; and PrPSc is the β-sheet rich misfolded form, which is infectious and forms neurotoxic species. Acidic pH induces the conversion of PrPC to PrPSc. We have performed molecular dynamics simulations of bovine PrP at various pH regimes. An acidic pH environment induced conformational changes that were not observed in neutral pH simulations. Putative misfolded structures, with nonnative β-strands formed in the flexible N-terminal domain, were found in acidic pH simulations. Two distinct pathways were observed for the formation of nonnative β-strands: at low pH, hydrophobic contacts with M129 nucleated the nonnative β-strand; at mid-pH, polar contacts involving Q168 and D178 facilitated the formation of a hairpin at the flexible N-terminus. These mid- and low pH simulations capture the process of nonnative β-strand formation, thereby improving our understanding of how PrPC misfolds into the β-sheet rich PrPSc and how pH factors into the process. [ABSTRACT FROM AUTHOR]

There are over 100 mutations in Cu/Zn superoxide dismutase (SOD1) that result in a subset of familial amyotrophic lateral sclerosis (fALS) cases. The hypothesis that dissociation of the dimer, misfolding of the monomer and subsequent aggregation of mutant SOD1 leads to fALS has been gaining support as an explanation for how these disparate missense mutations cause the same disease. These forms are only responsible for a fraction of the ALS cases; however, the rest are sporadic. Starting with a folded apo monomer, the species considered most likely to be involved in misfolding, we used high-temperature all-atom molecular dynamics simulations to explore the events of the wild-type protein unfolding through the denatured state. All simulations showed early loss of structure along the β5–β6 edge of the β-sandwich, supporting earlier findings of instability in this region. Transition state structures identified from the simulations are in good agreement with experiment, providing detailed, validated molecular models for this elusive state. Furthermore, we compare the process of thermal unfolding investigated here to that of the lethal A4V mutant-induced unfolding at physiological temperature and find that the pathways are very similar. [ABSTRACT FROM AUTHOR]

A thermostabilized variant (UVF) of the engrailed homeodomain (EnHD) was previously engineered by Mayo and co-workers. The melting temperature of the non-natural, designed protein is 50°C higher than the natural wild-type protein (>99 vs. 52°C), and the two proteins share 22% sequence identity. We have performed extensive (1 μs) all-atom, explicit solvent molecular dynamics simulations of the wild-type and engineered proteins to investigate their structural and dynamic properties at room temperature and at 100°C. Our simulations are in good agreement with nuclear magnetic resonance data available for the two proteins [nuclear Overhauser effect crosspeaks (NOEs), J-coupling constants and order parameters for EnHD; and NOEs for UVF], showing that we reproduce the backbone dynamics and side chain packing in the native state of both proteins. UVF was more dynamic at room temperature than EnHD, with respect to both its backbone and side chain motion. When the temperature was raised, the thermostable protein maintained this mobility while retaining its native conformation. EnHD, on the other hand, was unable to maintain its more rigid native structure at higher temperature and began to unfold. Heightened protein dynamics leading to promiscuous and dynamically interchangeable amino acid contacts makes UVF more tolerant to increasing temperature, providing a molecular explanation for heightened thermostability of this protein. [ABSTRACT FROM AUTHOR]

Membership in a protein domain database does not a domain make; a feature we realized when generating a consensus view of protein fold space with our consensus domain dictionary (CDD). This dictionary was used to select representative structures for characterization of the protein dynameome: the Dynameomics initiative. Through this endeavor we rejected a surprising 40% of the 1,695 folds in the CDD as being non-autonomous folding units. Although some of this was due to the challenges of grouping similar fold topologies, the dissonance between the cataloguing and structural qualification of protein domains remains surprising. Another potential factor is previously overlooked intrinsic disorder; predictions suggest that 40% of proteins have either local or global disorder. One thing is clear, filtering a structural database and ensuring a consistent definition for protein domains is crucial, and caution is prescribed when generalizations of globular domains are drawn from unfiltered protein domain datasets. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations of protein folding or unfolding, unlike most in vitro experimental methods, are performed on a single molecule. The effects of neighboring molecules on the unfolding/ folding pathway are largely ignored experimentally and simply not modeled computationally. Here, we present two all-atom, explicit solvent molecular dynamics simulations of 32 copies of the Engrailed homeodomain (EnHD), an ultrafast-folding and -unfolding protein for which the folding/unfolding pathway is well-characterized. These multimolecule simulations, in comparison with single-mole-cule simulations and experimental data, show that intermolecular interactions have little effect on the folding/unfolding pathway. EnHD unfolded by the same mechanism whether it was simulated in only water or also in the presence of other EnHD molecules. It populated the same native state, transition state, and folding intermediate in both simulation systems, and was in good agree-ment with experimental data available for each of the three states. Unfolding was slowed slightly by interactions with neighboring proteins, which were mostly hydrophobic in nature and ultimately caused the proteins to aggregate. Protein-water hydrogen bonds were also replaced with protein-protein hydrogen bonds, addition-ally contributing to aggregation. Despite the increase in protein-pro-tein interactions, the protein aggregates formed in simulation did not do so at the total exclusion of water. These simulations support the use of single-molecule techniques to study protein unfolding and also provide insight into the types of interactions that occur as pro-teins aggregate at high temperature at an atomic level. [ABSTRACT FROM AUTHOR]

High performance computing is leading to unprecedented volumes of data. Relational databases offer a robust and scalable model for storing and analyzing scientific data. However, these features do not come without a cost-significant design effort is required to build a functional and efficient repository. Modeling protein simulation data in a relational database presents several challenges: The data captured from individual simulations are large, multidimensional, and must integrate with both simulation software and external data sites. Here, we present the dimensional design and relational implementation of a comprehensive data warehouse for storing and analyzing molecular dynamics simulations using SQL Server. [ABSTRACT FROM AUTHOR]

Graphs are rapidly becoming a powerful and ubiquitous tool for the analysis of protein structure and for event detection in dynamical protein systems. Despite their rise in popularity, however, the graph representations employed to date have shared certain features and parameters that have not been thoroughly investigated. Here, we examine and compare variations on the construction of graph nodes and graph edges. We propose a graph representation based on chemical groups of similar atoms within a protein rather than residues or secondary structure and find that even very simple analyses using this representation form a powerful event detection system with significant advantages over residue-based graph representations. We additionally compare graph edges based on probability of contact to graph edges based on contact strength and analyses of the entire graph structure to an alternative and more computationally tractable node-based analysis. We develop the simplest useful technique for analyzing protein simulations based on these comparisons and use it to shed light on the speed with which static protein structures adjust to a solvated environment at room temperature in simulation. [ABSTRACT FROM AUTHOR]

As high-throughput molecular dynamics simulations of proteins become more common and the databases housing the results become larger and more prevalent, more sophisticated methods to quickly and accurately mine large numbers of trajectories for relevant information will have to be developed. One such method, which is only recently gaining popularity in molecular biology, is the continuous wavelet transform, which is especially well-suited for time course data such as molecular dynamics simulations. We describe techniques for the calculation and analysis of wavelet transforms of molecular dynamics trajectories in detail and present examples of how these techniques can be useful in data mining. We demonstrate that wavelets are sensitive to structural rearrangements in proteins and that they can be used to quickly detect physically relevant events. Finally, as an example of the use of this approach, we show how wavelet data mining has led to a novel hypothesis related to the mechanism of the protein γδ resolvase. [ABSTRACT FROM AUTHOR]

The prion diseases are a class of neurodegenerative diseases caused by the misfolding and aggregation of the prion protein (PrPC) into toxic and infectious oligomers (PrPSc). These oligomers are critical to understanding and combating these diseases. Differences in the sequence of PrP affect disease susceptibility, likely by shifting the tolerance of the protein for adaptation to PrPSc conformations and/or the recognition event between PrPSc and PrPC prior to conversion of the PrPC. We selected two sets of PrPSc-resistant mutant sequences for solvated atomistic molecular dynamics simulation to investigate the structural basis of resistance. The first group involved mutation in the X-loop (residues 164-171) resulting from selective breeding of sheep. The second group included eight mutants in mice identified by random mutagenesis targeting helix C followed by screening in cell cultures. Multiple simulations were performed of 14 different mutant and control constructs under different pH conditions for a total of 3.6 μs of simulation time. The X-loop formed a stable turn at neutral pH in wild-type PrP from both species. PrPSc-resistant mutations disrupted this turn even though only one of the mutants is in the X-loop. The X-loop is compact and buried in our previously described spiral models of PrPSc-like oligomers. On the basis of the findings presented here and in the context of the spiral oligomer model, we propose that expansion of the X-loop disrupts protofibril packing, providing a structural basis for resistance. [ABSTRACT FROM PUBLISHER]

Using molecular dynamics simulations and steady-state fluorescence spectroscopy, we have identified a conformational change in the active site of a thermophilic flavoenzyme, NADH oxidase from Thermus thermophilus HB8 (NOX). The enzyme's far-UV circular dichroism spectrum, intrinsic tryptophan fluorescence, and apparent molecular weight measured by dynamic light scattering varied little between 25 and 75°C. However, the fluorescence of the tightly bound FAD cofactor increased approximately fourfold over this temperature range. This effect appears not to be due to aggregation, unfolding, cofactor dissociation, or changes in quaternary structure. We therefore attribute the change in flavin fluorescence to a temperature-dependent conformational change involving the NOX active site. Molecular dynamics simulations and the effects of mutating aromatic residues near the flavin suggest that the change in fluorescence results from a decrease in quenching by electron transfer from tyrosine 137 to the flavin. Proteins 2012. © 2011 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

Members of the homeodomain superfamily are three-helix bundle proteins whose second and third helices form a helix-turn-helix motif (HTH). Their folding mechanism slides from the ultrafast, three-state framework mechanism for the engrailed homeodomain (EnHD), in which the HTH motif is independently stable, to an apparent two-state nucleation-condensation, model for family members with an unstable HTH motif. The folding intermediate of EnHO has nearly native HTH structure, but it is not docked with helixi. The determinant of whether twoor three-state folding was hypothesized to be the stability of the HTH substructure. Here, we describe a detailed 0-value analysis of the folding of the Piti homeodomain, which has similar ultrafast kinetics to that of EnHD. Formation of helixi was strongly coupled with formation of HTH, which was initially surprising because they are uncoupled in the EnHD folding intermediate. However, we found a key difference between Piti and EnHD: The isolated peptide corresponding to the HTH motif in Piti was not folded in the absence of Hi. Independent molecular dynamics simulations of Piti unfolding found an intermediate with Hi misfolded onto the HTH motif. The Piti folding pathway is the connection between that of EnHD and the slower folding homeodomains and provides a link in the transition of mechanisms from twoto three-state folding in this superfamily. The malleability of folding intermediates can lead to unstable substructures being stabilized by a variety of nonnative interactions, adding to the continuum of folding mechanisms. [ABSTRACT FROM AUTHOR]

We have recently completed systematic molecular dynamics simulations of 807 different proteins representing 95% of the known autonomous protein folds in an effort we refer to as Dynameomics. Here we focus on the analysis of side chain conformations and dynamics to create a dynamic rotamer library. Overall this library is derived from 31,000 occurrences of each of 86,217 different residues, or 2.7 × 10 rotamers. This dynamic library has 74% overlap of rotamer distributions with rotamer libraries derived from static high-resolution crystal structures. Seventy-five percent of the residues had an assignable primary conformation, and 68% of the residues had at least one significant alternate conformation. The average correlation time for switching between rotamers ranged from 22 ps for Met to over 8 ns for Cys; this time decreased 20-fold on the surface of the protein and modestly for dihedral angles further from the main chain. Side chain S axis order parameters were calculated and they correlated well with those derived from NMR relaxation experiments ( R = 0.9). Relationships relating the S axis order parameters to rotamer occupancy were derived. Overall the Dynameomics rotamer library offers a comprehensive depiction of side chain rotamer preferences and dynamics in solution, and more realistic distributions for dynamic proteins in solution at ambient temperature than libraries derived from crystal structures, in particular charged surface residues are better represented. Details of the rotamer library are presented here and the library itself can be downloaded at . [ABSTRACT FROM AUTHOR]

The classification of protein folds is necessarily based on the structural elements that distinguish domains. Classification of protein domains consists of two problems: the partition of structures into domains and the classification of domains into sets of similar structures (or folds). Although similar topologies may arise by convergent evolution, the similarity of their respective folding pathways is unknown. The discovery and the characterization of the majority of protein folds will be followed by a similar enumeration of available protein folding pathways. Consequently, understanding the intricacies of structural domains is necessary to understanding their collective folding pathways. We review the current state of the art in the field of protein domain classification and discuss methods for the systematic and comprehensive study of protein folding across protein fold space via atomistic molecular dynamics simulation. Finally, we discuss our large-scale Dynameomics project, which includes simulations of representatives of all autonomous protein folds. [ABSTRACT FROM AUTHOR]

Motivation: The discovery of new protein folds is a relatively rare occurrence even as the rate of protein structure determination increases. This rarity reinforces the concept of folds as reusable units of structure and function shared by diverse proteins. If the folding mechanism of proteins is largely determined by their topology, then the folding pathways of members of existing folds could encompass the full set used by globular protein domains.Results: We have used recent versions of three common protein domain dictionaries (SCOP, CATH and Dali) to generate a consensus domain dictionary (CDD). Surprisingly, 40% of the metafolds in the CDD are not composed of autonomous structural domains, i.e. they are not plausible independent folding units. This finding has serious ramifications for bioinformatics studies mining these domain dictionaries for globular protein properties. However, our main purpose in deriving this CDD was to generate an updated CDD to choose targets for MD simulation as part of our dynameomics effort, which aims to simulate the native and unfolding pathways of representatives of all globular protein consensus folds (metafolds). Consequently, we also compiled a list of representative protein targets of each metafold in the CDD.Availability and implementation: This domain dictionary is available at www.dynameomics.org.Contact: daggett@u.washington.eduSupplementary information: Supplementary data are available at Bioinformatics online. [ABSTRACT FROM PUBLISHER]

All currently known structures of proteins together define 'protein fold space'. To increase the general understanding of protein dynamics and protein folding, we selected a set of 807 proteins and protein domains that represent 95% of the currently known autonomous folded domains present in globular proteins. Native state and unfolding simulations of these representatives are now complete and accessible via a novel database containing over 11 000 simulations. Because protein folding is a microscopically reversible process, these simulations effectively sample protein folding across all of protein fold space. Here, we give an overview of how the representative proteins were selected and how the simulations were performed and validated. We then provide examples of different types of analyses that can be performed across our large set of simulations, made possible by the database approach. We further show how the unfolding simulations can be used to compare unfolding of structural elements in isolation and in different structural contexts, using as an example a short, triple stranded β-sheet that forms the WW domain and is present in several larger unrelated proteins. [ABSTRACT FROM AUTHOR]

A widely held hypothesis regarding the thermostability of thermophilic proteins states asserts that, at any given temperature, thermophilic proteins are more rigid than their mesophilic counterparts. Many experimental and computational studies have addressed this question with conflicting results. Here, we compare two homologous enzymes, one mesophilic (Escherichia coli FMN-dependent nitroreductase; NTR) and one thermophilic (Thermus thermophilus NADH oxidase; NOX), by multiple molecular dynamics simulations at temperatures from 5 to 100°C. We find that the global rigidity/flexibility of the two proteins, assessed by a variety of metrics, is similar on the time scale of our simulations. However, the thermophilic enzyme retains its native conformation to a much greater degree at high temperature than does the mesophilic enzyme, both globally and within the active site. The simulations identify the helix F–helix G ‘arm’ as the region with the greatest difference in loss of native contacts between the two proteins with increasing temperature. In particular, a network of electrostatic interactions holds helix F to the body of the protein in the thermophilic protein, and this network is absent in the mesophilic counterpart. [ABSTRACT FROM PUBLISHER]

It has been suggested that denatured proteins are predisposed toward the left-handed polyproline II (PII) conformation. One possible source of PII stability in the denatured state is water bridges. Water bridges are networks of water molecules that link nearby hydrogen bond acceptors and/or donors on proteins. On the basis of the proposed behavior of PII and water bridges, the propensity of a residue to participate in water bridges should be correlated with its PII propensity. To test this hypothesis, we analyzed the following data sets: 2351 high-resolution crystal structures, and the native and denatured states of 188 different proteins from all-atom, explicit-solvent molecular dynamics (MD) simulations, which are part of our Dynameomics effort. We found that water bridges do not explain the high frequency of PII in denatured states; such bridges are less frequent around PII than around other conformations. Thus, this analysis casts doubt on water bridges as a dominant factor determining the residue-based PII propensities. [ABSTRACT FROM PUBLISHER]

Protein L-isoaspartate O-methyltransferase (PIMT) repairs isoaspartate residues in damaged proteins, and it contains a Val–Ile polymorphismin in α5, ∼13 Å from its active site. Val119 has lower activity and thermal stability but increased affinity for endogenous substrates. Studies suggest that heterozygosity for Val/Ile favors efficient isoaspartate repair. We have performed multiple molecular dynamics simulations of 119I and 119V PIMT. Both V119 and I119 interact with the same residues throughout all of the simulations. However, the larger Ile altered the orientations of α5 and β5, both of which have co-substrate binding residues on their distal ends. I119 increases the flexibility of several residues, loosening up the S-adenosylmethionine (SAM)-binding site. These subtle changes are propagated towards the isoaspartate-docking site via residues common to both active sites. The increased mobility in 119I PIMT reorients α3, resulting in a salt-bridge network at the substrate-binding interface that disrupts several key side-chain interactions in the isoaspartate site. In contrast, 119V PIMT remains quite rigid with little change to the co-substrate binding site, which could hinder SAM's binding and release, accounting for the decreased activity. These results shed light on the molecular basis behind the decreased activity and increased specificity for endogenous substrates of 119V PIMT relative to the 119I variant. 119I PIMT catalyzes the methylation reaction but may have difficulties recognizing and orienting specific substrates due to its distorted substrate-binding site. Heterozygosity for both the Ile and Val alleles may provide the best of both worlds, allowing the fast and specific methylation of damaged proteins. [ABSTRACT FROM PUBLISHER]

Prion diseases, in which the conformational transition of the native prion protein (PrP) to a misfolded form causes aggregation and subsequent neurodegeneration, have fascinated the scientific community as this transmissible disease appears to be purely protein-based. Disease can arise due to genetic factors only. At least 30 single point mutations have been indicated to cause disease in humans. Somehow, these mutations must influence the stability, processing and/or cellular interactions of PrP, such that aggregation can occur and disease develops. In this review, the current evidence for such effects of single point mutations is discussed, indicating that PrP can be affected in many different ways, although questions remain about the mechanism by which mutations cause disease. [ABSTRACT FROM PUBLISHER]

The prion protein (PrP) is responsible for several fatal neurodegenerative diseases via conversion from its normal to disease-related isoform. The recombinant form of the protein is typically studied to investigate the conversion process. This constructs lacks the co- and post-translational modifications present in vivo, there the protein has two N-linked glycans and is bound to the outer leaflet of the plasma membrane via a glycosylphosphatidylinositol (GPI) anchor. The inherent flexibility and heterogeneity of the glycans, the plasticity of the GPI anchor, and the localization of the protein in a membrane make experimental structural characterization of biological constructs of cellular prion protein (PrPC) challenging. Yet this characterization is central in determining not only the suitability of recombinant (rec)-PrPC as a model for biological forms of the protein but also the potential role of co- and post-translational modifications on the disease process. Here, we present molecular dynamics simulations of three human prion protein constructs: (i) a protein-only construct modeling the recombinant form, (ii) a diglycosylated and soluble construct, and (iii) a diglycosylated and GPI-anchored construct bound to a lipid bilayer. We found that glycosylation and membrane anchoring do not significantly alter the structure or dynamics of PrPC, but they do appreciably modify the accessibility of the polypeptide surface PrPC. In addition, the simulations of membrane-bound PrPC revealed likely recognition domains for the disease-initiating PrPC:PrPSc (infectious and/or misfolded form of the prion protein) binding event and a potential mechanism for the observed inefficiency of conversion associated with differentially glycosylated PrP species. [ABSTRACT FROM AUTHOR]

Structure is only the first step in understanding the interactions and functions of proteins. In this paper, we explore the flexibility of proteins across a broad database of over 250 solvated protein molecular dynamics simulations in water for an aggregate simulation time of ∼6 μs. These simulations are from our Dynameomics project, and these proteins represent approximately 75% of all known protein structures. We employ principal component analysis of the atomic coordinates over time to determine the primary axis and magnitude of the flexibility of each atom in a simulation. This technique gives us both a database of flexibility for many protein fold families and a compact visual representation of a particular protein's native-state conformational space, neither of which are available using experimental methods alone. These tools allow us to better understand the nature of protein motion and to describe its relationship to other structural and dynamical characteristics. In addition to reporting general properties of protein flexibility and detailing many dynamic motifs, we characterize the relationship between protein native-state flexibility and early events in thermal unfolding and show that flexibility predicts how a protein will begin to unfold. We provide evidence that fold families have conserved flexibility patterns, and family members who deviate from the conserved patterns have very low sequence identity. Finally, we examine novel aspects of highly inflexible loops that are as important to structural integrity as conventional secondary structure. These loops, which are difficult if not impossible to locate without dynamic data, may constitute new structural motifs. [ABSTRACT FROM AUTHOR]

Here, we compare the distributions of main chain (ΦΨ) angles (i.e., Ramachandran maps) of the 20 naturally occurring amino acids in three contexts: (i) molecular dynamics (MD) simulations of Gly-Gly-X-Gly-Gly pentapeptides in water at 298 K with exhaustive sampling, where X = the amino acid in question; (ii) 188 independent protein simulations in water at 298 K from our Dynameomics Project; and (iii) static crystal and NMR structures from the Protein Data Bank. The GGXGG peptide series is often used as a model of the unstructured denatured state of proteins. The sampling in the peptide MD simulations is neither random nor uniform. Instead, individual amino acids show preferences for particular conformations, but the peptide is dynamic, and interconversion between conformers is facile. For a given amino acid, the (ΦΨ) distributions in the protein simulations and the Protein Data Bank are very similar and often distinct from those in the peptide simulations. Comparison between the peptide and protein simulations shows that packing constraints, solvation. and the tendency for particular amino acids to be used for specific structural motifs can overwhelm the "intrinsic propensities" of amino acids for particular (4,'I') conformations. We also compare our helical propensities with experimental consensus values using the host-guest method, which appear to be determined largely by context and not necessarily the intrinsic conformational propensities of the guest residues. These simulations represent an improved coil library free from contextual effects to better model intrinsic conformational propensities and provide a detailed view of conformations making up the "random coil" state. [ABSTRACT FROM AUTHOR]

The misfolding and aggregation of the prion protein (PrP) is the primary cause of a group of infectious neurodegenerative diseases including Creutzfeldt-Jacob disease in humans and Bovine Spongiform Encephalopathy in cows. A single disease can exhibit different infectious strains distinguishable by incubation time and morphology or distribution of the aggregates. Infected brain tissue from one species can be used to infect other species, but with different efficiencies, suggesting a spectrum of species compatibility. If PrP is, as widely believed, the sole component of infection, then the species and strain differences must be accounted for by the structure of the aggregates, likely influenced by each species’ PrP sequence. As there are no high-resolution data exploring this hypothesis, we performed molecular dynamics simulations of PrP for human, bovine, hamster, and D147N mutant hamster sequences at low pH to induce misfolding of the protein. We selected representative converted structures from each of the four sequences and, with the guidance of experimental data, constructed models of the infectious aggregates. Both hamster monomers showed high flexibility during conversion, suggesting hamster may more easily adopt altered conformations, which in turn may explain why it is more easily infected by some other species. Human and bovine aggregates were similar, with monomers docking in P31 symmetry to form a left-handed spiral. In contrast, hamster aggregates formed a P31 right-handed spiral. We detail the differences in the converted monomers that give rise to this difference and show that our results compare favorably with experimental data. [ABSTRACT FROM AUTHOR]