MS2 bacteriophage capsid studied using all-atom molecular dynamics

Exploring the Effects of Intersubunit Interface Mutations on Virus-Like Particle Structure and Stability

Virus-like particles (VLPs) from bacteriophage MS2 provide a platform to study protein self-assembly and create engineered systems for drug delivery. Here, we aim to understand the impact of intersubunit interface mutations on the local and global structure and function of MS2-based VLPs. In previous work, our lab identified locally supercharged double mutants [T71K/G73R] that concentrate positive charge at capsid pores, enhancing uptake into mammalian cells. To study the effects of particle size on cellular internalization, we combined these double mutants with a single point mutation [S37P] that was previously reported to switch particle geometry from T = 3 to T = 1 icosahedral symmetry. These new variants retained their enhanced cellular uptake activity and could deliver small-molecule drugs with efficacy levels similar to our first-generation capsids. Surprisingly, these engineered triple mutants exhibit increased thermostability and unexpected geometry, producing T = 3 particles instead of the anticipated T = 1 assemblies. Transmission electron microscopy revealed various capsid assembly states, including wild-type (T = 3), T = 1, and rod-like particles, that could be accessed using different combinations of these point mutations. Molecular dynamics experiments recapitulated the structural rationale in silico for the single point mutation [S37P] forming a T = 1 virus-like particle and showed that this assembly state was not favored when combined with mutations that favor rod-like architectures. Through this work, we investigated how interdimer interface dynamics influence VLP size and morphology and how these properties affect particle function in applications such as drug delivery.

The coagulation process for enveloped and non-enveloped virus removal in turbid water: Removal efficiencies, mechanisms and its application to SARS-CoV-2 Omicron BA.2

The coagulation process has a high potential as a treatment method that can handle pathogenic viruses including emerging enveloped viruses in drinking water treatment process which can lower infection risk through drinking water consumption. In this study, a surrogate enveloped virus, bacteriophage Փ6, and surrogate non-enveloped viruses, including bacteriophage MS-2, T4, ՓX174, were used to evaluate removal efficiencies and mechanisms by the conventional coagulation process with alum, poly‑aluminum chloride, and ferric chloride at pH 5, 7, and 9 in turbid water. Also, treatability of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a recent virus of global concern by coagulation was evaluated as SARS-CoV-2 can presence in drinking water sources. It was observed that an increase in the coagulant dose enhanced the removal efficiency of turbidity and viruses, and the condition that provided the highest removal efficiency of enveloped and non-enveloped viruses was 50 mg/L of coagulants at pH 5. In addition, the coagulation process was more effective for enveloped virus removal than for the non-enveloped viruses, and it demonstrated reduction of SARS-CoV-2 Omicron BA.2 over 0.83-log with alum. According to culture- and molecular-based assays (qPCR and CDDP-qPCR), the virus removal mechanisms were floc adsorption and coagulant inactivation. Through inactivation with coagulants, coagulants caused capsid destruction, followed by genome damage in non-enveloped viruses; however, damage to a lipid envelope is suggested to contribute to a great extend for enveloped virus inactivation. We demonstrated that conventional coagulation is a promising method for controlling emerging and re-emerging viruses in drinking water.

Complete Virion Simulated: All-Atom Model of an MS2 Bacteriophage with Native Genome

For the first time, a complete all-atom molecular dynamics (MD) model of a virus, bacteriophage MS2, in its entirety, including a protein outer shell, native genomic RNA with necessary divalent ions, and surrounding explicit aqueous solution with ions at physiological concentration, was built. The model is based on an experimentally measured cryo-EM structure, which was substantially augmented by reconstructing missing or low-resolution parts of the measured density (where the atomistic structure cannot be fit unambiguously). The model was tested by a quarter of a microsecond MD run, and various biophysical characteristics are obtained and analyzed. The developed methodology of building the model can be used for reconstructing other large biomolecular structures when experimental data are fragmented and/or of varying resolution, while the model itself can be used for studying the biology of MS2, including the dynamics of its interaction with the host bacteria.

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

The balance between fitness advantages and costs drives adaptation of bacteriophage Qβ to changes in host density at different temperatures

Introduction

Host density is one of the main factors affecting the infective capacity of viruses. When host density is low, it is more difficult for the virus to find a susceptible cell, which increases its probability of being damaged by the physicochemical agents of the environment. Nevertheless, viruses can adapt to variations in host density through different strategies that depend on the particular characteristics of the life cycle of each virus. In a previous work, using the bacteriophage Qβ as an experimental model, we found that when bacterial density was lower than optimal the virus increased its capacity to penetrate into the bacteria through a mutation in the minor capsid protein (A1) that is not described to interact with the cell receptor.

Results

