Molecular dynamics simulations reveal structural differences among wild-type NPC1 protein and its mutant forms

Abstract

NPC1 is a 25-exon gene located on the long arm of chromosome 18q11.2 and encodes NPC1, a transmembrane protein comprising 1278 amino acid residues. Mutations in the NPC1 gene can cause Niemann-Pick disease type C (NP-C), a rare autosomal-recessive neurovisceral disease. We assessed mutant protein folding using computer-based molecular dynamics (MD) simulations and molecular docking of the three most common NPC1 mutations, all of which result in changes in a cysteine-rich luminal loop region of the protein: a) I1061T is the most commonly detected variant in patients with NP-C worldwide; b) P1007A is the second most common variant, frequently detected in Portuguese, British and German patients; c) G992W occurs most often in patients of Acadian descent. Analyses of molecular structural information and related cellular physiological processes revealed that mutant NPC1 proteins exhibited altered function despite being far from the N-terminal domain cholesterol binding. MD simulations revealed that mutant I1061T protein shows remarkable instability in comparison the WT and also de other mutants, and interestingly this mutant has been identified as the most common variant. In the case of the mutant P1007A, it is presumed that this substitution promotes larger structural changes than proline due to their greater hydrophobic properties.

Structural changes related to the G992W mutation may affect the physicochemical space of G992W variant protein because tryptophan induces hydrophobic interactions. Cholesterol docking studies focused on binding recognition showed differences in the binding positions of variants versus the wild-type protein that go some way to explaining the molecular pathogenesis.

Communicated by Ramaswamy H. Sarma     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

In silico studying of the whole protein structure and dynamics of Dickkopf family members showed that N-terminal domain of Dickkopf 2 in contrary to other Dickkopfs facilitates its interaction with low density lipoprotein receptor related protein 5/6

Wnt (Wingless Int) signaling pathway has been known to be dysregulated in several human cancers, especially colorectal cancer (CRC). The Dickkopf (DKK) family which consists of four secreted proteins in vertebrates (DKK 1, 2, 3, 4) is one of the most critical antagonist families for Wnt signaling pathway. They typically antagonize Wnt/β-catenin signaling by binding and inhibiting Wnt co-receptors, LRP5/6 (low density lipoprotein receptor related protein 5/6). However, except for DKK1 (Dickkopf 1), details about structure and function of the members of this family are poorly defined. In this study, main Dickkopf family members were analyzed structurally, using protein structure prediction tools, molecular dynamics (MD), molecular docking and energy analyses. Three dimensional structure of whole DKKs was predicted and their interaction with LRP6 was investigated in detail. The results indicated that in DKK family members, a considerable diversity, in the case of structure, activity and physicochemical properties was seen. This diversity was more profound in DKK3 (Dickkopf3). Interestingly, the interaction mode of DKK2 (Dickkopf2) with its receptor, LRP6, was shown to be substantially different from other Dickkopf family members while N-terminal region of this ligand was also involved in the binding to the LRP6-P3P4. Moreover, the cysteine-rich domain 2 (CRD2) of DKK1 and DKK3 had a higher binding affinity to LRP6 in comparison with the whole protein structures.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics simulations of wild type and mutant of Pin1 peptidyl-prolyl isomerase

Pin1 catalyses the intrinsically slow process of cis-trans isomerisation and has been identified as a possible drug target in many diseases. Recently, the wild type (WT) and the Cys113Asp mutant of the Pin1 peptidyl-prolyl isomerase (PPIase) domain were determined by nuclear magnetic resonance. In this article, the WT and Cys113Asp mutant of PPIase domain are studied by molecular dynamics simulations. The structural stability analysis shows that the Cys113Asp mutation leads to the higher fluctuation of hydrophobic core in PPIase domain. The intrinsic correlated motions are important for the catalytic function of Pin1, whereas the Cys113Asp mutant system loses pivotal dynamical properties and develops wider conformational states than those in WT system. The intramolecular hydrogen bonds play crucial roles in the structural stability of PPIase domain. The mutated residue Asp113 attracts the side chain of His59 in the Cys113Asp system, which unbalances the internal interactions inside the catalytic tetrad. Meanwhile, the conformational changes of PPIase domain affect the side chain orientations of Lys63 and Arg69, which limit their binding with substrates. The Cys113Asp mutation destabilises the whole binding region of Pin1 PPIase domain, so the catalysis activity is severely reduced. These results are consistent with experimental studies and may help to understand the isomerisation mechanisms of Pin1.     

