Rate of translocation of bacteriophage T7 DNA across the membranes of Escherichia coli

Phage against the Machine: The SIE-ence of Superinfection Exclusion

Prophages can alter their bacterial hosts to prevent other phages from infecting the same cell, a mechanism known as superinfection exclusion (SIE). Such alterations are facilitated by phage interactions with critical bacterial components involved in motility, adhesion, biofilm production, conjugation, antimicrobial resistance, and immune evasion. Therefore, the impact of SIE extends beyond the immediate defense against superinfection, influencing the overall fitness and virulence of the bacteria. Evaluating the interactions between phages and their bacterial targets is critical for leading phage therapy candidates like Pseudomonas aeruginosa, a Gram-negative bacterium responsible for persistent and antibiotic-resistant opportunistic infections. However, comprehensive literature on the mechanisms underlying SIE remains scarce. Here, we provide a compilation of well-characterized and potential mechanisms employed by Pseudomonas phages to establish SIE. We hypothesize that the fitness costs imposed by SIE affect bacterial virulence, highlighting the potential role of this mechanism in the management of bacterial infections.

Coinfecting phages impede each other's entry into the cell

The developmental choice made by temperate phages, between cell death (lysis) and viral dormancy (lysogeny), is influenced by the relative abundance of viruses and hosts in the environment. The paradigm for this abundance-driven decision is phage lambda of E. coli, whose propensity to lysogenize increases with the number of viruses coinfecting the same bacterium. It is believed that lambda uses this number to infer whether phages or bacteria outnumber each other. However, this interpretation is premised on an accurate mapping between the extracellular phage-to-bacteria ratio and the intracellular multiplicity of infection (MOI). Here, we show this premise to be faulty. By simultaneously labeling phage capsids and genomes, we find that, while the number of phages landing on each cell reliably samples the population ratio, the number of phages entering the cell does not. Single-cell infections, performed in a microfluidic device and interpreted using a stochastic model, reveal that the probability and rate of phage entry decrease with the number of adsorbed phages. This decrease reflects an MOI-dependent perturbation to host physiology caused by phage attachment, as evidenced by compromised membrane integrity and loss of membrane potential. The dependence of entry dynamics on the surrounding medium results in a strong impact on the infection outcome, while the protracted entry of coinfecting phages increases the heterogeneity in infection outcome at a given MOI. Our findings in lambda, and similar results we obtained for phages T5 and P1, demonstrate the previously unappreciated role played by entry dynamics in determining the outcome of bacteriophage infection.

CO-INFECTING PHAGES IMPEDE EACH OTHER'S ENTRY INTO THE CELL

Bacteriophage lambda tunes its propensity to lysogenize based on the number of viral genome copies inside the infected cell. Viral self-counting is believed to serve as a way of inferring the abundance of available hosts in the environment. This interpretation is premised on an accurate mapping between the extracellular phage-to-bacteria ratio and the intracellular multiplicity of infection (MOI). However, here we show this premise to be untrue. By simultaneously labeling phage capsids and genomes, we find that, while the number of phages landing on each cell reliably samples the population ratio, the number of phages entering the cell does not. Single-cell infections, followed in a microfluidic device and interpreted using a stochastic model, reveal that the probability and rate of individual phage entries decrease with MOI. This decrease reflects an MOI-dependent perturbation to host physiology caused by phage landing, evidenced by compromised membrane integrity and loss of membrane potential. The dependence of phage entry dynamics on the surrounding medium is found to result in a strong impact of environmental conditions on the infection outcome, while the protracted entry of co-infecting phages increases the cell-to-cell variability in infection outcome at a given MOI. Our findings demonstrate the previously unappreciated role played by entry dynamics in determining the outcome of bacteriophage infection.

Engineering a Dynamic Controllable Infectivity Switch in Bacteriophage T7

Transcriptional repressors play an important role in regulating phage life cycle. Here, we examine how synthetic transcription repressors can be used in bacteriophage T7 to create a dynamic, controllable infectivity switch. We engineered T7 phage by replacing a large region of the early phage genome with different combinations of ligand-responsive promoters and ribosome binding sites (RBS) designed to control the phage RNA polymerase, gp1. Phages with engineered infectivity switch are fully viable at levels comparable to wildtype T7, when not repressed, indicating the phage can be engineered without loss of fitness. The most effective switch used a TetR-responsive promoter and an attenuated RBS, resulting in a 2-fold increase in latent period and a 10-fold decrease in phage titer when repressed. Phage activity can be further tuned using different inducer concentrations. Our study provides a proof of concept for how a simple synthetic circuit introduced into the phage genome enables user control over phage infectivity.

Isolation and characterization of a novel Escherichia coli Kayfunavirus phage DY1

Phage therapy has been revitalized since antibiotic resistance in bacteria is increasing. Compared with antibiotics, phages can target specific bacteria precisely, which requires more understanding of phage-host interactions by investigating different phages. Escherichia coli is a common pathogen with very high diversity. Based on the O antigens, E. coli can be classified into at least 183 serotypes and existing phages are far from being able to lyse all E. coli. Therefore, a novel phage specific to E. coli, named DY1, was identified and characterized to enhance our understanding of the phages of E. coli and expand the phage library. Phage DY1 belongs to the family Autographiviridae which is derived from Podoviridae. The genome of DY1 was determined to be 39,817 bp and comprises 54 putative open reading frames. Comparative genome and phylogenetic analysis demonstrated that DY1 was highly similar to phages belonging to the genus Kayfunavirus; however, the highest average nucleotide identity (ANI) values of DY1 with known phages was 0.82 suggesting that DY1 was a novel phage. Through stability tests, DY1 was very stable at temperatures ranging from 20 to 50 °C and pH levels from 5 to 11. Taken together, we report that phage DY1 is a novel Kayfunavirus phage with the potential for phage therapy.

Mutagenic Analysis of a DNA Translocating Tube's Interior Surface

