Exploring the multiscale signaling behavior of phototropin1 from Chlamydomonas reinhardtii using a full‐residue space kinetic Monte Carlo molecular dynamics technique

Devising analysis tools for elucidating the regulatory mechanism of complex enzymes has been a challenging task for many decades. It generally requires the determination of the structural‐dynamical information of protein solvent systems far from equilibrium over multiple length and time scales, which is still difficult both theoretically and experimentally. To cope with the problem, we introduce a full‐residue space multiscale simulation method based on a combination of the kinetic Monte Carlo and molecular dynamics techniques, in which the rates of the rate‐determining processes are evaluated from a biomolecular forcefield on the fly during the simulation run by taking into account the full space of residues. To demonstrate its reliability and efficiency, we explore the light‐induced functional behavior of the full‐length phototropin1 from Chlamydomonas reinhardtii (Cr‐phot1) and its various subdomains. Our results demonstrate that in the dark state the light oxygen voltage‐2‐Jα (LOV2‐Jα) photoswitch inhibits the enzymatic activity of the kinase, whereas the LOV1‐Jα photoswitch controls the dimerization with the LOV2 domain. This leads to the repulsion of the LOV1‐LOV2 linker out of the interface region between both LOV domains, which results in a positively charged surface suitable for cell–membrane interaction. By contrast, in the light state, we observe that the distance between both LOV domains is increased and the LOV1‐LOV2 linker forms a helix–turn–helix (HTH) motif, which enables gene control through nucleotide binding. Finally, we find that the kinase is activated through the disruption of the Jα‐helix from the LOV2 domain, which is followed by a stretching of the activation loop (A‐loop) and broadening of the catalytic cleft of the kinase. Proteins 2014; 82:2018–2040. © 2014 Wiley Periodicals, Inc.     

Identification of potential allosteric binding sites in cathepsin K based on intramolecular communication

Network theory methods and molecular dynamics (MD) simulations are accepted tools to study allosteric regulation. Indeed, dynamic networks built upon correlation analysis of MD trajectories provide detailed information about communication paths between distant sites. In this context, we aimed to understand whether the efficiency of intramolecular communication could be used to predict the allosteric potential of a given site. To this end, we performed MD simulations and network theory analyses in cathepsin K (catK), whose allosteric sites are well defined. To obtain a quantitative measure of the efficiency of communication, we designed a new protocol that enables the comparison between properties related to ensembles of communication paths obtained from different sites. Further, we applied our strategy to evaluate the allosteric potential of different catK cavities not yet considered for drug design. Our predictions of the allosteric potential based on intramolecular communication correlate well with previous catK experimental and theoretical data. We also discuss the possibility of applying our approach to other proteins from the same family.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

14-3-3蛋白研究进展

14-3-3蛋白是高度保守的、所有真核生物细胞中都普遍存在的、在大多数生物物种中由一个基因家族编码的一类蛋白调控家族。它几乎参与生命体所有的生理反应过程,人们在各种组织细胞中发现了各种不同的14-3-3蛋白。作为与磷酸丝氨酸／苏氨酸结合的第一信号分子,14-3-3蛋白在细胞的信号转导中起着至关重要的作用,尤其是它直接参与调节蛋白激酶和蛋白磷酸化酶的活性,被称为蛋白质与蛋白质相互作用的”桥梁蛋白”;它可以与转录因子结合形成复合体,调节相关基因的表达。一些研究表明,14-3-3蛋白调控机制的紊乱可以直接导致疾病的发生,在临床上14-3-3蛋白常常可以作为诊断的标志物。     

The interplay between the clustering of transmembrane proteins and coupling of anchored membrane proteins

