Energies and pressures in viruses: contribution of nonspecific electrostatic interactions

Implications of AAV affinity column reuse and vector stability on product quality attributes

In this work, the implications of AAV9 capsid design and column reuse on AAV9 vector product quality were assessed with POROS CaptureSelect (PCS) AAVX and AAV9 resins using sf9 insect cell-derived model AAV9 vectors with varying viral protein (VP) ratios. Chromatographic experiments with purified drug substance AAV9 model feeds indicated consistent vector elution profiles, independent of adeno-associated virus (AAV) VP ratio, or cycle number. In contrast, the presence of process impurities in the clarified lysate feeds resulted in clear changes in the elution patterns. This included increased aggregate content in the vector eluates over multiple cycles as well as clear differences in the performance of these affinity resin systems. The AAV9-serotype specific PCS AAV9 column, with lower vector elution pH, resulted in higher aggregate content over multiple cycles as compared to the serotype-independent PCS AAVX column. Further, the results with vectors of varying VP ratio indicated that while one vector type eluate displayed higher aggregation in both affinity columns over column reuse, the eluate with the other vector type did not exhibit changes in the aggregation profile. Interestingly, vector aggregates in the affinity eluates also contained double-stranded DNA impurities and histone proteins, with similar trends to the aggregate levels. This behavior upon column reuse indicates that these host cell impurities are likely carried over to subsequent runs due to incomplete clean-in-place (CIP). These results indicate that feed impurities, affinity resin characteristics, elution pH, column CIP, and vector stability can impact the reusability of AAV affinity columns and product quality.

The conformational phase diagram of charged polymers in the presence of attractive bridging crowders