Here we show that the adaptive pathway followed by Qβ in the face of similar variations in host density depends on environmental temperature. When the value for this parameter is lower than optimal (30°C), the mutation selected is the same as at the optimal temperature (37°C). However, when temperature increases to 43°C, the mutation selected is located in a different protein (A2), which is involved both in the interaction with the cell receptor and in the process of viral progeny release. The new mutation increases the entry of the phage into the bacteria at the three temperatures assayed. However, it also considerably increases the latent period at 30 and 37°C, which is probably the reason why it is not selected at these temperatures.

Conclusion

The conclusion is that the adaptive strategies followed by bacteriophage Qβ, and probably other viruses, in the face of variations in host density depend not only on their advantages at this selective pressure, but also on the fitness costs that particular mutations may present in function of the rest of environmental parameters that influence viral replication and stability.

Heteroaggregation of virions and microplastics reduces the number of active bacteriophages in aqueous environments

The objective of this study is to explore the effects of microplastics on the viability of the bacteriophages in an aqueous environment. Bacteriophages (phages), that is, viruses of bacteria, are essential in homeostasis. It is estimated that phages cause up to 40% of the death of all bacteria daily. Any factor affecting phage activity is vital for the whole food chain and the ecology of numerous niches. We hypothesize that the number of active phages decreases due to the virions' adsorption on microplastic particles or by the released leachables from additives used in the production of plastic, for example, stabilizers, plasticizers, colorants, and reinforcements. We exposed three diverse phages, namely, T4 (tailed), MS2 (icosahedral), and M13 (filamentous), to 1 mg/mL suspension of 12 industrial-grade plastics [acrylonitrile butadiene styrene, high-impact polystyrene, poly-ε-caproamide, polycarbonate, polyethylene, polyethylene terephthalate, poly(methyl methacrylate), polypropylene, polystyrene, polytetrafluoroethylene, polyurethane, and polyvinyl chloride] shredded to obtain microparticles of radius ranging from 2 to 50 μm. The effect of leachables was measured upon exposure of phages not to particles themselves but to the buffer preincubated with microplastics. A double-overlay plaque counting method was used to assess phage titers. We employed a classical linear regression model to verify which physicochemical parameters (65 variables were tested) govern the decrease of phage titers. The key finding is that adsorption mechanisms result in up to complete scavenging of virions, whereas leachables deactivate up to 50% of phages. This study reveals microplastic pollution's plausible and unforeseen ecotoxicological effect causing phage deactivation. Moreover, phage transmission through adsorption can alter the balance of the food chain in the new environment. The effect depends mainly on the zeta potentials of the polymers and the phage type.

Breathing and Tilting: Mesoscale Simulations Illuminate Influenza Glycoprotein Vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from convalescent human donor, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

Reconstruction and validation of entire virus model with complete genome from mixed resolution cryo-EM density

It is very difficult to reconstruct computationally a large biomolecular complex in its biological entirety from experimental data. The resulting atomistic model should not contain gaps structurally and it should yield stable dynamics. We, for the first time, reconstruct from the published incomplete cryo-EM density a complete MS2 virus at atomistic resolution, that is, the capsid with the genome, and validate the result by all-atom molecular dynamics with explicit water. The available experimental data includes a high resolution protein capsid and an inhomogeneously resolved genome map. For the genomic RNA, apart from 16 hairpins with atomistic resolution, the strands near the capsid's inner surface were resolved up to the nucleic backbone level, and the innermost density was completely unresolved. As a result, only 242 nucleotides (out of 3569) were positioned, while only a fragmented backbone was outlined for the rest of the genome, making a detailed model reconstruction necessary. For model reconstruction, in addition to the available atomistic structure information, we extensively used the predicted secondary structure of the genome (base pairing). The technique was based on semi-automatic building of relatively large strands of RNA with subsequent manual positioning over the traced backbone. The entire virus structure (capsid + genome) was validated by a molecular dynamics run in physiological solution with ions at standard conditions confirming the stability of the model.

Stability and conformational memory of electrosprayed and rehydrated bacteriophage MS2 virus coat proteins

Proteins are innately dynamic, which is important for their functions, but which also poses significant challenges when studying their structures. Gas-phase techniques can utilise separation and a range of sample manipulations to transcend some of the limitations of conventional techniques for structural biology in crystalline or solution phase, and isolate different states for separate interrogation. However, the transfer from solution to the gas phase risks affecting the structures, and it is unclear to what extent different conformations remain distinct in the gas phase, and if resolution in silico can recover the native conformations and their differences. Here, we use extensive molecular dynamics simulations to study the two distinct conformations of dimeric capsid protein of the MS2 bacteriophage. The protein undergoes notable restructuring of its peripheral parts in the gas phase, but subsequent simulation in solvent largely recovers the native structure. Our results suggest that despite some structural loss due to the experimental conditions, gas-phase structural biology techniques provide meaningful data that inform not only about the structures but also conformational dynamics of proteins.

