Exposure of Hydrophobic Core in Human Prion Protein Pathogenic Mutant H187R

Dysbiosis of the gut microbiota and its effect on α-synuclein and prion protein misfolding: consequences for neurodegeneration

Abnormal behavior of α-synuclein and prion proteins is the hallmark of Parkinson's disease (PD) and prion illnesses, respectively, being complex neurological disorders. A primary cause of protein aggregation, brain injury, and cognitive loss in prion illnesses is the misfolding of normal cellular prion proteins (PrPC) into an infectious form (PrPSc). Aggregation of α-synuclein causes disruptions in cellular processes in Parkinson's disease (PD), leading to loss of dopamine-producing neurons and motor symptoms. Alteration in the composition or activity of gut microbes may weaken the intestinal barrier and make it possible for prions to go from the gut to the brain. The gut-brain axis is linked to neuroinflammation; the metabolites produced by the gut microbiota affect the aggregation of α-synuclein, regulate inflammation and immunological responses, and may influence the course of the disease and neurotoxicity of proteins, even if their primary targets are distinct proteins. This thorough analysis explores the complex interactions that exist between the gut microbiota and neurodegenerative illnesses, particularly Parkinson's disease (PD) and prion disorders. The involvement of the gut microbiota, a complex collection of bacteria, archaea, fungi, viruses etc., in various neurological illnesses is becoming increasingly recognized. The gut microbiome influences neuroinflammation, neurotransmitter synthesis, mitochondrial function, and intestinal barrier integrity through the gut-brain axis, which contributes to the development and progression of disease. The review delves into the molecular mechanisms that underlie these relationships, emphasizing the effects of microbial metabolites such as bacterial lipopolysaccharides (LPS), and short-chain fatty acids (SCFAs) in regulating brain functioning. Additionally, it looks at how environmental influences and dietary decisions affect the gut microbiome and whether they could be risk factors for neurodegenerative illnesses. This study concludes by highlighting the critical role that the gut microbiota plays in the development of Parkinson's disease (PD) and prion disease. It also provides a promising direction for future research and possible treatment approaches. People afflicted by these difficult ailments may find hope in new preventive and therapeutic approaches if the role of the gut microbiota in these diseases is better understood.

Human prion protein: exploring the thermodynamic stability and structural dynamics of its pathogenic mutants

Human familial prion diseases are known to be associated with different single-point mutants of the gene coding for prion protein with a primary focus at several locations of the globular domain. We have identified 12 different single-point pathogenic mutants of human prion protein (HuPrP) with the help of extensive perturbations/mutation technique at multiple locations of HuPrP sequence related to potentiality towards conformational disorders. Among these, some of the mutants include pathogenic variants that corroborate well with the literature reported proteins while majority include some unique single-point mutants that are either not explicitly studied early or studied for variants with different residues at the specific position. Primarily, our study sheds light on the unfolding mechanism of the above mentioned mutants in depth. Besides, we could identify some mutants under investigation that demonstrates not only unfolding of the helical structures but also extension and generation of the β-sheet structures and or simultaneously have highly exposed hydrophobic surface which is assumed to be linked with the production of aggregate/fibril structures of the prion protein. Among the identified mutants, Q212E needs special attention due to its maximum exposure of hydrophobic core towards solvent and E200Q is found to be important due to its maximum extent of β-content. We are also able to identify different respective structural conformations of the proteins according to their degree of structural unfolding and those conformations can be extracted and further studied in detail. Communicated by Ramaswamy H. Sarma.

Unraveling the Molecular Mechanism of Prion H2 C-Terminus Misfolding by Metadynamics Simulations

Conformational transition from the normal cellular form of prion protein (PrP^C) to the pathogenic "scrapie" form (PrP^Sc) is considered to be a key event in the occurrence of prion disease. Additionally, the H2 C-terminus is widely considered to be a vital site for PrP conformational transition, which can be used as an important region to explore the potential mechanism of PrP misfolding. Therefore, to study the misfolding mechanism of PrP, 500 ns well-tempered metadynamics simulations were performed by focusing on the H2 C-terminus of PrP. For comparison, three systems were designed in total, including PrP in neutral and acidic conditions, as well as H187R mutant. The resulting free energy surfaces (FESs) obtained from metadynamics simulations reveal that acidic conditions and H187R mutation can facilitate PrP misfolding by decreasing free energy barriers for conformational transition and forming energy stable conformational states. Further analyses aimed at H2 C-terminus show that due to the increase of positive charge on residue 187 in both acidic and H187R systems, the electrostatic repulsion of residue 187 and R136/R156 increases greatly, which disrupts the electrostatic interaction network around H2 C-terminus and exposes the hydrophobic core to the solvent. Taken together, acidic conditions and H187R mutation can accelerate PrP misfolding mainly by forming more energetically stable metastable conformations with lower free energy barriers, and electrostatic network disruption involving residue 187 drives the initial misfolding of H2 C-terminus. This study provides quantitative insight into the related function of the H2 C-terminus in the PrP misfolding process, which may guide H2 C-terminus mediated drug design in the future.

