Conformational sampling with implicit solvent models: application to the PHF6 peptide in tau protein

Combination of structure-based virtual screening, molecular docking and molecular dynamics approaches for the discovery of anti-prion fibril flavonoids

Prion diseases are a group of rare neurodegenerative diseases caused by the structural conversion of cellular prion into Scrapie prion resulting aggregated fibrils. Therapy of prion diseases has been developed for several decades, especially drug designs based on the structure of prion monomers. Unfortunately, none of the designed anti-prion drugs function well clinically. To fight against prion fibrils, a drug design based on the precise structure of mammalian prion fibrils is highly required. Fortunately, based on the advantage of newly advanced cryo-electron microscopy (cryo-EM) in the deconvolution of large complexes, three prion fibril structures were resolved in the last 2 years. Based on the cryo-EM solved prion fibril structures, we are able to find some molecules fighting against prion fibrils. Quercetin, one flavonoid molecule in the polyphenol group, has been found to disaggregate the prion fibrils in vitro. In this study, we performed the molecular docking and molecular dynamics simulation on quercetin-like molecules possessing pharmacological properties to evaluate the anti-prion ability of tested molecules. As a result, four quercetin-like molecules interact with prion fibril and decrease the β-strand content by converting some β-strands into loop and helical structures to disintegrate the existing fibril structure. The results of this study are significant in the treatment of prion diseases, and the approaches used in this study are applicable to other amyloid diseases.

Comparison of the force fields on monomeric and fibrillar PHF6 of tau protein

The hexapeptide 306VQIVYK311 (PHF6) plays an important role in the aggregation of Tau protein, which is a hallmark of the Alzheimer's disease (AD). In this article, we systematically compare the effects of eight popular all-atom force fields on the monomeric and fibrillar PHF6 in the molecular dynamics (MD) simulations, which could be helpful in the computer-aided drug design against PHF6. We show that the fibrillar PHF6 prefers β-strand-like structures in all the force fields while the monomer has different structural preferences depending on the force fields. The interactions for stabilizing the fibril are further investigated. In the end, according to the interactions revealed by NMR and the stability of the fibril in the literature, we benchmark the force fields.

A Mechanism for the Inhibition of Tau Neurotoxicity: Studies with Artificial Membranes, Isolated Mitochondria, and Intact Cells

It is currently believed that molecular agents that specifically bind to and neutralize the toxic proteins/peptides, amyloid β (Aβ42), tau, and the tau-derived peptide PHF6, hold the key to attenuating the progression of Alzheimer's disease (AD). We thus tested our previously developed nonaggregating Aβ42 double mutant (Aβ42_DM) as a multispecific binder for three AD-associated molecules, wild-type Aβ42, the tau_K174Q mutant, and a synthetic PHF6 peptide. Aβ42_DM acted as a functional inhibitor of these molecules in in vitro assays and in neuronal cell-based models of AD. The double mutant bound both cytotoxic tau_K174Q and synthetic PHF6 and protected neuronal cells from the accumulation of tau in cell lysates and mitochondria. Aβ42_DM also reduced toxic intracellular levels of calcium and the overall cell toxicity induced by overexpressed tau, synthetic PHF6, Aβ42, or a combination of PHF6and Aβ42. Aβ42_DM inhibited PHF6-induced overall mitochondrial dysfunction: In particular, Aβ42_DM inhibited PHF6-induced damage to submitochondrial particles (SMPs) and suppressed PHF6-induced elevation of the ζ-potential of inverted SMPs (proxy for the inner mitochondrial membrane, IMM). PHF6 reduced the lipid fluidity of cardiolipin/DOPC vesicles (that mimic the IMM) but not DOPC (which mimics the outer mitochondrial membrane), and this effect was inhibited by Aβ42_DM. This inhibition may be explained by the conformational changes in PHF6 induced by Aβ42_DM in solution and in membrane mimetics. On this basis, the paper presents a mechanistic explanation for the inhibitory activity of Aβ42_DM against Aβ42- and tau-induced membrane permeability and cell toxicity and provides confirmatory evidence for its protective function in neuronal cells.

Molecular mechanism of amyloidogenicity and neurotoxicity of a pro-aggregated tau mutant in the presence of histidine tautomerism via replica-exchange simulation

Considering ΔK280 tau mutation, δε isomer with highest sheet content may accelerate aggregation; generating small compounds to inhibit this would help tp prevent tauopathies.

A droplet reactor on a super-hydrophobic surface allows control and characterization of amyloid fibril growth