Bacteriophage ϕX174 uses a decamer of DNA piloting proteins to penetrate its host. These proteins oligomerize into a cell wall-spanning tube, wide enough for genome passage. While the inner surface of the tube is primarily lined with inward-facing amino acid side chains containing amide and guanidinium groups, there is a 28 Å-long section near the tube’s C-terminus that does not exhibit this motif. The majority of the inward-facing residues in this region are conserved across the three ϕX174-like clades, suggesting that they play an important role during genome delivery. To test this hypothesis, and explore the general function of the tube’s inner surface, non-glutamine residues within this region were mutated to glutamine, while existing glutamine residues were changed to serine. Four of the resulting mutants had temperature-dependent phenotypes. Virion assembly, host attachment, and virion eclipse, defined as the cell’s ability to inactivate the virus, were not affected. Genome delivery, however, was inhibited. The results support a model in which a balance of forces governs genome delivery: potential energy provided by the densely packaged viral genome and/or an osmotic gradient move the genome into the cell, while the tube’s inward facing glutamine residues exert a frictional force, or drag, that controls genome release.

Hitchhiking, collapse, and contingency in phage infections of migrating bacterial populations

Natural bacterial populations are subjected to constant predation pressure by bacteriophages. Bacteria use a variety of molecular mechanisms to defend themselves from phage predation. However, since phages are nonmotile, perhaps the simplest defense against phage is for bacteria to move faster than phages. In particular, chemotaxis, the active migration of bacteria up attractant gradients, may help the bacteria escape slowly diffusing phages. Here we study phage infection dynamics in migrating bacterial populations driven by chemotaxis through low viscosity agar plates. We find that expanding phage-bacteria populations supports two moving fronts, an outermost bacterial front driven by nutrient uptake and chemotaxis and an inner phage front at which the bacterial population collapses due to phage predation. We show that with increasing adsorption rate and initial phage population, the speed of the moving phage front increases, eventually overtaking the bacterial front and driving the system across a transition from a regime where bacterial front speed exceeds that of the phage front to one where bacteria must evolve phage resistance to survive. Our data support the claim that this process requires phage to hitchhike with moving bacteria. A deterministic model recapitulates the transition under the assumption that phage virulence declines with host growth rate which we confirm experimentally. Finally, near the transition between regimes we observe macroscopic fluctuations in bacterial densities at the phage front. Our work opens a new, spatio-temporal, line of investigation into the eco-evolutionary struggle between bacteria and phage.

Transcript degradation and codon usage regulate gene expression in a lytic phage

Many viral genomes are small, containing only single- or double-digit numbers of genes and relatively few regulatory elements. Yet viruses successfully execute complex regulatory programs as they take over their host cells. Here, we propose that some viruses regulate gene expression via a carefully balanced interplay between transcription, translation, and transcript degradation. As our model system, we employ bacteriophage T7, whose genome of approximately sixty genes is well annotated and for which there is a long history of computational models of gene regulation. We expand upon prior modeling work by implementing a stochastic gene expression simulator that tracks individual transcripts, polymerases, ribosomes, and ribonucleases participating in the transcription, translation, and transcript-degradation processes occurring during a T7 infection. By combining this detailed mechanistic modeling of a phage infection with high-throughput gene expression measurements of several strains of bacteriophage T7, evolved and engineered, we can show that both the dynamic interplay between transcription and transcript degradation, and between these two processes and translation, appear to be critical components of T7 gene regulation. Our results point to targeted degradation as a generic gene regulation strategy that may have evolved in many other viruses. Further, our results suggest that detailed mechanistic modeling may uncover the biological mechanisms at work in both evolved and engineered virus variants.

Time-resolved DNA release from an O-antigen-specific Salmonella bacteriophage with a contractile tail

Myoviruses, bacteriophages with T4-like architecture, must contract their tails prior to DNA release. However, quantitative kinetic data on myovirus particle opening are lacking, although they are promising tools in bacteriophage-based antimicrobial strategies directed against Gram-negative hosts. For the first time, we show time-resolved DNA ejection from a bacteriophage with a contractile tail, the multi-O-antigen-specific Salmonella myovirus Det7. DNA release from Det7 was triggered by lipopolysaccharide (LPS) O-antigen receptors and notably slower than in noncontractile-tailed siphoviruses. Det7 showed two individual kinetic steps for tail contraction and particle opening. Our in vitro studies showed that highly specialized tailspike proteins (TSPs) are necessary to attach the particle to LPS. A P22-like TSP confers specificity for the Salmonella Typhimurium O-antigen. Moreover, crystal structure analysis at 1.63 Å resolution confirmed that Det7 recognized the Salmonella Anatum O-antigen via an ϵ15-like TSP, DettilonTSP. DNA ejection triggered by LPS from either host showed similar velocities, so particle opening is thus a process independent of O-antigen composition and the recognizing TSP. In Det7, at permissive temperatures TSPs mediate O-antigen cleavage and couple cell surface binding with DNA ejection, but no irreversible adsorption occurred at low temperatures. This finding was in contrast to short-tailed Salmonella podoviruses, illustrating that tailed phages use common particle-opening mechanisms but have specialized into different infection niches.

Identification of a bacteriophage from an environmental multidrug-resistant E. coli isolate and its function in horizontal transfer of ARGs

Horizontal transfer of ARGs was generally considered to be mediated by three methods - transformation, conjugation and transduction through phages - during which the contribution of bacteriophages to gene transfer in the environment is unclear or even questioned. In this study, a multiple-antibiotic-resistant Escherichia coli strain and its phage (YZ1) were isolated from a municipal wastewater treatment system. The results of the morphological and genomic analyses of phage YZ1 showed that it is a member of the T7 viral genus in the subfamily Autographivirinae. Its genome is similar to that of the E. coli phage K1F in both organization and sequence and does not encode ARGs. However, 28 paired reads in the raw sequencing data aligned to ARGs, including those promoting β-lactam, aminoglycoside, and fluoroquinolone resistance, among others. Quantitative PCR showed that ARGs were present in bacteriophage DNA (approximately 10³ copies/mL) and were also detected in the bacterial host DNA. The results suggested that while infrequent, some ARG-carrying transducing phages were presumably generated by erroneous packaging during infection of antibiotic-resistant bacteria, which may create the possibility of horizontal transfer of ARGs.

Predicting Evolution of the Transcription Regulatory Network in a Bacteriophage