Clustering of membrane proteins plays an important role in many cellular activities such as protein sorting and signal transduction. In this study, we used dissipative particle dynamics simulation method to investigate the clustering of anchored membrane proteins (AMPs) in the presence of transmembrane proteins (TMPs). First, our simulation results show that clustering of AMPs and that of TMPs are in fact interdependent, and depending on their hydrophobic length, both protein mixing and protein demixing are observed. Especially, the protein demixing occurs only when the hydrophobic mismatch of TMPs is negative while that of AMPs is positive. Second, our simulation results indicate that the clustering of TMPs also modulates the coupling of the clustering of AMPs in both leaflets. On the one hand, the coupling between AMPs in different leaflets will be strongly restrained if TMPs form protein mixing with AMPs in one leaflet and protein demixing with AMPs in the other leaflet. On the other hand, the coupling between AMPs can be enhanced or mediated by TMPs when TMPs mix with AMPs in both leaflets. Our results may have some implications on our understanding of how different types of membrane proteins cluster and provide a possible explanation of how TMPs participate in signal transduction across cellular membranes.     

Seasonal patterns of dehydrins and 70-kDa heat-shock proteins in bark tissues of eight species of woody plants

Although considerable effort has been directed at identifying and understanding the function and regulation of stress-induced proteins in herbaceous plants, reports concerning woody plants are limited. Studies with herbaceous crops have revealed similarities in the types of proteins that accumulate in response to a wide array of abiotic stresses and hormonal cues such as the accumulation of abscisic acid. Many of the identified proteins appear to be related to dehydrins (the D-11 subgroup of late-embryogenesis-abundant proteins). The objective of the present study was to determine if seasonal induction of dehydrins is a common feature in woody plants and to see if seasonal patterns existed for other stress-induced proteins. Bark tissues from eight species of woody plants were collected monthly for a period of 1.5 years. The species included: peach (Prunus persica) cv. Loring; apple (Malus domestica) cv. Golden Delicious; thornless blackberry (Rubus sp.) cv. Chester; hybrid poplar (Populus nigra); weeping willow (Salix babylonica); flowering dogwood (Cornus florida); sassafras (Sassafras albidum); and black locust (Robinia pseudo-acacia). Immunoblots of bark proteins were probed with a polyclonal antibody recognizing a conserved region of dehydrin proteins, and monoclonal antibodies directed against members of the HS70 family of heat-shock proteins. Some proteins, immunologically related to dehydrins, appeared to be constitutive; however, distinct seasonal patterns associated with winter acclimation were also observed in all species. The molecular masses of these proteins varied widely, although similarities were observed in related species (willow and poplar). Identification of proteins using the monoclonal antibodies (HSP70, HSC70, BiP) was more definitive because of their inherent specificity, but seasonal patterns were more variable among the eight species examined. This study represents only a precursory examination of several proteins reported to be stress related in herbaceous plants, but the results indicate that these proteins are also common to woody plants and that further research to characterize their regulation and function in relation to stress adaptation and the perennial life cycle of woody plants is warranted.     

Reduction of protein disulfide isomerase results in open conformations and stimulates dynamic exchange between structural ensembles

Human protein disulfide isomerase (PDI) is an essential redox-regulated enzyme required for oxidative protein folding. It comprises four thioredoxin domains, two catalytically active (a, a’) and two inactive (b, b’), organized to form a flexible abb’a’ U-shape. Snapshots of unbound oxidized and reduced PDI have been obtained by X-ray crystallography. Yet, how PDI’s structure changes in response to the redox environment and inhibitor binding remains controversial. Here, we used multiparameter confocal single-molecule FRET to track the movements of the two catalytic domains with high temporal resolution. We found that at equilibrium, PDI visits three structurally distinct conformational ensembles, two “open” (O₁ and O₂) and one “closed” (C). We show that the redox environment dictates the time spent in each ensemble and the rate at which they exchange. While oxidized PDI samples O₁, O₂, and C more evenly and in a slower fashion, reduced PDI predominantly populates O₁ and O₂ and exchanges between them more rapidly, on the submillisecond timescale. These findings were not expected based on crystallographic data. Using mutational analyses, we further demonstrate that the R300-W396 cation-π interaction and active site cysteines dictate, in unexpected ways, how the catalytic domains relocate. Finally, we show that irreversible inhibitors targeting the active sites of reduced PDI did not abolish these protein dynamics but rather shifted the equilibrium toward the closed ensemble. This work introduces a new structural framework that challenges current views of PDI dynamics, helps rationalize its multifaceted role in biology, and should be considered when designing PDI-targeted therapeutics.     

A universal allosteric mechanism for G protein activation

