Regulation of the Hsp104 Middle Domain Activity Is Critical for Yeast Prion Propagation

Molecular chaperones play a significant role in preventing protein misfolding and aggregation. Indeed, some protein conformational disorders have been linked to changes in the chaperone network. Curiously, in yeast, chaperones also play a role in promoting prion maintenance and propagation. While many amyloidogenic proteins are associated with disease in mammals, yeast prion proteins, and their ability to undergo conformational conversion into a prion state, are proposed to play a functional role in yeast biology. The chaperone Hsp104, a AAA+ ATPase, is essential for yeast prion propagation. Hsp104 fragments large prion aggregates to generate a population of smaller oligomers that can more readily convert soluble monomer and be transmitted to daughter cells. Here, we show that the middle (M) domain of Hsp104, and its mobility, plays an integral part in prion propagation. We generated and characterized mutations in the M-domain of Hsp104 that are predicted to stabilize either a repressed or de-repressed conformation of the M-domain (by analogy to ClpB in bacteria). We show that the predicted stabilization of the repressed conformation inhibits general chaperone activity. Mutation to the de-repressed conformation, however, has differential effects on ATP hydrolysis and disaggregation, suggesting that the M-domain is involved in coupling these two activities. Interestingly, we show that changes in the M-domain differentially affect the propagation of different variants of the [PSI+] and [RNQ+] prions, which indicates that some prion variants are more sensitive to changes in the M-domain mobility than others. Thus, we provide evidence that regulation of the M-domain of Hsp104 is critical for efficient prion propagation. This shows the importance of elucidating the function of the M-domain in order to understand the role of Hsp104 in the propagation of different prions and prion variants.     

Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants

Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.     

A Translocation Motif in Relaxase TrwC Specifically Affects Recruitment by Its Conjugative Type IV Secretion System

Type IV secretion system (T4SS) substrates are recruited through a translocation signal that is poorly defined for conjugative relaxases. The relaxase TrwC of plasmid R388 is translocated by its cognate conjugative T4SS, and it can also be translocated by the VirB/D4 T4SS of Bartonella henselae, causing DNA transfer to human cells. In this work, we constructed a series of TrwC variants and assayed them for DNA transfer to bacteria and human cells to compare recruitment requirements by both T4SSs. Comparison with other reported relaxase translocation signals allowed us to determine two putative translocation sequence (TS) motifs, TS1 and TS2. Mutations affecting TS1 drastically affected conjugation frequencies, while mutations affecting either motif had only a mild effect on DNA transfer rates through the VirB/D4 T4SS of B. henselae. These results indicate that a single substrate can be recruited by two different T4SSs through different signals. The C terminus affected DNA transfer rates through both T4SSs tested, but no specific sequence requirement was detected. The addition of a Bartonella intracellular delivery (BID) domain, the translocation signal for the Bartonella VirB/D4 T4SS, increased DNA transfer up to 4% of infected human cells, providing an excellent tool for DNA delivery to specific cell types. We show that the R388 coupling protein TrwB is also required for this high-efficiency TrwC-BID translocation. Other elements apart from the coupling protein may also be involved in substrate recognition by T4SSs.     

Signal propagation in feedforward neuronal networks with unreliable synapses

In this paper, we systematically investigate both the synfire propagation and firing rate propagation in feedforward neuronal network coupled in an all-to-all fashion. In contrast to most earlier work, where only reliable synaptic connections are considered, we mainly examine the effects of unreliable synapses on both types of neural activity propagation in this work. We first study networks composed of purely excitatory neurons. Our results show that both the successful transmission probability and excitatory synaptic strength largely influence the propagation of these two types of neural activities, and better tuning of these synaptic parameters makes the considered network support stable signal propagation. It is also found that noise has significant but different impacts on these two types of propagation. The additive Gaussian white noise has the tendency to reduce the precision of the synfire activity, whereas noise with appropriate intensity can enhance the performance of firing rate propagation. Further simulations indicate that the propagation dynamics of the considered neuronal network is not simply determined by the average amount of received neurotransmitter for each neuron in a time instant, but also largely influenced by the stochastic effect of neurotransmitter release. Second, we compare our results with those obtained in corresponding feedforward neuronal networks connected with reliable synapses but in a random coupling fashion. We confirm that some differences can be observed in these two different feedforward neuronal network models. Finally, we study the signal propagation in feedforward neuronal networks consisting of both excitatory and inhibitory neurons, and demonstrate that inhibition also plays an important role in signal propagation in the considered networks.     

