A rate distortion approach to protein symmetry

A spontaneous symmetry breaking argument is applied to the problem of protein folding, via a rate distortion analysis of the relation between genome coding and the final condensation of the protein molten globule that is, in spirit, analogous to Tlusty’s (2007) exploration of the evolution of the genetic code. In the ‘energy’ picture, the average distortion between codon message and final protein structure, under constraints driven by evolutionary selection, serves as a temperature analog, so that low values limit the possible distribution of protein forms, producing the canonical folding funnel. A dual ‘developmental’ perspective sees the rate distortion function itself as the temperature analog, and permits incorporation of chaperons or toxic exposures as catalysts, driving the system to different possible outcomes or affecting the rate of convergence. The rate distortion function appears constrained by the availability of metabolic free energy, with implications for prebiotic evolution, and a nonequilibrium empirical Onsager treatment provides an adaptable statistical model that can be fitted to data, in the same manner as a regression equation. In sum, mechanistic models of protein folding fail to account for the observed spectrum of protein folding and aggregation disorders, suggesting that a biologically based cognitive paradigm describing folding will be needed for understanding the etiology, prevention, and treatment of these diseases. The developmental formalism introduced here may contribute substantially to such a paradigm.     

Protein folding disorders: Toward a basic biological paradigm

Mechanistic ‘physics’ models of protein folding fail to account for the observed spectrum and rate of protein folding and aggregation disorders in human populations, showing that more appropriately in vivo paradigms reflecting biological and other embedding contexts are needed for understanding the etiology, prevention, and treatment of these diseases. Here, a topological rate distortion analysis is applied to the problem that is analogous to Tlusty (2007) elegant exploration of the genetic code. A ‘developmental’ perspective sees the rate distortion function as a temperature analog in a spontaneous symmetry breaking argument, and permits incorporation of external factors as catalysts, driving the system to different possible outcomes via a nonequilibrium empirical Onsager treatment, viewed as a kind of dynamic regression equation. The formalism produces large-scale, quasi-equilibrium ‘resilience’ states representing normal and pathological protein folding. Generalization to long times produces diffusion models of protein folding disorders in which epigenetic or life history factors determine the rate of onset of dysfunction.     

Matching Biochemical Reaction Kinetics to the Timescales of Life: Structural Determinants That Influence the Autodephosphorylation Rate of Response Regulator Proteins

In two-component regulatory systems, covalent phosphorylation typically activates the response regulator signaling protein, and hydrolysis of the phosphoryl group reestablishes the inactive state. Despite highly conserved three-dimensional structures and active-site features, the rates of catalytic autodephosphorylation for different response regulators vary by a factor of almost 10⁶. Previous studies identified two variable active-site residues, corresponding to Escherichia coli CheY residues 59 and 89, that modulate response regulator autodephosphorylation rates about 100-fold. Here, a set of five CheY mutants, which match other “model” response regulators (ArcA, CusR, DctD, FixJ, PhoB, or Spo0F) at variable active-site positions corresponding to CheY residues 14, 59, and 89, were characterized functionally and structurally in an attempt to identify mechanisms that modulate autodephosphorylation rate. As expected, the autodephosphorylation rates of the CheY mutants were reduced 6- to 40-fold relative to wild-type CheY, but all still autodephosphorylated 12- to 80-fold faster than their respective model response regulators. Comparison of X-ray crystal structures of the five CheY mutants (complexed with the phosphoryl group analogue BeF₃⁻) to wild-type CheY or corresponding model response regulator structures gave strong evidence for steric obstruction of the phosphoryl group from the attacking water molecule as one mechanism to enhance phosphoryl group stability. Structural data also suggested that impeding the change of a response regulator from the active to the inactive conformation might retard the autodephosphorylation reaction if the two processes are coupled, and that the residue at position ‘58’ may contribute to rate modulation. A given combination of amino acids at positions ‘14’, ‘59’, and ‘89’ adopted similar conformations regardless of protein context (CheY or model response regulator), suggesting that knowledge of residue identity may be sufficient to predict autodephosphorylation rate, and hence the kinetics of the signaling response, in the response regulator family of proteins.     