Methods to produce protein amyloid fibrils, in vitro, and in situ structure characterization, are of primary importance in biology, medicine, and pharmacology. We first demonstrated the droplet on a super-hydrophobic substrate as the reactor to produce protein amyloid fibrils with real-time monitoring of the growth process by using combined light-sheet microscopy and thermal imaging. The molecular structures were characterized by Raman spectroscopy, X-ray diffraction and X-ray scattering. We demonstrated that the convective flow induced by the temperature gradient of the sample is the main driving force in the growth of well-ordered protein fibrils. Particular attention was devoted to PHF6 peptide and full-length Tau441 protein to form amyloid fibrils. By a combined experimental with the molecular dynamics simulations, the conformational polymorphism of these amyloid fibrils were characterized. The study provided a feasible procedure to optimize the amyloid fibrils formation and characterizations of other types of proteins in future studies.

Zhang et al present an integrated real-time imaging and flow field control platform based on water droplet evaporation on super-hydrophobic substrate (SHS) to enable amyloid fibril aggregation. They apply this methodology to observe structural polymorphism in PHF6 peptide and full length Tau441.

What makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces

Molecular dynamics simulations of the aqueous biomolecule–graphene interface have predicted the free energy of adsorption of amino acids and the structure of peptides.

Mechanism of amyloid-β fibril elongation

Amyloid-β is an intrinsically disordered protein that forms fibrils in the brains of patients with Alzheimer's disease. To explore factors that affect the process of fibril growth, we computed the free energy associated with disordered amyloid-β monomers being added to growing amyloid fibrils using extensive molecular dynamics simulations coupled with umbrella sampling. We find that the mechanisms of Aβ40 and Aβ42 fibril elongation have many features in common, including the formation of an obligate on-pathway β-hairpin intermediate that hydrogen bonds to the fibril core. In addition, our data lead to new hypotheses for how fibrils may serve as secondary nucleation sites that can catalyze the formation of soluble oligomers, a finding in agreement with recent experimental observations. These data provide a detailed mechanistic description of amyloid-β fibril elongation and a structural link between the disordered free monomer and the growth of amyloid fibrils and soluble oligomers.

Variational Optimization of an All-Atom Implicit Solvent Force Field to Match Explicit Solvent Simulation Data

The development of accurate implicit solvation models with low computational cost is essential for addressing many large-scale biophysical problems. Here, we present an efficient solvation term based on a Gaussian solvent-exclusion model (EEF1) for simulations of proteins in aqueous environment, with the primary aim of having a good overlap with explicit solvent simulations, particularly for unfolded and disordered states - as would be needed for multiscale applications. In order to achieve this, we have used a recently proposed coarse-graining procedure based on minimization of an entropy-related objective function to train the model to reproduce the equilibrium distribution obtained from explicit water simulations. Via this methodology, we have optimized both a charge screening parameter and a backbone torsion term against explicit solvent simulations of an α-helical and a β-stranded peptide. The performance of the resulting effective energy function, termed EEF1-SB, is tested with respect to the properties of folded proteins, the folding of small peptides or fast-folding proteins, and NMR data for intrinsically disordered proteins. The results show that EEF1-SB provides a reasonable description of a wide range of systems, but its key advantage over other methods tested is that it captures very well the structure and dimension of disordered or weakly structured peptides. EEF1-SB is thus a computationally inexpensive (~ 10 times faster than Generalized-Born methods) and transferable approximation for treating solvent effects.

The role of electrostatic interactions on klentaq1 insight for domain separation

We investigated the relationship between the thermostability of Klentaq1 and factors stabilizing interdomain interactions. When thermal adaptation of Klentaq1 was analyzed at the atomic level, the protein was stable at 300 and 350 K. It gradually unfolded at 373 K and almost spontaneously unfolded at 400 K. Domain separation was induced by disrupting electrostatic interactions in two salt bridges formed by Lys354-Glu445 and Asp371-Arg435 on the interface domain. The role of these interactions in protein stability was evaluated by comparing free energy solvation (ΔΔG_solv) between wild type and mutants. Substitution of Asp371 by Glu or Asn, and also Glu445 by Asn resulted in a positive value of ΔΔG_solv, suggesting that mutations destabilized the protein structure. Nevertheless, substitution of Glu445 by Asp gave a negative value to ΔΔG_solv reflecting increasing protein stability. Our results demonstrate that interactions at the interface domains of Klentaq1 are essential factors correlated with the Klentaq1 thermostability.

Finding order within disorder: elucidating the structure of proteins associated with neurodegenerative disease