Directed evolution to re-adapt a co-evolved network within an enzyme

We have previously used targeted active-site saturation mutagenesis to identify a number of transketolase single mutants that improved activity towards either glycolaldehyde (GA), or the non-natural substrate propionaldehyde (PA). Here, all attempts to recombine the singles into double mutants led to unexpected losses of specific activity towards both substrates. A typical trade-off occurred between soluble expression levels and specific activity for all single mutants, but many double mutants decreased both properties more severely suggesting a critical loss of protein stability or native folding. Statistical coupling analysis (SCA) of a large multiple sequence alignment revealed a network of nine co-evolved residues that affected all but one double mutant. Such networks maintain important functional properties such as activity, specificity, folding, stability, and solubility and may be rapidly disrupted by introducing one or more non-naturally occurring mutations. To identify variants of this network that would accept and improve upon our best D469 mutants for activity towards PA, we created a library of random single, double and triple mutants across seven of the co-evolved residues, combining our D469 variants with only naturally occurring mutations at the remaining sites. A triple mutant cluster at D469, E498 and R520 was found to behave synergistically for the specific activity towards PA. Protein expression was severely reduced by E498D and improved by R520Q, yet variants containing both mutations led to improved specific activity and enzyme expression, but with loss of solubility and the formation of inclusion bodies. D469S and R520Q combined synergistically to improve k_cat 20-fold for PA, more than for any previous transketolase mutant. R520Q also doubled the specific activity of the previously identified D469T to create our most active transketolase mutant to date. Our results show that recombining active-site mutants obtained by saturation mutagenesis can rapidly destabilise critical networks of co-evolved residues, whereas beneficial single mutants can be retained and improved upon by randomly recombining them with natural variants at other positions in the network.     

Functional Diversification of Hsp40: Distinct J-Protein Functional Requirements for Two Prions Allow for Chaperone-Dependent Prion Selection

Yeast prions are heritable amyloid aggregates of functional yeast proteins; their propagation to subsequent cell generations is dependent upon fragmentation of prion protein aggregates by molecular chaperone proteins. Mounting evidence indicates the J-protein Sis1 may act as an amyloid specificity factor, recognizing prion and other amyloid aggregates and enabling Ssa and Hsp104 to act in prion fragmentation. Chaperone interactions with prions, however, can be affected by variations in amyloid-core structure resulting in distinct prion variants or ‘strains’. Our genetic analysis revealed that Sis1 domain requirements by distinct variants of [PSI ⁺] are strongly dependent upon overall variant stability. Notably, multiple strong [PSI ⁺] variants can be maintained by a minimal construct of Sis1 consisting of only the J-domain and glycine/phenylalanine-rich (G/F) region that was previously shown to be sufficient for cell viability and [RNQ ⁺] prion propagation. In contrast, weak [PSI ⁺] variants are lost under the same conditions but maintained by the expression of an Sis1 construct that lacks only the G/F region and cannot support [RNQ ⁺] propagation, revealing mutually exclusive requirements for Sis1 function between these two prions. Prion loss is not due to [PSI ⁺]-dependent toxicity or dependent upon a particular yeast genetic background. These observations necessitate that Sis1 must have at least two distinct functional roles that individual prions differentially require for propagation and which are localized to the glycine-rich domains of the Sis1. Based on these distinctions, Sis1 plasmid-shuffling in a [PSI ⁺]/[RNQ ⁺] strain permitted J-protein-dependent prion selection for either prion. We also found that, despite an initial report to the contrary, the human homolog of Sis1, Hdj1, is capable of [PSI ⁺] prion propagation in place of Sis1. This conservation of function is also prion-variant dependent, indicating that only one of the two Sis1-prion functions may have been maintained in eukaryotic chaperone evolution.     