A dual expression platform to optimize the soluble production of heterologous proteins in the periplasm of <Emphasis Type="Italic">Escherichia coli</Emphasis>

The functional analysis of individual proteins or of multiprotein complexes—since the completion of several genome sequencing projects—is in focus of current scientific work. Many heterologous proteins contain disulfide-bonds, required for their correct folding and activity, and therefore, need to be transported to the periplasm. The production of soluble and functional protein in the periplasm often needs target-specific regulatory genetic elements, leader peptides, and folding regimes. Usually, the optimization of periplasmic expression is a step-wise and time-consuming procedure. To overcome this problem we developed a dual expression system, containing a degP-promoter-based reporter system and a highly versatile plasmid set. This combines the differential protein expression with the selection of a target-specific expression plasmid. For the validation of this expression tool, two different molecular formats of a recombinant antibody directed to the human epidermal growth factor receptor and human 11β-hydroxysteroid dehydrogenase type 2 (11β-HSD2) were used. By application of this expression system we demonstrated that the amount of functional protein is inversely proportional to the on-line luciferase signal. We showed that this technology offers a simple tool to evaluate and improve the yield of functionally expressed proteins in the periplasm, which depends on the used regulatory elements and folding strategies.     

Salt-stress induced changes in the leaf proteome of diploid and tetraploid mandarins with contrasting Na and Cl accumulation behaviour

To understand the genotypic variation of citrus to mild salt stress, a proteomic approach has been carried out in parallel on two citrus genotypes (‘Cleopatra’ and ‘Willow leaf’ mandarins), which differ for Na⁺ and Cl⁻ accumulation, and their cognate autotetraploids (4×). Using two-dimensional electrophoresis approximately 910 protein spots were reproducibly detected in control and salt-stressed leaves of all genotypes. Among them, 44 protein spots showing significant variations at least in one genotype were subjected to mass spectrometry analysis for identification. Salt-responsive proteins were involved in several functions, including photosynthetic processes, ROS scavenging, stress defence, and signalling. Genotype factors affect the salt-responsive pattern, especially that of carbon metabolism. The no ion accumulator ‘Cleopatra’ mandarin genotype showed the highest number of salt-responsive proteins, and up-regulation of Calvin cycle-related proteins. Conversely the ion accumulator ‘Willow leaf’ mandarin showed high levels of several photorespiration-related enzymes. A common set of proteins (twelve spots) displayed higher levels in salt-stressed leaves of 2× and 4× ‘Cleopatra’ and 4× ‘Willow leaf’ mandarin. Interestingly, antioxidant enzymes and heat shock proteins showed higher constitutive levels in 4× ‘Cleopatra’ mandarin and 4× ‘Willow leaf’ mandarin compared with the cognate 2× genotype. This work provides for the first time information on the effect of 8 weeks of salt stress on citrus genotypes contrasting for ion accumulation and their cognate autotetraploids. Results underline that genetic factors have a predominant effect on the salt response, although a common stress response independent from genotype was also found.     

Quantum biology at the cellular level—Elements of the research program

Quantum biology is emerging as a new field at the intersection between fundamental physics and biology, promising novel insights into the nature and origin of biological order. We discuss several elements of QBCL (quantum biology at cellular level) – a research program designed to extend the reach of quantum concepts to higher than molecular levels of biological organization. We propose a new general way to address the issue of environmentally induced decoherence and macroscopic superpositions in biological systems, emphasizing the ‘basis-dependent’ nature of these concepts. We introduce the notion of ‘formal superposition’ and distinguish it from that of Schroedinger's cat (i.e., a superposition of macroscopically distinct states). Whereas the latter notion presents a genuine foundational problem, the former one contradicts neither common sense nor observation, and may be used to describe cellular ‘decision-making’ and adaptation. We stress that the interpretation of the notion of ‘formal superposition’ should involve non-classical correlations between molecular events in a cell. Further, we describe how better understanding of the physics of Life can shed new light on the mechanism driving evolutionary adaptation (viz., ‘Basis-Dependent Selection’, BDS). Experimental tests of BDS and the potential role of synthetic biology in closing the ‘evolvability mechanism’ loophole are also discussed.     