Prediction of evolutionary trajectories has been an elusive goal, requiring a deep knowledge of underlying mechanisms that relate genotype to phenotype plus understanding how phenotype impacts organismal fitness. We tested our ability to predict molecular regulatory evolution in a bacteriophage (T7) whose RNA polymerase (RNAP) was altered to recognize a heterologous promoter differing by three nucleotides from the wild-type promoter. A mutant of wild-type T7 lacking its RNAP gene was passaged on a bacterial strain providing the novel RNAP in trans. Higher fitness rapidly evolved. Predicting the evolutionary trajectory of this adaptation used measured in vitro transcription rates of the novel RNAP on the six promoter sequences capturing all possible one-step pathways between the wild-type and the heterologous promoter sequences. The predictions captured some of the regulatory evolution but failed both in explaining 1) a set of T7 promoters that consistently failed to evolve and 2) some promoter evolution that fell outside the expected one-step pathways. Had a more comprehensive set of transcription assays been undertaken initially, all promoter evolution would have fallen within predicted bounds, but the lack of evolution in some promoters is unresolved. Overall, this study points toward the increasing feasibility of predicting evolution in well-characterized, simple systems.

Role of the vacuolar ATPase in the Alphavirus replication cycle

We have shown that Alphaviruses can enter cells by direct penetration at the plasma membrane (R. Vancini, G. Wang, D. Ferreira, R. Hernandez, and D. Brown, J Virol, 87:4352–4359, 2013). Direct penetration removes the requirement for receptor-mediated endocytosis exposure to low pH and membrane fusion in the process of RNA entry. Endosomal pH as well as the pH of the cell cytoplasm is maintained by the activity of the vacuolar ATPase (V-ATPase). Bafilomycin is a specific inhibitor of V-ATPase. To characterize the roll of the V-ATPase in viral replication we generated a Bafilomycin A1(BAF) resistant mutant of Sindbis virus (BRSV). BRSV produced mature virus and virus RNA in greater amounts than parent virus in BAF-treated cells. Sequence analysis revealed mutations in the E2 glycoprotein, T15I/Y18H, were responsible for the phenotype. These results show that a functional V-ATPase is required for efficient virus RNA synthesis and virus maturation in Alphavirus infection.

CRISPR-based herd immunity can limit phage epidemics in bacterial populations

Herd immunity, a process in which resistant individuals limit the spread of a pathogen among susceptible hosts has been extensively studied in eukaryotes. Even though bacteria have evolved multiple immune systems against their phage pathogens, herd immunity in bacteria remains unexplored. Here we experimentally demonstrate that herd immunity arises during phage epidemics in structured and unstructured Escherichia coli populations consisting of differing frequencies of susceptible and resistant cells harboring CRISPR immunity. In addition, we develop a mathematical model that quantifies how herd immunity is affected by spatial population structure, bacterial growth rate, and phage replication rate. Using our model we infer a general epidemiological rule describing the relative speed of an epidemic in partially resistant spatially structured populations. Our experimental and theoretical findings indicate that herd immunity may be important in bacterial communities, allowing for stable coexistence of bacteria and their phages and the maintenance of polymorphism in bacterial immunity.

Late-Arriving Signals Contribute Less to Cell-Fate Decisions

Gene regulatory networks are largely responsible for cellular decision-making. These networks sense diverse external signals and respond by adjusting gene expression, enabling cells to reach environment-dependent decisions crucial for their survival or reproduction. However, information-carrying signals may arrive at variable times. Besides the intrinsic strength of these signals, their arrival time (timing) may also carry information about the environment and can influence cellular decision-making in ways that are poorly understood. For example, it is unclear how the timing of individual phage infections affects the lysis-lysogeny decision of bacteriophage λ despite variable infection times being likely in the wild and even in laboratory conditions. In this work, we combine mathematical modeling with experimentation to address this question. We develop an experimentally testable theory, which reveals that late-infecting phages contribute less to cellular decision-making. This implies that infection delays lower the probability of lysogeny compared to simultaneous infections. Furthermore, we show that infection delays reduce lysogenization by providing insufficient CII for threshold crossing during the critical decision-making period. We find evidence for a cutoff time after which subsequent infections cannot influence the cellular decision. We derive an intuitive formula that approximates the probability of lysogeny for variable infection times by a time-weighted average of probabilities for simultaneous infections. We validate these theoretical predictions experimentally. Similar concepts and simplifying modeling approaches may help elucidate the mechanisms underlying other cellular decisions.

A non-invasive method for studying viral DNA delivery to bacteria reveals key requirements for phage SPP1 DNA entry in Bacillus subtilis cells

Bacteriophages use most frequently a tail apparatus to create a channel across the entire bacterial cell envelope to transfer the viral genome to the host cell cytoplasm, initiating infection. Characterization of this critical step remains a major challenge due to the difficulty to monitor DNA entry in the bacterium and its requirements. In this work we developed a new method to study phage DNA entry that has the potential to be extended to many tailed phages. Its application to study genome delivery of bacteriophage SPP1 into Bacillus subtilis disclosed a key role of the host cell membrane potential in the DNA entry process. An energized B. subtilis membrane and a millimolar concentration of calcium ions are shown to be major requirements for SPP1 DNA entry following the irreversible binding of phage particles to the receptor YueB.

The T7 ejection nanomachine components gp15-gp16 form a spiral ring complex that binds DNA and a lipid membrane

Bacteriophage T7 initiates infection by ejecting several internal capsid proteins into the host cell; these proteins then assemble into a nanomachine that translocates the viral genome from the phage head into the cytoplasm. The ejected proteins are thought to partially unfold as they pass through the lumen of the portal and the short stubby T7 tail during their entry into the cell. In vivo, the internal proteins gp15 and gp16 assemble into a tubular structure that spans the periplasm and cytoplasmic membrane. We show here that purified gp15 and gp16 can refold from a partially denatured state in vitro, and that gp15 interacts with gp16 to form a spiral ring structure. Purified gp15 binds to DNA, whereas gp16 binds protein-free liposomes; the gp15-gp16 complex binds both DNA and liposomes. Limited proteolysis of the liposome-bound gp16 reveals that its C-terminal region is protected, suggesting a partial membrane insertion of the protein.

Bacteriophage P22 ejects all of its internal proteins before its genome

Double-stranded DNA bacteriophages are highly pressurized, providing a force driving ejection of a significant fraction of the genome from its capsid. In P22-like Podoviridae, internal proteins ("E proteins") are packaged into the capsid along with the genome, and without them the virus is not infectious. However, little is known about how and when these proteins come out of the virus. We employed an in vitro osmotic suppression system with high-molecular-weight polyethylene glycol to study P22 E protein release. While slow ejection of the DNA can be triggered by lipopolysaccharide (LPS), the rate is significantly enhanced by the membrane protein OmpA from Salmonella. In contrast, E proteins are not ejected unless both OmpA and LPS are present and their ejection when OmpA is present is largely complete before any genome is ejected, suggesting that E proteins play a key role in the early stage of transferring P22 DNA into the host.