Discovery of intramolecular signal transduction network based on a new protein dynamics model of energy dissipation

A novel approach to reveal intramolecular signal transduction network is proposed in this work. To this end, a new algorithm of network construction is developed, which is based on a new protein dynamics model of energy dissipation. A key feature of this approach is that direction information is specified after inferring protein residue-residue interaction network involved in the process of signal transduction. This enables fundamental analysis of the regulation hierarchy and identification of regulation hubs of the signaling network. A well-studied allosteric enzyme, E. coli aspartokinase III, is used as a model system to demonstrate the new method. Comparison with experimental results shows that the new approach is able to predict all the sites that have been experimentally proved to desensitize allosteric regulation of the enzyme. In addition, the signal transduction network shows a clear preference for specific structural regions, secondary structural types and residue conservation. Occurrence of super-hubs in the network indicates that allosteric regulation tends to gather residues with high connection ability to collectively facilitate the signaling process. Furthermore, a new parameter of propagation coefficient is defined to determine the propagation capability of residues within a signal transduction network. In conclusion, the new approach is useful for fundamental understanding of the process of intramolecular signal transduction and thus has significant impact on rational design of novel allosteric proteins.     

Influence of Electromagnetic Signal of Antibiotics Excited by Low-Frequency Pulsed Electromagnetic Fields on Growth of Escherichia coli

Energy medicine (EM) provides a new medical choice for patients, and its advantages are the noninvasive detection and nondrug treatment. An electromagnetic signal, a kind of EM, induced from antibiotic coupling with weak, extremely low-frequency pulsed electromagnetic fields (PEMFs) is utilized for investigating the growth speed of Escherichia coli (E. coli). PEMFs are produced by solenoidal coils for coupling the electromagnetic signal of antibiotics (penicillin). The growth retardation rate (GRR) of E. coli is used to investigate the efficacy of the electromagnetic signal of antibiotics. The E. coli is cultivated in the exposure of PEMFs coupling with the electromagnetic signal of antibiotics. The maximum GRR of PEMFs with and without the electromagnetic signal of antibiotics on the growth of E. coli cells in the logarithmic is 17.4 and 9.08 %, respectively. The electromagnetic signal of antibiotics is successfully coupled by the electromagnetic signal coupling instrument to affect the growth of E. coli. In addition, the retardation effect on E. coli growth can be improved of by changing the carrier frequency of PEMFs coupling with the electromagnetic signal of antibiotics. GRR caused by the electromagnetic signal of antibiotics can be fixed by a different carrier frequency in a different phase of E. coli growth.     

Directed co-evolution of interacting protein–peptide pairs by compartmentalized two-hybrid replication (C2HR)

Directed evolution methodologies benefit from read-outs quantitatively linking genotype to phenotype. We therefore devised a method that couples protein–peptide interactions to the dynamic read-out provided by an engineered DNA polymerase. Fusion of a processivity clamp protein to a thermostable nucleic acid polymerase enables polymerase activity and DNA amplification in otherwise prohibitive high-salt buffers. Here, we recapitulate this phenotype by indirectly coupling the Sso7d processivity clamp to Taq DNA polymerase via respective fusion to a high affinity and thermostable interacting protein–peptide pair. Escherichia coli cells co-expressing protein–peptide pairs can directly be used in polymerase chain reactions to determine relative interaction strengths by the measurement of amplicon yields. Conditional polymerase activity is further used to link genotype to phenotype of interacting protein–peptide pairs co-expressed in E. coli using the compartmentalized self-replication directed evolution platform. We validate this approach, termed compartmentalized two-hybrid replication, by selecting for high-affinity peptides that bind two model protein partners: SpyCatcher and the large fragment of NanoLuc luciferase. We further demonstrate directed co-evolution by randomizing both protein and peptide components of the SpyCatcher–SpyTag pair and co-selecting for functionally interacting variants.     