Using extensive molecular dynamics simulations, we obtain the conformational phase diagram of a charged polymer in the presence of oppositely charged counterions and neutral attractive crowders for monovalent, divalent, and trivalent counterion valencies. We demonstrate that the charged polymer can exist in three phases: (1) an extended phase for low charge densities and weak polymer-crowder attractive interactions [Charged Extended (CE)]; (2) a collapsed phase for high charge densities and weak polymer-crowder attractive interactions, primarily driven by counterion condensation [Charged Collapsed due to Intra-polymer interactions [(CCI)]; and (3) a collapsed phase for strong polymer-crowder attractive interactions, irrespective of the charge density, driven by crowders acting as bridges or cross-links [Charged Collapsed due to Bridging interactions [(CCB)]. Importantly, simulations reveal that the interaction with crowders can induce collapse, despite the presence of strong repulsive electrostatic interactions, and can replace condensed counterions to facilitate a direct transition from the CCI and CE phases to the CCB phase.

Evolutionary changes in the number of dissociable amino acids on spike proteins and nucleoproteins of SARS-CoV-2 variants

The spike protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for target recognition, cellular entry, and endosomal escape of the virus. At the same time, it is the part of the virus that exhibits the greatest sequence variation across the many variants which have emerged during its evolution. Recent studies have indicated that with progressive lineage emergence, the positive charge on the spike protein has been increasing, with certain positively charged amino acids (AAs) improving the binding of the spike protein to cell receptors. We have performed a detailed analysis of dissociable AAs of more than 1400 different SARS-CoV-2 lineages, which confirms these observations while suggesting that this progression has reached a plateau with Omicron and its subvariants and that the positive charge is not increasing further. Analysis of the nucleocapsid protein shows no similar increase in positive charge with novel variants, which further indicates that positive charge of the spike protein is being evolutionarily selected for. Furthermore, comparison with the spike proteins of known coronaviruses shows that already the wild-type SARS-CoV-2 spike protein carries an unusually large amount of positively charged AAs when compared to most other betacoronaviruses. Our study sheds light on the evolutionary changes in the number of dissociable AAs on the spike protein of SARS-CoV-2, complementing existing studies and providing a stepping stone towards a better understanding of the relationship between the spike protein charge and viral infectivity and transmissibility.

Hidden symmetry of the flavivirus protein shell and pH-controlled reconstruction of the viral surface

Using recent Zika virus structural data we reveal a hidden symmetry of protein order in immature and mature flavivirus shells, violating the Caspar-Klug paradigmatic model of capsid structures. We show that proteins of the outer immature shell layer exhibit trihexagonal tiling, while proteins from inner and outer layers conjointly form a double-shelled close-packed structure, based on a common triangular spherical lattice. Within the proposed structural model, we furthermore rationalize the structural organization of misassembled non-infectious subviral particles that have no inner capsid. We consider a pH-controlled structural reconstruction of the outer shell from the trimeric to the dimeric state, and demonstrate that this transition, occurring during the virus maturation, can be induced by changes in protein charges at lower pH, leading to a decrease in the electrostatic interaction free energy. This transition could also be assisted by electrostatic attraction of shell proteins to the interposed lipid membrane substrate separating the shells.

Spanning tree model and the assembly kinetics of RNA viruses

Single-stranded RNA (ssRNA) viruses self-assemble spontaneously in solutions that contain the viral RNA genome molecules and viral capsid proteins. The self-assembly of empty capsids can be understood on the basis of free energy minimization. However, during the self-assembly of complete viral particles in the cytoplasm of an infected cell, the viral genome molecules must be selected from a large pool of very similar host messenger RNA molecules and it is not known whether this also can be understood by free energy minimization. We address this question using a simple mathematical model, the spanning tree model, that was recently proposed for the assembly of small ssRNA viruses. We present a statistical physics analysis of the properties of this model. RNA selection takes place via a kinetic mechanism that operates during the formation of the nucleation complex and that is related to Hopfield kinetic proofreading.

Biophysical Modeling of SARS-CoV-2 Assembly: Genome Condensation and Budding

The COVID-19 pandemic caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spurred unprecedented and concerted worldwide research to curtail and eradicate this pathogen. SARS-CoV-2 has four structural proteins: Envelope (E), Membrane (M), Nucleocapsid (N), and Spike (S), which self-assemble along with its RNA into the infectious virus by budding from intracellular lipid membranes. In this paper, we develop a model to explore the mechanisms of RNA condensation by structural proteins, protein oligomerization and cellular membrane-protein interactions that control the budding process and the ultimate virus structure. Using molecular dynamics simulations, we have deciphered how the positively charged N proteins interact and condense the very long genomic RNA resulting in its packaging by a lipid envelope decorated with structural proteins inside a host cell. Furthermore, considering the length of RNA and the size of the virus, we find that the intrinsic curvature of M proteins is essential for virus budding. While most current research has focused on the S protein, which is responsible for viral entry, and it has been motivated by the need to develop efficacious vaccines, the development of resistance through mutations in this crucial protein makes it essential to elucidate the details of the viral life cycle to identify other drug targets for future therapy. Our simulations will provide insight into the viral life cycle through the assembly of viral particles de novo and potentially identify therapeutic targets for future drug development.

Packaging contests between viral RNA molecules and kinetic selectivity

The paper presents a statistical-mechanics model for the kinetic selection of viral RNA molecules by packaging signals during the nucleation stage of the assembly of small RNA viruses. The effects of the RNA secondary structure and folding geometry of the packaging signals on the assembly activation energy barrier are encoded by a pair of characteristics: the wrapping number and the maximum ladder distance. Kinetic selection is found to be optimal when assembly takes place under conditions of supersaturation and also when the concentration ratio of capsid protein and viral RNA concentrations equals the stoichiometric ratio of assembled viral particles. As a function of the height of the activation energy barrier, there is a form of order-disorder transition such that for sufficiently low activation energy barriers, kinetic selectivity is erased by entropic effects associated with the number of assembly pathways.

The effect of explicit polarity on the conformational behavior of a single polyelectrolyte chain

In this work using dissipative particle dynamics simulations with explicit treatment of polar species we demonstrate that the molecular nature of dielectric media has a significant impact on swelling and collapse of a polyelectrolyte chain in a dilute solution. We show that the small-scale effects related to the presence of polar species lead to the intensification of the electrostatic interactions when the charges are close to each other and/or their density is high enough. As a result, the electrostatic strength , usually regarded as the main parameter governing the polyelectrolyte chain collapse, does not have a universal meaning: the value of λ at which the coil-to-globule transition occurs is found to be dependent on the specific fixed value of the solvent bulk permittivity ε while varying the monomer unit charge Q and vice versa. This effect is observed even when the backbone and the counterions have the same polarity as the solvent beads, i.e. no dielectric mismatch is present. The reason for such behavior is rationalized in terms of the "effective" dielectric permittivity ε^eff which depends on the volume fraction φ of charged units inside the polymer chain volume; using ε^eff instead of ε collapses all data onto one master curve describing the chain shrinking with λ. Furthermore, it is shown that a polar chain adopts less swollen conformations in the polyelectrolyte regime and collapses more easily compared to a non-polar chain.

Mechanisms for destabilisation of RNA viruses at air-water and liquid-liquid interfaces

Understanding the interactions between viruses and surfaces or interfaces is important, as they provide the principles underpinning the cleaning and disinfection of contaminated surfaces. Yet, the physics of such interactions is currently poorly understood. For instance, there are longstanding experimental observations suggesting that the presence of air-water interfaces can generically inactivate and kill viruses, yet the mechanism underlying this phenomenon remains unknown. Here we use theory and simulations to show that electrostatics may provide one such mechanism, and that this is very general. Thus, we predict that the electrostatic free energy of an RNA virus should increase by several thousands of k_BT as the virion breaches an air-water interface. We also show that the fate of a virus approaching a generic liquid-liquid interface depends strongly on the detailed balance between interfacial and electrostatic forces, which can be tuned, for instance, by choosing different media to contact a virus-laden respiratory droplet. Tunability arises because both the electrostatic and interfacial forces scale similarly with viral size. We propose that these results can be used to design effective strategies for surface disinfection.

We know that air-water interfaces can generically inactivate viruses, but the mechanisms behind this observation are unclear. Here the authors use simulations to uncover those mechanisms and find that the electrostatic repulsive free energy of an RNA virus increases by several thousands of kBT as it approaches an air-water interface, providing a mechanism for viral destabilization which may induce inactivation.

Kinetics of charged polymer collapse in poor solvents

Extensive molecular dynamics simulations, using simple charged polymer models, have been employed to probe the collapse kinetics of a single flexible polyelectrolyte (PE) chain under implicit poor solvent conditions. We investigate the role of the charged nature of PE chain (A), valency of counterions (Z) on the kinetics of such PE collapse. Our study shows that the collapse kinetics of charged polymers are significantly different from those of the neutral polymer and that the finite-size scaling behavior of PE collapse times does not follow the Rouse scaling as observed in the case of neutral polymers. The critical exponent for charged PE chains is found to be less than that of neutral polymers and also exhibits dependence on counterion valency. The coarsening of clusters along the PE chain suggests a multi-stage collapse and exhibits opposite behavior of exponents compared to neutral polymers: faster in the early stages and slower in the later stages of collapse.

Effects of surface charge and environmental factors on the electrostatic interaction of fiber with virus-like particle: A case of coronavirus

We propose a theoretical model to elucidate intermolecular electrostatic interactions between a virus and a substrate. Our model treats the virus as a homogeneous particle having surface charge and the polymer fiber of the respirator as a charged plane. Electric potentials surrounding the virus and fiber are influenced by the surface charge distribution of the virus. We use Poisson-Boltzmann equations to calculate electric potentials. Then, Derjaguin's approximation and a linear superposition of the potential function are extended to determine the electrostatic force. In this work, we apply this model for coronavirus or SARS-CoV-2 case and numerical results quantitatively agree with prior simulation. We find that the influence of fiber's potential on the surface charge of the virus is important and is considered in interaction calculations to obtain better accuracy. The electrostatic interaction significantly decays with increasing separation distance, and this curve becomes steeper when adding more salt. Although the interaction force increases with heating, one can observe the repulsive-attractive transition when the environment is acidic.

Surfactants - Compounds for inactivation of SARS-CoV-2 and other enveloped viruses

We provide here a general view on the interactions of surfactants with viruses, with a particular emphasis on how such interactions can be controlled and employed for inhibiting the infectivity of enveloped viruses, including coronaviruses. The aim is to provide to interested scientists from different fields, including chemistry, physics, biochemistry, and medicine, an overview of the basic properties of surfactants and (corona)viruses, which are relevant to understanding the interactions between the two. Various types of interactions between surfactant and virus are important, and they act on different components of a virus such as the lipid envelope, membrane (envelope) proteins and nucleocapsid proteins. Accordingly, this cannot be a detailed account of all relevant aspects but instead a summary that bridges between the different disciplines. We describe concepts and cover a selection of the relevant literature as an incentive for diving deeper into the relevant material. Our focus is on more recent developments around the COVID-19 pandemic caused by SARS-CoV-2, applications of surfactants against the virus, and on the potential future use of surfactants for pandemic relief. We also cover the most important aspects of the historical development of using surfactants in combatting virus infections. We conclude that surfactants are already playing very important roles in various directions of defence against viruses, either directly, as in disinfection, or as carrier components of drug delivery systems for prophylaxis or treatment. By designing tailor-made surfactants, and consequently, advanced formulations, one can expect more and more effective use of surfactants, either directly as antiviral compounds or as part of more complex formulations.

Relationships between RNA topology and nucleocapsid structure in a model icosahedral virus

The process of genome packaging in most of viruses is poorly understood, notably the role of the genome itself in the nucleocapsid structure. For simple icosahedral single-stranded RNA viruses, the branched topology due to the RNA secondary structure is thought to lower the free energy required to complete a virion. We investigate the structure of nucleocapsids packaging RNA segments with various degrees of compactness by small-angle x-ray scattering and cryotransmission electron microscopy. The structural differences are mild even though compact RNA segments lead on average to better-ordered and more uniform particles across the sample. Numerical calculations confirm that the free energy is lowered for the RNA segments displaying the larger number of branch points. The effect is, however, opposite with synthetic polyelectrolytes, in which a star topology gives rise to more disorder in the capsids than a linear topology. If RNA compactness and size account in part for the proper assembly of the nucleocapsid and the genome selectivity, other factors most likely related to the host cell environment during viral assembly must come into play as well.

Relative humidity in droplet and airborne transmission of disease

A large number of infectious diseases are transmitted by respiratory droplets. How long these droplets persist in the air, how far they can travel, and how long the pathogens they might carry survive are all decisive factors for the spread of droplet-borne diseases. The subject is extremely multifaceted and its aspects range across different disciplines, yet most of them have only seldom been considered in the physics community. In this review, we discuss the physical principles that govern the fate of respiratory droplets and any viruses trapped inside them, with a focus on the role of relative humidity. Importantly, low relative humidity-as encountered, for instance, indoors during winter and inside aircraft-facilitates evaporation and keeps even initially large droplets suspended in air as aerosol for extended periods of time. What is more, relative humidity affects the stability of viruses in aerosol through several physical mechanisms such as efflorescence and inactivation at the air-water interface, whose role in virus inactivation nonetheless remains poorly understood. Elucidating the role of relative humidity in the droplet spread of disease would permit us to design preventive measures that could aid in reducing the chance of transmission, particularly in indoor environment.

Virus Mechanics under Molecular Crowding

Viruses avoid exposure of the viral genome to harmful agents with the help of a protective protein shell known as the capsid. A secondary effect of this protective barrier is that macromolecules that may be in high concentration on the outside cannot freely diffuse across it. Therefore, inside the cell and possibly even outside, the intact virus is generally under a state of osmotic stress. Viruses deal with this type of stress in various ways. In some cases, they might harness it for infection. However, the magnitude and influence of osmotic stress on virus physical properties remains virtually unexplored for single-stranded RNA viruses-the most abundant class of viruses. Here, we report on how a model system for the positive-sense RNA icosahedral viruses, brome mosaic virus (BMV), responds to osmotic pressure. Specifically, we study the mechanical properties and structural stability of BMV under controlled molecular crowding conditions. We show that BMV is mechanically reinforced under a small external osmotic pressure but starts to yield after a threshold pressure is reached. We explain this mechanochemical behavior as an effect of the molecular crowding on the entropy of the "breathing" fluctuation modes of the virus shell. The experimental results are consistent with the viral RNA imposing a small negative internal osmotic pressure that prestresses the capsid. Our findings add a new line of inquiry to be considered when addressing the mechanisms of viral disassembly inside the crowded environment of the cell.

Theory of Inhomogeneous Rod-like Coulomb Fluids

A field theoretic representation of the classical partition function is derived for a system composed of a mixture of anisotropic and isotropic mobile charges that interact via long range Coulomb and short range nematic interactions. The field theory is then solved on a saddle-point approximation level, leading to a coupled system of Poisson–Boltzmann and Maier–Saupe equations. Explicit solutions are finally obtained for a rod-like counterion-only system in proximity to a charged planar wall. The nematic order parameter profile, the counterion density profile and the electrostatic potential profile are interpreted within the framework of a nematic–isotropic wetting phase with a Donnan potential difference.

Buckling transitions and soft-phase invasion of two-component icosahedral shells

What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.

Virus Isoelectric Point Estimation: Theories and Methods

Much of virus fate, both in the environment and in physical/chemical treatment, is dependent on electrostatic interactions. Developing an accurate means of predicting virion isoelectric point (pI) would help to understand and anticipate virus fate and transport, especially for viruses that are not readily propagated in the lab. One simple approach to predicting pI estimates the pH at which the sum of charges from ionizable amino acids in capsid proteins approaches zero. However, predicted pIs based on capsid charges frequently deviate by several pH units from empirically measured pIs. Recently, the discrepancy between empirical and predicted pI was attributed to the electrostatic neutralization of predictable polynucleotide-binding regions (PBRs) of the capsid interior. In this paper, we review models presupposing (i) the influence of the viral polynucleotide on surface charge or (ii) the contribution of only exterior residues to surface charge. We then compare these models to the approach of excluding only PBRs and hypothesize a conceptual electrostatic model that aligns with this approach. The PBR exclusion method outperformed methods based on three-dimensional (3D) structure and accounted for major discrepancies in predicted pIs without adversely affecting pI prediction for a diverse range of viruses. In addition, the PBR exclusion method was determined to be the best available method for predicting virus pI, since (i) PBRs are predicted independently of the impact on pI, (ii) PBR prediction relies on proteome sequences rather than detailed structural models, and (iii) PBR exclusion was successfully demonstrated on a diverse set of viruses. These models apply to nonenveloped viruses only. A similar model for enveloped viruses is complicated by a lack of data on enveloped virus pI, as well as uncertainties regarding the influence of the phospholipid envelope on charge and ion gradients.

Effect of the charge distribution of virus coat proteins on the length of packaged RNAs

Single-stranded RNA viruses efficiently encapsulate their genome into a protein shell called the capsid. Electrostatic interactions between the positive charges in the capsid protein's N-terminal tail and the negatively charged genome have been postulated as the main driving force for virus assembly. Recent experimental results indicate that the N-terminal tail with the same number of charges and same lengths packages different amounts of RNA, which reveals that electrostatics alone cannot explain all the observed outcomes of the RNA self-assembly experiments. Using a mean-field theory, we show that the combined effect of genome configurational entropy and electrostatics can explain to some extent the amount of packaged RNA with mutant proteins where the location and number of charges on the tails are altered. Understanding the factors contributing to the virus assembly could promote the attempt to block viral infections or to build capsids for gene therapy applications.

Genome Modularization Reveals Overlapped Gene Topology Is Necessary for Efficient Viral Reproduction

Sequence overlap between two genes is common across all genomes, with viruses having high proportions of these gene overlaps. Genome modularization and refactoring is the process of disrupting natural gene overlaps to separate coding sequences to enable their individual manipulation. The biological function and fitness effects of gene overlaps are not fully understood, and their effects on gene cluster and genome-level refactoring are unknown. The bacteriophage φX174 genome has ∼26% of nucleotides involved in encoding more than one gene. In this study we use an engineered φX174 phage containing a genome with all gene overlaps removed to show that gene overlap is critical to maintaining optimal viral fecundity. Through detailed phenotypic measurements we reveal that genome modularization in φX174 causes virion replication, stability, and attachment deficiencies. Quantitation of the complete phage proteome across an infection cycle reveals 30% of proteins display abnormal expression patterns. Taken together, we have for the first time comprehensively demonstrated that gene modularization severely perturbs the coordinated functioning of a bacteriophage replication cycle. This work highlights the biological importance of gene overlap in natural genomes and that reducing gene overlap disruption should be an integral part of future genome engineering projects.

Making sense of virus size and the tradeoffs shaping viral fitness

Viruses span an impressive size range, with genome length varying a thousandfold and virion volume nearly a millionfold. For cellular organisms the scaling of traits with size is a pervasive influence on ecological processes, but whether size plays a central role in viral ecology is unknown. Here, we focus on viruses of aquatic unicellular organisms, which exhibit the greatest known range of virus size. We outline hypotheses within a quantitative framework, and analyse data where available, to consider how size affects the primary components of viral fitness. We argue that larger viruses have fewer offspring per infection and slower contact rates with host cells, but a larger genome tends to increase infection efficiency, broaden host range, and potentially increase attachment success and decrease decay rate. These countervailing selective pressures may explain why a breadth of sizes exist and even coexist when infecting the same host populations. Oligotrophic ecosystems may be enriched in "giant" viruses, because environments with resource-limited phagotrophs at low concentrations may select for broader host range, better control of host metabolism, lower decay rate and a physical size that mimics bacterial prey. Finally, we describe where further research is needed to understand the ecology and evolution of viral size diversity.

Theoretical Modeling of Chemical Equilibrium in Weak Polyelectrolyte Layers on Curved Nanosystems

Surface functionalization with end-tethered weak polyelectrolytes (PE) is a versatile way to modify and control surface properties, given their ability to alter their degree of charge depending on external cues like pH and salt concentration. Weak PEs find usage in a wide range of applications, from colloidal stabilization, lubrication, adhesion, wetting to biomedical applications such as drug delivery and theranostics applications. They are also ubiquitous in many biological systems. Here, we present an overview of some of the main theoretical methods that we consider key in the field of weak PE at interfaces. Several applications involving engineered nanoparticles, synthetic and biological nanopores, as well as biological macromolecules are discussed to illustrate the salient features of systems involving weak PE near an interface or under (nano)confinement. The key feature is that by confining weak PEs near an interface the degree of charge is different from what would be expected in solution. This is the result of the strong coupling between structural organization of weak PE and its chemical state. The responsiveness of engineered and biological nanomaterials comprising weak PE combined with an adequate level of modeling can provide the keys to a rational design of smart nanosystems.

Soft matter science and the COVID-19 pandemic

Much of the science underpinning the global response to the COVID-19 pandemic lies in the soft matter domain. Coronaviruses are composite particles with a core of nucleic acids complexed to proteins surrounded by a protein-studded lipid bilayer shell. A dominant route for transmission is via air-borne aerosols and droplets. Viral interaction with polymeric body fluids, particularly mucus, and cell membranes controls their infectivity, while their interaction with skin and artificial surfaces underpins cleaning and disinfection and the efficacy of masks and other personal protective equipment. The global response to COVID-19 has highlighted gaps in the soft matter knowledge base. We survey these gaps, especially as pertaining to the transmission of the disease, and suggest questions that can (and need to) be tackled, both in response to COVID-19 and to better prepare for future viral pandemics.

Topological-Defect-Induced Surface Charge Heterogeneities in Nematic Electrolytes

We show that topological defects in an ion-doped nematic liquid crystal can be used to manipulate the surface charge distribution on chemically homogeneous, charge-regulating external surfaces, using a minimal theoretical model. In particular, the location and type of the defect encodes the precise distribution of surface charges and the effect is enhanced when the liquid crystal is flexoelectric. We demonstrate the principle for patterned surfaces and charged colloidal spheres. More generally, our results indicate an interesting approach to control surface charges on external surfaces without changing the surface chemistry.

Influence of slip velocity at the core of a diffuse soft particle and ion partition effects on mobility

Nonlinear effects on the electrophoresis of a soft particle, consisting of a rigid hydrophobic core coated with a diffuse polymer layer (PEL) suspended in an electrolyte medium, are studied. The impact of the ion partitioning effect arising due to the Born energy difference between the PEL and the electrolyte is approximated based on the equilibrium Boltzmann equation, with which the ion distribution and hence, the charge density is modified. The equations describing the electrokinetic transport comprising the Darcy-Brinkman extended Navier-Stokes equations which includes the ion partitioning effect coupled with the modified Nernst-Planck equations and Poisson equations for electric field are solved numerically. The present numerical model for the soft particle compares well with the existing theoretical solutions and experimental results in the limiting cases. A deviation from existing simplified models based on the Boltzmann distribution of ions occurs when the Debye layer polarization, relaxation and the electroosmosis induced by the PEL immobile charge become significant. The hydrophobicity of the inner core strongly influences the nonlinear electrokinetic effects by modifying the Debye layer, electroosmotic flow in the PEL and surface conduction. The results indicate that the ion partitioning can significantly increase the electrophoretic mobility of the soft particle by attenuating the shielding effect. When the Debye layer is in the order of the particle size the hydrophobicity of the core surface and the ion partitioning effect manifest the surface conduction, which implies that the Boltzmann distribution of ions is no longer valid. The core hydrophobicity and ion partitioning effect have influence on the condensation of the PEL immobile charge, which creates a significant impact on the mobility.

Icosadeltahedral Geometry of Geodesic Domes, Fullerenes and Viruses: A Tutorial on the T-Number

The Caspar–Klug (CK) classification of viruses is discussed by parallel examination of geometry of icosahedral geodesic domes, fullerenes, and viruses. The underlying symmetry of all structures is explained and thoroughly visually represented. Euler’s theorem on polyhedra is used to calculate the number of vertices, edges, and faces in domes, number of atoms, bonds, and pentagonal and hexagonal rings in fullerenes, and number of proteins and protein–protein contacts in viruses. The T-number, the characteristic for the CK classification, is defined and discussed. The superposition of fullerene and dome designs is used to obtain a representation of a CK virus with all the proteins indicated. Some modifications of the CK classifications are sketched, including elongation of the CK blueprint, fusion of two CK blueprints, dodecahedral view of the CK shapes, and generalized CK designs without a clearly visible geometry of the icosahedron. These are compared to cases of existing viruses.

Analytic Expression for Electrophoretic Mobility of Soft Particles with a Hydrophobic Inner Core at Different Electrostatic Conditions

This paper presents a simplified model for the electrophoresis of a soft particle with a nonwettable rigid core with charged polyelectrolyte corona under a weak-field and low-charge density consideration. We have derived a closed form solution for the mobility, which reduces to the well-known expressions for mobility as derived by Ohshima for limiting cases such as a hydrophilic charged core coated with an uncharged polymer (Ohshima, H. J. Colloid Interface Sci. 2002, 252, 119-125) or an uncharged no-slip core coated with a polyelectrolyte layer (Ohshima, H. Electrophoresis 2006, 27, 526-533). The generalized mobility expression reduces to the existing expression for mobility of a rigid hydrophobic colloid as the soft layer shrinks to zero. The general form of the mobility expression involves elliptic integrals, which can be computed easily through a software like Mathematica. We have derived analytical solutions for mobility pertaining to several particular cases. The occurrence of mobility reversal when the core and polyelectrolyte layer has a charge of opposite polarity is demonstrated in this paper.

Osmotic stress and pore nucleation in charged biological nanoshells and capsids

A model system is proposed to investigate the chemical equilibrium and mechanical stability of biological spherical-like nanoshells in contact with an aqueous solution with added dissociated electrolyte of a given concentration. The ionic chemical equilibrium across the permeable shell is investigated in the framework of an accurate Density Functional Theory (DFT) that incorporates electrostatic and hardcore correlations beyond the traditional mean-field (e.g., Poisson-Boltzmann) limit. The accuracy of the theory is tested by a direct comparison with Monte Carlo (MC) simulations. A simple analytical expression is then deduced which clearly highlights the entropic, electrostatic, and self-energy contributions to the osmotic stress over the shell in terms of the calculated ionic profiles. By invoking a continuum mean-field elastic approach to account for the shell surface stress upon osmotic stretching, the mechanical equilibrium properties of the shell under a wide variety of ionic strengths and surface charges are investigated. The model is further coupled to a continuum mechanical approach similar in structure to a Classical Nucleation Theory (CNT) to address the question of mechanical stability of the shells against a pore nucleation. This allows us to construct a phase diagram which delimits the mechanical stability of capsids for different ionic strengths and shell surface charges.

On the Response of Nanoelectrode Impedance Spectroscopy Measures to Plant, Animal, and Human Viruses

A simplified lumped geometrical and electrical model for the high-frequency impedance spectroscopy (HFIS) response of nanoelectrodes to capsids and full viruses is developed starting from atomistic descriptions, in order to test the theoretical response of a realistic HFIS CMOS biosensor platform to different viruses. Capacitance spectra are computed for plant (cowpea chlorotic mottle virus), animal (rabbit haemorrhagic disease virus), and human (hepatitis A virus) viruses. A few common features of the spectra are highlighted, and the role of virus charge, pH, and ionic strength on the expected signal is discussed. They suggest that the frequency of highest sensitivity at nearly physiological concentrations (100 mM) is within reach of existing HFIS platform designs.

Spherical structure factor and classification of hyperuniformity on the sphere

Understanding how particles are arranged on the surface of a sphere is not only central to numerous physical, biological, soft matter, and materials systems but also finds applications in computational problems, approximation theory, and analysis of geophysical and meteorological measurements. Objects that lie on a sphere experience constraints that are not present in Euclidean (flat) space and that influence both how the particles can be arranged as well as their statistical properties. These constraints, coupled with the curved geometry, require a careful extension of quantities used for the analysis of particle distributions in Euclidean space to distributions confined to the surface of a sphere. Here, we introduce a framework designed to analyze and classify structural order and disorder in particle distributions constrained to the sphere. The classification is based on the concept of hyperuniformity, which was first introduced 15 years ago and since then studied extensively in Euclidean space, yet has only very recently been considered also for spherical surfaces. We employ a generalization of the structure factor on the sphere, related to the power spectrum of the corresponding multipole expansion of particle density distribution. The spherical structure factor is then shown to couple with cap number variance, a measure of density variations at different scales, allowing us to analytically derive different forms of the variance pertaining to different types of distributions. Based on these forms, we construct a classification of hyperuniformity for scale-free particle distributions on the sphere and show how it can be extended to include other distribution types as well. We demonstrate that hyperuniformity on the sphere can be defined either through a vanishing spherical structure factor at low multipole numbers or through a scaling of the cap number variance-in both cases extending the Euclidean definition, while at the same time pointing out crucial differences. Our work thus provides a comprehensive tool for detecting global, long-range order on spheres and for the analysis of spherical computational meshes, biological and synthetic spherical assemblies, and ordering phase transitions in spherically distributed particles.

Charged nanorods at heterogeneously charged surfaces

We study the spatial and orientational distribution of charged nanorods (rodlike counterions) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The nanorods are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude.

Electrostatics-Driven Inflation of Elastic Icosahedral Shells as a Model for Swelling of Viruses

We develop a clear theoretical description of radial swelling in virus-like particles that delineates the importance of electrostatic contributions to swelling in the absence of any conformational changes. The model couples the elastic parameters of the capsid-represented as a continuous elastic shell-to the electrostatic pressure acting on it. We show that different modifications of the electrostatic interactions brought about by, for instance, changes in pH or solution ionic strength are often sufficient to achieve the experimentally observed swelling (∼10% of the capsid radius). Additionally, we derive analytical expressions for the electrostatics-driven radial swelling of virus-like particles that enable one to quickly estimate the magnitudes of physical quantities involved.

Systematic analysis of biological roles of charged amino acid residues located throughout the structured inner wall of a virus capsid

Structure-based mutational analysis of viruses is providing many insights into the relationship between structure and biological function of macromolecular complexes. We have systematically investigated the individual biological roles of charged residues located throughout the structured capsid inner wall (outside disordered peptide segments) of a model spherical virus, the minute virus of mice (MVM). The functional effects of point mutations that altered the electrical charge at 16 different positions at the capsid inner wall were analyzed. The results revealed that MVM capsid self-assembly is rather tolerant to point mutations that alter the number and distribution of charged residues at the capsid inner wall. However, mutations that either increased or decreased the number of positive charges around capsid-bound DNA segments reduced the thermal resistance of the virion. Moreover, mutations that either removed or changed the positions of negatively charged carboxylates in rings of acidic residues around capsid pores were deleterious by precluding a capsid conformational transition associated to through-pore translocation events. The results suggest that number, distribution and specific position of electrically charged residues across the inner wall of a spherical virus may have been selected through evolution as a compromise between several different biological requirements.

Adsorption and encapsulation of flexible polyelectrolytes in charged spherical vesicles

We present a theory of adsorption of flexible polyelectrolytes on the interior and exterior surfaces of a charged vesicle in an electrolyte solution. The criteria for adsorption and the density profiles of the adsorbed polymer chain are derived in terms of various characteristics of the polymer, vesicle, and medium, such as the charge density and length of the polymer, charge density and size of the vesicle, electrolyte concentration and dielectric constant of the medium. For adsorption inside the vesicle, the competition between the loss of conformational entropy and gain in adsorption energy results in two kinds of encapsulated states, depending on the strength of the polymer-vesicle interaction. By considering also the adsorption from outside the vesicle, we derive the entropic and energy contributions to the free energy change to transfer an adsorbed chain in the interior to an adsorbed chain on the exterior. In this paper, we have used the Wentzel-Kramers-Brillouin (WKB) method to solve the equation for the probability distribution function of the chain. The present WKB results are compared with the previous results based on variational methods. The WKB and variational results are in good agreement for both the interior and exterior states of adsorption, except in the zero-salt limit for adsorption in the exterior region. The adsorption criteria and density profiles for both the interior and exterior states are presented in terms of various experimentally controllable variables. Calculation of the dependencies of free energy change to transfer an adsorbed chain from the interior to the exterior surface on salt concentration and vesicle radius shows that the free energy penalty to expel a chain from a vesicle is only of the order of thermal energy.

Compactness of viral genomes: effect of disperse and localized random mutations

Genomes of single-stranded RNA viruses have evolved to optimize several concurrent properties. One of them is the architecture of their genomic folds, which must not only feature precise structural elements at specific positions, but also allow for overall spatial compactness. The latter was shown to be disrupted by random synonymous mutations, a disruption which can consequently negatively affect genome encapsidation. In this study, we use three mutation schemes with different degrees of locality to mutate the genomes of phage MS2 and Brome Mosaic virus in order to understand the observed sensitivity of the global compactness of their folds. We find that mutating local stretches of their genomes' sequence or structure is less disruptive to their compactness compared to inducing randomly-distributed mutations. Our findings are indicative of a mechanism for the conservation of compactness acting on a global scale of the genomes, and have several implications for understanding the interplay between local and global architecture of viral RNA genomes.

From discrete to continuous description of spherical surface charge distributions

The importance of electrostatic interactions in soft matter and biological systems can often be traced to non-uniform charge effects, which are commonly described using a multipole expansion of the corresponding charge distribution. The standard approach when extracting the charge distribution of a given system is to treat the constituent charges as points. This can, however, lead to an overestimation of multipole moments of high order, such as dipole, quadrupole, and higher moments. Focusing on distributions of charges located on a spherical surface - characteristic of numerous biological macromolecules, such as globular proteins and viral capsids, as well as of inverse patchy colloids - we develop a novel way of representing spherical surface charge distributions based on the von Mises-Fisher distribution. This approach takes into account the finite spatial extension of individual charges, and leads to a simple yet powerful way of describing surface charge distributions and their multipole expansions. In this manner, we analyze charge distributions and the derived multipole moments of a number of different spherical configurations of identical charges with various degrees of symmetry. We show how the number of charges, their size, and the geometry of their configuration influence the behavior and relative importance of multipole magnitudes of different order. Importantly, we clearly demonstrate how neglecting the effect of charge size leads to an overestimation of high-order multipoles. The results of our work can be applied to construct analytical models of electrostatic interactions and multipole expansion of charged particles in diverse soft matter and biological systems.

Varieties of charge distributions in coat proteins of ssRNA+ viruses

A major part of the interactions involved in the assembly and stability of icosahedral, positive-sense single-stranded RNA (ssRNA+) viruses is electrostatic in nature, as can be inferred from the strong pH- and salt-dependence of their assembly phase diagrams. Electrostatic interactions do not act only between the capsid coat proteins (CPs), but just as often provide a significant contribution to the interactions of the CPs with the genomic RNA, mediated to a large extent by positively charged, flexible N-terminal tails of the CPs. In this work, we provide two clear and complementary definitions of an N-terminal tail of a protein, and use them to extract the tail sequences of a large number of CPs of ssRNA+  viruses. We examine the pH-dependent interplay of charge on both tails and CPs alike, and show that-in contrast to the charge on the CPs-the net positive charge on the N-tails persists even to very basic pH values. In addition, we note a limit to the length of the wild-type genomes of those viruses which utilize positively charged tails, when compared to viruses without charged tails and similar capsid size. At the same time, we observe no clear connection between the charge on the N-tails and the genome lengths of the viruses included in our study.

Charge regulation of nonpolar colloids

Individual colloids often carry a charge as a result of the dissociation (or adsorption) of weakly-ionized surface groups. The magnitude depends on the precise chemical environment surrounding a particle, which in a concentrated dispersion is a function of the colloid packing fraction η. Theoretical studies have suggested that the effective charge Z_eff in regulated systems could, in general, decrease with increasing η. We test this hypothesis for nonpolar dispersions by determining Z_eff(η) over a wide range of packing fractions (10^-5 ≤ η ≤ 0.3) using a combination of small-angle X-ray scattering and electrophoretic mobility measurements. All dispersions remain entirely in the fluid phase regime. We find a complex dependence of the particle charge as a function of the packing fraction, with Z_eff initially decreasing at low concentrations before finally increasing at high η. We attribute the non-monotonic density dependence to a crossover from concentration-independent screening at low η, to a high packing fraction regime in which counterions outnumber salt ions and electrostatic screening becomes η-dependent. The efficiency of charge stabilization at high concentrations may explain the unusually high stability of concentrated nanoparticle dispersions which has been reported.

50th Anniversary Perspective: A Perspective on Polyelectrolyte Solutions

From the beginning of life with the information-containing polymers until the present era of a plethora of water-based materials in health care industry and biotechnology, polyelectrolytes are ubiquitous with a broad range of structural and functional properties. The main attribute of polyelectrolyte solutions is that all molecules are strongly correlated both topologically and electrostatically in their neutralizing background of charged ions in highly polarizable solvent. These strong correlations and the necessary use of numerous variables in experiments on polyelectrolytes have presented immense challenges toward fundamental understanding of the various behaviors of charged polymeric systems. This Perspective presents the author's subjective summary of several conceptual advances and the remaining persistent challenges in the contexts of charge and size of polymers, structures in homogeneous solutions, thermodynamic instability and phase transitions, structural evolution with oppositely charged polymers, dynamics in polyelectrolyte solutions, kinetics of phase separation, mobility of charged macromolecules between compartments, and implications to biological systems.

pH Dependence of Charge Multipole Moments in Proteins

Electrostatic interactions play a fundamental role in the structure and function of proteins. Due to ionizable amino acid residues present on the solvent-exposed surfaces of proteins, the protein charge is not constant but varies with the changes in the environment-most notably, the pH of the surrounding solution. We study the effects of pH on the charge of four globular proteins by expanding their surface charge distributions in terms of multipoles. The detailed representation of the charges on the proteins is in this way replaced by the magnitudes and orientations of the multipole moments of varying order. Focusing on the three lowest-order multipoles-the total charge, dipole, and quadrupole moment-we show that the value of pH influences not only their magnitudes, but more notably and importantly also the spatial orientation of their principal axes. Our findings imply important consequences for the study of protein-protein interactions and the assembly of both proteinaceous shells and patchy colloids with dissociable charge groups.

Role of Condensing Particles in Polymer Confinement: A Model for Virus-Packed "Minichromosomes"

Confined mixtures of a polymer and nonspecifically binding particles (condensers) are studied as models for viruses containing double-stranded DNA (polymer) and condensing proteins (particles). We explore a model in which all interactions between the packed content (polymer and particles) and its confinement are purely repulsive, with only a short-range attraction between the condensers and polymer to simulate binding. In the range of physical parameters applicable to viruses, the model predicts reduction of pressure in the system effected by the condensers, despite the reduction in free volume. Condensers are found to be interspersed throughout the spherical confinement and only partially wrapped in the polymer, which acts as an effective medium for the condenser interactions. Crowding of the viral interior influences the DNA and protein organization, producing a picture inconsistent with a chromatin-like, beads-on-a-string structure. The model predicts an organization of the confined interior compatible with experimental data on unperturbed adenoviruses and polyomaviruses, at the same time providing insight into the role of condensing proteins in the viral infectious cycles of related viral families.

Impact of Hydrogen Bonding in the Binding Site between Capsid Protein and MS2 Bacteriophage ssRNA

MS2 presents a well-studied example of a single-stranded RNA virus for which the genomic RNA plays a pivotal role in the virus assembly process based on the packaging signal-mediated mechanism. Packaging signals (PSs) are multiple dispersed RNA sequence/structure motifs varying around a central recognition motif that interact in a specific way with the capsid protein in the assembly process. Although the discovery and identification of these PSs was based on bioinformatics and geometric approaches, in tandem with sophisticated experimental protocols, we approach this problem using large-scale ab initio computation centered on critical aspects of the consensus protein-RNA interactions recognition motif. DFT calculations are carried out on two nucleoprotein complexes: wild-type and mutated (PDB IDs: 1ZDH and 5MSF ). The calculated partial charge distribution of residues and the strength of hydrogen bonding (HB) between them enabled us to locate the exact binding sites with the strongest HBs, identified to be LYS43-A^-4, ARG49-C^-13, TYR85-C^-5, and LYS61-C^-5, due to the change in the sequence of the mutated RNA.