Estimation of Nanoparticle's Surface Electrostatic Potential in Solution Using Acid-Base Molecular Probes. III. Experimental Hydrophobicity/Hydrophilicity and Charge Distribution of MS2 Virus Surface

MS2 bacteriophage is often used as a model for evaluating pathogenic viruses' behavior in aqueous solution. However, the questions of the virus surface's hydrophilic/hydrophobic balance, the charge distribution, and the binding mechanism are open. Using the dynamic light scattering method and laser Doppler electrophoresis, the hydrodynamic diameter and the ζ-potential of the virus particles were measured at their concentration of 5 × 10¹¹ particles per mL and ionic strength 0.03 M. The values were found to be 30 nm and -29 or -34 mV (by Smoluchowski or Ohshima approximations), respectively. The MS2 bacteriophage surface was also investigated using a series of acid-base indicator dyes of various charge type, size, and structure. Their spectral and acid-base properties (pK_a) are very sensitive to the microenvironment in aqueous solution, including containing nanoparticles. The electrostatic potential of the surface Ψ was estimated using the common formula: Ψ = 59 × (pK_aⁱ - pK_a) in mV at 25 °C. The Ψ values were -50 and +10 mV, respectively, which indicate the "mosaic" way of the charge distribution on the surface. These data are in good agreement with the obtained ζ-potential values and provide even more information about the virus surface. It was found that the surface of the MS2 virus is hydrophilic in solution in contrast to the commonly accepted hypothesis of the hydrophobicity of virus particles. No hydrophobic interactions between various molecular probes and the capsid were observed.

Engineered Riboswitch Nanocarriers as a Possible Disease-Modifying Treatment for Metabolic Disorders

Both DNA- and RNA-based nanotechnologies are remarkably useful for the engineering of molecular devices in vitro and are applied in a vast collection of applications. Yet, the ability to integrate functional nucleic acid nanostructures in applications outside of the lab requires overcoming their inherent degradation sensitivity and subsequent loss of function. Viruses are minimalistic yet sophisticated supramolecular assemblies, capable of shielding their nucleic acid content in nuclease-rich environments. Inspired by this natural ability, we engineered RNA-virus-like particles (VLPs) nanocarriers (NCs). We showed that the VLPs can function as an exceptional protective shell against nuclease-mediated degradation. We then harnessed biological recognition elements and demonstrated how engineered riboswitch NCs can act as a possible disease-modifying treatment for genetic metabolic disorders. The functional riboswitch is capable of selectively and specifically binding metabolites and preventing their self-assembly process and its downstream effects. When applying the riboswitch nanocarriers to an in vivo yeast model of adenine accumulation and self-assembly, significant inhibition of the sensitivity to adenine feeding was observed. In addition, using an amyloid-specific dye, we proved the riboswitch nanocarriers' ability to reduce the level of intracellular amyloid-like metabolite cytotoxic structures. The potential of this RNA therapeutic technology does not apply only to metabolic disorders, as it can be easily fine-tuned to be applied to other conditions and diseases.

Breathing and tilting: mesoscale simulations illuminate influenza glycoprotein vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from human convalescent plasma, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

How Pore Architecture Regulates the Function of Nanoscale Protein Compartments

Self-assembling proteins can form porous compartments that adopt well-defined architectures at the nanoscale. In nature, protein compartments act as semipermeable barriers to enable spatial separation and organization of complex biochemical processes. The compartment pores play a key role in their overall function by selectively controlling the influx and efflux of important biomolecular species. By engineering the pores, the functionality of compartments can be tuned to facilitate non-native applications, such as artificial nanoreactors for catalysis. In this review, we analyze how protein structure determines the porosity and impacts the function of both native and engineered compartments, highlighting the wealth of structural data recently obtained by cryo-EM and X-ray crystallography. Through this analysis, we offer perspectives on how current structural insights can inform future studies into the design of artificial protein compartments as nanoreactors with tunable porosity and function.

All-atom molecular dynamics study of hepatitis B virus containing pregenome RNA in solution

Immature hepatitis B virus (HBV) captures nucleotides in its capsid for reverse transcription. The nucleotides and nucleotide analog drugs, which are triphosphorylated and negatively charged in the cell, approach the capsid via diffusion and are absorbed into it. In this study, we performed a long-time molecular dynamics calculation of the entire HBV capsid containing pregenome RNA to investigate the interactions between the capsid and negatively charged substances. Electric field analysis demonstrated that negatively charged substances can approach the HBV capsid by thermal motion, avoiding spikes. The substances then migrate all over the floor of the HBV capsid. Finally, they find pores through which they can pass through the HBV capsid shell. Free energy profiles were calculated along these pores for small ions to understand their permeability through the pores. Anions (Cl^-) showed higher free energy barriers than cations (Na⁺ and K⁺) through all pores, and the permeation rate of Cl^- was eight times slower than that of K⁺ or Na⁺. Furthermore, the ions were more stable in the capsid than in the bulk water. Thus, the HBV capsid exerts ion selectivity for uptake and provides an environment for ions, such as nucleotides and nucleotide analog drugs, to be stabilized within the capsid.