Morphing structures and signal transduction in Mimosa pudica L. induced by localized thermal stress

Leaf movements in Mimosa pudica, are in response to thermal stress, touch, and light or darkness, appear to be regulated by electrical, hydrodynamical, and chemical signal transduction. The pulvinus of the M. pudica shows elastic properties. We have found that the movements of the petiole, or pinnules, are accompanied by a change of the pulvinus morphing structures. After brief flaming of a pinna, the volume of the lower part of the pulvinus decreases and the volume of the upper part increases due to the redistribution of electrolytes between these parts of the pulvinus; as a result of these changes the petiole falls. During the relaxation of the petiole, the process goes in the opposite direction. Ion and water channel blockers, uncouplers as well as anesthetic agents diethyl ether or chloroform decrease the speed of alert wave propagation along the plant. Brief flaming of a pinna induces bidirectional propagation of electrical signal in pulvini. Transduction of electrical signals along a pulvinus induces generation of an action potential in perpendicular direction between extensor and flexor sides of a pulvinus. Inhibition of signal transduction and mechanical responses in M. pudica by volatile anesthetic agents chloroform or by blockers of voltage gated ion channels shows that the generation and propagation of electrical signals is a primary effect responsible for turgor change and propagation of an excitation. There is an electrical coupling in a pulvinus similar to the electrical synapse in the animal nerves.     

Multiple barriers to recombination between divergent HIV-1 variants revealed by a dual-marker recombination assay

Recombination is a major force for generating human immunodeficiency virus type 1 (HIV-1) diversity and produces numerous recombinants circulating in the human population. We previously established a cell-based system using green fluorescent protein gene (gfp) as a reporter to study the mechanisms of HIV-1 recombination. We now report an improved system capable of detecting recombination using authentic viral sequences. Frameshift mutations were introduced into the gag gene so that parental viruses do not express full-length Gag; however, recombination can generate a progeny virus that expresses a functional Gag. We demonstrate that this Gag reconstitution assay can be used to detect recombination between two group M HIV-1 variants of the same or of different subtypes. Using both gfp and gag assays, we found that, similar to group M viruses, group O viruses also recombine frequently. When recombination between a group M virus and a group O virus was examined, we found three distinct barriers for intergroup recombination. First, similar to recombination within group M viruses, intergroup recombination is affected by the identity of the dimerization initiation signal (DIS); variants with the same DIS recombined at a higher rate than those with different DIS. Second, using the gfp recombination assay, we showed that intergroup recombination occurs much less frequently than intragroup recombination, even though the gfp target sequence is identical in all viruses. Finally, Gag reconstitution between variants from different groups is further reduced compared with green fluorescent protein, indicating that sequence divergence interferes with recombination efficiency in the gag gene. Compared with identical sequences, we estimate that recombination rates are reduced by 3-fold and by 10- to 13-fold when the target regions in gag contain 91% and 72-73% sequence identities, respectively. These results show that there are at least three distinct mechanisms preventing exchange of genetic information between divergent HIV-1 variants from different groups.     

Analysis of Mid-Infrared Surface Plasmon Modes in a Graphene-Based Cylindrical Hybrid Waveguide

A graphene-based cylindrical hybrid surface plasmon polariton waveguide, composed of a silicon nanowire core surrounded by a silica layer and then a graphene layer, is investigated using the finite-difference time-domain method. The analytical solutions and the numerical simulation show that an ultra-small mode area and a large propagation length can be achieved with this waveguide. Utilizing the perturbation theory of coupled mode, we demonstrate that the six lowest-order coupling modes originate from the coupling of the three lowest-order single-waveguide modes, and the m?=?1 order yy-coupling mode possesses the maximum coupling length and the minimum crosstalk. This waveguide can be used for photonic integrated circuits in the mid-infrared range.     

Signal propagation of the MAPK cascade in Xenopus oocytes: role of bistability and ultrasensitivity for a mixed problem