Molecular dynamics simulations and free energy analyses on the dimer formation of an amyloidogenic heptapeptide from human beta2-microglobulin: implication for the protofibril structure

Amyloid formation is associated with many neurodegenerative diseases. Recent findings suggest that early oligomeric aggregates could be major sources of toxicity. We present a computational investigation of the first step of amyloid initiation-dimer formation of a seven residue peptide (NHVTLSQ) from human beta2-microglobulin at pH 2.0, which renders +2.0 units charges to each peptide. A total of over 1.2 micros of simulations with explicit solvent and 1.0 micros of simulations with implicit solvent were conducted. Main-chain conformational restraint was applied to facilitate the formation of ordered dimers. An antiparallel beta-sheet with six main-chain hydrogen bonds was dominant in the implicit solvent simulations. In contrast, no stable dimers were observed in the two negative controls, the mouse heptapeptide (KHDSMAE, +3.0 units charges) and the scrambled human heptapeptide (QVLHTSN). Explicit solvent simulations presented a more complex scenario. The wild-type human heptapeptide formed predominantly antiparallel beta-sheets ( approximately 38%) although parallel ones ( approximately 12%) were also observed. Hydrophobic contacts preceded hydrogen bond saturation in the majority of the association events in the explicit solvent simulations, highlighting the important role of hydrophobic interaction in amyloid initiation. The fact that the mouse dimer dissociated immediately after the removal of conformational restraint suggests that the higher conformational entropy barrier, along with the stronger charge repulsion and weaker hydrophobic interaction, contributed to its inability to form amyloid fibril. The closeness of positive charge pairs in the dimers of the scrambled human heptapeptide may prohibit further beta-sheet extension and fibril growth. Combining the results from simulations and free energy analyses, we propose that the building block for this amyloid fibril is an antiparallel dimer with a two-residue register shift and six main-chain hydrogen bonds. A double-layer protofibril structure is also proposed in which two antiparallel beta-sheets face each other and are held together by hydrophobic staples and hydrogen bonds of the polar side-chains.     

Reconstruction of the src-SH3 protein domain transition state ensemble using multiscale molecular dynamics simulations

We use an integrated computational approach to reconstruct accurately the transition state ensemble (TSE) for folding of the src-SH3 protein domain. We first identify putative TSE conformations from free energy surfaces generated by importance sampling molecular dynamics for a fully atomic, solvated model of the src-SH3 protein domain. These putative TSE conformations are then subjected to a folding analysis using a coarse-grained representation of the protein and rapid discrete molecular dynamics simulations. Those conformations that fold to the native conformation with a probability (P(fold)) of approximately 0.5, constitute the true transition state. Approximately 20% of the putative TSE structures were found to have a P(fold) near 0.5, indicating that, although correct TSE conformations are populated at the free energy barrier, there is a critical need to refine this ensemble. Our simulations indicate that the true TSE conformations are compact, with a well-defined central beta sheet, in good agreement with previous experimental and theoretical studies. A structured central beta sheet was found to be present in a number of pre-TSE conformations, however, indicating that this element, although required in the transition state, does not define it uniquely. An additional tight cluster of contacts between highly conserved residues belonging to the diverging turn and second beta-sheet of the protein emerged as being critical elements of the folding nucleus. A number of commonly used order parameters to identify the transition state for folding were investigated, with the number of native Cbeta contacts displaying the most satisfactory correlation with P(fold) values.     

Molecular dynamics simulation,binding free energy calculation and molecular docking of human D-amino acid oxidase (DAAO) with its inhibitors

