Odd viscosity-induced Hall-like transport of an active chiral fluid.

Broken time-reversal and parity symmetries in active spinner fluids imply a nondissipative "odd viscosity," engendering phenomena unattainable in traditional passive or active fluids. Here we show that the odd viscosity itself can lead to a Hall-like transport when the active chiral fluid flows through a quenched matrix of obstacles, reminiscent of the anomalous Hall effect. The Hall-like velocity depends significantly on the spinner activity and longitudinal flow due to the interplay between odd viscosity and spinner-obstacle collisions. Our findings underscore the importance of odd viscosity in active chiral matter and elucidate its essential role in the anomalous Hall-like effect.

Erratum: Spontaneous Domain Formation in Spherically Confined Elastic Filaments [Phys. Rev. Lett. 123, 047801 (2019)].

This corrects the article DOI: 10.1103/PhysRevLett.123.047801.

On the nature of screening in charge-regulated macroion solutions.

We present a derivation of the screening length for a solution containing a charge-regulated macroion, e.g. protein, with its counterions. We show that it can be obtained directly from the second derivatives of the total free energy by taking recourse to the "uncertainty relation" of the Legendre transform, which connects the Hessians or the local curvatures of the free energy as a function of density and its Legendre transform, i.e., osmotic pressure, as a function of chemical potentials. Based on the Fowler-Guggenheim-Frumkin model of charge regulation, we then analyze the "screening resonance" and the "overscreening" of the screening properties of the charge-regulated macroion solution.

Close packings of identical proteins in small spherical capsids and similar proteinaceous shells.

Understanding the principles governing protein arrangement in viral capsids and structurally similar protein shells can enable the development of new antiviral strategies and the design of artificial protein cages for various applications. We study these principles within the context of the close packing problem, by analyzing dozens of small spherical shells assembled from a single type of protein. First, we use icosahedral spherical close packings containing 60T identical disks, where T ≤ 4, to rationalize the protein arrangement in twenty real icosahedral shells both satisfying and violating the paradigmatic Caspar-Klug model. We uncover a striking correspondence between the protein mass centers in the considered shells and the centers of disks in the close packings. To generalize the packing model, we consider proteins with a weak shape anisotropy and propose an interaction energy, minimization of which allows us to obtain spherical dense packings of slightly anisotropic structural units. In the case of strong anisotropy, we model the proteins as sequences of overlapping discs of different sizes, with minimum energy configuration not only resulting in packings, accurately reproducing locations and orientations of individual proteins, but also revealing that icosahedral packings that display the handedness of real capsids are energetically more favorable. Finally, by introducing effective disc charges, we rationalize the formation of inter-protein bonds in protein shells.

Electrostatic interactions between charge regulated spherical macroions.

We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte-thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions.

Curvature effects in charge-regulated lipid bilayers.

We formulate a theory of electrostatic interactions in lipid bilayer membranes where both monolayer leaflets contain dissociable moieties that are subject to charge regulation. We specifically investigate the coupling between membrane curvature and charge regulation of a lipid bilayer vesicle using both the linear Debye-Hückel (DH) and the non-linear Poisson-Boltzmann (PB) theory. We find that charge regulation of an otherwise symmetric bilayer membrane can induce charge symmetry breaking, non-linear flexoelectricity and anomalous curvature dependence of free energy. The pH effects investigated go beyond the paradigm of electrostatic renormalization of the mechano-elastic properties of membranes.

On the nature of screening in Voorn-Overbeek type theories.

By using a recently formulated Legendre transform approach to the thermodynamics of charged systems, we explore the general form of the screening length in the Voorn-Overbeek-type theories, which remains valid also in the cases where the entropy of the charged component(s) is not given by the ideal gas form as in the Debye-Hückel theory. The screening length consistent with the non-electrostatic terms in the free energy ansatz for the Flory-Huggins and Voorn-Overbeek type theories, derived from the local curvature properties of the Legendre transform, has distinctly different behavior than the often invoked standard Debye screening length, though it reduces to it in some special cases.

Mutations of Omicron Variant at the Interface of the Receptor Domain Motif and Human Angiotensin-Converting Enzyme-2.