The MAPK signaling cascade is nowadays understood as a network module highly conserved across species. Its main function is to transfer a signal arriving at the plasma membrane to the cellular interior. Current understanding of ‘how’ this is achieved involves the notions of ultrasensitivity and bistability which relate to the nonlinear dynamics of the biochemical network, ignoring spatial aspects. Much less, indeed, is so far known about the propagation of the signal through the cytoplasm. In this work we formulate, starting from a Michaelis–Menten model for the MAPK cascade in Xenopus oocytes, a reaction-diffusion model of the cascade. We study this model in one space dimension. Basing ourselves on previous general results on reaction diffusion models, we particularly study for our model the conditions for signal propagation. We show that the existence of a propagating front depends sensitively on the initial and boundary conditions at the plasma membrane. Possible biological consequences of this finding are discussed.     

In silico study on the effect of surface lysines and arginines on the electrostatic interactions and protein stability

Charged amino acids are mostly exposed on a protein surface, thereby forming a network of interactions with the surrounding amino acids as well as with water. In particular, positively charged arginine and lysine have different side chain geometries and provide a different number of potential electrostatic interactions. This study reports a comparative analysis of the difference in the number of two representative electrostatic interactions, such as salt-bridges and hydrogen bonds, contributed by surface arginine and lysine, as well as their effect on protein stability using molecular modeling and dynamics simulation techniques. Two in silico variants, the R variant with all arginines and the K variant with all lysines on the protein surface, were modeled by mutating all the surface lysines to arginines and the surface arginines to lysines, respectively, for each of the 10 model proteins. A structural comparison of the respective two variants showed that the majority of R variants possessed more salt-bridges and hydrogen bond interactions than the K variants, indicating that arginine provides a higher probability of electrostatic interactions than lysine owing to its side chain geometry. Molecular dynamics simulations of these variants revealed the R variants to be more stable than the K variants at room temperature but this effect was not prominent under protein denaturating conditions, such as 353 and 333 K with 8 M urea. These results suggest that the arginine residues on a protein surface contribute to the protein stability slightly more than lysine by enhancing the electrostatic interactions.     

Human Chondrocytes Respond Discordantly to the Protein Encoded by the Osteoarthritis Susceptibility Gene GDF5

A genetic deficit mediated by SNP rs143383 that leads to reduced expression of GDF5 is strongly associated with large-joint osteoarthritis. We speculated that this deficit could be attenuated by the application of exogenous GDF5 protein and as a first step we have assessed what effect such application has on primary osteoarthritis chondrocyte gene expression. Chondrocytes harvested from cartilage of osteoarthritic patients who had undergone joint replacement were cultured with wildtype recombinant mouse and human GDF5 protein. We also studied variants of GDF5, one that has a higher affinity for the receptor BMPR-IA and one that is insensitive to the GDF5 antagonist noggin. As a positive control, chondrocytes were treated with TGF-β1. Chondrocytes were cultured in monolayer and micromass and the expression of genes coding for catabolic and anabolic proteins of cartilage were measured by quantitative PCR. The expression of the GDF5 receptor genes and the presence of their protein products was confirmed and the ability of GDF5 signal to translocate to the nucleus was demonstrated by the activation of a luciferase reporter construct. The capacity of GDF5 to elicit an intracellular signal in chondrocytes was demonstrated by the phosphorylation of intracellular Smads. Chondrocytes cultured with TGF-β1 demonstrated a consistent down regulation of MMP1, MMP13 and a consistent upregulation of TIMP1 and COL2A1 with both culture techniques. In contrast, chondrocytes cultured with wildtype GDF5, or its variants, did not show any consistent response, irrespective of the culture technique used. Our results show that osteoarthritis chondrocytes do not respond in a predictable manner to culture with exogenous GDF5. This may be a cause or a consequence of the osteoarthritis disease process and will need to be surmounted if treatment with exogenous GDF5 is to be advanced as a potential means to overcome the genetic deficit conferring osteoarthritis susceptibility at this gene.     

Noise Propagation in Gene Regulation Networks Involving Interlinked Positive and Negative Feedback Loops