Schizophrenia is a mental illness; most affected people live in developing countries, and neither appropriate treatment nor commercial drugs are currently available. One possibility is to inhibit human-d-amino acid oxidase (h-DAAO). In this study, molecular dynamic simulations of the monomer, dimer and tetramer forms of h-DAAO complexed with the inhibitor 3-hydroxyquinolin-2(1H)-one(2) were performed. Seven residues, Leu51, Gln53, Leu215, Tyr228, Ile230, Arg283 and Gly313, were identified as essential for interacting with the inhibitor. Molecular docking of h-DAAO with pyrrole, quinoline and kojic acid derivatives, representing 69 known or potential h-DAAO inhibitors, was also performed. The results indicated that the activity of the inhibitor can be improved by modifying the compounds to have a substituent group capable of interacting with the side chain of Tyr228. Van der Waals interactions of the inhibitor with the hydrophobic pocket of h-DAAO and electrostatic interactions or H-bonds with Arg283 and Gly313 were important elements in determining the efficiency of the inhibitor. These results provide information on the interaction between h-DAAO and its inhibitors at the molecular level and can aid in the design of novel inhibitors against h-DAAO for new drug development in the treatment of schizophrenia.     

Molecular dynamics simulations of the thermal stability of tteRBP and ecRBP

Molecular dynamics simulations were performed for investigating the thermal stability of the extremely thermophilic Thermoanaerobacter tengcongensis ribose binding protein (tteRBP) and the mesophilic homologous Escherichia coli ribose binding protein (ecRBP). The simulations for the two proteins were carried out under the room temperature (300?K) and the optimal activity temperature (tteRBP 375?K and ecRBP 329?K), respectively. The comparative analyses of the trajectories show that the two proteins have stable overall structures at the two temperatures; further analyses indicate that they both have strong side-chain interactions and different backbone flexibilities at the different temperatures. The tteRBP 375?K and ecRBP 329?K have stronger internal motion and higher flexibility than tteRBP 300?K and ecRBP 300?K, respectively, it is noted that the flexibility of tteRBP is much higher than that of ecRBP at the two temperatures. Therefore, tteRBP 375?K can adapt to high temperature due to its higher flexibility of backbone. Combining with the researches by Cuneo et al., it is concluded that the side-chain interactions and flexibility of backbone are both the key factors to maintain thermal stability of the two proteins.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:22     

The structure of the major transition state for folding of an FF domain from experiment and simulation

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.     

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.     

Molecular dynamics simulations of transitions for ECD epidermal growth factor receptors show key differences between human and drosophila forms of the receptors

Recent X‐ray structural work on the Drosophila epidermal growth factor receptor (EFGR) has suggested an asymmetric dimer that rationalizes binding affinity measurements that go back decades (Alvarado et al., Cell 2010;142:568–579; Dawson et al., Structure 2007;15:942–954; Lemmon et al., Embo J 1997;16:281–294; Mattoon et al., Proc Natl Acad Sci USA 2004;101:923–928; Mayawala et al., Febs Lett 2005;579:3043–3047; Ozcan et al., Proc Natl Acad Sci USA 2006;103:5735–5740). This type of asymmetric structure has not been seen for the human EGF receptor family and it may or may not be important for function in that realm. We hypothesize that conformational changes in the Drosophila system have been optimized for the transition, whereas the barrier for the same transition is much higher in the human forms. To address our hypothesis we perform dynamic importance sampling (DIMS) (Perilla et al., J Comput Chem 2010;32:196–209) for barrier crossing transitions in both Drosophila and human EFGRs. For each set of transitions, we work from the hypothesis, based on results from the AdK system, that salt‐bridge pairs making and breaking connections are central to the conformational change. To evaluate the effectiveness of the salt‐bridges as drivers for the conformational change, we use the effective transfer entropy based on stable state MD calculations (Kamberaj and Der Vaart, Biophys J 2009;97:1747–1755) to define a reduced subset of degrees of freedom that seem to be important for driving the transition (Perilla and Woolf, J Chem Phys 2012;136:164101). Our results suggest that salt‐bridge making and breaking is not the dominant factor in driving the symmetric to asymmetric transition, but that instead it is a result of more concerted and correlated functional motions within a subset of the dimer structures. Furthermore, the analysis suggests that the set of residues involved in the transitions from the Drosophila relative to the human forms differs and that this difference in substate distributions relates to why the asymmetric form may be more common to Drosophila than to the human forms. We close with a discussion about the residues that may be changed in the human and the Drosophila forms to potentially shift the kinetics of the symmetric to asymmetric transition. Proteins 2013; 81:1113–1126. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulation and automated docking of the pro-apoptotic Bax protein and its complex with a peptide designed from the Bax-binding domain of anti-apoptotic Ku70