Synergy as a rationale for phage therapy using phage cocktails

Where phages are used to treat bacterial contaminations and infections, multiple phages are typically applied at once as a cocktail. When two or more phages in the cocktail attack the same bacterium, the combination may produce better killing than any single phage (synergy) or the combination may be worse than the best single phage (interference). Synergy is of obvious utility, especially if it can be predicted a priori, but it remains poorly documented with few examples known. This study addresses synergy in which one phage improves adsorption by a second phage. It first presents evidence of synergy from an experimental system of two phages and a mucoid E. coli host. The synergy likely stems from a tailspike enzyme produced by one of the phages. We then offer mathematical models and simulations to understand the dynamics of synergy and the enhanced magnitude of bacterial control possible. The models and observations complement each other and suggest that synergy may be of widespread utility and may be predictable from easily observed phenotypes.

Challenges in predicting the evolutionary maintenance of a phage transgene

BACKGROUND

In prior work, a phage engineered with a biofilm-degrading enzyme (dispersin B) cleared artificial, short-term biofilms more fully than the phage lacking the enzyme. An unresolved question is whether the transgene will be lost or maintained during phage growth - its loss would limit the utility of the engineering. Broadly supported evolutionary theory suggests that transgenes will be lost through a 'tragedy of the commons' mechanism unless the ecology of growth in biofilms meets specific requirements. We test that theory here.

RESULTS

Functional properties of the transgenic phage were identified. Consistent with the previous study, the dispersin phage was superior to unmodified phage at clearing short term biofilms grown in broth, shown here to be an effect attributable to free enzyme. However, the dispersin phage was only marginally better than control phages on short term biofilms in minimal media and was no better than control phages in clearing long term biofilms. There was little empirical support for the tragedy of the commons framework despite a strong theoretical foundation for its supposed relevance. The framework requires that the transgene imposes an intrinsic cost, yet the transgene was intrinsically neutral or beneficial when expressed from one part of the phage genome. Expressed from a different part of the genome, the transgene did behave as if intrinsically costly, but its maintenance did not benefit from spatially structured growth per se - violating the tragedy framework.

CONCLUSIONS

Overall, the transgene was beneficial under many conditions, but no insight to its maintenance was attributable to the established evolutionary framework. The failure likely resides in system details that would be used to parameterize the models. Our study cautions against naive applications of evolutionary theory to synthetic biology, even qualitatively.

The tripartite virions of the brome mosaic virus have distinct physical properties that affect the timing of the infection process

The three subsets of virions that comprise the Brome mosaic virus (BMV) were previously thought to be indistinguishable. This work tested the hypothesis that distinct capsid-RNA interactions in the BMV virions allow different rates of viral RNA release. Several results support distinct interactions between the capsid and the BMV genomic RNAs. First, the deletion of the first eight residues of the BMV coat protein (CP) resulted in the RNA1-containing particles having altered morphologies, while those containing RNA2 were unaffected. Second, subsets of the BMV particles separated by density gradients into a pool enriched for RNA1 (B1) and for RNA2 and RNA3/4 (B2.3/4) were found to have different physiochemical properties. Compared to the B2.3/4 particles, the B1 particles were more sensitive to protease digestion and had greater resistivity to nanoindentation by atomic force microscopy and increased susceptibility to nuclease digestion. Mapping studies showed that portions of the arginine-rich N-terminal tail of the CP could interact with RNA1. Mutational analysis in the putative RNA1-contacting residues severely reduced encapsidation of BMV RNA1 without affecting the encapsidation of RNA2. Finally, during infection of plants, the more easily released RNA1 accumulated to higher levels early in the infection.

IMPORTANCE

Viruses with genomes packaged in distinct virions could theoretically release the genomes at different times to regulate the timing of gene expression. Using an RNA virus composed of three particles, we demonstrated that the RNA in one of the virions is released more easily than the other two in vitro. The differential RNA release is due to distinct interactions between the viral capsid protein and the RNAs. The ease of RNA release is also correlated with the more rapid accumulation of that RNA in infected plants. Our study identified a novel role for capsid-RNA interactions in the regulation of a viral infection.

Langevin dynamics simulation of DNA ejection from a phage

We have performed Langevin dynamics simulations of a coarse-grained model of ejection of dsDNA from Φ29 phage. Our simulation results show significant variations in the local ejection speed, consistent with experimental observations reported in the literature for both in vivo and in vitro systems. In efforts to understand the origin of such variations in the local speed of ejection, we have investigated the correlations between the local ejection kinetics and the packaged structures created at various motor forces and chain flexibility. At lower motor forces, the packaged DNA length is shorter with better organization. On the other hand, at higher motor forces typical of realistic situations, the DNA organization inside the capsid suffers from significant orientational disorder, but yet with long orientational correlation times. This in turn leads to lack of registry between the direction of the DNA segments just to be ejected and the direction of exit. As a result, a significant amount of momentum transfer is required locally for successful exit. Consequently, the DNA ejection temporarily slows down exhibiting pauses. This slowing down occurs at random times during the ejection process, completely determined by the particular starting conformation created by prescribed motor forces. In order to augment our inference, we have additionally investigated the ejection of chains with deliberately changed persistence length. For less inflexible chains, the demand on the occurrence of large momentum transfer for successful ejection is weaker, resulting in more uniform ejection kinetics. While being consistent with experimental observations, our results show the nonergodic nature of the ejection kinetics and call for better theoretical models to portray the kinetics of genome ejection from phages.

A non-bacterial transcription factor inhibits bacterial transcription by a multipronged mechanism

The process of transcription initiation is the major target for regulation of gene expression in bacteria and is performed by a multi-subunit RNA polymerase enzyme (RNAp). A complex network of regulatory elements controls the activity of the RNAp to fine-tune transcriptional output. Thus, RNAp is a nexus for controlling bacterial gene expression at the transcription level. Many bacteriophages, viruses that infect bacteria, encode transcription factors that specifically target and modulate the activity of the host RNAp and, thereby, facilitate the acquisition of the host bacteria by the phage. Here, we describe the modus operandi of a T7 bacteriophage-encoded small protein called Gp2 and define Gp2 as a non-bacterial regulator of bacterial transcription.

