Search

The limits of evolution have long fascinated biologists. However, the causes of evolutionary constraint have remained elusive due to a poor mechanistic understanding of studied phenotypes. Recently, a range of innovative approaches have leveraged mechanistic information on regulatory networks and cellular biology. These methods combine systems biology models with population and single-cell quantification and with new genetic tools, and they have been applied to a range of complex cellular functions and engineered networks. In this article, we review these developments, which are revealing the mechanistic causes of epistasis at different levels of biological organization¤mdash¤in molecular recognition, within a single regulatory network, and between different networks¤mdash¤providing first indications of predictable features of evolutionary constraint.

Funder: Nederlandse Organisatie voor Wetenschappelijk Onderzoek (Netherlands Organisation for Scientific Research); doi: https://doi.org/10.13039/501100003246, Many proteins form dynamic complexes with DNA, RNA, and other proteins, which often involves protein conformational changes that are key to function. Yet, methods to probe these critical dynamics are scarce. Here we combine optical tweezers with fluorescence imaging to simultaneously monitor the conformation of individual proteins and their binding to partner proteins. Central is a protein–DNA coupling strategy, which uses exonuclease digestion and partial re-synthesis to generate DNA overhangs of different lengths, and ligation to oligo-labeled proteins. It provides up to 40 times higher coupling yields than existing protocols and enables new fluorescence-tweezers assays, which require particularly long and strong DNA handles. We demonstrate the approach by detecting the emission of a tethered fluorescent protein and of a molecular chaperone (trigger factor) complexed with its client. We conjecture that our strategy will be an important tool to study conformational dynamics within larger biomolecular complexes.

Protein folding is well known to be supervised by a dedicated class of proteins called chaperones. However, the core mode of action of these molecular machines has remained elusive due to several reasons including the promiscuous nature of the interactions between chaperones and their many clients, as well as the dynamics and heterogeneity of chaperone conformations and the folding process itself. While troublesome for traditional bulk techniques, these properties make an excellent case for the use of single-molecule approaches. In this review, we will discuss how force spectroscopy, fluorescence microscopy, FCS, and FRET methods are starting to zoom in on this intriguing and diverse molecular toolbox that is of direct importance for protein quality control in cells, as well as numerous degenerative conditions that depend on it.

Hsp70 binds unfolded protein segments in its groove, but can also bind and stabilize folded protein structures, owing to its moveable lid, with ATP hydrolysis and co-chaperones allowing control of these contrasting effects. The protein-chaperone system centred on Hsp70 performs a variety of cellular control tasks, including folding assistance, protection against aggregation, trafficking and regulation of enzyme activity, a versatility that has been hard to reconcile with structural data, which suggest that Hsp70 only binds extended polypeptide segments. Now, using laser molecular tweezers, Sander Tans and colleagues show that the bacterial homolog of Hsp70, known as DnaK, relies on its 'groove' to bind unfolded proteins, but can also bind folded structures, thanks to its 'lid', with control of ATP hydrolysis by co-chaperones allowing regulation of such contrasting effects. Contrary to known stabilization mechanisms, through precise structural fit, Hsp70 can stabilize a vast repertoire of client proteins, through a clamp-like, ATP-driven conformational change. The Hsp70 system is a central hub of chaperone activity in all domains of life. Hsp70 performs a plethora of tasks, including folding assistance, protection against aggregation, protein trafficking, and enzyme activity regulation1,2,3,4,5, and interacts with non-folded chains, as well as near-native, misfolded, and aggregated proteins6,7,8,9,10. Hsp70 is thought to achieve its many physiological roles by binding peptide segments that extend from these different protein conformers within a groove that can be covered by an ATP-driven helical lid11,12,13,14,15. However, it has been difficult to test directly how Hsp70 interacts with protein substrates in different stages of folding and how it affects their structure. Moreover, recent indications of diverse lid conformations in Hsp70–substrate complexes raise the possibility of additional interaction mechanisms15,16,17,18. Addressing these issues is technically challenging, given the conformational dynamics of both chaperone and client, the transient nature of their interaction, and the involvement of co-chaperones and the ATP hydrolysis cycle19. Here, using optical tweezers, we show that the bacterial Hsp70 homologue (DnaK) binds and stabilizes not only extended peptide segments, but also partially folded and near-native protein structures. The Hsp70 lid and groove act synergistically when stabilizing folded structures: stabilization is abolished when the lid is truncated and less efficient when the groove is mutated. The diversity of binding modes has important consequences: Hsp70 can both stabilize and destabilize folded structures, in a nucleotide-regulated manner; like Hsp90 and GroEL, Hsp70 can affect the late stages of protein folding; and Hsp70 can suppress aggregation by protecting partially folded structures as well as unfolded protein chains. Overall, these findings in the DnaK system indicate an extension of the Hsp70 canonical model that potentially affects a wide range of physiological roles of the Hsp70 system.