Bax, a multi-domain protein belonging to the large family of Bcl-2 proteins, has a pivotal role for the initiation of the cytochrome c-mediated apoptosis, a vital physiologic process to eliminate damaged or unwanted cells. In response to specific stimuli Bax translocates from cytosol to mitochondria outer membrane where a process of oligomerization occurs with pore formation through which cytochrome c and other death molecules escape. The pro-death action of Bax is regulated by the interaction with other pro-survival proteins. However, the conformational changes and the structural details necessary for homo and hetero interaction with other regulating proteins are largely unknown. This article reports a combined investigation of molecular dynamics (MD) simulation and automated docking that evidence the molecular regions of Bax involved in the binding with anti-apoptotic exapeptide (Bip) designed from Ku70, a subunit of the protein complex essential for non-homologous DNA repair but that inhibits also the Bax translocation to mitochondria. Since Bip suppresses apoptosis induced by several anti-cancer drugs, it appears relevant to achieve a better understanding of the structural and dynamical aspects that characterize the Bip-Bax complex in view of potential therapeutic implications. The present results show that the Bax region with the highest affinity for Bip is located in proximity of BH3 homology domain of Bax and also involves the alpha-helices 1 and 8. Moreover, the comparison of essential motions of Bax at 300 and 400 K before and after the formation of the complex with Bip evidences how the binding with the exa-peptide affects the collective motions of specific molecular districts of Bax considered to have functional relevance.     

Molecular dynamics simulations and functional characterization of the interactions of the PAR2 ectodomain with factor VIIa

Signaling of the tissue factor‐FVIIa complex regulates angiogenesis, tumor growth, and inflammation. TF‐FVIIa triggers cell signaling events by cleavage of protease activated receptor (PAR2) at the Arg36‐Ser37 scissile bond. The recognition of PAR2 by the FVIIa protease domain is poorly understood. We perform molecular modeling and dynamics simulations to derive the PAR2‐FVIIa interactions. Docking of the PAR2 Arg36‐Ser37 scissile bond to the S1 site and subsequent molecular dynamics leads to interactions of the PAR2 ectodomain with P and P′ sites of the FVIIa catalytic cleft as well as to electrostatic interactions between a stably folded region of PAR2 and a cluster of basic residues remote from the catalytic cleft of FVIIa. To address the functional significance of this interaction for PAR2 cleavage, we employed two antibodies with epitopes previously mapped to this cluster of basic residues. Although these antibodies do not block the catalytic cleft, both antibodies completely abrogated PAR2 activation by TF‐FVIIa. Our simulations indicate a conformation of the PAR2 ectodomain that limits the cleavage site to no more than 33 Å from its membrane proximal residue. Since the active site of FVIIa in the TF‐FVIIa complex is ～75 Å above the membrane, cleavage of the folded conformation of PAR2 would require tilting of the TF‐FVIIa complex toward the membrane, indicating that additional cellular factors may be required to properly align the scissile bond of PAR2 with TF‐FVIIa. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Effects of PI(4,5)P2 concentration on the F-BAR domain membrane binding as revealed by coarse-grained simulations

Bin/Amphyphysin/Rvs (BAR) domain proteins form a key link between membrane remodeling and cytoskeleton dynamics. They are dimers that bind to membranes via electrostatic interactions with different preferences toward negatively charged lipids. In the present article, we examine the interactions of the F-BAR domain of nervous wreck (Nwk) with phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂)-containing membranes using coarse-grained molecular dynamics. We demonstrated PI(4,5)P₂ concentration effects, identified the sequence of events that underlies the protein binding and identified amino acids involved in protein–lipid interactions. Our simulations point out the primary role of the basic stretch at the tips of the dimer, which anchors the protein to the membrane and initiates the binding process. When the PI(4,5)P₂ concentration is high, the protein stably associates with the membrane by its concave surface or by the opposite side. At low PI(4,5)P₂ concentration, the former orientation becomes more favorable; also a state with only one tip bound is observed, due to the weaker attachment and more pronounced association/dissociation events. Our results provide a theoretical model that describes the lipid-binding behavior of Nwk observed in vitro.     