A number of neurodegenerative disorders, such as Alzheimer's disease and Parkinson's disease, involve the formation of protein aggregates. The primary constituent of these aggregates belongs to a unique class of heteropolymers known as intrinsically disordered proteins (IDPs). While many proteins fold to a unique conformation that is determined by their amino acid sequence, IDPs do not adopt a single well-defined conformation in solution. Instead, they populate a heterogeneous set of conformers under physiological conditions. Despite this intrinsic propensity for disorder, a number of these proteins can form ordered aggregates both in vitro and in vivo. As the formation of these aggregates may play an important role in disease pathogenesis, a detailed structural characterization of these proteins and their mechanism of aggregation is of critical importance. However, new methods are needed to understand the diversity of structures that make up the unfolded ensemble of these systems. In this review, we discuss recent advances in the structural analysis and modeling of IDPs involved in neurodegenerative diseases. While there are challenges in both the experimental characterization and the modeling of such proteins, a comprehensive understanding of the structure of IDPs will likely facilitate the development of effective therapies for a number of neurodegenerative diseases.

Structure and dynamics of human vimentin intermediate filament dimer and tetramer in explicit and implicit solvent models

Intermediate filaments, in addition to microtubules and microfilaments, are one of the three major components of the cytoskeleton in eukaryotic cells, and play an important role in mechanotransduction as well as in providing mechanical stability to cells at large stretch. The molecular structures, mechanical and dynamical properties of the intermediate filament basic building blocks, the dimer and the tetramer, however, have remained elusive due to persistent experimental challenges owing to the large size and fibrillar geometry of this protein. We have recently reported an atomistic-level model of the human vimentin dimer and tetramer, obtained through a bottom-up approach based on structural optimization via molecular simulation based on an implicit solvent model (Qin et al. in PLoS ONE 2009 4(10):e7294, 9). Here we present extensive simulations and structural analyses of the model based on ultra large-scale atomistic-level simulations in an explicit solvent model, with system sizes exceeding 500,000 atoms and simulations carried out at 20 ns time-scales. We report a detailed comparison of the structural and dynamical behavior of this large biomolecular model with implicit and explicit solvent models. Our simulations confirm the stability of the molecular model and provide insight into the dynamical properties of the dimer and tetramer. Specifically, our simulations reveal a heterogeneous distribution of the bending stiffness along the molecular axis with the formation of rather soft and highly flexible hinge-like regions defined by non-alpha-helical linker domains. We report a comparison of Ramachandran maps and the solvent accessible surface area between implicit and explicit solvent models, and compute the persistence length of the dimer and tetramer structure of vimentin intermediate filaments for various subdomains of the protein. Our simulations provide detailed insight into the dynamical properties of the vimentin dimer and tetramer intermediate filament building blocks, which may guide the development of novel coarse-grained models of intermediate filaments, and could also help in understanding assembly mechanisms.

Residual structure within the disordered C-terminal segment of p21(Waf1/Cip1/Sdi1) and its implications for molecular recognition

Probably the most unusual class of proteins in nature is the intrinsically unstructured proteins (IUPs), because they are not structured yet play essential roles in protein-protein signaling. Many IUPs can bind different proteins, and in many cases, adopt different bound conformations. The p21 protein is a small IUP (164 residues) that is ubiquitous in cellular signaling, for example, cell cycle control, apoptosis, transcription, differentiation, and so forth; it binds to approximately 25 targets. How does this small, unstructured protein recognize each of these targets with high affinity? Here, we characterize residual structural elements of the C-terminal segment of p21 encompassing residues 145-164 using a combination of NMR measurements and molecular dynamics simulations. The N-terminal half of the peptide has a significant helical propensity which is recognized by calmodulin while the C-terminal half of the peptide prefers extended conformations that facilitate binding to the proliferating cell nuclear antigen (PCNA). Our results suggest that the final bound conformations of p21 (145-164) pre-exist in the free peptide even without its binding partners. While the conformational flexibility of the p21 peptide is essential for adapting to diverse binding environments, the intrinsic structural preferences of the free peptide enable promiscuous yet high affinity binding to a diverse array of molecular targets.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

PROTDES: CHARMM toolbox for computational protein design

We present an open-source software able to automatically mutate any residue positions and find the best aminoacids in an arbitrary protein structure without requiring pairwise approximations. Our software, PROTDES, is based on CHARMM and it searches automatically for mutations optimizing a protein folding free energy. PROTDES allows the integration of molecular dynamics within the protein design. We have implemented an heuristic optimization algorithm that iteratively searches the best aminoacids and their conformations for an arbitrary set of positions within a structure. Our software allows CHARMM users to perform protein design calculations and to create their own procedures for protein design using their own energy functions. We show this by implementing three different energy functions based on different solvent treatments: surface area accessibility, generalized Born using molecular volume and an effective energy function. PROTDES, a tutorial, parameter sets, configuration tools and examples are freely available at http://soft.synth-bio.org/protdes.html.

Calculation of protein heat capacity from replica-exchange molecular dynamics simulations with different implicit solvent models