Genetic variation and drought response in two Populus × euramericana genotypes through 2-DE proteomic analysis of leaves from field and glasshouse cultivated plants

Genotype and water deficit effects on leaf 2-DE protein profiles of two Populus deltoides × Populus nigra, cv. ‘Agathe_F’ and ‘Cima’, were analysed over a short-term period of 18 days in glasshouse using 4-month-old rooted cuttings and over a long-lasting period of 86 days in open field using 4-year-old rooted cuttings. Leaf proteomes were analyzed using two-dimensional gel electrophoresis, and proteins were identified after database searching from MS peptide spectra.A reliable genotype effect was observed in the leaf proteome over experiment locations, water regimes and sampling dates. Quantitative differences between genotypes were found. Most of them corresponded to proteins matching isoforms or post-translational modification variants. However, ‘Cima’ displayed the highest abundance of antioxidant enzymes.In response to water deficit, about 10% of the reproducible spots significantly varied regardless of the experiment location, among which about 25% also displayed genotype-dependent variations. As a whole, while ‘Cima’ differed from ‘Agathe_F’ by increased abundance of enzymes involved in photorespiration and in oxidative stress, ‘Agathe_F’ was mainly differentiated by increased abundance of enzymes involved in photosynthesis.     

Kinetics of the Multistep Rupture of Fibrin ‘A-a’ Polymerization Interactions Measured Using Atomic Force Microscopy

Fibrin, the structural scaffold of blood clots, spontaneously polymerizes through the formation of ‘A-a’ knob-hole bonds. When subjected to external force, the dissociation of this bond is accompanied by two to four abrupt changes in molecular dimension observable as rupture events in a force curve. Herein, the configuration, molecular extension, and kinetic parameters of each rupture event are examined. The increases in contour length indicate that the D region of fibrinogen can lengthen by ∼50% of the length of a fibrin monomer before rupture of the ‘A-a’ interaction. The dependence of the dissociation rate on applied force was obtained using probability distributions of rupture forces collected at different pull-off velocities. These distributions were fit using a model in which the effects of the shape of the binding potential are used to quantify the kinetic parameters of forced dissociation. We found that the weak initial rupture (i.e., event 1) was not well approximated by these models. The ruptured bonds comprising the strongest ruptures, events 2 and 3, had kinetic parameters similar to those commonly found for the mechanical unfolding of globular proteins. The bonds ruptured in event 4 were well described by these analyses, but were more loosely bound than the bonds in events 2 and 3. We propose that the first event represents the rupture of an unknown interaction parallel to the ‘A-a’ bond, events 2 and 3 represent unfolding of structures in the D region of fibrinogen, and event 4 is the rupture of the ‘A-a’ knob-hole bond weakened by prior structural unfolding. Comparison of the activation energy obtained via force spectroscopy measurements with the thermodynamic free energy of ‘A-a’ bond dissociation indicates that the ‘A-a’ bond may be more resistant to rupture by applied force than to rupture by thermal dissociation.     

Acceleration of oxidative protein folding by curcumin through novel non-redox chemistry