Molecular dynamics simulations of peptides containing an unnatural amino acid: dimerization, folding, and protein binding

We have performed molecular dynamics (MD) simulations to study the dimerization, folding, and binding to a protein of peptides containing an unnatural amino acid. NMR studies have shown that the substitution of one residue in a tripeptide beta-strand by the unnatural amino acid Hao (5-HO2CCONH-2-MeO-C6H3-CO-NHNH2) modifies the conformational flexibility of the beta-strand and the hydrogen-bonding properties of its two edges: The number of hydrogen-bond donors and acceptors increases at one edge, whereas at the other, they are sterically hindered. In simulations in chloroform, the Hao-containing peptide 9 (i-PrCO-Phe-Hao-Val-NHBu) forms a beta-sheet-like hydrogen-bonded dimer, in good agreement with the available experimental data. Addition of methanol to the solution induces instability of this beta-sheet, as confirmed by the experiments. MD simulations also reproduce the folding of the synthetic peptide 1a (i-PrCO-Hao-Ut-Phe-Ile-Leu-NHMe) into a beta-hairpin-like structure in chloroform. Finally, the Hao-containing peptide, Ac-Ala-Hao-Ala-NHMe, is shown to form a stable complex with the Ras analogue, Rap1 A, in water at room temperature. Together with the available experimental data, these simulation studies indicate that Hao-containing peptides may serve as inhibitors of beta-sheet interactions between proteins.     

Effect of altered glycosylation on the structure of the I-like domain of beta1 integrin: a molecular dynamics study

Glycosylation plays an important role in the regulation of integrin function. Molecular mechanisms underlying the effects of altered glycosylation on beta1 integrin structure and function are still largely unknown. In this study, we used a molecular modeling approach to study the effects of altered glycosylation, with alpha2-6 sialic acid and without alpha2-6 sialic acid, on the structure of the I-like domain of the beta1 integrin. Our results demonstrated that altered glycosylation affected the interactions between oligosaccharides and the I-like domain, which in turn changed the accessibility of the specificity-determining loop for ligand binding. Altered glycosylation caused significant conformational changes for most of the key functional regions of the I-like domain of beta1 integrin, including the metal ion-dependent adhesion site that contains a DLSYS motif, and other critical residues for ligand binding (Asn-224, Glu-229, Asp-233, Asp-267, and Asp-295). In addition, altered glycosylation caused significant movement of the alpha1 and alpha7 helices, which are important for the activation of beta1 integrin. The results from this study offered molecular mechanisms for the experimental observations that variant glycosylation regulates integrin function.     

Characterization of differences in substrate specificity among CYP1A1, CYP1A2 and CYP1B1: an integrated approach employing molecular docking and molecular dynamics simulations

Recent trends in new drug discovery of anticancer drugs have made oncologists more aware of the fact that the new drug discovery must target the developing mechanism of tumorigenesis to improve the therapeutic efficacy of antineoplastic drugs. The drugs designed are expected to have high affinity towards the novel targets selectively. Current research highlights overexpression of CYP450s, particularly cytochrome P450 1A1 (CYP1A1), in tumour cells, representing a novel target for anticancer therapy. However, the CYP1 family is identified as posing significant problems in selectivity of anticancer molecules towards CYP1A1. Three members have been identified in the human CYP1 family: CYP1A1, CYP1A2 and CYP1B1. Although sequences of the three isoform have high sequence identity, they have distinct substrate specificities. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics, design novel antitumour compounds that could be specifically metabolized by only CYP1A1 to mediate their antitumour activity and elucidate the reasons for differences in substrate specificity profile among the three proteins. In the present study, we employed a combination of computational methodologies: molecular docking and molecular dynamics simulations. We utilized eight substrates for elucidating the difference in substrate specificity of the three isoforms. Lastly, we conclude that the substrate specificity of a particular substrate depends upon the type of the active site residues, the dynamic motions in the protein structure upon ligand binding and the physico‐chemical characteristics of a particular ligand. Copyright © 2016 John Wiley & Sons, Ltd.