Effects of pulling forces, osmotic pressure, condensing agents and viscosity on the thermodynamics and kinetics of DNA ejection from bacteriophages to bacterial cells: a computational study

In this work, we report on simulations of double-stranded DNA (dsDNA) ejection from bacteriophage φ29 into a bacterial cell. The ejection was studied with a coarse-grained model, in which viral dsDNA was represented by beads on a torsion-less string. The bacteriophage's capsid and the bacterial cell were defined by sets of spherical constraints. To account for the effects of the viscous medium inside the bacterial cell, the simulations were carried out using a Langevin dynamics protocol. Our simplest simulations (involving constant viscosity and no external biasing forces) produced results compatible with the push-pull model of DNA ejection, with an ejection rate significantly higher in the first part of ejection than in the latter parts. Additionally, we performed more complicated simulations, in which we included additional factors such as external forces, osmotic pressure, condensing agents and ejection-dependent viscosity. The effects of these factors (independently and in combination) on the thermodynamics and kinetics of DNA ejection were studied. We found that, in general, the dependence of ejection forces and ejection rates on the amount of DNA ejected becomes more complex if the ejection is modeled with a broader, more realistic set of parameters and influences (such as variation in the solvent's viscosity and the application of an external force). However, certain combinations of factors and numerical parameters led to the opposition of some ejection-driving and ejection-inhibiting influences, ultimately causing an apparent simplification of the ejection profiles.

Popping the cork: mechanisms of phage genome ejection

Sixty years after Hershey and Chase showed that nucleic acid is the major component of phage particles that is ejected into cells, we still do not fully understand how the process occurs. Advances in electron microscopy have revealed the structure of the condensed DNA confined in a phage capsid, and the mechanisms and energetics of packaging a phage genome are beginning to be better understood. Condensing DNA subjects it to high osmotic pressure, which has been suggested to provide the driving force for its ejection during infection. However, forces internal to a phage capsid cannot, alone, cause complete genome ejection into cells. Here, we describe the structure of the DNA inside mature phages and summarize the current models of genome ejection, both in vitro and in vivo.

The bacteriophage t7 virion undergoes extensive structural remodeling during infection

Adsorption and genome ejection are fundamental to the bacteriophage life cycle, yet their molecular mechanisms are not well understood. We used cryo-electron tomography to capture T7 virions at successive stages of infection of Escherichia coli minicells at ~4-nm resolution. The six phage tail fibers were folded against the capsid, extending and orienting symmetrically only after productive adsorption to the host cell surface. Receptor binding by the tail triggered conformational changes resulting in the insertion of an extended tail, which functions as the DNA ejection conduit into the cell cytoplasm. After ejection, the extended phage tail collapsed or disassembled, which allowed resealing of the infected cell membrane. These structural studies provide a detailed series of intermediates during phage infection.

Engineered cell-cell communication via DNA messaging

BACKGROUND

Evolution has selected for organisms that benefit from genetically encoded cell-cell communication. Engineers have begun to repurpose elements of natural communication systems to realize programmed pattern formation and coordinate other population-level behaviors. However, existing engineered systems rely on system-specific small molecules to send molecular messages among cells. Thus, the information transmission capacity of current engineered biological communication systems is physically limited by specific biomolecules that are capable of sending only a single message, typically "regulate transcription."

RESULTS

We have engineered a cell-cell communication platform using bacteriophage M13 gene products to autonomously package and deliver heterologous DNA messages of varying lengths and encoded functions. We demonstrate the decoupling of messages from a common communication channel via the autonomous transmission of various arbitrary genetic messages. Further, we increase the range of engineered DNA messaging across semisolid media by linking message transmission or receipt to active cellular chemotaxis.

CONCLUSIONS

We demonstrate decoupling of a communication channel from message transmission within engineered biological systems via the autonomous targeted transduction of user-specified heterologous DNA messages. We also demonstrate that bacteriophage M13 particle production and message transduction occurs among chemotactic bacteria. We use chemotaxis to improve the range of DNA messaging, increasing both transmission distance and communication bit rates relative to existing small molecule-based communication systems. We postulate that integration of different engineered cell-cell communication platforms will allow for more complex spatial programming of dynamic cellular consortia.

A single-molecule Hershey-Chase experiment

Ever since Hershey and Chase used phages to establish DNA as the carrier of genetic information in 1952, the precise mechanisms of phage DNA translocation have been a mystery. Although bulk measurements have set a timescale for in vivo DNA translocation during bacteriophage infection, measurements of DNA ejection by single bacteriophages have only been made in vitro. Here, we present direct visualization of single bacteriophages infecting individual Escherichia coli cells. For bacteriophage λ, we establish a mean ejection time of roughly 5 min with significant cell-to-cell variability, including pausing events. In contrast, corresponding in vitro single-molecule ejections are more uniform and finish within 10 s. Our data reveal that when plotted against the amount of DNA ejected, the velocity of ejection for two different genome lengths collapses onto a single curve. This suggests that in vivo ejections are controlled by the amount of DNA ejected. In contrast, in vitro DNA ejections are governed by the amount of DNA left inside the capsid. This analysis provides evidence against a purely intrastrand repulsion-based mechanism and suggests that cell-internal processes dominate. This provides a picture of the early stages of phage infection and sheds light on the problem of polymer translocation.

Lipid-containing viruses: bacteriophage PRD1 assembly

PRD1 is a tailless icosahedrally symmetric virus containing an internal lipid membrane beneath the protein capsid. Its linear dsDNA genome and covalently attached terminal proteins are delivered into the cell where replication occurs via a protein-primed mechanism. Extensive studies have been carried out to decipher the roles of the 37 viral proteins in PRD1 assembly, their association in virus particles and lately, especially the functioning of the unique packaging machinery that translocates the genome into the procapsid. These issues will be addressed in this chapter especially in the context of the structure of PRD1. We will also discuss the major challenges still to be addressed in PRD1 assembly.

Bacteriophage protein-protein interactions