Organoids are self-organizing 3D structures grown from stem cells that recapitulate essential aspects of organ structure and function. Here, we describe a method to establish long-term-expanding human airway organoids from broncho-alveolar resections or lavage material. The pseudostratified airway organoids consist of basal cells, functional multi-ciliated cells, mucus-producing secretory cells, and CC10-secreting club cells. Airway organoids derived from cystic fibrosis (CF) patients allow assessment of CFTR function in an organoid swelling assay. Organoids established from lung cancer resections and metastasis biopsies retain tumor histopathology as well as cancer gene mutations and are amenable to drug screening. Respiratory syncytial virus (RSV) infection recapitulates central disease features, dramatically increases organoid cell motility via the non-structural viral NS2 protein, and preferentially recruits neutrophils upon co-culturing. We conclude that human airway organoids represent versatile models for the in vitro study of hereditary, malignant, and infectious pulmonary disease.

Chaperones assist in protein folding, but what this common phrase means in concrete terms has remained surprisingly poorly understood. We can readily measure chaperone binding to unfolded proteins, but how they bind and affect proteins along folding trajectories has remained obscure. Here we review recent efforts by our labs and others that are beginning to pry into this issue, with a focus on the chaperones trigger factor and Hsp70. Single-molecule methods are central, as they allow the stepwise process of folding to be followed directly. First results have already revealed contrasts with long-standing paradigms: rather than acting only "early" by stabilizing unfolded chain segments, these chaperones can bind and stabilize partially folded structures as they grow to their native state. The findings suggest a fundamental redefinition of the protein folding problem and a more extensive functional repertoire of chaperones than previously assumed.

Elucidating the role of molecular stochasticity in cellular growth is central to understanding phenotypic heterogeneity and the stability of cellular proliferation. The inherent stochasticity of metabolic reaction events should have negligible effect, because of averaging over the many reaction events contributing to growth. Indeed, metabolism and growth are often considered to be constant for fixed conditions. Stochastic fluctuations in the expression level of metabolic enzymes could produce variations in the reactions they catalyse. However, whether such molecular fluctuations can affect growth is unclear, given the various stabilizing regulatory mechanisms, the slow adjustment of key cellular components such as ribosomes, and the secretion and buffering of excess metabolites. Here we use time-lapse microscopy to measure fluctuations in the instantaneous growth rate of single cells of Escherichia coli, and quantify time-resolved cross-correlations with the expression of lac genes and enzymes in central metabolism. We show that expression fluctuations of catabolically active enzymes can propagate and cause growth fluctuations, with transmission depending on the limitation of the enzyme to growth. Conversely, growth fluctuations propagate back to perturb expression. Accordingly, enzymes were found to transmit noise to other unrelated genes via growth. Homeostasis is promoted by a noise-cancelling mechanism that exploits fluctuations in the dilution of proteins by cell-volume expansion. The results indicate that molecular noise is propagated not only by regulatory proteins but also by metabolic reactions. They also suggest that cellular metabolism is inherently stochastic, and a generic source of phenotypic heterogeneity.