Gold-Polyoxoborates Nanocomposite Prohibits Adsorption of Bacteriophages on Inner Surfaces of Polypropylene Labware and Protects Samples from Bacterial and Yeast Infections

Bacteriophages (phages) are a specific type of viruses that infect bacteria. Because of growing antibiotic resistance among bacterial strains, phage-based therapies are becoming more and more attractive. The critical problem is the storage of bacteriophages. Recently, it was found that bacteriophages might adsorb on the surfaces of plastic containers, effectively decreasing the titer of phage suspensions. Here, we showed that a BOA nanocomposite (gold nanoparticles embedded in polyoxoborate matrix) deposited onto the inner walls of the containers stabilizes phage suspensions against uncontrolled adsorption and titer decrease. Additionally, BOA provides antibacterial and antifungal protection. The application of BOA assures safe and sterile means for the storage of bacteriophages.

Fighting viruses with computers, right now

The synergistic conjunction of various technological revolutions with the accumulated knowledge and workflows is rapidly transforming several scientific fields. Particularly, Virology can now feed from accurate physical models, polished computational tools, and massive computational power to readily integrate high-resolution structures into biological representations of unprecedented detail. That preparedness allows for the first time to get crucial information for vaccine and drug design from in-silico experiments against emerging pathogens of worldwide concern at relevant action windows. The present work reviews some of the main milestones leading to these breakthroughs in Computational Virology, providing an outlook for future developments in capacity building and accessibility to computational resources.

Interaction of the maturation protein of the bacteriophage MS2 and the sex pilus of the Escherichia coli F plasmid

One promising strategy to combat antimicrobial resistance is to use bacteriophages that attach to the sex pili produced by transmissible antimicrobial resistance (AMR) plasmids, infect AMR bacteria and select for loss of the AMR plasmids, prolonging the life of existing antimicrobials. The maturation protein of the bacteriophage MS2 attaches to the pili produced by Incompatibility group F plasmid-containing bacteria. This interaction initiates delivery of the viral genetic material into the bacteria. Using protein-protein docking we constructed a model of the F pilus comprising a trimer of subunits binding to the maturation protein. Interactions between the maturation protein and the F pilus were investigated using molecular dynamics simulations. In silico alanine scanning and in silico single-point mutations were explored, with the longer term aim of increasing the affinity of the maturation protein to other Incompatibility group pili, without reducing the strength of binding to F pilin. We report our computational findings on which residues are required for the maturation protein and F pilin to interact, those which had no effect on the interaction and the mutations which led to a stronger interaction.

Multiscale modelling and simulation of viruses

In recent years, advances in structural biology, integrative modelling, and simulation approaches have allowed us to gain unprecedented insights into viral structure and dynamics. In this article we survey recent studies utilizing this wealth of structural information to build computational models of partial or complete viruses and to elucidate mechanisms of viral function. Additionally, the close interplay of viral pathogens with host factors - such as cellular and intracellular membranes, receptors, antibodies, and other host proteins - makes accurate models of viral interactions and dynamics essential. As viruses continue to pose severe challenges in prevention and treatment, enhancing our mechanistic understanding of viral infection is vital to enable the development of novel therapeutic strategies.

Coarse-graining molecular dynamics: stochastic models with non-Gaussian force distributions

Incorporating atomistic and molecular information into models of cellular behaviour is challenging because of a vast separation of spatial and temporal scales between processes happening at the atomic and cellular levels. Multiscale or multi-resolution methodologies address this difficulty by using molecular dynamics (MD) and coarse-grained models in different parts of the cell. Their applicability depends on the accuracy and properties of the coarse-grained model which approximates the detailed MD description. A family of stochastic coarse-grained (SCG) models, written as relatively low-dimensional systems of nonlinear stochastic differential equations, is presented. The nonlinear SCG model incorporates the non-Gaussian force distribution which is observed in MD simulations and which cannot be described by linear models. It is shown that the nonlinearities can be chosen in such a way that they do not complicate parametrization of the SCG description by detailed MD simulations. The solution of the SCG model is found in terms of gamma functions.

Multi-resolution simulations of intracellular processes

This is an introduction to the Thematic Issue of Interface Focus containing papers by speakers of the Theo Murphy International Scientific Meeting on ‘Multi-resolution simulations of intracellular processes’