The heat capacity has played a major role in relating microscopic and macroscopic properties of proteins and their disorder-order phase transition of folding. Its calculation by atomistic simulation methods remains a significant challenge due to the complex and dynamic nature of protein structures, their solvent environment, and configurational averaging. To better understand these factors on calculating a protein heat capacity, we provide a comparative analysis of simulation models that differ in their implicit solvent description and force-field resolution. Our model protein system is the src Homology 3 (SH3) domain of alpha-spectrin, and we report a series of 10 ns replica-exchange molecular dynamics simulations performed at temperatures ranging from 298 to 550 K, starting from the SH3 native structure. We apply the all-atom CHARMM22 force field with different modified analytical generalized Born solvent models (GBSW and GBMV2) and compare these simulation models with the distance-dependent dielectric screening of charge-charge interactions. A further comparison is provided with the united-atom CHARMM19 plus a pairwise GB model. Unfolding-folding transition temperatures of SH3 were estimated from the temperature-dependent profiles of the heat capacity, root-mean-square distance from the native structure, and the fraction of native contacts, each calculated from the density of states by using the weighted histogram analysis method. We observed that, for CHARMM22, the unfolding transition and energy probability density were quite sensitive to the implicit solvent description, in particular, the treatment of the protein-solvent dielectric boundary in GB models and their surface-area-based hydrophobic term. Among the solvent models tested, the calculated melting temperature varied in the range 353-438 K and was higher than the experimental value near 340 K. A reformulated GBMV2 model of employing a smoother molecular-volume dielectric interface was the most accurate in reproducing the native conformation and a two-state folding landscape, although the melting transition temperature did not show the smallest deviation from experiment. For the lower-resolution CHARMM19/GB model, the simulations failed to yield a bimodal energy distribution, yet the melting temperature was observed to be a good estimate of higher-resolution simulation models. We also demonstrate that a careful analysis of a relatively long simulation is necessary to avoid trapping in local minima and to find a true thermodynamic transition temperature.

The effect of a DeltaK280 mutation on the unfolded state of a microtubule-binding repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimer's disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a DeltaK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and DeltaK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of beta-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates.

Structural and Pathway Complexity of β-Strand Reorganization within Aggregates of Human Transthyretin(105−115) Peptide

Interstrand conformational rearrangements of human transthyretin peptide (TTR(105-115)) within dimeric aggregates were simulated by means of molecular dynamics (MD) with implicit solvation model for a total length of 48 micros. The conformations sampled in the MD simulations were clustered to identify free energy minima without any projections of free energy surface. A connected graph was constructed with nodes (=clusters) and edges corresponding to free energy minima and transitions between nodes, respectively. This connected graph which reflects the complexity of the free energy surface was used to extract the transition disconnectivity graph, which reflects the overall free energy barriers between pairs of free energy minima but does not contain information on transition paths. The routes of transitions between important free energy minima were obtained by further processing the original graph and the MD data. We have found that both parallel and antiparallel aggregates are populated. The parallel aggregates with different alignment patterns are separated by nonnegligible free energy barriers. Multiroutes exist in the interstrand conformational reorganization. Most visited routes do not dominant the kinetics, while less visited routes contribute a little each but they are numerous and their total contributions are actually dominant. There are various kinds of reptation motions, including those through a beta-bulge, side-chain aided reptation, and flipping or rotation of a hairpin formed by one strand.

Structural malleability of plasticins: Preorganized conformations in solution and relevance for antimicrobial activity

Plasticins (23 long-residue glycine-leucine-rich dermaseptin-related peptides produced by the skin of South American hylids) have very similar amino acid sequences, hydrophobicities, and amphipathicities, but differ in their membrane-damaging properties and structurations (i.e. destabilized helix states, beta-hairpin, beta-sheet, and disordered states) at anionic and zwitterionic membrane interfaces. Structural malleability of plasticins in aqueous solutions together with parameters that may govern their ability to fold within beta-hairpin like structures were analyzed through circular dichroism and FTIR spectroscopic studies completed by molecular dynamics simulations in polar mimetic media. The goal of this study was to probe to which extent pre-existent peptide conformations, i.e. intrinsic "conformational landscape", may be responsible for variability in bioactive conformation and antimicrobial/hemolytic mechanisms of action of these peptides in relation with their various membrane disturbing properties. All plasticins present a turn region that does not always result in folding into a beta-hairpin shaped conformation. Residue at position 8 plays a major role in initiating the folding, while position 12 is not critical. Conformational stability has no major impact on antimicrobial efficacy. However, preformed beta-hairpin in solution may act as a conformational lock that prevents switch to alpha-helical structure. This lock lowers the antimicrobial efficiency and explains subtle differences in potencies of the most active antimicrobial plasticins.