The folding of complex proteins can be dramatically affected by misfolding transitions. Directly observing misfolding and distinguishing it from aggregation is challenging. Experiments with optical tweezers revealed transitions between the folded states of a single protein in the absence of mechanical tension. Nonfolded chains of the multidomain protein luciferase folded within seconds to different partially folded states, one of which was stable over several minutes and was more resistant to forced unfolding than other partially folded states. Luciferase monomers can thus adopt a stable misfolded state and can do so without interacting with aggregation partners. This result supports the notion that luciferase misfolding is the cause of the low refolding yields and aggregation observed with this protein. This approach could be used to study misfolding transitions in other large proteins, as well as the factors that affect misfolding.

Inferring the directionality of interactions between cellular processes is a major challenge in systems biology. Time-lagged correlations allow to discriminate between alternative models, but they still rely on assumed underlying interactions. Here, we use the transfer entropy (TE), an information-theoretic quantity that quantifies the directional influence between fluctuating variables in a model-free way. We present a theoretical approach to compute the transfer entropy, even when the noise has an extrinsic component or in the presence of feedback. We re-analyze the experimental data from Kiviet et al. (2014) where fluctuations in gene expression of metabolic enzymes and growth rate have been measured in single cells of E. coli. We confirm the formerly detected modes between growth and gene expression, while prescribing more stringent conditions on the structure of noise sources. We furthermore point out practical requirements in terms of length of time series and sampling time which must be satisfied in order to infer optimally transfer entropy from times series of fluctuations., 24 pages, 7 figures

Microorganisms have developed an elaborate spectrum of mechanisms to respond and adapt to environmental stress conditions. Among these is the expression of dps , coding for the D NA-binding p rotein from s tarved cells. Dps becomes the dominant nucleoid-organizing protein in stationary-phase Escherichia coli cells and is required for robust survival under stress conditions, including carbon or nitrogen starvation, oxidative stress, metal exposure, and irradiation. To study the complex regulation of Dps in E. coli , we utilized time-lapse fluorescence microscopy imaging to examine the kinetics, input encoding, and variability of the Dps response in single cells. In the presence of an oxidative stressor, we observed a single pulse of activation of Dps production. Increased concentrations of H 2 O 2 led to increased intensity and duration of the pulse. While lower concentrations of H 2 O 2 robustly activated the Dps response with little effect on the growth rate, higher concentrations of H 2 O 2 resulted in dramatically lower and highly varied growth rates. A comparison of cells exposed to the same concentration of H 2 O 2 revealed that increased levels of Dps expression did not confer a growth advantage, indicating that recovery from stress may rely primarily upon variation in the amount of damage caused to individual cells. IMPORTANCE We show for the first time the response of the DNA-binding protein from starved cells (Dps) to oxidative stress in single cells of E. coli . Through time-lapse fluorescence microscopy, a single pulse of Dps production is observed in cells exposed to H 2 O 2 , with a duration and intensity of induction proportional to the concentration of the applied stress. More intense Dps expression did not provide a growth benefit to the bacteria, suggesting that healing from oxidative stress may largely depend upon the amount of damage in each individual cell.

How chaperone interactions affect protein folding pathways is a central problem in biology. With the use of optical tweezers and all-atom molecular dynamics simulations, we studied the effect of chaperone SecB on the folding and unfolding pathways of maltose binding protein (MBP) at the single-molecule level. In the absence of SecB, we find that the MBP polypeptide first collapses into a molten globulelike compacted state and then folds into a stable core structure onto which several α helices are finally wrapped. Interactions with SecB completely prevent stable tertiary contacts in the core structure but have no detectable effect on the folding of the external α helices. It appears that SecB only binds to the extended or molten globulelike structure and retains MBP in this latter state. Thus during MBP translocation, no energy is required to disrupt stable tertiary interactions.