Structural insight into conformational change in prion protein by breakage of electrostatic network around H187 due to its protonation

A conformational change from normal prion protein(PrP^C) to abnormal prion protein(PrP^SC) induces fatal neurodegenerative diseases. Acidic pH is well-known factors involved in the conformational change. Because the protonation of H187 is strongly linked to the change in PrP stability, we examined the charged residues R156, E196, and D202 around H187. Interestingly, there have been reports on pathological mutants, such as H187R, E196A, and D202N. In this study, we focused on how an acidic pH and pathological mutants disrupt this electrostatic network and how this broken network destabilizes PrP structure. To do so, we performed a temperature-based replica-exchange molecular dynamics (T-REMD) simulation using a cumulative 252 μs simulation time. We measured the distance between amino acids comprising four salt bridges (R156–E196/D202 and H187–E196/D202). Our results showed that the spatial configuration of the electrostatic network was significantly altered by an acidic pH and mutations. The structural alteration in the electrostatic network increased the RMSF value around the first helix (H1). Thus, the structural stability of H1, which is anchored to the H2–H3 bundle, was decreased. It induces separation of R156 from the electrostatic network. Analysis of the anchoring energy also shows that two salt-bridges (R156-E196/D202) are critical for PrP stability.

Structural mechanisms of oligomer and amyloid fibril formation by the prion protein

Misfolding and aggregation of the prion protein is responsible for multiple neurodegenerative diseases. Works from several laboratories on folding of both the WT and multiple pathogenic mutant variants of the prion protein have identified several structurally dissimilar intermediates, which might be potential precursors to misfolding and aggregation. The misfolded aggregates themselves are morphologically distinct, critically dependent on the solution conditions under which they are prepared, but always β-sheet rich. Despite the lack of an atomic resolution structure of the infectious pathogenic agent in prion diseases, several low resolution models have identified the β-sheet rich core of the aggregates formed in vitro, to lie in the α2-α3 subdomain of the prion protein, albeit with local stabilities that vary with the type of aggregate. This feature article describes recent advances in the investigation of in vitro prion protein aggregation using multiple spectroscopic probes, with particular focus on (1) identifying aggregation-prone conformations of the monomeric protein, (2) conditions which trigger misfolding and oligomerization, (3) the mechanism of misfolding and aggregation, and (4) the structure of the misfolded intermediates and final aggregates.

Characterization of mutations in PRNP (prion) gene and their possible roles in neurodegenerative diseases

Abnormal prion proteins are responsible for several fatal neurodegenerative diseases in humans and in animals, including Creutzfeldt-Jakob disease (CJD), Gerstmann-Sträussler-Scheinker disease, and fatal familial insomnia. Genetics is important in prion diseases, but in the most cases, cause of diseases remained unknown. Several mutations were found to be causative for prion disorders, and the effect of mutations may be heterogeneous. In addition, different prion mutations were suggested to play a possible role in additional phenotypes, such as Alzheimer's type pathology, spongiform encephalopathy, or frontotemporal dementia. Pathogenic nature of several prion mutations remained unclear, such as M129V and E219K. These two polymorphic sites were suggested as either risk factors for different disorders, such as Alzheimer's disease (AD), variant CJD, or protease-sensitive prionopathy, and they can also be disease-modifying factors. Pathological overlap may also be possible with AD or progressive dementia, and several patients with prion mutations were initially diagnosed with AD. This review also introduces briefly the diagnosis of prion diseases and the issues with their diagnosis. Since prion diseases have quite heterogeneous phenotypes, a complex analysis, a combination of genetic screening, cerebrospinal fluid biomarker analysis and imaging technologies could improve the early disease diagnosis.

Peculiarities of copper binding to alpha-synuclein