Bacteriophages T7, λ, P22, and P2/P4 (from Escherichia coli), as well as ϕ29 (from Bacillus subtilis), are among the best-studied bacterial viruses. This chapter summarizes published protein interaction data of intraviral protein interactions, as well as known phage-host protein interactions of these phages retrieved from the literature. We also review the published results of comprehensive protein interaction analyses of Pneumococcus phages Dp-1 and Cp-1, as well as coliphages λ and T7. For example, the ≈55 proteins encoded by the T7 genome are connected by ≈43 interactions with another ≈15 between the phage and its host. The chapter compiles published interactions for the well-studied phages λ (33 intra-phage/22 phage-host), P22 (38/9), P2/P4 (14/3), and ϕ29 (20/2). We discuss whether different interaction patterns reflect different phage lifestyles or whether they may be artifacts of sampling. Phages that infect the same host can interact with different host target proteins, as exemplified by E. coli phage λ and T7. Despite decades of intensive investigation, only a fraction of these phage interactomes are known. Technical limitations and a lack of depth in many studies explain the gaps in our knowledge. Strategies to complete current interactome maps are described. Although limited space precludes detailed overviews of phage molecular biology, this compilation will allow future studies to put interaction data into the context of phage biology.

Dynamics of bacteriophage genome ejection in vitro and in vivo

Bacteriophages, phages for short, are viruses of bacteria. The majority of phages contain a double-stranded DNA genome packaged in a capsid at a density of ∼500 mg ml(-1). This high density requires substantial compression of the normal B-form helix, leading to the conjecture that DNA in mature phage virions is under significant pressure, and that pressure is used to eject the DNA during infection. A large number of theoretical, computer simulation and in vitro experimental studies surrounding this conjecture have revealed many--though often isolated and/or contradictory--aspects of packaged DNA. This prompts us to present a unified view of the statistical physics and thermodynamics of DNA packaged in phage capsids. We argue that the DNA in a mature phage is in a (meta)stable state, wherein electrostatic self-repulsion is balanced by curvature stress due to confinement in the capsid. We show that in addition to the osmotic pressure associated with the packaged DNA and its counterions, there are four different pressures within the capsid: pressure on the DNA, hydrostatic pressure, the pressure experienced by the capsid and the pressure associated with the chemical potential of DNA ejection. Significantly, we analyze the mechanism of force transmission in the packaged DNA and demonstrate that the pressure on DNA is not important for ejection. We derive equations showing a strong hydrostatic pressure difference across the capsid shell. We propose that when a phage is triggered to eject by interaction with its receptor in vitro, the (thermodynamic) incentive of water molecules to enter the phage capsid flushes the DNA out of the capsid. In vivo, the difference between the osmotic pressures in the bacterial cell cytoplasm and the culture medium similarly results in a water flow that drags the DNA out of the capsid and into the bacterial cell.

The role of the T7 Gp2 inhibitor of host RNA polymerase in phage development

Bacteriophage T7 relies on its own RNA polymerase (RNAp) to transcribe its middle and late genes. Early genes, which include the viral RNAp gene, are transcribed by the host RNAp from three closely spaced strong promoters-A1, A2, and A3. One middle T7 gene product, gp2, is a strong inhibitor of the host RNAp. Gp2 is essential and is required late in infection, during phage DNA packaging. Here, we explore the role of gp2 in controlling host RNAp transcription during T7 infection. We demonstrate that in the absence of gp2, early viral transcripts continue to accumulate throughout the infection. Decreasing transcription from early promoter A3 is sufficient to make gp2 dispensable for phage infection. Gp2 also becomes dispensable when an antiterminating element boxA, located downstream of early promoters, is deleted. The results thus suggest that antiterminated transcription by host RNAp from the A3 promoter is interfering with phage development and that the only essential role for gp2 is to prevent this transcription.

Intracellular kinetics of a growing virus: a genetically structured simulation for bacteriophage T7

Viruses have evolved to efficiently direct the resources of their hosts toward their own reproduction. A quantitative understanding of viral growth will help researchers develop antiviral strategies, design metabolic pathways, construct vectors for gene therapy, and engineer molecular systems that self-assemble. As a model system we examine here the growth of bacteriophage T7 in Escherichia coli using a chemical-kinetic framework. Data published over the last three decades on the genetics, physiology, and biophysics of phage T7 are incorporated into a genetically structured simulation that accounts for entry of the T7 genome into its host, expression of T7 genes, replication of T7 DNA, assembly of T7 procapsids, and packaging of T7 DNA to finally produce intact T7 progeny. Good agreement is found between the simulated behavior and experimental observations for the shift in transcription capacity from the host to the phage, the initiation times of phage protein synthesis, and the intracellular assembly of both wild-type phage and a fast-growing deletion mutant. The simulation is utilized to predict the effect of antisense molecules targeted to different T7 mRNA. Further, a postulated mechanism for the down regulation of T7 transcription in vivo is quantitatively examined and shown to agree with available data. The simulation is found to be a useful tool for exploring and understanding the dynamics of virus growth at the molecular level. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55: 375-389, 1997.

Gp15 and gp16 cooperate in translocating bacteriophage T7 DNA into the infected cell

Loss of up to four amino acids from the C terminus of the 1318 residue bacteriophage T7 gp16 allows plaque formation at normal efficiencies. Loss of five residues results in non-infective virions, and loss of twelve prevents assembly of stable particles. However, replacing the C-terminal seven with nineteen non-native residues allows assembly of non-infective virions. The latter adsorb and eject internal core proteins into the cell envelope but no phage DNA enters the cytoplasm. Extragenic suppressors of the defective gene 16 lie in gene 15; the mutant gp15 proteins not only re-establish infectivity, they fully restore the kinetics of genome internalization to those exhibited by wild-type phage. After ejection from the infecting particle, gp15 and gp16 thus function together in ratcheting the leading end of the T7 genome into the cytoplasm of the infected cell.

The force acting on a polymer partially confined in a tube

We consider the force acting on a polymer part of whose length is configurationally confined in a tube and the rest of which is free. This situation arises in many different physical contexts, including a flexible synthetic polymer partially confined in a nanopore and a stiff viral genome partially ejected from its capsid. In both cases the force acting to pull the chain molecule out of its confinement is argued to be constant once a few persistence lengths are "free"/"outside". We present Brownian dynamics simulations that confirm the constancy of the force for different chain lengths and illustrate the dependence of the force on the strength of tube confinement. Experimental results are reported for genome ejection from viral capsids, from which we estimate the pulling force to be a few tenths of a piconewton.

The role of interactions between phage and bacterial proteins within the infected cell: a diverse and puzzling interactome