1. Download : Download high-res image (229KB)
2. Download : Download full-size image

     

Prediction of helix-helix contacts and interacting helices in polytopic membrane proteins using neural networks

Despite rapidly increasing numbers of available 3D structures, membrane proteins still account for less than 1% of all structures in the Protein Data Bank. Recent high-resolution structures indicate a clearly broader structural diversity of membrane proteins than initially anticipated, motivating the development of reliable structure prediction methods specifically tailored for this class of molecules. One important prediction target capturing all major aspects of a protein's 3D structure is its contact map. Our analysis shows that computational methods trained to predict residue contacts in globular proteins perform poorly when applied to membrane proteins. We have recently published a method to identify interacting alpha-helices in membrane proteins based on the analysis of coevolving residues in predicted transmembrane regions. Here, we present a substantially improved algorithm for the same problem, which uses a newly developed neural network approach to predict helix-helix contacts. In addition to the input features commonly used for contact prediction of soluble proteins, such as windowed residue profiles and residue distance in the sequence, our network also incorporates features that apply to membrane proteins only, such as residue position within the transmembrane segment and its orientation toward the lipophilic environment. The obtained neural network can predict contacts between residues in transmembrane segments with nearly 26% accuracy. It is therefore the first published contact predictor developed specifically for membrane proteins performing with equal accuracy to state-of-the-art contact predictors available for soluble proteins. The predicted helix-helix contacts were employed in a second step to identify interacting helices. For our dataset consisting of 62 membrane proteins of solved structure, we gained an accuracy of 78.1%. Because the reliable prediction of helix interaction patterns is an important step in the classification and prediction of membrane protein folds, our method will be a helpful tool in compiling a structural census of membrane proteins.     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

Clarifying allosteric control of flap conformations in the 1TW7 crystal structure of HIV-1 protease

The 1TW7 crystal structure of HIV-1 protease shows the flaps placed wider and more open than what is seen in other examples of the semi-open, apo form. It has been proposed that this might be experimental evidence of allosteric control, because crystal packing creates contacts to the "elbow region" of the protease, which may cause deformation of the flaps. Recent dynamics simulations have shown that the conformation seen in 1TW7 relaxes into the typical semi-open conformation in the absence of the crystal contacts, definitively showing that the crystal contacts cause the deformation (Layten et al., J Am Chem Soc 2006;128:13360-13361). However, this does not prove or disprove allosteric modulation at the elbow. In this study, we have conducted additional simulations, supplemented with experimental testing, to further probe the possibility of 1TW7 providing an example of allosteric control of the flap region. We show that the contacts are unstable and do not restrict the conformational sampling of the flaps. The deformation seen in the 1TW7 crystal structure is simply opportunistic crystal packing and not allosteric control.     

Synergy and antagonism between Notch and BMP receptor signaling pathways in endothelial cells

Notch and bone morphogenetic protein signaling pathways are important for cellular differentiation, and both have been implicated in vascular development. In many cases the two pathways act similarly, but antagonistic effects have also been reported. The underlying mechanisms and whether this is caused by an interplay between Notch and BMP signaling is unknown. Here we report that expression of the Notch target gene, Herp2, is synergistically induced upon activation of Notch and BMP receptor signaling pathways in endothelial cells. The synergy is mediated via RBP-Jkappa/CBF-1 and GC-rich palindromic sites in the Herp2 promoter, as well as via interactions between the Notch intracellular domain and Smad that are stabilized by p/CAF. Activated Notch and its downstream effector Herp2 were found to inhibit endothelial cell (EC) migration. In contrast, BMP via upregulation of Id1 expression has been reported to promote EC migration. Interestingly, Herp2 was found to antagonize BMP receptor/Id1-induced migration by inhibiting Id1 expression. Our results support the notion that Herp2 functions as a critical switch downstream of Notch and BMP receptor signaling pathways in ECs.     

Conformational equilibria in allosteric control of Hsp70 chaperones

1. Download : Download high-res image (219KB)
2. Download : Download full-size image

     

Reaction path and free energy calculations of the transition between alternate conformations of HIV-1 protease