The most recent Omicron variant of SARS-CoV-2 has caused global concern and anxiety. The only thing certain about this strain, with a large number of mutations in the spike protein, is that it spreads quickly, seems to evade immune defense, and mitigates the benefits of existing vaccines. Based on the ultra-large-scale ab initio computational modeling of the receptor binding motif (RBM) and the human angiotensin-converting enzyme-2 (ACE2) interface, we provide the details of the effect of Omicron mutations at the fundamental atomic scale level. In-depth analysis anchored in the novel concept of amino acid-amino acid bond pair units (AABPU) indicates that mutations in the Omicron variant are connected with (i) significant changes in the shape and structure of AABPU components, together with (ii) significant increase in the positive partial charge, which facilitates the interaction with ACE2. We have identified changes in bonding due to mutations in the RBM. The calculated bond order, based on AABPU, reveals that the Omicron mutations increase the binding strength of RBM to ACE2. Our findings correlate with and are instrumental to explain the current observations and can contribute to the prediction of next potential new variant of concern.

Effect of Delta and Omicron Mutations on the RBD-SD1 Domain of the Spike Protein in SARS-CoV-2 and the Omicron Mutations on RBD-ACE2 Interface Complex.

The receptor-binding domain (RBD) is the essential part in the Spike-protein (S-protein) of SARS-CoV-2 virus that directly binds to the human ACE2 receptor, making it a key target for many vaccines and therapies. Therefore, any mutations at this domain could affect the efficacy of these treatments as well as the viral-cell entry mechanism. We introduce ab initio DFT-based computational study that mainly focuses on two parts: (1) Mutations effects of both Delta and Omicron variants in the RBD-SD1 domain. (2) Impact of Omicron RBD mutations on the structure and properties of the RBD-ACE2 interface system. The in-depth analysis is based on the novel concept of amino acid-amino acid bond pair units (AABPU) that reveal the differences between the Delta and/or Omicron mutations and its corresponding wild-type strain in terms of the role played by non-local amino acid interactions, their 3D shapes and sizes, as well as contribution to hydrogen bonding and partial charge distributions. Our results also show that the interaction of Omicron RBD with ACE2 significantly increased its bonding between amino acids at the interface providing information on the implications of penetration of S-protein into ACE2, and thus offering a possible explanation for its high infectivity. Our findings enable us to present, in more conspicuous atomic level detail, the effect of specific mutations that may help in predicting and/or mitigating the next variant of concern.

Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors.

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a -1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.

Electrostatic interactions between the SARS-CoV-2 virus and a charged electret fibre.

While almost any kind of face mask offers some protection against particles and pathogens of different sizes, the most efficient ones make use of a layered structure where one or more layers are electrically charged. These electret layers are essential for the efficient filtration of difficult-to-capture small particles, yet the exact nature of electrostatic capture with respect to the charge on both the particles and the electret fibres as well as the effect of the immediate environment remain unclear. Here, we explore in detail the electrostatic interactions between the surface of a single charged electret fibre and a model of the SARS-CoV-2 virus. Using Poisson-Boltzmann electrostatics coupled to a detailed spike protein charge regulation model, we show how pH and salt concentration drastically change both the scale and the sign of the interaction. Furthermore, the configuration of the few spike proteins closest to the electret fibre turns out to be as important for the strength of the interaction as their total number on the virus envelope, a direct consequence of spike protein charge regulation. The results of our work elucidate the details of virus electrostatics and contribute to the general understanding of efficient virus filtration mechanisms.

Key Interacting Residues between RBD of SARS-CoV-2 and ACE2 Receptor: Combination of Molecular Dynamics Simulation and Density Functional Calculation.