Interactions between bacteriophage proteins and bacterial proteins are important for efficient infection of the host cell. The phage proteins involved in these bacteriophage-host interactions are often produced immediately after infection. A survey of the available set of published bacteriophage-host interactions reveals the targeted host proteins are inhibited, activated or functionally redirected by the phage protein. These interactions protect the bacteriophage from bacterial defence mechanisms or adapt the host-cell metabolism to establish an efficient infection cycle. Regrettably, a large majority of bacteriophage early proteins lack any identified function. Recent research into the antibacterial potential of bacteriophage-host interactions indicates that phage early proteins seem to target a wide variety of processes in the host cell - many of them non-essential. Since a clear understanding of such interactions may become important for regulations involving phage therapy and in biotechnological applications, increased scientific emphasis on the biological elucidation of such proteins is warranted.

Transport of phage P22 DNA across the cytoplasmic membrane

Although a great deal is known about the life cycle of bacteriophage P22, the mechanism of phage DNA transport into Salmonella is poorly understood. P22 DNA is initially ejected into the periplasmic space and subsequently transported into the host cytoplasm. Three phage-encoded proteins (gp16, gp20, and gp7) are coejected with the DNA. To test the hypothesis that one or more of these proteins mediate transport of the DNA across the cytoplasmic membrane, we purified gp16, gp20, and gp7 and analyzed their ability to associate with membranes and to facilitate DNA uptake into membrane vesicles in vitro. Membrane association experiments revealed that gp16 partitioned into the membrane fraction, while gp20 and gp7 remained in the soluble fraction. Moreover, the addition of gp16, but not gp7 or gp20, to liposomes preloaded with a fluorescent dye promoted release of the dye. Transport of (32)P-labeled DNA into liposomes occurred only in the presence of gp16 and an artificially created membrane potential. Taken together, these results suggest that gp16 partitions into the cytoplasmic membrane and mediates the active transport of P22 DNA across the cytoplasmic membrane of Salmonella.

TABASCO: A single molecule, base-pair resolved gene expression simulator

Background

Experimental studies of gene expression have identified some of the individual molecular components and elementary reactions that comprise and control cellular behavior. Given our current understanding of gene expression, and the goals of biotechnology research, both scientists and engineers would benefit from detailed simulators that can explicitly compute genome-wide expression levels as a function of individual molecular events, including the activities and interactions of molecules on DNA at single base pair resolution. However, for practical reasons including computational tractability, available simulators have not been able to represent genome-scale models of gene expression at this level of detail.

Results

Here we develop a simulator, TABASCO , which enables the precise representation of individual molecules and events in gene expression for genome-scale systems. We use a single molecule computational engine to track individual molecules interacting with and along nucleic acid polymers at single base resolution. Tabasco uses logical rules to automatically update and delimit the set of species and reactions that comprise a system during simulation, thereby avoiding the need for a priori specification of all possible combinations of molecules and reaction events. We confirm that single molecule, base-pair resolved simulation using TABASCO (Tabasco) can accurately compute gene expression dynamics and, moving beyond previous simulators, provide for the direct representation of intermolecular events such as polymerase collisions and promoter occlusion. We demonstrate the computational capacity of Tabasco by simulating the entirety of gene expression during bacteriophage T7 infection; for reference, the 39,937 base pair T7 genome encodes 56 genes that are transcribed by two types of RNA polymerases active across 22 promoters.

Conclusion

Tabasco enables genome-scale simulation of transcription and translation at individual molecule and single base-pair resolution. By directly representing the position and activity of individual molecules on DNA, Tabasco can directly test the effects of detailed molecular processes on system-wide gene expression. Tabasco would also be useful for studying the complex regulatory mechanisms controlling eukaryotic gene expression. The computational engine underlying Tabasco could also be adapted to represent other types of processive systems in which individual reaction events are organized across a single spatial dimension (e.g., polysaccharide synthesis).

Real-time observations of single bacteriophage lambda DNA ejections in vitro

The physical, chemical, and structural features of bacteriophage genome release have been the subject of much recent attention. Many theoretical and experimental studies have centered on the internal forces driving the ejection process. Recently, Mangenot et al. [Mangenot S, Hochrein M, Rädler J, Letellier L (2005) Curr Biol 15:430-435.] reported fluorescence microscopy of phage T5 ejections, which proceeded stepwise between DNA nicks, reaching a translocation speed of 75 kbp/s or higher. It is still unknown how high the speed actually is. This paper reports real-time measurements of ejection from phage lambda, revealing how the speed depends on key physical parameters such as genome length and ionic state of the buffer. Except for a pause before DNA is finally released, the entire 48.5-kbp genome is translocated in approximately 1.5 s without interruption, reaching a speed of 60 kbp/s. The process gives insights particularly into the effects of two parameters: a shorter genome length results in lower speed but a shorter total time, and the presence of divalent magnesium ions (replacing sodium) reduces the pressure, increasing ejection time to 8-11 s. Pressure caused by DNA-DNA interactions within the head affects the initiation of ejection, but the close packing is also the dominant source of friction: more tightly packed phages initiate ejection earlier, but with a lower initial speed. The details of ejection revealed in this study are probably generic features of DNA translocation in bacteriophages and have implications for the dynamics of DNA in other biological systems.

Viral capsids: mechanical characteristics, genome packaging and delivery mechanisms

The main functions of viral capsids are to protect, transport and deliver their genome. The mechanical properties of capsids are supposed to be adapted to these tasks. Bacteriophage capsids also need to withstand the high pressures the DNA is exerting onto it as a result of the DNA packaging and its consequent confinement within the capsid. It is proposed that this pressure helps driving the genome into the host, but other mechanisms also seem to play an important role in ejection. DNA packaging and ejection strategies are obviously dependent on the mechanical properties of the capsid. This review focuses on the mechanical properties of viral capsids in general and the elucidation of the biophysical aspects of genome packaging mechanisms and genome delivery processes of double-stranded DNA bacteriophages in particular.

Is phage DNA 'injected' into cells--biologists and physicists can agree

The double-stranded DNA inside bacteriophages is packaged at a density of approximately 500 mg/ml and exerts an osmotic pressure of tens of atmospheres. This pressure is commonly assumed to cause genome ejection during infection. Indeed, by the addition of their natural receptors, some phages can be induced in vitro to completely expel their genome from the virion. However, the osmotic pressure of the bacterial cytoplasm exerts an opposing force, making it impossible for the pressure of packaged DNA to cause complete genome ejection in vivo. Various processes for complete genome ejection are discussed, but we focus on a novel proposal suggesting that the osmotic gradient between the extracellular environment and the cytoplasm results in fluid flow through the phage virion at the initiation of infection. The phage genome is thereby sucked into the cell by hydrodynamic drag.