Curcumin, the major constituent of turmeric is a known antioxidant. We have examined the oxidative folding of the model four-disulfide-bond-containing protein bovine pancreatic ribonuclease A (RNase A) in its presence; results indicate that RNase A regeneration rate increases in a curcumin-dependent manner. Examination of the native tendency of the fully-reduced polypeptide and the stability of key folding intermediates suggests that the increased oxidative folding rate can be attributed to native-like elements induced within the fully-reduced polypeptide and the stabilization of native-like species by this non-redox-active natural product. Our results provide a template for the design of curcuminoid-based synthetic small-molecule fold catalysts that accelerate the folding of ER-processed proteins; this assumes significance given that nitrosative stress and dysfunction of the ER-resident oxidoreductase protein disulfide isomerise due to S-nitrosylation are factors associated with the pathogenesis of Alzheimer’s and Parkinson’s diseases.     

A genetic screening strategy identifies novel regulators of the proteostasis network

A hallmark of diseases of protein conformation and aging is the appearance of protein aggregates associated with cellular toxicity. We posit that the functional properties of the proteostasis network (PN) protect the proteome from misfolding and combat the proteotoxic events leading to cellular pathology. In this study, we have identified new components of the proteostasis network that can suppress aggregation and proteotoxicity, by performing RNA interference (RNAi) genetic screens for multiple unrelated conformationally challenged cytoplasmic proteins expressed in Caenorhabditis elegans. We identified 88 suppressors of polyglutamine (polyQ) aggregation, of which 63 modifiers also suppressed aggregation of mutant SOD1(G93A). Of these, only 23 gene-modifiers suppressed aggregation and restored animal motility, revealing that aggregation and toxicity can be genetically uncoupled. Nine of these modifiers were shown to be effective in restoring the folding and function of multiple endogenous temperature-sensitive (TS) mutant proteins, of which five improved folding in a HSF-1-dependent manner, by inducing cytoplasmic chaperones. This triage screening strategy also identified a novel set of PN regulatory components that, by altering metabolic and RNA processing functions, establish alternate cellular environments not generally dependent on stress response activation and that are broadly protective against misfolded and aggregation-prone proteins.     

Protein folding rates and thermodynamic stability are key determinants for interaction with the Hsp70 chaperone system

The Hsp70 family of molecular chaperones participates in vital cellular processes including the heat shock response and protein homeostasis. E. coli''s Hsp70, known as DnaK, works in concert with the DnaJ and GrpE co-chaperones (K/J/E chaperone system), and mediates cotranslational and post-translational protein folding in the cytoplasm. While the role of the K/J/E chaperones is well understood in the presence of large substrates unable to fold independently, it is not known if and how K/J/E modulates the folding of smaller proteins able to fold even in the absence of chaperones. Here, we combine experiments and computation to evaluate the significance of kinetic partitioning as a model to describe the interplay between protein folding and binding to the K/J/E chaperone system. First, we target three nonobligatory substrates, that is, proteins that do not require chaperones to fold. The experimentally observed chaperone association of these client proteins during folding is entirely consistent with predictions from kinetic partitioning. Next, we develop and validate a computational model (CHAMP70) that assumes kinetic partitioning of substrates between folding and interaction with K/J/E. CHAMP70 quantitatively predicts the experimentally measured interaction of RNase H^D as it refolds in the presence of various chaperones. CHAMP70 shows that substrates are posed to interact with K/J/E only if they are slow-folding proteins with a folding rate constant k_f <50 s⁻¹, and/or thermodynamically unstable proteins with a folding free energy ΔG⁰_UN ≥−2 kcal mol⁻¹. Hence, the K/J/E system is tuned to use specific protein folding rates and thermodynamic stabilities as substrate selection criteria.     

Certain heptapeptide and large sequences representing an entire helix, strand or coil conformation in proteins are associated as chameleon sequences