It is well known that noise is inevitable in gene regulatory networks due to the low-copy numbers of molecules and local environmental fluctuations. The prediction of noise effects is a key issue in ensuring reliable transmission of information. Interlinked positive and negative feedback loops are essential signal transduction motifs in biological networks. Positive feedback loops are generally believed to induce a switch-like behavior, whereas negative feedback loops are thought to suppress noise effects. Here, by using the signal sensitivity (susceptibility) and noise amplification to quantify noise propagation, we analyze an abstract model of the Myc/E2F/MiR-17-92 network that is composed of a coupling between the E2F/Myc positive feedback loop and the E2F/Myc/miR-17-92 negative feedback loop. The role of the feedback loop on noise effects is found to depend on the dynamic properties of the system. When the system is in monostability or bistability with high protein concentrations, noise is consistently suppressed. However, the negative feedback loop reduces this suppression ability (or improves the noise propagation) and enhances signal sensitivity. In the case of excitability, bistability, or monostability, noise is enhanced at low protein concentrations. The negative feedback loop reduces this noise enhancement as well as the signal sensitivity. In all cases, the positive feedback loop acts contrary to the negative feedback loop. We also found that increasing the time scale of the protein module or decreasing the noise autocorrelation time can enhance noise suppression; however, the systems sensitivity remains unchanged. Taken together, our results suggest that the negative/positive feedback mechanisms in coupled feedback loop dynamically buffer noise effects rather than only suppressing or amplifying the noise.     

Petri net-based method for the analysis of the dynamics of signal propagation in signaling pathways

MOTIVATION: Cellular signaling networks are dynamic systems that propagate and process information, and, ultimately, cause phenotypical responses. Understanding the circuitry of the information flow in cells is one of the keys to understanding complex cellular processes. The development of computational quantitative models is a promising avenue for attaining this goal. Not only does the analysis of the simulation data based on the concentration variations of biological compounds yields information about systemic state changes, but it is also very helpful for obtaining information about the dynamics of signal propagation. RESULTS: This article introduces a new method for analyzing the dynamics of signal propagation in signaling pathways using Petri net theory. The method is demonstrated with the Ca(2+)/calmodulin-dependent protein kinase II (CaMKII) regulation network. The results constitute temporal information about signal propagation in the network, a simplified graphical representation of the network and of the signal propagation dynamics and a characterization of some signaling routes as regulation motifs.     

Electrostatic effects in a network of polar and ionizable groups in staphylococcal nuclease

His121 and His124 are embedded in a network of polar and ionizable groups on the surface of staphylococcal nuclease. To examine how membership in a network affects the electrostatic properties of ionizable groups, the tautomeric state and the pK_a values of these histidines were measured with NMR spectroscopy in the wild-type nuclease and in 13 variants designed to disrupt the network. In the background protein, His121 and His124 titrate with pK_a values of 5.2 and 5.6, respectively. In the variants, where the network was disrupted, the pK_a values range from 4.03 to 6.46 for His121, and 5.04 to 5.99 for His124. The largest decrease in a pK_a was observed when the favorable Coulomb interaction between His121 and Glu75 was eliminated; the largest increase was observed when Tyr91 or Tyr93 was substituted with Ala or Phe. In all variants, the dominant tautomeric state at neutral pH was the N^ε2 state. At one level the network behaves as a rigid unit that does not readily reorganize when disrupted: crystal structures of the E75A or E75Q variants show that even when the pivotal Glu75 is removed, the overall configuration of the network was unaffected. On the other hand, a few key hydrogen bonds appear to govern the conformation of the network, and when these bonds are disrupted the network reorganizes. Coulomb interactions within the network report an effective dielectric constant of 20, whereas a dielectric constant of 80 is more consistent with the magnitude of medium to long-range Coulomb interactions in this protein. The data demonstrate that when structures are treated as static, rigid bodies, structure-based pK_a calculations with continuum electrostatics method are not useful to treat ionizable groups in cases where pK_a values are governed by short-range polar and Coulomb interactions.