How chaperones interact with protein chains to assist in their folding is a central open question in biology. Obtaining atomistic insight is challenging in particular, given the transient nature of the chaperone-substrate complexes and the large system sizes. Recent single-molecule experiments have shown that the chaperone Trigger Factor (TF) not only binds unfolded protein chains, but can also guide protein chains to their native state by interacting with partially folded structures. Here, we used all-atom MD simulations to provide atomistic insights into how Trigger Factor achieves this chaperone function. Our results indicate a crucial role for the tips of the finger-like appendages of TF in the early interactions with both unfolded chains and partially folded structures. Unfolded chains are kinetically trapped when bound to TF, which suppresses the formation of transient, non-native end-to-end contacts. Mechanical flexibility allows TF to hold partially folded structures with two tips (in a pinching configuration), and to stabilize them by wrapping around its appendages. This encapsulation mechanism is distinct from that of chaperones such as GroEL, and allows folded structures of diverse size and composition to be protected from aggregation and misfolding interactions. The results suggest that an ATP cycle is not required to enable both encapsulation and liberation., Author Summary Trigger Factor (TF) is an ATP-independent chaperone protein that assists in folding and prevents misfolding. Up to now, it is a general unsolved question how chaperones assist in the folding of protein chains. Experimental methods that can probe at the length and timescales of inter-residue interactions are scarce, while the systems are too large—and the folding process too long—to be studied by computer simulations. To overcome these obstacles, the authors performed molecular dynamics simulations at key moments along the folding pathway, and address the changes in the folding and unfolding dynamics of protein chains while in contact with TF. This study provides the first detailed view on a chaperone-protein complex in different stages of folding and offers an explanation for the ability of TF to guide chains to their native state. Moreover, the results demonstrates the role of TF’s flexibility in interacting with a wide range of client states. Overall, it explains how TF can interact with many types of substrates in various stages of folding, without the need for an ATP cycle to switch between encapsulation and liberation of client proteins.

As part of the viral infection cycle, viruses must package their newly replicated genomes for delivery to other host cells. Bacteriophage straight phi29 packages its 6.6-microm long, double-stranded DNA into a 42 x 54 nm capsid by means of a portal complex that hydrolyses ATP. This process is remarkable because entropic, electrostatic and bending energies of the DNA must be overcome to package the DNA to near-crystalline density. Here we use optical tweezers to pull on single DNA molecules as they are packaged, thus demonstrating that the portal complex is a force-generating motor. This motor can work against loads of up to 57 pN on average, making it one of the strongest molecular motors reported to date. Movements of over 5 microm are observed, indicating high processivity. Pauses and slips also occur, particularly at higher forces. We establish the force-velocity relationship of the motor and find that the rate-limiting step of the motor's cycle is force dependent even at low loads. Notably, the packaging rate decreases as the prohead is filled, indicating that an internal force builds up to approximately 50 pN owing to DNA confinement. Our data suggest that this force may be available for initiating the ejection of the DNA from the capsid during infection.

The epistatic interactions that underlie evolutionary constraint have mainly been studied for constant external conditions. However, environmental changes may modulate epistasis and hence affect genetic constraints. Here we investigate genetic constraints in the adaptive evolution of a novel regulatory function in variable environments, using the lac repressor, LacI, as a model system. We have systematically reconstructed mutational trajectories from wild type LacI to three different variants that each exhibit an inverse response to the inducing ligand IPTG, and analyzed the higher-order interactions between genetic and environmental changes. We find epistasis to depend strongly on the environment. As a result, mutational steps essential to inversion but inaccessible by positive selection in one environment, become accessible in another. We present a graphical method to analyze the observed complex higher-order interactions between multiple mutations and environmental change, and show how the interactions can be explained by a combination of mutational effects on allostery and thermodynamic stability. This dependency of genetic constraint on the environment should fundamentally affect evolutionary dynamics and affects the interpretation of phylogenetic data., Author Summary Epistatic interactions limit the number of adaptive trajectories to peaks on evolutionary fitness landscapes, and may therefore hamper the progress of evolution. Recent research has focused on adaptive landscapes in one constant environment. However, adaptive evolution is generally known to occur in variable, heterogeneous environments. Here, we have constructed fitness landscapes of three inverse lac repressor variants in two contrasting environments. We find that the epistatic interactions between the pairs of mutations are profoundly altered upon an environmental change. We develop a new graphical method to analyze the underlying higher-order interactions between genetic changes and the environment, and explain the complex environmental dependencies in terms of simple molecular mechanisms. Our results show that the information about epistatic interactions acquired in one environment does not inform on the true limitations of adaptive evolution. We argue that this dependency of genetic constraints on the environment will have important effects on the progress of adaptation in heterogeneous environments, and will affect our ability to establish realistic genealogies from the phylogenic record.