Heavy metals have been implicated as the causative agents for the pathogenesis of the most prevalent neurodegenerative disease. Various mechanisms have been proposed to explain the toxic effects of metals ranging from metal-induced oxidation of protein to metal-induced changes in the protein conformation. Aggregation of a-synuclein is implicated in Parkinson's disease (PD), and various metals, including copper, constitute a prominent group of alpha-synuclein aggregation enhancers. In this study, we have systematically characterized the a-synuclein-Cu21 binding sites and analyzed the possible role of metal binding in a-synuclein fibrillation using a set of biophysical techniques, such as electron paramagnetic resonance (EPR), electron spin-echo envelope modulation (ESEEM), circular dichroism (CD), and size exclusion chromatography (SEC). Our analyses indicated that a-synuclein possesses at least two binding sites for Cu21. We have been able to locate one of the binding sites in the N-terminal region. Furthermore, based on the EPR studies of model peptides and Beta-synuclein, we concluded that the suspected His residue did not appear to participate in strong Cu21 binding.

Unraveling the Molecular Mechanism of pH-Induced Misfolding and Oligomerization of the Prion Protein

The misfolding of the prion protein (PrP) to aggregated forms is linked to several neurodegenerative diseases. Misfolded oligomeric forms of PrP are associated with neurotoxicity and/or infectivity, but the molecular mechanism by which they form is still poorly understood. A reduction in pH is known to be a key factor that triggers misfolded oligomer formation by PrP, but the residues whose protonation is linked with misfolding remain unidentified. The structural consequences of the protonation of these residues also remain to be determined. In the current study, amino acid residues whose protonation is critical for PrP misfolding and oligomerization have been identified using site-directed mutagenesis and misfolding/oligomerization assays. It is shown that the protonation of either H186 or D201, which mimics the effects of pathogenic mutations (H186R and D201N) at both residue sites, is critically linked to the stability, misfolding and oligomerization of PrP. Hydrogen-deuterium exchange studies coupled with mass spectrometry show that the protonation of either H186 or D201 leads to the same common structural change: increased structural dynamics in helix 1 and that in the loop between helix 1 and β-strand 2. It is shown that the protonation of either of these residues is sufficient for accelerating misfolded oligomer formation, most likely because the protonation of either residue causes the same structural perturbation. Hence, the increased structural dynamics in helix 1 and that in the loop between helix 1 and β-strand 2 appear to play an early critical role in acid-induced misfolding of PrP.

Molecular Mechanism of the Misfolding and Oligomerization of the Prion Protein: Current Understanding and Its Implications

Prion diseases, also known as transmissible spongiform encephalopathies, make up a group of fatal neurodegenerative disorders linked with the misfolding and aggregation of the prion protein (PrP). Although it is not yet understood how the misfolding of PrP induces neurodegeneration, it is widely accepted that the formation of misfolded prion protein (termed PrP(Sc)) is both the triggering event in the disease and the main component of the infectious agent responsible for disease transmission. Despite the clear involvement of PrP(Sc) in prion diseases, the exact composition of PrP(Sc) is not yet well-known. Recent studies show that misfolded oligomers of PrP could, however, be responsible for neurotoxicity and/or infectivity in the prion diseases. Hence, understanding the molecular mechanism of formation of the misfolded oligomers of PrP is critical for developing an understanding about the prion diseases and for developing anti-prion therapeutics. This review discusses recent advances in understanding the molecular mechanism of misfolded oligomer formation by PrP and its implications for the development of anti-prion therapeutics.

Molecular dynamics studies on the buffalo prion protein

It was reported that buffalo is a low susceptibility species resisting to transmissible spongiform encephalopathies (TSEs) (same as rabbits, horses, and dogs). TSEs, also called prion diseases, are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of species (except for rabbits, dogs, horses, and buffalo), manifesting as scrapie in sheep and goats; bovine spongiform encephalopathy (BSE or "mad-cow" disease) in cattle; chronic wasting disease in deer and elk; and Creutzfeldt-Jakob diseases, Gerstmann-Sträussler-Scheinker syndrome, fatal familial insomnia, and Kulu in humans etc. In molecular structures, these neurodegenerative diseases are caused by the conversion from a soluble normal cellular prion protein (PrP(C)), predominantly with α-helices, into insoluble abnormally folded infectious prions (PrP(Sc)), rich in β-sheets. In this article, we studied the molecular structure and structural dynamics of buffalo PrP(C) (BufPrP(C)), in order to understand the reason why buffalo is resistant to prion diseases. We first did molecular modeling of a homology structure constructed by one mutation at residue 143 from the NMR structure of bovine and cattle PrP(124-227); immediately we found that for BufPrP(C)(124-227), there are five hydrogen bonds (HBs) at Asn143, but at this position, bovine/cattle do not have such HBs. Same as that of rabbits, dogs, or horses, our molecular dynamics studies also revealed there is a strong salt bridge (SB) ASP178-ARG164 (O-N) keeping the β2-α2 loop linked in buffalo. We also found there is a very strong HB SER170-TYR218 linking this loop with the C-terminal end of α-helix H3. Other information, such as (i) there is a very strong SB HIS187-ARG156 (N-O) linking α-helices H2 and H1 (if mutation H187R is made at position 187, then the hydrophobic core of PrP(C) will be exposed (L.H. Zhong (2010). Exposure of hydrophobic core in human prion protein pathogenic mutant H187R. Journal of Biomolecular Structure and Dynamics 28(3), 355-361)), (ii) at D178, there is a HB Y169-D178 and a polar contact R164-D178 for BufPrP(C) instead of a polar contact Q168-D178 for bovine PrP(C) (C.J. Cheng, & V. Daggett. (2014). Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH. Biomolecules 4(1), 181-201), (iii) BufPrP(C) owns three 310 helices at 125-127, 152-156, and in the β2-α2 loop, respectively, and (iv) in the β2-α2 loop, there is a strong π-π stacking and a strong π-cation F175-Y169-R164.(N)NH2, has been discovered.