Two different structures of ligand-free HIV protease have been determined by X-ray crystallography. These structures differ in the position of two 12 residue, β-hairpin regions (or “flaps”) which cap the active site. The movements of the flaps must be involved in the binding of substrates since, in either conformation, the flaps block the binding site. One of these structures is similar to structures of the ligand-bound enzyme; however, the importance of both structures to enzyme function is unclear. This transformation takes place on a time scale too long for conventional molecular dynamics simulations, so the process was studied by first identifying a reaction path between the two structures and then calculating the free energy along this path using umbrella sampling. For the ligand-free enzyme, it is found that the two structures are nearly equally stable, with the ligand-bound-type structure being less stable, consistent with X-ray crystallography data. The more stable open structure does not have a lower potential energy, but is stabilized by entropy. The transition occurs through a collapse and reformation of the β-sheet structure of the conformationally flexible, glycine-rich flap ends. Additionally, some problems in studying conformational changes in proteins through the use of a single reaction path are addressed. Proteins 32:7–16, 1998. © 1998 Wiley-Liss, Inc.     

Divergent requirements for fibroblast growth factor signaling in zebrafish maxillary barbel and caudal fin regeneration

The zebrafish maxillary barbel is an integumentary organ containing skin, glands, pigment cells, taste buds, nerves, and endothelial vessels. The maxillary barbel can regenerate (LeClair & Topczewski 2010); however, little is known about its molecular regulation. We have studied fibroblast growth factor (FGF) pathway molecules during barbel regeneration, comparing this system to a well‐known regenerating appendage, the zebrafish caudal fin. Multiple FGF ligands (fgf20a, fgf24), receptors (fgfr1‐4) and downstream targets (pea3, il17d) are expressed in normal and regenerating barbel tissue, confirming FGF activation. To test if specific FGF pathways were required for barbel regeneration, we performed simultaneous barbel and caudal fin amputations in two temperature‐dependent zebrafish lines. Zebrafish homozygous for a point mutation in fgf20a, a factor essential for caudal fin blastema formation, regrew maxillary barbels normally, indicating that the requirement for this ligand is appendage‐specific. Global overexpression of a dominant negative FGF receptor, Tg(hsp70l:dn‐fgfr1:EGFP)^pd1 completely blocked fin outgrowth but only partially inhibited barbel outgrowth, suggesting reduced requirements for FGFs in barbel tissue. Maxillary barbels expressing dn‐fgfr1 regenerated peripheral nerves, dermal connective tissue, endothelial tubes, and a glandular epithelium; in contrast to a recent report in which dn‐fgfr1 overexpression blocks pharyngeal taste bud formation in zebrafish larvae (Kapsimali et al. 2011), we observed robust formation of calretinin‐positive tastebuds. These are the first experiments to explore the molecular mechanisms of maxillary barbel regeneration. Our results suggest heterogeneous requirements for FGF signaling in the regeneration of different zebrafish appendages (caudal fin versus maxillary barbel) and taste buds of different embryonic origin (pharyngeal endoderm versus barbel ectoderm).     

Sorghum bicolor SbHSP110 has an elongated shape and is able of protecting against aggregation and replacing human HSPH1/HSP110 in refolding and disaggregation assays

Perturbations in the native structure, often caused by stressing cellular conditions, not only impair protein function but also lead to the formation of aggregates, which can accumulate in the cell leading to harmful effects. Some organisms, such as plants, express the molecular chaperone HSP100 (homologous to HSP104 from yeast), which has the remarkable capacity to disaggregate and reactivate proteins. Recently, studies with animal cells, which lack a canonical HSP100, have identified the involvement of a distinct system composed of HSP70/HSP40 that needs the assistance of HSP110 to efficiently perform protein breakdown. As sessile plants experience stressful conditions more severe than those experienced by animals, we asked whether a plant HSP110 could also play a role in collaborating with HSP70/HSP40 in a system that increases the efficiency of disaggregation. Thus, the gene for a putative HSP110 from the cereal Sorghum bicolor was cloned and the protein, named SbHSP110, purified. For comparison purposes, human HsHSP110 (HSPH1/HSP105) was also purified and investigated in parallel. First, a combination of spectroscopic and hydrodynamic techniques was used for the characterization of the conformation and stability of recombinant SbHSP110, which was produced folded. Second, small-angle X-ray scattering and combined predictors of protein structure indicated that SbHSP110 and HsHSP110 have similar conformations. Then, the chaperone activities, which included protection against aggregation, refolding, and reactivation, were investigated, showing that SbHSP110 and HsHSP110 have similar functional activities. Altogether, the results add to the structure/function relationship study of HSP110s and support the hypothesis that plants have multiple strategies to act upon the reactivation of protein aggregates.     