Helices, strands and coils in proteins of known three-dimensional structure, corresponding to heptapeptide and large sequences (‘probe’ peptides), were scanned against peptide sequences of variable length, comprising seven or more residues that correspond to a different conformation (‘target’ peptides) in protein crystal structures available from the Protein Data Bank (PDB). Where the ‘probe’ and ‘target’ peptide sequences exactly match, they correspond to ‘chameleon’ sequences in protein structures. We observed ∼548 heptapeptide and large chameleon sequences that included peptides in the coil conformation from 53,794 PDB files that were analyzed. However, after excluding several chameleon peptides based on the quality of protein structure data, redundancy and peptides associated with cloning artifacts, such as, histidine-tags, we observed only ten chameleon peptides in structurally different proteins and the maximum length comprised seven amino acid residues. Our analysis suggests that the quality of protein structure data is important for identifying possibly, the ‘true chameleons’ in PDB. Majority of the chameleon sequences correspond to an entire strand in one protein that is observed as part of helix sequence in another protein. The heptapeptide chameleons are characterized with a high propensity of alanine, leucine and valine amino acid residues. The total hydropathy values range between −11.2 and 22.9, the difference in solvent accessibility between 2.0 Å² and 373 Å² units and the difference in total number of residue neighbor contacts between 0 and 7 residues. Our work identifies for the first time heptapeptide and large sequences that correspond to a single complete helix, strand or coil, which adopt entirely different secondary structures in another protein.     

Genetic diversity of a relict plant species, Ligularia sibirica (L.) Cass. (Asteraceae)

Rare plant species can be divided into naturally, ‘old rare’ species and anthropogenically, ‘new rare’ species. Many recent studies explored genetic diversity of ‘new rare’ species. Less is, however, known about genetic diversity of ‘old rare’ species. We examined isozyme genetic variability of 20 populations of an ‘old rare’ plant species, Ligularia sibirica (Asteraceae) in the Czech and Slovak Republic. It is a long-lived perennial herb with mixed-mating breeding system, widely distributed from East Asia to European Russia, with few isolated relict populations in the remaining part of Europe.The results showed high genetic diversity within populations (80.8%) and a low level of genetic differentiation (F_ST = 0.179). Genetic distance between populations correlated significantly with geographic distance. There was also a significant positive correlation between genetic diversity and population size. This is probably caused by destruction of habitats in last centuries and subsequent decrease of population size. Patterns of genetic diversity suggest that the recent distribution is a result of stepwise postglacial migration of the species and subsequent natural fragmentation.We conclude that L. sibirica populations preserve high levels of genetic diversity and are not yet threatened by genetic factors. However, this may change if changes in habitat conditions continue.     

A new formal perspective on ‘Cambrian explosions’

The ‘Cambrian explosion’ 500 Myr ago saw a relatively sudden proliferation of organism Bauplan and ecosystem niche structure that continues to haunt evolutionary biology. Here, adapting standard methods from information theory and statistical mechanics, we model the phenomenon as a noise-driven phase transition, in the context of deep-time relaxation of current path-dependent evolutionary constraints. The result is analogous to recent suggestions that multiple ‘explosions’ of increasing complexity in the genetic code were driven by rising intensities of available metabolic free energy. In the absence of severe path-dependent lock-in, ‘Cambrian explosions’ are standard features of blind evolutionary process, representing outliers in the ongoing routine of evolutionary punctuated equilibrium.     

The implications of gene heterozygosity for protein folding and protein turnover

The offspring of closely related parents often suffer from inbreeding depression, sometimes resulting in a slower growth rate for inbred offspring relative to non-inbred offspring. Previous research has shown that some of the slower growth rate of inbred organisms can be attributed to the inbred organisms’ increased levels of protein turnover. This paper attempts to show that the higher levels of protein turnover among inbred organisms can be attributed to accumulations of misfolded and aggregated proteins that require degradation by the inbred organisms’ protein quality control systems. The accumulation of misfolded and aggregated proteins within inbred organisms are the result of more negative free energies of folding for proteins encoded at homozygous gene loci and higher concentrations of potentially aggregating non-native protein species within the cell. The theory presented here makes several quantitative predictions that suggest a connection between protein misfolding/aggregation and polyploidy that can be tested by future research.