Many applications in biosensing, biomaterial engineering and single molecule biophysics require multiple non-covalent linkages between DNA, protein molecules, and surfaces that are specific yet strong. Here, we present a novel method to join proteins and dsDNA molecule at their ends, in an efficient, rapid and specific manner, based on the recently developed linkage between the protein StrepTactin (STN) and the peptide StrepTag II (ST). We introduce a two-step approach, in which we first construct a hybrid between DNA and a tandem of two STs peptides (tST). In a second step, this hybrid is linked to polystyrene bead surfaces and Maltose Binding Protein (MBP) using STN. Furthermore, we show the STN-tST linkage is more stable against forces applied by optical tweezers than the commonly used biotin-Streptavidin (STV) linkage. It can be used in conjunction with Neutravidin (NTV)-biotin linkages to form DNA tethers that can sustain applied forces above 65 pN for tens of minutes in a quarter of the cases. The method is general and can be applied to construct other surface-DNA and protein-DNA hybrids. The reversibility, high mechanical stability and specificity provided by this linking procedure make it highly suitable for single molecule mechanical studies, as well as biosensing and lab on chip applications.

Trigger factor (TF) is a chaperone, found in bacterial cells and chloroplasts, that interacts with nascent polypeptide chains to suppress aggregation. While its crystal structure has been resolved, the solution structure and dynamics are largely unknown. We performed multiple molecular dynamics simulations on Trigger factor in solution, and show that its tertiary domains display collective motions hinged about inter-domain linkers with minimal or no loss in secondary structure. Moreover, we find that isolated TF typically adopts a collapsed state, with the formation of domain pairs. This collapse of TF in solution is induced by hydrophobic interactions and stabilised by hydrophilic contacts. To determine the nature of the domain interactions, we analysed the hydrophobicity of the domain surfaces by using the hydrophobic probe method of Acharya et al. [1], [2], as the standard hydrophobicity scales predictions are limited due to the complex environment. We find that the formation of domain pairs changes the hydrophobic map of TF, making the N-terminal and arm2 domain pair more hydrophilic and the head and arm1 domain pair more hydrophobic. These insights into the dynamics and interactions of the TF domains are important to eventually understand chaperone-substrate interactions and chaperone function.

Stochasticity in gene regulation has been characterized extensively, but how it affects cellular growth and fitness is less clear. We study the growth of E. coli cells as they shift from glucose to lactose metabolism, which is characterized by an obligatory growth arrest in bulk experiments that is termed the lag phase. Here, we follow the growth dynamics of individual cells at minute-resolution using a single-cell assay in a microfluidic device during this shift, while also monitoring lac expression. Mirroring the bulk results, the majority of cells displays a growth arrest upon glucose exhaustion, and resume when triggered by stochastic lac expression events. However, a significant fraction of cells maintains a high rate of elongation and displays no detectable growth lag during the shift. This ability to suppress the growth lag should provide important selective advantages when nutrients are scarce. Trajectories of individual cells display a highly non-linear relation between lac expression and growth, with only a fraction of fully induced levels being sufficient for achieving near maximal growth. A stochastic molecular model together with measured dependencies between nutrient concentration, lac expression level, and growth accurately reproduces the observed switching distributions. The results show that a growth arrest is not obligatory in the classic diauxic shift, and underscore that regulatory stochasticity ought to be considered in terms of its impact on growth and survival.