The spike protein of SARS-CoV-2 binds to the ACE2 receptor via its receptor-binding domain (RBD), with the RBD-ACE2 complex presenting an essential molecular target for vaccine development to stall the virus infection proliferation. The computational analyses at molecular, amino acid (AA), and atomic levels have been performed systematically to identify the key interacting AAs in the formation of the RBD-ACE2 complex for SARS-CoV and SARS-CoV-2 with its Alpha and Beta variants. Our study uses the molecular dynamics (MD) simulations with the molecular mechanics generalized Born surface area (MM-GBSA) method to predict the binding free energy (BFE) and to determine the actual interacting AAs, as well as two ab initio quantum chemical protocols based on the density functional theory (DFT) implementation. Based on MD results, Q493, Y505, Q498, N501, T500, N487, Y449, F486, K417, Y489, F456, Y495, and L455 have been identified as hotspots in SARS-CoV-2 RBD, while those in ACE2 are K353, K31, D30, D355, H34, D38, Q24, T27, Y83, Y41, and E35. RBD with Alpha and Beta variants has slightly different interacting AAs due to N501Y mutation. Both the electrostatic and hydrophobic interactions are the main driving force to form the AA-AA binding pairs. We confirm that Q493, Q498, N501, F486, K417, and F456 in RBD are the key residues responsible for the tight binding of SARS-CoV-2 with ACE2 compared to SARS-CoV. RBD with the Alpha variant binds with ACE2 stronger than the wild-type RBD or Beta. In the Beta variant, K417N reduces the binding, E484K slightly enhances it, and N501Y significantly increases it as in Alpha. The DFT results reveal that N487, Q493, Y449, T500, G496, G446, and G502 in RBD of SARS2 form pairs via specific hydrogen bonding with Q24, H34, E35, D38, Y41, Q42, and K353 in ACE2.

Contribution of dipolar bridging to phospholipid membrane interactions: A mean-field analysis.

We develop a model of interacting zwitterionic membranes with rotating surface dipoles immersed in a monovalent salt and implement it in a field theoretic formalism. In the mean-field regime of monovalent salt, the electrostatic forces between the membranes are characterized by a non-uniform trend: at large membrane separations, the interfacial dipoles on the opposing sides behave as like-charge cations and give rise to repulsive membrane interactions; at short membrane separations, the anionic field induced by the dipolar phosphate groups sets the behavior in the intermembrane region. The attraction of the cationic nitrogens in the dipolar lipid headgroups leads to the adhesion of the membrane surfaces via dipolar bridging. The underlying competition between the opposing field components of the individual dipolar charges leads to the non-uniform salt ion affinity of the zwitterionic membrane with respect to the separation distance; large inter-membrane separations imply anionic excess, while small nanometer-sized separations favor cationic excess. This complex ionic selectivity of zwitterionic membranes may have relevant repercussions on nanofiltration and nanofluidic transport techniques.

Irreversible and reversible morphological changes in the φ6 capsid and similar viral shells: symmetry and micromechanics.

Understanding the physicochemical processes occurring in viruses during their maturation is of fundamental importance since only mature viruses can infect host cells. Here we consider the irreversible and reversible morphological changes that occur with the dodecahedral φ6 procapsid during the sequential packaging of 3 RNA segments forming the viral genome. It is shown that the dodecahedral shape of all the four observed capsid states is perfectly reproduced by a sphere radially deformed by only two irreducible spherical harmonics with icosahedral symmetry and wave numbers l = 6 and l = 10. The rotation of proteins around the 3-fold axes at the Procapsid → Intermediate 1 irreversible transformation is in fact also well described with the shear field containing only two irreducible harmonics with the same two wave numbers. The high stability of the Intermediate 1 state is discussed and the shapes of the Intermediate 2 state and Capsid (reversibly transforming back to the Intermediate 1 state) are shown to be mainly due to the isotropic pressure that the encapsidated RNA segments exert on the shell walls. The hidden symmetry of the capsid and the physicochemical features of the in vitro genome extraction from the viral shell are also elucidated.

Intra- and intermolecular atomic-scale interactions in the receptor binding domain of SARS-CoV-2 spike protein: implication for ACE2 receptor binding.

The COVID-19 pandemic poses a severe threat to human health with unprecedented social and economic disruption. Spike (S) glycoprotein in the SARS-CoV-2 virus is pivotal in understanding the virus anatomy, since it initiates the early contact with the ACE2 receptor in the human cell. The subunit S1 in chain A of S-protein has four structural domains: the receptor binding domain (RBD), the n-terminal domain (NTD) and two subdomains (SD1, SD2). We report details of the intra- and inter-molecular binding mechanism of RBD using density functional theory, including electronic structure, interatomic bonding and partial charge distribution. We identify five strong hydrogen bonds and analyze their roles in binding. This provides a pathway to a quantum-chemical understanding of the interaction between the S-protein and the ACE2 receptor with insights into the function of conserved features in the ACE2 receptor binding domain that could inform vaccine and drug development.