In silico analysis of prion protein mutants: a comparative study by molecular dynamics approach

Polymorphisms in the human prion proteins lead to amino acid substitutions by the conversion of PrPC to PrPSc and amyloid formation, resulting in prion diseases such as familial Creutzfeldt-Jakob disease, Gerstmann-Straussler-Scheinker disease and fatal familial insomnia. Cation-π interaction is a non-covalent binding force that plays a significant role in protein stability. Here, we employ a novel approach by combining various in silico tools along with molecular dynamics simulation to provide structural and functional insight into the effect of mutation on the stability and activity of mutant prion proteins. We have investigated impressions of prevalent mutations including 1E1S, 1E1P, 1E1U, 1E1P, 1FKC and 2K1D on the human prion proteins and compared them with wild type. Structural analyses of the models were performed with the aid of molecular dynamics simulation methods. According to our results, frequently occurred mutations were observed in conserved sequences of human prion proteins and the most fluctuation values appear in the 2K1D mutant model at around helix 4 with residues ranging from 190 to 194. Our observations in this study could help to further understand the structural stability of prion proteins.

Two misfolding routes for the prion protein around pH 4.5

Using molecular dynamics simulations, we show that the prion protein (PrP) exhibits a dual behavior, with two possible transition routes, upon protonation of H187 around pH 4.5, which mimics specific conditions encountered in endosomes. Our results suggest a picture in which the protonated imidazole ring of H187 experiences an electrostatic repulsion with the nearby guanidinium group of R136, to which the system responds by pushing either H187 or R136 sidechains away from their native cavities. The regions to which H187 and R136 are linked, namely the C-terminal part of H2 and the loop connecting S1 to H1, respectively, are affected in a different manner depending on which pathway is taken. Specific in vivo or in vitro conditions, such as the presence of molecular chaperones or a particular experimental setup, may favor one transition pathway over the other, which can result in very different [Formula: see text] monomers. This has some possible connections with the observation of various fibril morphologies and the outcome of prion strains. In addition, the finding that the interaction of H187 with R136 is a weak point in mammalian PrP is supported by the absence of the [Formula: see text] residue pair in non-mammalian species that are known to be resistant to prion diseases.

Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies

Human pancreatic cholesterol esterase (hCEase) is one of the lipases found to involve in the digestion of large and broad spectrum of substrates including triglycerides, phospholipids, cholesteryl esters, etc. The presence of bile salts is found to be very important for the activation of hCEase. Molecular dynamic simulations were performed for the apoform and bile salt complexed form of hCEase using the co-ordinates of two bile salts from bovine CEase. The stability of the systems throughout the simulation time was checked and two representative structures from the highly populated regions were selected using cluster analysis. These two representative structures were used in pharmacophore model generation. The generated pharmacophore models were validated and used in database screening. The screened hits were refined for their drug-like properties based on Lipinski's rule of five and ADMET properties. The drug-like compounds were further refined by molecular docking simulation using GOLD program based on the GOLD fitness score, mode of binding, and molecular interactions with the active site amino acids. Finally, three hits of novel scaffolds were selected as potential leads to be used in novel and potent hCEase inhibitor design. The stability of binding modes and molecular interactions of these final hits were re-assured by molecular dynamics simulations.

Exploring the mechanism of a regulatory SNP of KLK3 by molecular dynamics simulation