Single molecule studies on membrane proteins embedded in their native environment are hampered by the intrinsic difficulty of immobilizing elastic and sensitive biological membranes without interfering with protein activity. Here, we present hydrogels composed of nano-scaled fibers as a generally applicable tool to immobilize biological membrane vesicles of various size and lipid composition. Importantly, membrane proteins immobilized in the hydrogel as well as soluble proteins are fully active. The triggered opening of the mechanosensitive channel of large conductance (MscL) reconstituted in giant unilamellar vesicles (GUVs) was followed in time on single GUVs. Thus, kinetic studies of vectorial transport processes across biological membranes can be assessed on single, hydrogel immobilized, GUVs. Furthermore, protein translocation activity by the membrane embedded protein conducting channel of bacteria, SecYEG, in association with the soluble motor protein SecA was quantitatively assessed in bulk and at the single vesicle level in the hydrogel. This technique provides a new way to investigate membrane proteins in their native environment at the single molecule level by means of fluorescence microscopy.

SecB is a molecular chaperone in Gram-negative bacteria dedicated to the post-translational translocation of proteins across the cytoplasmic membrane. The entire surface of this chaperone is used for both of its native functions in protein targeting and unfolding. Single molecule studies revealed how SecB affects the folding pathway of proteins and how it prevents the tertiary structure formation and aggregation to support protein translocation.

The molecular chaperone SecB binds to hydrophobic sections of unfolded secretory proteins and thereby prevents their premature folding prior to secretion by the translocase of Escherichia coli. Here, we have Investigated the effect of the single-residue mutation of leucine 42 to arginine (L42R) centrally positioned in the polypeptide binding pocket of SecB on its chaperonin function. The mutant retains its tetrameric structure and SecA targeting function but is defective in its holdase activity. Isothermal titration calorimetry and single-molecule optical tweezer studies suggest that the SecB(L42R) mutant exhibits a reduced polypeptide binding affinity allowing for partial folding of the bound polypeptide chain rendering it translocation-incompetent.

The concept of epistasis has since long been used to denote non-additive fitness effects of genetic changes and has played a central role in understanding the evolution of biological systems. Owing to an array of novel experimental methodologies, it has become possible to experimentally determine epistatic interactions as well as more elaborate genotype-fitness maps. These data have opened up the investigation of a host of long-standing questions in evolutionary biology, such as the ruggedness of fitness landscapes and the accessibility of mutational trajectories, the evolution of sex, and the origin of robustness and modularity. Here we review this recent and timely marriage between systems biology and evolutionary biology, which holds the promise to understand evolutionary dynamics in a more mechanistic and predictive manner.

The defining feature of bacterial phase variation is a stochastic 'all-or-nothing' switching in gene expression. However, direct observations of these rare switching events have so far been lacking, obscuring possible correlations between switching events themselves, and between switching and other cellular events, such as division and DNA replication. We monitored the phase variation of type 1 fimbriae in individual Escherichia coli in real time and simultaneously tracked the chromosome replication process. We observed distinctive patterns of fim (fimbriae) expression in multiple genealogically related lineages. These patterns could be explained by a model that combines a single switching event with chromosomal fim replication, as well as the epigenetic inheritance of expressed fim protein and RNA, and their dilution by growth. Analysis of the moment of switching at sub-cell-cycle resolution revealed a correlation between fim switching and cell age, which challenges the traditional idea of phase variation as a random Poissonian phenomenon.

In Escherichia coli, secretory proteins (preproteins) are translocated across the cytoplasmic membrane by the Sec system composed of a protein-conducting channel, SecYEG, and an ATP-dependent motor protein, SecA. After binding of the preprotein to SecYEG-bound SecA, cycles of ATP binding and hydrolysis by SecA are thought to drive the stepwise translocation of the preprotein across the membrane. To address how the length of a preprotein substrate affects the SecA-driven translocation process, we constructed derivatives of the precursor of the outer membrane protein A ( proOmpA) with 2, 4, 6, and 8 in-tandem repeats of the periplasmic domain. With increasing polypeptide length, an increasing delay in the time before full-length translocation was observed, but the translocation rate expressed as amino acid translocation per minute remained constant. These data indicate that in the ATP-dependent reaction, SecA drives a constant rate of preprotein translocation consistent with a stepping mechanism of translocation.