Critical behavior of charge-regulated macro-ions.

Based on a collective description of electrolytes composed of charge-regulated macro-ions and simple salt ions, we analyze their equilibrium charge state in the bulk and their behavior in the vicinity of an external electrified surface. The mean-field formulation of mobile macro-ions in an electrolyte bathing solution is extended to include interactions between association and dissociation sites. We demonstrate that above a critical concentration of salt, similar to the critical micelle concentration, a non-trivial distribution of charge states sets in. Such a charge state can eventually lead to a liquid-liquid phase separation based on charge regulation.

Ordering of adsorbed rigid rods mediated by the Boussinesq interaction on a soft substrate.

Orientational ordering driven by mechanical distortion of soft substrates plays a major role in material transformation processes such as elastocapillarity and surface anchoring. We present a theoretical model of the orientational response of anisotropic rods deposited onto a surface of a soft, elastic substrate of finite thickness. We show that anisotropic rods exhibit a continuous isotropic-nematic phase transition, driven by orientational interactions between surface deposited rods. This interaction is mediated by the deformation of the underlying elastic substrate and is quantified by the Boussinesq solution adapted to the case of slender, surface deposited rods. From the microscopic rod-rod interactions, we derive the appropriate Maier-Saupe mean-field description, which includes the Boussinesq elastic free energy contribution due to the substrate elasticity, derive the conditions for the existence of a continuous orientational ordering transition, and discuss the implication of results in the soft (bio)system context.

Spontaneous Domain Formation in Spherically Confined Elastic Filaments.

Although the free energy of a genome packing into a virus is dominated by DNA-DNA interactions, ordering of the DNA inside the capsid is elasticity driven, suggesting general solutions with DNA organized into spool-like domains. Using analytical calculations and computer simulations of a long elastic filament confined to a spherical container, we show that the ground state is not a single spool as assumed hitherto, but an ordering mosaic of multiple homogeneously ordered domains. At low densities, we observe concentric spools, while at higher densities, other morphologies emerge, which resemble topological links. We discuss our results in the context of metallic wires, viral DNA, and flexible polymers.

Casimir-like interactions and surface anchoring duality in bookshelf geometry of smectic-A liquid crystals.

We analyze the transverse intersubstrate Casimir-like force, arising as a result of thermal fluctuations of the liquid crystalline layers of a smectic-A film confined between two planar substrates in a bookshelf geometry, in which the equidistant smectic layers are placed perpendicular to the bounding surfaces. We discuss the variation of the interaction force as a function of the intersubstrate separation in the presence of surface anchoring to the substrates, showing that the force induced by confined fluctuations is attractive and depends on the penetration length as well as the layer spacing. The strongest effect occurs for tightly confined fluctuations, in which the surface anchoring energy is set to infinity, where the force per area scales linearly with the thermal energy and inversely with the fourth power of the intersubstrate separation. By reducing the strength of the surface anchoring energy, the force first becomes weaker in magnitude but then increases in magnitude as the surface anchoring strength is further reduced down to zero, in which case the force tends to that obtained in the limit of strong anchoring.

Hidden symmetry of the anomalous bluetongue virus capsid and its role in the infection process.

Clear understanding of the principles that control the arrangement of proteins and their self-assembly into viral shells is very important for the development of antiviral strategies. Here we consider the structural peculiarities and hidden symmetry of the anomalous bluetongue virus (BTV) capsid. Each of its three concentric shells violates the paradigmatic geometrical model of Caspar and Klug, which is otherwise well suited to describe most of the known icosahedral viral shells. As we show, three icosahedral spherical lattices, which are commensurate with each other and possess locally hexagonal (primitive or honeycomb) order, underlie the proteinaceous shells of the BTV capsid. This interpretation of the multishelled envelope allows us to discuss the so-called "symmetry mismatch" between its layers. We also analyze the structural stability of the considered spherical lattices on the basis of the classical theory of spherical packing and relate the proximity of the outer spherical lattice to destabilization with the fact that during infection of the cell VP2 trimers are detached from the surface of the BTV capsid. An electrostatic mechanism that can assist in this detachment is discussed in detail.