The SNP -158G>A of KLK3 has been validated as a regulatory SNP (rSNP) by molecular biology assays, but the mechanism of how it affects the binding of an androgen receptor (AR) homodimer with DNA is unclear. In the current study, molecular dynamics simulation was adopted to explain its inner cause. Based on a recent review), three types of intermolecular forces were analyzed, and the differences among them were compared between complexes containing -158 A:T and -158 G:C. Extra hydrophobic contacts caused by the methyl group on the mutated thymine were the most crucial factor to the regulatory effect of this rSNP. Further analysis concerning the relative motion of the two recognition helixes of the AR homodimer indicated that the hydrophobic interactions between the recognition helix B and the major groove containing -158 A:T changed that helix's motion greatly from swaying in a plane at free state to vibrating slightly around an equilibrium position. A relatively full explanation on the occurrence of rSNP -158G>A is presented here.

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (-1.89 kcal/mol) and the hexamer (-1.97 kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.

Role of ligand binding in structural organization of add A-riboswitch aptamer: a molecular dynamics simulation

The specific binding of ligands is the first step of gene expression or translation regulation by riboswitches. However, understanding the mechanism of the specific binding is still difficult because the tertiary structures of the riboswitch aptamers are available almost only for ligand-bound state at present. In this paper we hope to give some insights into this problem through the studies of the role of ligand-aptamer interaction in the structural organization of add A-riboswitch aptamer, based on the crystal structure of the ligand-bound aptamer. We use all-atom molecular dynamics to simulate the behaviors of the aptamer in ligand-bound, free and mutated states by Amber force field. The results show that the correct paring of the ligand adenine with the nucleotide U74 in the binding pocket is crucial to stabilizing the conformations of the ligand-bound aptamer, especially the helix P1 connecting the expression platform. Our results also suggest that both the nucleotide U74 and U51 may be the key sites of the ligand recognition but the former has much higher probability as the initial docking site. This is in agreement with previous experimental results.

Mechanism of antisense oligonucleotide interaction with natural RNAs

Oligonucleotides find several numbers of applications: as diagnostic probes, RT and PCR primers and antisense agents due to their ability of forming specific interactions with complementary nucleotide sequences within nucleic acids. These interactions are strongly affected by accessibility of the target sequence in the RNA structure. In the present work the mechanism of invasion of RNA structure by oligonucleotide was investigated using a model system: yeast tRNA(Phe) and oligonucleotides complementary to the 3'-part of this molecule. Kinetics of interaction of oligonucleotides with in vitro transcript of yeast tRNAPhe was studied using stopped-flow technique with fluorescence quenching detection, 5'-DABCYL labeled oligonucleotide was hybridized with 3'-fluorescein labeled tRNA(Phe). The results evidence for a four-step invasion process of the oligonucleotide-RNA complex formation. The process is initiated by formation of transition complexes with nucleotides in the T-loop and ACCA sequence. This complex formation is followed by RNA unfolding and formation of an extended heteroduplex with the oligonucleotide via strand displacement process. Computer modeling of oligonucleotide-tRNA(Phe) interaction revealed potential factors that could favor transition complexes formation and confirmed the proposed mechanism, showing the oligonucleotide to be a molecular "wedge". Our data evidence that oligonucleotide invasion into structured RNA is initiated by loop-single strand interactions, similar to the initial step of the antisense RNA-RNA interactions. The obtained results can be used for choosing efficient oligonucleotide probes.

Predicting sumoylation site by feature selection method

The small ubiquitin-like modifier (SUMO) proteins are a kind of proteins that can be attached to a series of proteins. The sumoylation of protein is an important posttranslational modification. Thus, the prediction of the sumoylation site of a given protein is significant. Here we employed a combined method to perform this task. We predicted the sumoylation site of a protein by a two-staged procedure. At the first stage, whether a protein would be sumoylated was predicted; whereas at the second stage, the sumoylation sites of the protein were predicted if it was determined to be modified by SUMO at the first stage. At the first stage, we encoded a protein with protein families (PFAM) and trained the predictor with nearest network algorithm (NNA); at the second stage, we encoded nonapeptides (peptides that contain nine residues) of the protein containing the lysine residues, with Amino Acid Index, and trained the predictor with NNA. The predictor was tested by the k-fold cross-validation method. The highest accuracy of the second-staged predictor was 99.55% when 12 features were incorporated in the predictor. The corresponding Matthews Correlation Coefficient was 0.7952. These results indicate that the method is a promising tool to predict the sumoylation site of a protein. At last, the features used in the predictor are discussed. The software is available at request.

Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes

The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics, ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.