Equilibrium transitions between side‐chain conformations in leucine and isoleucine

Despite recent improvements in computational methods for protein design, we still lack a quantitative, predictive understanding of the intrinsic probabilities for amino acids to adopt particular side‐chain conformations. Surprisingly, this question has remained unsettled for many years, in part because of inconsistent results from different experimental approaches. To explicitly determine the relative populations of different side‐chain dihedral angles, we performed all‐atom hard‐sphere Langevin Dynamics simulations of leucine (Leu) and isoleucine (Ile) dipeptide mimetics with stereo‐chemical constraints and repulsive‐only steric interactions between non‐bonded atoms. We determine the relative populations of the different χ₁ and χ₂ dihedral angle combinations as a function of the backbone dihedral angles ? and ψ. We also propose, and test, a mechanism for inter‐conversion between the different side‐chain conformations. Specifically, we discover that some of the transitions between side‐chain dihedral angle combinations are very frequent, whereas others are orders of magnitude less frequent, because they require rare coordinated motions to avoid steric clashes. For example, to transition between different values of χ₂, the Leu side‐chain bond angles κ₁ and κ₂ must increase, whereas to transition in χ₁, the Ile bond angles λ₁ and λ₂ must increase. These results emphasize the importance of computational approaches in stimulating further experimental studies of the conformations of side‐chains in proteins. Moreover, our studies emphasize the power of simple steric models to inform our understanding of protein structure, dynamics, and design. Proteins 2015; 83:1488–1499. © 2015 Wiley Periodicals, Inc.     

Protein-protein recognition: an experimental and computational study of the R89K mutation in Raf and its effect on Ras binding

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras-Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [gamma-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9+/-1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein-protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 A away, namely residues 37-41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2-3 specific water molecules. This analysis will provide insights into the large experimental database of Ras-Raf mutations.     

分子伴侣的功能和应用

本文综述了分子伴侣的分类、功能、作用机理、研究现状及应用前景。分子伴侣是在生物大分子的折叠、组装、转运及降解等过程中起协助作用,参与协助抗原的呈递和遗传物质的复制、转录及构象的确立,但自身并不发生任何变化的一大类广泛存在于生物体内的蛋白质分子。随着对分子伴侣的进一步研究和相关知识的不断深入,分子伴侣在生物产品开发、物种改良、抗衰老,疾病预防、诊断和治疗以及环境监测方面具有广阔的前景。     

Investigation of the molecular similarity in closely related protein systems: The PrP case study

The amyloid conversion is a massive detrimental modification affecting several proteins upon specific physical or chemical stimuli characterizing a plethora of diseases. In many cases, the amyloidogenic stimuli induce specific structural features to the protein conferring the propensity to misfold and form amyloid deposits. The investigation of mutants, structurally similar to their native isoform but inherently prone to amyloid conversion, may be a viable strategy to elucidate the structural features connected with amyloidogenesis. In this article, we present a computational protocol based on the combination of molecular dynamics (MD) and grid‐based approaches suited for the pairwise comparison of closely related protein structures. This method was applied on the cellular prion protein (PrP^C) as a case study and, in particular, addressed to the quali/quantification of the structural features conferred by either E200K mutations and treatment with CaCl₂, both able to induce the scrapie conversion of PrP. Several schemes of comparison were developed and applied to this case study, and made up suitable of application to other protein systems. At this purpose an in‐house python codes has been implemented that, together with the parallelization of the GRID force fields program, will spread the applicability of the proposed computational procedure. Proteins 2015; 83:1751–1765. © 2015 Wiley Periodicals, Inc.