Protection of antiterminator RNA by the transcript elongation complex

Nascent transcripts encoded by the putL and putR sites of phage HK022 bind the transcript elongation complex and suppress termination at downstream transcription terminators. We report here that the chemical stability of putL RNA is considerably greater than that of the typical Escherichia coli message because the elongation complex protects this RNA from degradation. When binding to the elongation complex was prevented by mutation of either putL or RNA polymerase, RNA stability decreased more than 50-fold. The functional modification conferred by putL RNA on the elongation complex is also long-lived: the efficiency of terminator suppression remained high for at least 10 kb from the putL site. We find that RNase III rapidly and efficiently cleaved the transcript just downstream of the putL sequences, but such cleavage changed neither the stability of putL RNA nor the efficiency of antitermination. These results argue that the continuity of the RNA that connects put sequences to the growing point is not required for persistence of the antiterminating modification in vivo.

Multicopy plasmid modification with phage lambda Red recombineering

Recombineering, in vivo genetic engineering using the bacteriophage lambda Red generalized recombination system, was used to create various modifications of a multicopy plasmid derived from pBR322. All genetic modifications possible on the Escherichia coli chromosome and on bacterial artificial chromosomes (BACs) are also possible on multicopy plasmids and are obtained with similar frequencies to their chromosomal counterparts, including creation of point mutations (5-10% unselected frequency), deletions and substitutions. Parental and recombinant plasmids are nearly always present as a mixture following recombination, and circular multimeric plasmid molecules are often generated during the recombineering.

Biological consequences of tightly bent DNA: the other life of a macromolecular celebrity

The mechanical properties of DNA play a critical role in many biological functions. For example, DNA packing in viruses involves confining the viral genome in a volume (the viral capsid) with dimensions that are comparable to the DNA persistence length. Similarly, eukaryotic DNA is packed in DNA-protein complexes (nucleosomes), in which DNA is tightly bent around protein spools. DNA is also tightly bent by many proteins that regulate transcription, resulting in a variation in gene expression that is amenable to quantitative analysis. In these cases, DNA loops are formed with lengths that are comparable to or smaller than the DNA persistence length. The aim of this review is to describe the physical forces associated with tightly bent DNA in all of these settings and to explore the biological consequences of such bending, as increasingly accessible by single-molecule techniques.

Model-based design of growth-attenuated viruses

Live-virus vaccines activate both humoral and cell-mediated immunity, require only a single boosting, and generally provide longer immune protection than killed or subunit vaccines. However, growth of live-virus vaccines must be attenuated to minimize their potential pathogenic effects, and mechanisms of attenuation by conventional serial-transfer viral adaptation are not well-understood. New methods of attenuation based on rational engineering of viral genomes may offer a potentially greater control if one can link defined genetic modifications to changes in virus growth. To begin to establish such links between genotype and growth phenotype, we developed a computer model for the intracellular growth of vesicular stomatitis virus (VSV), a well-studied, nonsegmented, negative-stranded RNA virus. Our model incorporated established regulatory mechanisms of VSV while integrating key wild-type infection steps: hijacking of host resources, transcription, translation, and replication, followed by assembly and release of progeny VSV particles. Generalization of the wild-type model to allow for genome rearrangements matched the experimentally observed attenuation ranking for recombinant VSV strains that altered the genome position of their nucleocapsid gene. Finally, our simulations captured previously reported experimental results showing how altering the positions of other VSV genes has the potential to attenuate the VSV growth while overexpressing the immunogenic VSV surface glycoprotein. Such models will facilitate the engineering of new live-virus vaccines by linking genomic manipulations to controlled changes in virus gene-expression and growth.

The effect of genome length on ejection forces in bacteriophage lambda

A variety of viruses tightly pack their genetic material into protein capsids that are barely large enough to enclose the genome. In particular, in bacteriophages, forces as high as 60 pN are encountered during packaging and ejection, produced by DNA bending elasticity and self-interactions. The high forces are believed to be important for the ejection process, though the extent of their involvement is not yet clear. As a result, there is a need for quantitative models and experiments that reveal the nature of the forces relevant to DNA ejection. Here, we report measurements of the ejection forces for two different mutants of bacteriophage lambda, lambdab221cI26 and lambdacI60, which differ in genome length by approximately 30%. As expected for a force-driven ejection mechanism, the osmotic pressure at which DNA release is completely inhibited varies with the genome length: we find inhibition pressures of 15 atm and 25 atm, for the short and long genomes, respectively, values that are in agreement with our theoretical calculations.

The MS2 Coat Protein Shell is Likely Assembled Under Tension: A Novel Role for the MS2 Bacteriophage A Protein as Revealed by Small-angle Neutron Scattering

Recombinant forms of the bacteriophage MS2 and its RNA-free (empty) MS2 capsid were analyzed in solution to determine if RNA content and/or the A (or maturation) protein play a role in the global arrangement of the virus protein shell. Analysis of the (coat) protein shell of recombinant versions of MS2 that lack the A protein revealed dramatic differences compared to wild-type MS2 in solution. Specifically, A protein-deficient virus particles form a protein shell of between 31(+/-1) A and 37(+/-1) A. This is considerably thicker than the protein shell formed by either the wild-type MS2 or the RNA-free MS2 capsid, whose protein shells have a thickness of 21(+/-1) A and 25(+/-1) A, respectively. Since the A protein is known to separate from the intact MS2 protein shell after infection, the thin shell form of MS2 represents the pre-infection state, while the post-infection state is thick. Interestingly, these A protein-dependent differences in the virus protein shell are not seen using crystallography, as the crystallization process seems to artificially compact the wild-type MS2 virion. Furthermore, when the A protein is absent from the virus shell (post-infection), the process of crystallization exerts sufficient force to convert the protein shell from the post-infection (thick) state to the pre-infection (thin) conformation. In summary, the data are consistent with the idea that RNA content or amount does not affect the structure of the MS2 virus shell. Rather, the A protein influences the global arrangement of the virus coat dramatically, possibly by mediating the storage of energy or tension within the protein shell during virus assembly. This tension may later be used to eject the MS2 genomic RNA and A protein fragments into the host during infection.