Coupling finite and boundary element methods to solve the Poisson-Boltzmann equation for electrostatics in molecular solvation

The Poisson-Boltzmann equation is widely used to model electrostatics in molecular systems. Available software packages solve it using finite difference, finite element, and boundary element methods, where the latter is attractive due to the accurate representation of the molecular surface and partial charges, and exact enforcement of the boundary conditions at infinity. However, the boundary element method is limited to linear equations and piecewise constant variations of the material properties. In this work, we present a scheme that couples finite and boundary elements for the linearised Poisson-Boltzmann equation, where the finite element method is applied in a confined solute region and the boundary element method in the external solvent region. As a proof-of-concept exercise, we use the simplest methods available: Johnson-Nédélec coupling with mass matrix and diagonal preconditioning, implemented using the Bempp-cl and FEniCSx libraries via their Python interfaces. We showcase our implementation by computing the polar component of the solvation free energy of a set of molecules using a constant and a Gaussian-varying permittivity. As validation, we compare against well-established finite difference solvers for an extensive binding energy data set, and with the finite difference code APBS (to 0.5%) for Gaussian permittivities. We also show scaling results from protein G B1 (955 atoms) up to immunoglobulin G (20,148 atoms). For small problems, the coupled method was efficient, outperforming a purely boundary integral approach. For Gaussian-varying permittivities, which are beyond the applicability of boundary elements alone, we were able to run medium to large-sized problems on a single workstation. The development of better preconditioning techniques and the use of distributed memory parallelism for larger systems remains an area for future work. We hope this work will serve as inspiration for future developments that consider space-varying field parameters, and mixed linear-nonlinear schemes for molecular electrostatics with implicit solvent models.

Minor electrostatic changes robustly increase VP40 membrane binding, assembly, and budding of Ebola virus matrix protein derived virus-like particles

Ebola virus (EBOV) is a filamentous negative-sense RNA virus, which causes severe hemorrhagic fever. There are limited vaccines or therapeutics for prevention and treatment of EBOV, so it is important to get a detailed understanding of the virus lifecycle to illuminate new drug targets. EBOV encodes for the matrix protein, VP40, which regulates assembly and budding of new virions from the inner leaflet of the host cell plasma membrane (PM). In this work, we determine the effects of VP40 mutations altering electrostatics on PM interactions and subsequent budding. VP40 mutations that modify surface electrostatics affect viral assembly and budding by altering VP40 membrane-binding capabilities. Mutations that increase VP40 net positive charge by one (e.g., Gly to Arg or Asp to Ala) increase VP40 affinity for phosphatidylserine and phosphatidylinositol 4,5-bisphosphate in the host cell PM. This increased affinity enhances PM association and budding efficiency leading to more effective formation of virus-like particles. In contrast, mutations that decrease net positive charge by one (e.g., Gly to Asp) lead to a decrease in assembly and budding because of decreased interactions with the anionic PM. Taken together, our results highlight the sensitivity of slight electrostatic changes on the VP40 surface for assembly and budding. Understanding the effects of single amino acid substitutions on viral budding and assembly will be useful for explaining changes in the infectivity and virulence of different EBOV strains, VP40 variants that occur in nature, and for long-term drug discovery endeavors aimed at EBOV assembly and budding.

Bridging Eulerian and Lagrangian Poisson-Boltzmann solvers by ESES

The Poisson-Boltzmann (PB) model is a widely used electrostatic model for biomolecular solvation analysis. Formulated as an elliptic interface problem, the PB model can be numerically solved on either Eulerian meshes using finite difference/finite element methods or Lagrangian meshes using boundary element methods. Molecular surface generators, which produce the discretized dielectric interfaces between solutes and solvents, are critical factors in determining the accuracy and efficiency of the PB solvers. In this work, we investigate the utility of the Eulerian Solvent Excluded Surface (ESES) software for rendering conjugated Eulerian and Lagrangian surface representations, which enables us to numerically validate and compare the quality of Eulerian PB solvers, such as the MIBPB solver, and the Lagrangian PB solvers, such as the TABI-PB solver. Furthermore, with the ESES software and its associated PB solvers, we are able to numerically validate an interesting and useful but often neglected source-target symmetric property associated with the linearized PB model.

A Universal Approximation for Conductance Blockade in Thin Nanopore Membranes

Nanopore-based sensing platforms have transformed single-molecule detection and analysis. The foundation of nanopore translocation experiments lies in conductance measurements, yet existing models, which are largely phenomenological, are inaccurate in critical experimental conditions such as thin and tightly fitting pores. Of the two components of the conductance blockade, channel and access resistance, the access resistance is poorly modeled. We present a comprehensive investigation of the access resistance and associated conductance blockade in thin nanopore membranes. By combining a first-principles approach, multiscale modeling, and experimental validation, we propose a unified theoretical modeling framework. The analytical model derived as a result surpasses current approaches across a broad parameter range. Beyond advancing our theoretical understanding, our framework's versatility enables analyte size inference and predictive insights into conductance blockade behavior. Our results will facilitate the design and optimization of nanopore devices for diverse applications, including nanopore base calling and data storage.

Ionic strength alters crosslinker-driven self-organization of microtubules

The microtubule cytoskeleton is a major structural element inside cells that directs self-organization using microtubule-associated proteins and motors. It has been shown that finite-sized, spindle-like microtubule organizations, called "tactoids," can form in vitro spontaneously from mixtures of tubulin and the antiparallel crosslinker, MAP65, from the MAP65/PRC1/Ase family. Here, we probe the ability of MAP65 to form tactoids as a function of the ionic strength of the buffer to attempt to break the electrostatic interactions binding MAP65 to microtubules and inter-MAP65 binding. We observe that, with increasing monovalent salts, the organizations change from finite tactoids to unbounded length bundles, yet the MAP65 binding and crosslinking appear to stay intact. We further explore the effects of ionic strength on the dissociation constant of MAP65 using both microtubule pelleting and single-molecule binding assays. We find that salt can reduce the binding, yet salt never negates it. Instead, we believe that the salt is affecting the ability of the MAP65 to form phase-separated droplets, which cause the nucleation and growth of tactoids, as recently demonstrated.

Tuning Electrostatic Interactions To Control Orientation of GFP Protein Adsorption on Silica Surface

The adsorption of green fluorescent protein (GFP) on silica surfaces has been the subject of growing interest due to its potential applications in various fields, including biotechnology and biomedicine. In this study, we used all-atom molecular dynamics simulations to investigate the charge-driven adsorption of wild type GFP and its supercharged variants on silica surfaces. The results showed that the positively charged variant of GFP adsorbed on the negatively charged silica surface with minimal loss in its secondary structure. Further studies were conducted to understand the role of surface charge distribution on two other positively charged variants of GFP, and the results showed that the orientation of GFP on silica can be easily tuned by careful mutations of the charged amino acid residues on the GFP. This study provides valuable molecular insights into the role of electrostatic-driven adsorption of GFP and highlights the importance of charge interactions in the adsorption process.

Maximizing Performance and Stability of Organic Solar Cells at Low Driving Force for Charge Separation

Thanks to the development of novel electron acceptor materials, the power conversion efficiencies (PCE) of organic photovoltaic (OPV) devices are now approaching 20%. Further improvement of PCE is complicated by the need for a driving force to split strongly bound excitons into free charges, causing voltage losses. This review discusses recent approaches to finding efficient OPV systems with minimal driving force, combining near unity quantum efficiency (maximum short circuit currents) with optimal energy efficiency (maximum open circuit voltages). The authors discuss apparently contradicting results on the amount of exciton binding in recent literature, and approaches to harmonize the findings. A comprehensive view is then presented on motifs providing a driving force for charge separation, namely hybridization at the donor:acceptor interface and polarization effects in the bulk, of which quadrupole moments (electrostatics) play a leading role. Apart from controlling the energies of the involved states, these motifs also control the dynamics of recombination processes, which are essential to avoid voltage and fill factor losses. Importantly, all motifs are shown to depend on both molecular structure and process conditions. The resulting high dimensional search space advocates for high throughput (HT) workflows. The final part of the review presents recent HT studies finding consolidated structure-property relationships in OPV films and devices from various deposition methods, from research to industrial upscaling.

Minor changes in electrostatics robustly increase VP40 membrane binding, assembly, and budding of Ebola virus matrix protein derived virus-like particles

Ebola virus (EBOV) is a filamentous negative-sense RNA virus which causes severe hemorrhagic fever. There are limited vaccines or therapeutics for prevention and treatment of EBOV, so it is important to get a detailed understanding of the virus lifecycle to illuminate new drug targets. EBOV encodes for the matrix protein, VP40, which regulates assembly and budding of new virions from the inner leaflet of the host cell plasma membrane (PM). In this work we determine the effects of VP40 mutations altering electrostatics on PM interactions and subsequent budding. VP40 mutations that modify surface electrostatics affect viral assembly and budding by altering VP40 membrane binding capabilities. Mutations that increase VP40 net positive charge by one (e.g., Gly to Arg or Asp to Ala) increase VP40 affinity for phosphatidylserine (PS) and PI(4,5)P2 in the host cell PM. This increased affinity enhances PM association and budding efficiency leading to more effective formation of virus-like particles (VLPs). In contrast, mutations that decrease net positive charge by one (e.g., Gly to Asp) lead to a decrease in assembly and budding because of decreased interactions with the anionic PM. Taken together our results highlight the sensitivity of slight electrostatic changes on the VP40 surface for assembly and budding. Understanding the effects of single amino acid substitutions on viral budding and assembly will be useful for explaining changes in the infectivity and virulence of different EBOV strains, VP40 variants that occur in nature, and for long-term drug discovery endeavors aimed at EBOV assembly and budding.

Electrostatically Doped Junctionless Graphene Nanoribbon Tunnel Field-Effect Transistor for High-Performance Gas Sensing Applications: Leveraging Doping Gates for Multi-Gas Detection

In this paper, a new junctionless graphene nanoribbon tunnel field-effect transistor (JLGNR TFET) is proposed as a multi-gas nanosensor. The nanosensor has been computationally assessed using a quantum simulation based on the self-consistent solutions of the mode space non-equilibrium Green's function (NEGF) formalism coupled with the Poisson's equation considering ballistic transport conditions. The proposed multi-gas nanosensor is endowed with two top gates ensuring both reservoirs' doping and multi-gas sensing. The investigations have included the IDS-VGS transfer characteristics, the gas-induced electrostatic modulations, subthreshold swing, and sensitivity. The order of change in drain current has been considered as a sensitivity metric. The underlying physics of the proposed JLGNR TFET-based multi-gas nanosensor has also been studied through the analysis of the band diagrams behavior and the energy-position-resolved current spectrum. It has been found that the gas-induced work function modulation of the source (drain) gate affects the n-type (p-type) conduction branch by modulating the band-to-band tunneling (BTBT) while the p-type (n-type) conduction branch still unaffected forming a kind of high selectivity from operating regime point of view. The high sensitivity has been recorded in subthermionic subthreshold swing (SS < 60 mV/dec) regime considering small gas-induced gate work function modulation. In addition, advanced simulations have been performed for the detection of two different types of gases separately and simultaneously, where high-performance has been recorded in terms of sensitivity, selectivity, and electrical behavior. The proposed detection approach, which is viable, innovative, simple, and efficient, can be applied using other types of junctionless tunneling field-effect transistors with emerging channel nanomaterials such as the transition metal dichalcogenides materials. The proposed JLGNRTFET-based multi-gas nanosensor is not limited to two specific gases but can also detect other gases by employing appropriate gate materials in terms of selectivity.

Electrostatic properties of human germlines and biodistribution of small biologics

Off-target biodistribution of biologics bears important toxicological consequences. Antibody fragments intended for use as vectors of cytotoxic payloads (e.g. antibody-drug conjugates, radiotherapy) can accumulate at clearance organs like kidneys and liver, where they can cause dose-limiting toxicities. Renal and hepatic uptakes are known to be affected by protein electrostatics, which promote protein internalization through pinocytosis. Using minibodies as a model of an antibody fragment lacking FcRn recycling, we compared the biodistributions of leads with different degrees of accumulation at the kidney and liver. We identified a positive electrostatic patch highly conserved in a germline family very commonly used in the humanization of approved biologics. Neutralization of this patch led to a drastic reduction in the kidney uptake, leading to a biodistribution more favorable to the delivery of highly cytotoxic payloads. Next, we conducted a high throughput study of the electrostatic properties for all combinations of VH and VL germlines. This analysis shows how different VH/VL combinations exhibit varying tendencies to create electrostatic patches, resulting in Fv variants with different isoelectric points. Our work emphasizes the importance of carefully selecting germlines for humanization with optimal electrostatic properties in order to control the unspecific tissue uptake of low molecular weight biologics.

A conserved electrostatic membrane-binding surface in synaptotagmin-like proteins revealed using molecular phylogenetic analysis and homology modeling

Protein structure prediction has emerged as a core technology for understanding biomolecules and their interactions. Here, we combine homology-based structure prediction with molecular phylogenetic analysis to study the evolution of electrostatic membrane binding among the vertebrate synaptotagmin-like protein (Slp) family. Slp family proteins play key roles in the membrane trafficking of large dense-core secretory vesicles. Our previous experimental and computational study found that the C2A domain of Slp-4 (also called granuphilin) binds with high affinity to anionic phospholipids in the cytoplasmic leaflet of the plasma membrane through a large positively charged protein surface centered on a cluster of phosphoinositide-binding lysine residues. Because this surface contributes greatly to Slp-4 C2A domain membrane binding, we hypothesized that the net charge on the surface might be evolutionarily conserved. To test this hypothesis, the known C2A sequences of Slp-4 among vertebrates were organized by class (from mammalia to pisces) using molecular phylogenetic analysis. Consensus sequences for each class were then identified and used to generate homology structures, from which Poisson-Boltzmann electrostatic potentials were calculated. For comparison, homology structures and electrostatic potentials were also calculated for the five human Slp protein family members. The results demonstrate that the charge on the membrane-binding surface is highly conserved throughout the evolution of Slp-4, and more highly conserved than many individual residues among the human Slp family paralogs. Such molecular phylogenetic-driven computational analysis can help to describe the evolution of electrostatic interactions between proteins and membranes which are crucial for their function.

A Conserved Electrostatic Membrane-Binding Surface in Synaptotagmin-Like Proteins Revealed Using Molecular Phylogenetic Analysis and Homology Modeling

Protein structure prediction has emerged as a core technology for understanding biomolecules and their interactions. Here, we combine homology-based structure prediction with molecular phylogenetic analysis to study the evolution of electrostatic membrane binding among vertebrate synaptotagmin-like proteins (Slps). Slp family proteins play key roles in the membrane trafficking of large dense-core secretory vesicles. Our previous experimental and computational study found that the C2A domain of Slp-4 (also called granuphilin) binds with high affinity to anionic phospholipids in the cytoplasmic leaflet of the plasma membrane through a large positively charged protein surface centered on a cluster of phosphoinositide-binding lysine residues. Because this surface contributes greatly to Slp-4 C2A domain membrane binding, we hypothesized that the net charge on the surface might be evolutionarily conserved. To test this hypothesis, the known C2A sequences of Slp-4 among vertebrates were organized by class (from mammalia to pisces) using molecular phylogenetic analysis. Consensus sequences for each class were then identified and used to generate homology structures, from which Poisson-Boltzmann electrostatic potentials were calculated. For comparison, homology structures and electrostatic potentials were also calculated for the five human Slp protein family members. The results demonstrate that the charge on the membrane-binding surface is highly conserved throughout the evolution of Slp-4, and more highly conserved than many individual residues among the human Slp family paralogs. Such molecular phylogenetic-driven computational analysis can help to describe the evolution of electrostatic interactions between proteins and membranes which are crucial for their function.

Ticks' attraction to electrically charged hosts

Ticks are blood-feeding parasites with limited locomotion, known for transmitting multiple pathogens to vertebrates. England et al. suggest that ticks can be easily pulled, via electrostatic induction, toward charged hosts with fluffy coats that are prone to accumulate higher electrostatic potentials. Thus, static electricity may influence ticks' ecology and management.

Dipolar Ligands Tune Plasmonic Properties of Tin-Doped Indium Oxide Nanocrystals

Surface functionalization with dipolar molecules is known to tune the electronic band alignment in semiconductor films and colloidal quantum dots. Yet, the influence of surface modification on plasmonic nanocrystals and their properties remains little explored. Here, we functionalize tin-doped indium oxide nanocrystals (ITO NCs) via ligand exchange with a series of cinnamic acids with different electron-withdrawing and -donating dipolar characters. Consistent with previous reports on semiconductors, we find that withdrawing (donating) ligands increase (decrease) the work function caused by an electrostatic potential shift across the molecular layer. Quantitative analyses of the plasmonic extinction spectra reveal that varying the ligand molecular dipole affects the near-surface depletion layer, with an anticorrelated trend between the electron concentration and electronic volume fraction, factors that are positively correlated in as-synthesized NCs. Electronic structure engineering through surface modification provides access to distinctive combinations of plasmonic properties that could enable optoelectronic applications, sensing, and hot electron-driven processes.

PKAD-2: New entries and expansion of functionalities of the database of experimentally measured pKa's of proteins

Almost all biological reactions are pH dependent and understanding the origin of pH dependence requires knowledge of the pKa's of ionizable groups. Here we report a new edition of PKAD, the PKAD-2, which is a database of experimentally measured pKa's of proteins, both wild type and mutant proteins. The new additions include 117 wild type and 54 mutant pKa values, resulting in total 1742 experimentally measured pKa's. The new edition of PKAD-2 includes 8 new wild type and 12 new mutant proteins, resulting in total of 220 proteins. This new edition incorporates a visual 3D image of the highlighted residue of interest within the corresponding protein or protein complex. Hydrogen bonds were identified, counted, and implemented as a search feature. Other new search features include the number of neighboring residues <4A from the heaviest atom of the side chain of a given amino acid. Here, we present PKAD-2 with the intention to continuously incorporate novel features and current data with the goal to be used as benchmark for computational methods.

Building polarization into protein-inhibitor binding dynamics in rational drug design for rheumatoid arthritis

Standard force field-based simulations to accomplish structure-based evaluations of lead molecules is a powerful tool. Combining protein fragmentation into tractable sub-systems with continuum solvation method is envisaged to enable quantum mechanics-based electronic structure calculations of macromolecules in their realistic environment. This along with incorporation of many-body polarization effect in molecular dynamics simulations may augment an accurate description of electrostatics of protein-inhibitor systems for effective drug design. Rheumatoid arthritis (RA) is a complex autoimmune disorder plagued by the ceiling effect of current targeted therapies, encouraging identification of new druggable targets and corresponding drug design to tackle the refractory form of disease. In this study, polarization-inclusive force field approach has been used to model protein solvation and ligand binding for 'Mitogen-activated protein kinase' (MAP3K8), a regulatory node of notable pharmacological relevance in RA synovial biology. For MAP3K8 inhibitors belonging to different scaffold series, the calculations illustrated differential electrostatic contribution to their relative binding affinities and successfully explained examples from available structure-activity relationship studies. Results from this study exemplified i) the advantage of this approach in reliably ranking inhibitors having close nanomolar range activities for the same target; and ii) its prospective application in lead molecule identification aiding drug discovery efforts in RA.Communicated by Ramaswamy H. Sarma.

Static electricity passively attracts ticks onto hosts

Most terrestrial animals naturally accumulate electrostatic charges, meaning that they will generate electric forces that interact with other charges in their environment, including those on or within other organisms. However, how this naturally occurring static electricity influences the ecology and life history of organisms remains largely unknown.¹ Mammals, birds, and reptiles are known to carry appreciable net electrostatic charges, equivalent to surface potentials on the order of hundreds to tens of thousands of volts.^{1,2,3,4,5,6,7} Therefore, we hypothesize that their parasites, such as ticks, are passively attracted onto their surfaces by electrostatic forces acting across air gaps. This biophysical mechanism is proposed by us to assist these ectoparasites in making contact with their hosts, increasing their effective "reach" because they are otherwise incapable of jumping. Herein, experimental and theoretical evidence show that the tick Ixodes ricinus (Figure 1A) can close the gap to their hosts using ecologically relevant electric fields. We also find that this electrostatic interaction is not significantly influenced by the polarity of the electric field, revealing that the mechanism of attraction relies upon induction of an electrical polarization within the tick, as opposed to a static charge on its surface. These findings open a new dimension to our understanding of how ticks, and possibly many other terrestrial organisms, find and attach to their hosts or vectors. Furthermore, this discovery may inspire novel solutions for mitigating the notable and often devastating economic, social, and public health impacts of ticks on humans and livestock.^{8,9,10,11,12,13,14,15}.

Discrete Element Method Evaluation of Triboelectric Charging Due to Powder Handling in the Capsule of a DPI

The generation and accumulation of an electrostatic charge from handling pharmaceutical powders is a well-known phenomenon, given the insulating nature of most APIs (Active Pharmaceutical Ingredients) and excipients. In capsule-based DPIs (Dry Powder Inhalers), the formulation is stored in a gelatine capsule placed in the inhaler just before inhalation. The action of capsule filling, as well as tumbling or vibration effects during the capsule life cycle, implies a consistent amount of particle-particle and particle-wall contacts. A significant contact-induced electrostatic charging can then take place, potentially affecting the inhaler's efficiency. DEM (Discrete Element Method) simulations were performed on a carrier-based DPI formulation (salbutamol-lactose) to evaluate such effects. After performing a comparison with the experimental data on a carrier-only system under similar conditions, a detailed analysis was conducted on two carrier-API configurations with different API loadings per carrier particle. The charge acquired by the two solid phases was tracked in both the initial particle settling and the capsule shaking process. Alternating positive-negative charging was observed. Particle charging was then investigated in relation to the collision statistics, tracking the particle-particle and particle-wall events for the carrier and API. Finally, an analysis of the relative importance of electrostatic, cohesive/adhesive, and inertial forces allowed the importance of each term in determining the trajectory of the powder particles to be estimated.

Analyzing the Geometry and Dynamics of Viral Structures: A Review of Computational Approaches Based on Alpha Shape Theory, Normal Mode Analysis, and Poisson-Boltzmann Theories

The current SARS-CoV-2 pandemic highlights our fragility when we are exposed to emergent viruses either directly or through zoonotic diseases. Fortunately, our knowledge of the biology of those viruses is improving. In particular, we have more and more structural information on virions, i.e., the infective form of a virus that includes its genomic material and surrounding protective capsid, and on their gene products. It is important to have methods that enable the analyses of structural information on such large macromolecular systems. We review some of those methods in this paper. We focus on understanding the geometry of virions and viral structural proteins, their dynamics, and their energetics, with the ambition that this understanding can help design antiviral agents. We discuss those methods in light of the specificities of those structures, mainly that they are huge. We focus on three of our own methods based on the alpha shape theory for computing geometry, normal mode analyses to study dynamics, and modified Poisson-Boltzmann theories to study the organization of ions and co-solvent and solvent molecules around biomacromolecules. The corresponding software has computing times that are compatible with the use of regular desktop computers. We show examples of their applications on some outer shells and structural proteins of the West Nile Virus.

Electrode potential driven dissociation of N-heterocycle-BF3 adducts: a possible manifestation of the electro-inductive effect

Recently, non-Faradaic effects were used to modify the electronic structure and reactivity of electrode-bound species. We hypothesize that these electrostatic perturbations could influence the chemical reactivity of electrolyte species near an electrode in the absence of Faradaic electron transfer. A prime example of non-Faradaic effects is acid-base dissociation near an interface. Here, we probed the near-electrode dissociation of N-heterocycle-BF3 Lewis adducts upon electrode polarization, well outside of the redox potential window of the adducts. Using scanning electrochemical microscopy and confocal fluorescence spectroscopy, we detected a potential-dependent depletion of the adduct near the electrode. We propose an electro-inductive effect where a more positive potential leads to electron withdrawal on the N-heterocycle. This study takes a step forward in the use of electrostatics at electrochemical interfaces for field-driven electrocatalytic and electrosynthetic processes.

How Do Quantum Effects Influence the Capacitance and Carrier Density of Monolayer MoS2 Transistors?

When transistor gate insulators have nanometer-scale equivalent oxide thickness (EOT), the gate capacitance (C_G) becomes smaller than the oxide capacitance (C_ox) due to the quantum capacitance and charge centroid capacitance of the channel. Here, we study the capacitance of monolayer MoS₂ as a prototypical two-dimensional (2D) channel while considering spatial variations in the potential, charge density, and density of states. At 0.5 nm EOT, the monolayer MoS₂ capacitance is smaller than its quantum capacitance, limiting the single-gated C_G of an n-type channel to between 63% and 78% of C_ox, for gate overdrive voltages between 0.5 and 1 V. Despite these limitations, for dual-gated devices, the on-state C_G of monolayer MoS₂ is 50% greater than that of silicon at 0.5 nm EOT and more than three times that of InGaAs at 1 nm EOT, indicating that such 2D semiconductors are promising for improved gate control of nanoscale transistors at future technology nodes.

Zwitterions fine-tune interactions in electrolyte solutions

Cellular organisms regulate electrolyte composition in the cytosol to optimize intracellular molecular interactions at the same time as balancing external osmotic pressure. While osmotic pressure can be tuned using multiple ionic, zwitterionic, and nonionic solutes, interactions between proteins and other macromolecules are sensitive to the precise composition of the medium. Nonetheless, the roles of individual ions and nonionic solutes in mediating cellular interactions remain relatively unexplored, and standard buffer solutions used in laboratory studies often contain only a few simple salts. Here, we report on model experiments investigating the combined effect of ionic and zwitterionic solutes on interaction forces across electrolytes, revealing a clear role for zwitterions in modifying interactions compared to simple salt solutions. First, we find that zwitterions act to disrupt water layering at interfaces, leading to smoothed interaction potentials. Second, we find that zwitterions strengthen electrostatic repulsions by enhancing effective surface charge. Third, zwitterions enhance the effective dielectric permittivity of the solution, and this "dielectricizer" effect extends the range of electrostatic repulsions compared to solutions without zwitterion present. The latter two effects are likely important in stabilizing proteins and other macromolecules when external osmotic and mechanical pressure are very high and simple ionic solutes alone would lead to collapse.

Substituent Effect and Its Halogen-Atom Dependence of Halogen Bonding Viewed through Electron Density Changes

Elucidating how the halogen-bonding ability and strength are controlled by the substituent effect and how this control depends on halogen atom will be essential for finely-tuned design of functionally important molecules. Here, this problem is tackled by analyzing the electron density differences/changes for variously substituted halobenzenes. It is shown that the anisotropy of the electron distribution around the halogen atom, which is an important factor for halogen-bonding ability, is not much affected by the substituent effect and rather simply depends on the halogen atom, while the partial charge on the halogen atom, which is related to the bond dipole of the C-X bond, is significantly modulated by the substituent effect and gives rise to enhancement of the electrostatic potential on the line extended from the C-X bond. The properties related to the polarization effect are also discussed.

Surface electrostatics dictate RNA-binding protein CAPRIN1 condensate concentration and hydrodynamic properties

Biomolecular condensates concentrate proteins, nucleic acids, and small molecules and play an essential role in many biological processes. Their formation is tuned by a balance between energetically favorable and unfavorable contacts, with charge-charge interactions playing a central role in some systems. The positively charged intrinsically disordered carboxy-terminal region of the RNA-binding protein CAPRIN1 is one such example, phase separating upon addition of negatively charged ATP or high concentrations of sodium chloride (NaCl). Using solution NMR spectroscopy, we measured residue-specific near-surface electrostatic potentials (ϕENS) of CAPRIN1 along its NaCl-induced phase separation trajectory to compare with those obtained using ATP. In both cases, electrostatic shielding decreases ϕENS values, yet surface potentials of CAPRIN1 in the two condensates can be different, depending on the amount of NaCl or ATP added. Our results establish that even small differences in ϕENS can significantly affect the level of protein enrichment and the mechanical properties of the condensed phase, leading, potentially, to the regulation of biological processes.

ES-Screen: A Novel Electrostatics-Driven Method for Drug Discovery Virtual Screening

Electrostatic interactions drive biomolecular interactions and associations. Computational modeling of electrostatics in biomolecular systems, such as protein-ligand, protein-protein, and protein-DNA, has provided atomistic insights into the binding process. In drug discovery, finding biologically plausible ligand-protein target interactions is challenging as current virtual screening and adjuvant techniques such as docking methods do not provide optimal treatment of electrostatic interactions. This study describes a novel electrostatics-driven virtual screening method called 'ES-Screen' that performs well across diverse protein target systems. ES-Screen provides a unique treatment of electrostatic interaction energies independent of total electrostatic free energy, typically employed by current software. Importantly, ES-Screen uses initial ligand pose input obtained from a receptor-based pharmacophore, thus independent of molecular docking. ES-Screen integrates individual polar and nonpolar replacement energies, which are the energy costs of replacing the cognate ligand for a target with a query ligand from the screening. This uniquely optimizes thermodynamic stability in electrostatic and nonpolar interactions relative to an experimentally determined stable binding state. ES-Screen also integrates chemometrics through shape and other physicochemical properties to prioritize query ligands with the greatest physicochemical similarities to the cognate ligand. The applicability of ES-Screen is demonstrated with in vitro experiments by identifying novel targets for many drugs. The present version includes a combination of many other descriptor components that, in a future version, will be purely based on electrostatics. Therefore, ES-Screen is a first-in-class unique electrostatics-driven virtual screening method with a unique implementation of replacement electrostatic interaction energies with broad applicability in drug discovery.

Swelling and Collapse of Cylindrical Polyelectrolyte Microgels

In this study, we propose computer simulations of charged cylindrical microgels. The effects of cross-linking density, aspect ratio, and fraction of charged groups on the microgel swelling and collapse with a variation in the solvent quality were studied. The results were compared with those obtained for equivalent neutral cylindrical microgels. The study demonstrated that microgels' degree of swelling strongly depends on the fraction of charged groups. Polyelectrolyte microgels under adequate solvent conditions are characterized by a larger length and thickness than their neutral analogues: the higher the fraction of charged groups, the longer their length and greater their thickness. Microgels' collapse upon solvent quality decline is characterized by a decrease in length and non-monotonous behavior of its thickness. First, the thickness decreases due to the attraction of monomer units (beads) upon collapse. The further thickness increase is related to the surface tension, which tends to reduce the anisotropy of collapsed objects (the minimum surface energy is known to be achieved for the spherical objects). This reduction is opposed by the network elasticity. The microgels with a low cross-linking density and/or a low enough aspect ratio reveal a cylinder-to-sphere collapse. Otherwise, the cylindrical shape is preserved in the course of the collapse. Aspect ratio as a function of the solvent quality (interaction parameter) demonstrates the maximum, which is solely due to the electrostatics. Finally, we plotted radial concentration profiles for network segments, their charged groups, and counterions.

Synthetic fertilizers alter floral biophysical cues and bumblebee foraging behavior

The use of agrochemicals is increasingly recognized as interfering with pollination services due to its detrimental effects on pollinators. Compared to the relatively well-studied chemical toxicity of agrochemicals, little is known on how they influence various biophysical floral cues that are used by pollinating insects to identify floral rewards. Here, we show that widely used horticultural and agricultural synthetic fertilizers affect bumblebee foraging behavior by altering a complex set of interlinked biophysical properties of the flower. We provide empirical and model-based evidence that synthetic fertilizers recurrently alter the magnitude and dynamics of floral electrical cues, and that similar responses can be observed with the neonicotinoid pesticide imidacloprid. We show that biophysical responses interact in modifying floral electric fields and that such changes reduce bumblebee foraging, reflecting a perturbation in the sensory events experienced by bees during flower visitation. This unveils a previously unappreciated anthropogenic interference elicited by agrochemicals within the electric landscape that is likely relevant for a wide range of chemicals and organisms that rely on naturally occurring electric fields.

Measuring Local Electrostatic Potentials Around Nucleic Acids by Paramagnetic NMR Spectroscopy

Electrostatic potentials around macromolecules in the presence of mobile charges are difficult to assess especially for highly charged systems. Here, we report measurements of local electrostatic potentials around DNA by paramagnetic NMR. Through quantitative analysis of NMR paramagnetic relaxation enhancement arising from positively charged or neutral paramagnetic cosolutes, we were able to determine local electrostatic potentials around 1H nuclei at >100 sites in major and minor grooves of 13C,15N-labeled 15-bp DNA at 100 mM NaCl. Our experimental electrostatic potential data directly confirmed the Coulombic end effects of DNA. The effective near-surface electrostatic potentials from the NMR data were in good agreement with the theoretical predictions with the Poisson-Boltzmann equation. This NMR method allows for unprecedented experimental investigations into the electrostatic properties of nucleic acids.

Protein-imprinted particles for coronavirus capture from solution

Molecular imprinting is a promising strategy to selectively adsorb viruses, but it requires discerning and validating epitopes that serve as effective imprinting templates. In this work, glycoprotein-imprinted particles were synthesized for coronavirus capture. Adsorption was maximized at pH 6 (the glycoprotein isoelectric point) where the glycoprotein-imprinted particles outperformed non-imprinted particles, adsorbing 4.96 × 10⁶  ± 3.33 × 10³ versus 3.54 × 10⁶  ± 1.39 × 10⁶ median tissue culture infectious dose/mg of the target coronavirus, human coronavirus - organ culture 43, within the first 30 min (p = 0.012). During competitive adsorption, with pH adjustment (pH 6), the glycoprotein-imprinted particles adsorbed more target virus than non-target coronavirus (human coronavirus - Netherland 63) with 2.34 versus 1.94 log removal in 90 min (p < 0.01). In contrast, the non-imprinted particles showed no significant difference in target versus non-target virus removal. Electrostatic potential calculation shows that the human coronavirus - organ culture 43 glycoprotein has positively charged pockets at pH 6, which may facilitate adsorption at lower pH values. Therefore, tuning the target virus glycoprotein charge via pH adjustment enhanced adsorption by minimizing repulsive electrostatic interactions with the particles. Overall, these results highlight the effective use of glycoprotein-imprinted particles for coronavirus capture and discern the merits and limitations of glycoprotein imprinting for the capture of enveloped viruses.

Electrostatics in Computational Biophysics and Its Implications for Disease Effects

This review outlines the role of electrostatics in computational molecular biophysics and its implication in altering wild-type characteristics of biological macromolecules, and thus the contribution of electrostatics to disease mechanisms. The work is not intended to review existing computational approaches or to propose further developments. Instead, it summarizes the outcomes of relevant studies and provides a generalized classification of major mechanisms that involve electrostatic effects in both wild-type and mutant biological macromolecules. It emphasizes the complex role of electrostatics in molecular biophysics, such that the long range of electrostatic interactions causes them to dominate all other forces at distances larger than several Angstroms, while at the same time, the alteration of short-range wild-type electrostatic pairwise interactions can have pronounced effects as well. Because of this dual nature of electrostatic interactions, being dominant at long-range and being very specific at short-range, their implications for wild-type structure and function are quite pronounced. Therefore, any disruption of the complex electrostatic network of interactions may abolish wild-type functionality and could be the dominant factor contributing to pathogenicity. However, we also outline that due to the plasticity of biological macromolecules, the effect of amino acid mutation may be reduced, and thus a charge deletion or insertion may not necessarily be deleterious.

RNA Electrostatics: How Ribozymes Engineer Active Sites to Enable Catalysis

Electrostatic interactions are fundamental to RNA structure and function, and intimately influenced by solvation and the ion atmosphere. RNA enzymes, or ribozymes, are catalytic RNAs that are able to enhance reaction rates over a million-fold, despite having only a limited repertoire of building blocks and available set of chemical functional groups. Ribozyme active sites usually occur at junctions where negatively charged helices come together, and in many cases leverage this strained electrostatic environment to recruit metal ions in solution that can assist in catalysis. Similar strategies have been implicated in related artificially engineered DNA enzymes. Herein, we apply Poisson-Boltzmann, 3D-RISM, and molecular simulations to study a set of metal-dependent small self-cleaving ribozymes (hammerhead, pistol, and Varkud satellite) as well as an artificially engineered DNAzyme (8-17) to examine electrostatic features and their relation to the recruitment of monovalent and divalent metal ions important for activity. We examine several fundamental roles for these ions that include: (1) structural integrity of the catalytically active state, (2) pKa tuning of residues involved in acid-base catalysis, and (3) direct electrostatic stabilization of the transition state via Lewis acid catalysis. Taken together, these examples demonstrate how RNA electrostatics orchestrates the site-specific and territorial binding of metal ions to play important roles in catalysis.

Gaussian Accelerated Molecular Dynamics in OpenMM

Gaussian accelerated molecular dynamics (GaMD) is a computational technique that provides both unconstrained enhanced sampling and free energy calculations of biomolecules. Here, we present the implementation of GaMD in the OpenMM simulation package and validate it on model systems of alanine dipeptide and RNA folding. For alanine dipeptide, 30 ns GaMD production simulations reproduced free energy profiles of 1000 ns conventional molecular dynamics (cMD) simulations. In addition, GaMD simulations captured the folding pathways of three hyperstable RNA tetraloops (UUCG, GCAA, and CUUG) and binding of the rbt203 ligand to the HIV-1 Tar RNA, both of which involved critical electrostatic interactions such as hydrogen bonding and base stacking. Together with previous implementations, GaMD in OpenMM will allow for wider applications in simulations of proteins, RNA, and other biomolecules.

Electronic Polarization at the Interface between the p53 Transactivation Domain and Two Binding Partners

Intrinsically disordered proteins (IDPs) are an abundant class of highly charged proteins that participate in numerous crucial biological processes, often in regulatory roles. IDPs do not have one major free energy minimum with a dominant structure, instead existing as conformational ensembles of multiple semistable conformations. p53 is a prototypical protein with disordered regions and binds to many structurally diverse partners, making it a useful model for exploring the role of electrostatic interactions at IDP binding interfaces. In this study, we used the Drude-2019 force field to simulate the p53 transactivation domain with two protein partners to probe the role of electrostatic interactions in IDP protein-protein interactions. We found that the Drude-2019 polarizable force field reasonably reproduced experimental chemical shifts of the p53 transactivation domain (TAD) in one complex for which these data are available. We also found that the proteins in these complexes displayed dipole response at specific residues of each protein and that residues primarily involved in binding showed a large percent change in dipole moment between the unbound and complexed states. Probing the role of electrostatic interactions in IDP binding can allow us greater fundamental understanding of these interactions and may help with targeting p53 or its partners for drug design.

Electrostatic Interaction with the Bacterial Cell Envelope Tunes the Lytic Activity of Two Novel Peptidoglycan Hydrolases

Peptidoglycan (PG) hydrolases, due to their crucial role in the metabolism of the bacterial cell wall (CW), are increasingly being considered suitable targets for therapies, and a potent alternative to conventional antibiotics. In the light of contradictory data reported, detailed mechanism of regulation of enzymes activity based on electrostatic interactions between hydrolase molecule and bacterial CW surface remains unknown. Here, we report a comprehensive study on this phenomenon using as a model two novel PG hydrolases, SpM23_A, and SpM23_B, which although share the same bacterial host, similarities in sequence conservation, domain architecture, and structure, display surprisingly distinct net charges (in 2D electrophoresis, pI 6.8, and pI 9.7, respectively). We demonstrate a strong correlation between hydrolases surface net charge and the enzymes activity by modulating the charge of both, enzyme molecule and bacterial cell surface. Teichoic acids, anionic polymers present in the bacterial CW, are shown to be involved in the mechanism of enzymes activity regulation by the electrostatics-based interplay between charged bacterial envelope and PG hydrolases. These data serve as a hint for the future development of chimeric PG hydrolases of desired antimicrobial specificity. IMPORTANCE This study shows direct relationship between the surface charge of two recently described enzymes, SpM23_A and SpM23_B, and bacterial cell walls. We demonstrate that by (i) surface charge probing of bacterial strains collection, (ii) reduction of the net charge of the positively charged enzyme, and (iii) altering the net charge of the bacterial surface by modifying the content and composition of teichoic acids. In all cases, we observed that lytic activity and binding strength of SpM23 enzymes, are regulated by electrostatic interactions with the bacterial cell envelope and that this interaction contributes to the determination of the spectrum of susceptible bacterial species. Moreover, we revealed the regulatory role of charged cell wall components, namely, teichoic and lipoteichoic acids, over the SpM23 enzymes. We believe that our findings make an important contribution to understand the means of hydrolases activity regulation in the complex environment of the bacterial cell wall.

TUPÃ: Electric field analyses for molecular simulations

We introduce TUPÃ, a Python-based algorithm to calculate and analyze electric fields in molecular simulations. To demonstrate the features in TUPÃ, we present three test cases in which the orientation and magnitude of the electric field exerted by biomolecules help explain biological phenomena or observed kinetics. As part of TUPÃ, we also provide a PyMOL plugin to help researchers visualize how electric fields are organized within the simulation system. The code is freely available and can be obtained at https://mdpoleto.github.io/tupa/.

Negatively Charged Disordered Regions are Prevalent and Functionally Important Across Proteomes

Intrinsically disordered regions (IDRs) of proteins are often characterized by a high fraction of charged residues, but differ in their overall net charge and in the organization of the charged residues. The function-encoding information stored via IDR charge composition and organization remains elusive. Here, we aim to decipher the sequence-function relationship in IDRs by presenting a comprehensive bioinformatic analysis of the charge properties of IDRs in the human, mouse, and yeast proteomes. About 50% of the proteins comprise at least a single IDR, which is either positively or negatively charged. Highly negatively charged IDRs are longer and possess greater net charge per residue compared with highly positively charged IDRs. A striking difference between positively and negatively charged IDRs is the characteristics of the repeated units, specifically, of consecutive Lys or Arg residues (K/R repeats) and Asp or Glu (D/E repeats) residues. D/E repeats are found to be about five times longer than K/R repeats, with the longest found containing 49 residues. Long stretches of consecutive D and E are found to be more prevalent in nucleic acid-related proteins. They are less common in prokaryotes, and in eukaryotes their abundance increases with genome size. The functional role of D/E repeats and the profound differences between them and K/R repeats are discussed.

Electrostatic Contribution to the Photo-Assisted Piezoresponse Force Microscopy by Photo-Induced Surface Charge

The surging interest in manipulating the polarization of piezo/ferroelectric materials by means of light has driven an increasing number of studies toward their light-polarization interaction. One way to investigate such interaction is by performing piezoresponse force microscopy (PFM) while/after the sample is exposed to light illumination. However, caution must be exercised when analyzing and interpreting the data, as demonstrated in this paper, because sizeable photo-response observed in the PFM amplitude image of the sample is shown to be caused by the electrostatic interaction between the photo-induced surface charge and tip. Through photo-assisted Kelvin probe force microscopy (KPFM), positive surface potential is found to be developed near the sample's surface under 405 nm light illumination, whose effects on the measured PFM signal is revealed by the comparative studies on its amplitude curves that are obtained using PFM spectroscopy mode with/without illumination. This work exemplifies the need for complementary use of KPFM, PFM imaging mode, and PFM spectroscopy mode in order to distinguish real behavior from artifacts.

Bicyclo-DNA mimics with enhanced protein binding affinities: insights from molecular dynamics simulations

DNA-protein interactions occur at all levels of DNA expression and replication and are crucial determinants for the survival of a cell. Several modified nucleotides have been utilized to manipulate these interactions and have implications in drug discovery. In the present article, we evaluated the binding of bicyclo-nucleotides (generated by forming a methylene bridge between C1' and C5' in sugar, leading to a bicyclo system with C2' axis of symmetry at the nucleotide level) to proteins. We utilized four ssDNA-protein complexes with experimentally known binding free energies and investigated the binding of modified nucleotides to proteins via all-atom explicit solvent molecular dynamics (MD) simulations (200 ns), and compared the binding with control ssDNA-protein systems. The modified ssDNA displayed enhanced binding to proteins as compared to the control ssDNA, as seen by means of MD simulations followed by MM-PBSA calculations. Further, the Delphi-based electrostatic estimation revealed that the high binding of modified ssDNA to protein might be related to the enhanced electrostatic complementarity displayed by the modified ssDNA molecules in all the four systems considered for the study. The improved binding achieved with modified nucleotides can be utilized to design and develop anticancer/antisense molecules capable of targeting proteins or ssRNAs.Communicated by Ramaswamy H. Sarma.

Towards optimal boundary integral formulations of the Poisson-Boltzmann equation for molecular electrostatics

The Poisson-Boltzmann equation offers an efficient way to study electrostatics in molecular settings. Its numerical solution with the boundary element method is widely used, as the complicated molecular surface is accurately represented by the mesh, and the point charges are accounted for explicitly. In fact, there are several well-known boundary integral formulations available in the literature. This work presents a generalized expression of the boundary integral representation of the implicit solvent model, giving rise to new forms to compute the electrostatic potential. Moreover, it proposes a strategy to build efficient preconditioners for any of the resulting systems, improving the convergence of the linear solver. We perform systematic benchmarking of a set of formulations and preconditioners, focusing on the time to solution, matrix conditioning, and eigenvalue spectrum. We see that the eigenvalue clustering is a good indicator of the matrix conditioning, and show that they can be easily manipulated by scaling the preconditioner. Our results suggest that the optimal choice is problem-size dependent, where a simpler direct formulation is the fastest for small molecules, but more involved second-kind equations are better for larger problems. We also present a fast Calderón preconditioner for first-kind formulations, which shows promising behavior for future analysis. This work sets the basis towards choosing the most convenient boundary integral formulation of the Poisson-Boltzmann equation for a given problem.

Electrostatic influence on IL-1 transport through the GSDMD pore

A variety of signals, including inflammasome activation, trigger the formation of large transmembrane pores by gasdermin D (GSDMD). There are primarily two functions of the GSDMD pore, to drive lytic cell death, known as pyroptosis, and to permit the release of leaderless interleukin-1 (IL-1) family cytokines, a process that does not require pyroptosis. We are interested in the mechanism by which the GSDMD pore channels IL-1 release from living cells. Recent studies revealed that electrostatic interaction, in addition to cargo size, plays a critical role in GSDMD-dependent protein release. Here, we determined computationally that to enable electrostatic filtering against pro-IL-1β, acidic lipids in the membrane need to effectively neutralize positive charges in the membrane-facing patches of the GSDMD pore. In addition, we predicted that salt has an attenuating effect on electrostatic filtering and then validated this prediction using a liposome leakage assay. A calibrated electrostatic screening factor is necessary to account for the experimental observations, suggesting that ion distribution within the pore may be different from the bulk solution. Our findings corroborate the electrostatic influence of IL-1 transport exerted by the GSDMD pore and reveal extrinsic factors, including lipid and salt, that affect the electrostatic environment.

Functionally Additive Fixed Positive and Negative Charges in the CFTR Channel Pore Control Anion Binding and Conductance

Ion channels use charged amino-acid residues to attract oppositely charged permeant ions into the channel pore. In the cystic fibrosis transmembrane conductance regulator (CFTR) Cl⁻ channel, a number of arginine and lysine residues have been shown to be important for Cl⁻ permeation. Among these, two in close proximity in the pore—Lys⁹⁵ and Arg¹³⁴—are indispensable for anion binding and high Cl⁻ conductance, suggesting that high positive charge density is required for pore function. Here we used mutagenesis and functional characterization to show that a nearby pore-lining negatively charged residue (Glu⁹²) plays a functionally additive role with these two positive charges. While neutralization of this negative charge had little effect on anion binding or Cl⁻ conductance, such neutralization was able to reverse the detrimental effects of removing the positive charge at either Lys⁹⁵ or Arg¹³⁴, as well as the similar effects of introducing a negative charge at a neighboring residue (Ser¹¹⁴¹). Furthermore, neutralization of Glu⁹² greatly increased the susceptibility of the channel to blockage by divalent S₂O₃²⁻ anions, mimicking the effect of introducing additional positive charge in this region; this effect was reversed by concurrent neutralization of either Lys⁹⁵ or Arg¹³⁴. Across a panel of mutant channels that introduced or removed fixed charges at these four positions, we found that many pore properties are dependent on the overall charge or charge density. We propose that the CFTR pore uses a combination of positively and negatively charged residues to optimize the anion binding and Cl⁻ conductance properties of the channel.

Electrostatic Discovery Atomic Force Microscopy

While offering high resolution atomic and electronic structure, scanning probe microscopy techniques have found greater challenges in providing reliable electrostatic characterization on the same scale. In this work, we offer electrostatic discovery atomic force microscopy, a machine learning based method which provides immediate maps of the electrostatic potential directly from atomic force microscopy images with functionalized tips. We apply this to characterize the electrostatic properties of a variety of molecular systems and compare directly to reference simulations, demonstrating good agreement. This approach offers reliable atomic scale electrostatic maps on any system with minimal computational overhead.

Electrostatic pull-in application in flexible devices: A review

The electrostatic pull-in effect is a common phenomenon and a key parameter in the design of microscale and nanoscale devices. Flexible electronic devices based on the pull-in effect have attracted increasing attention due to their unique ductility. This review summarizes nanoelectromechanical switches made by flexible materials and classifies and discusses their applications in, among others, radio frequency systems, microfluidic systems, and electrostatic discharge protection. It is supposed to give researchers a more comprehensive understanding of the pull-in phenomenon and the development of its applications. Also, the review is meant to provide a reference for engineers to design and optimize devices.

Enhanced attraction between drops carrying fluctuating charge distributions

The electrostatic force between conductive spheres is always attractive at small separations, irrespective of their mean charge, when the charge on the spheres is constant. In many situations, the charge may not be fixed, such as for water drops in the natural atmosphere which vary in size and charge. We show that the attractive force between charged conductive spheres increases with increasing charge variance. The importance of this unrecognized electrostatic effect between water drops is evaluated for its potential to enhance rain formation.

Free Energy Analyses of Cell-Penetrating Peptides Using the Weighted Ensemble Method

Cell-penetrating peptides (CPPs) have been widely used for drug-delivery agents; however, it has not been fully understood how they translocate across cell membranes. The Weighted Ensemble (WE) method, one of the most powerful and flexible path sampling techniques, can be helpful to reveal translocation paths and free energy barriers along those paths. Within the WE approach we show how Arg9 (nona-arginine) and Tat interact with a DOPC/DOPG(4:1) model membrane, and we present free energy (or potential mean of forces, PMFs) profiles of penetration, although a translocation across the membrane has not been observed in the current simulations. Two different compositions of lipid molecules were also tried and compared. Our approach can be applied to any CPPs interacting with various model membranes, and it will provide useful information regarding the transport mechanisms of CPPs.

The effect of charged residue substitutions on the thermodynamics of protein-surface interactions

The interactions of proteins with surfaces are important in both biological processes and biotechnologies. In contrast to decades of study regarding the biophysics of proteins in bulk solution, however, our mechanistic understanding of the biophysics of proteins interacting with surfaces remains largely qualitative. In response, we have set to explore quantitatively the thermodynamics of protein-surface interactions. In this work, we explore systematically the role of electrostatics in modulating the interaction between proteins and charged surfaces. In particular, we use electrochemistry to explore the extent to which a macroscopic, hydroxyl-coated surface held at a slightly negative potential affects the folding thermodynamics of surface-attached protein variants with different composition of charged amino acids. Doing so, we find that attachment to the surface generally leads to a net stabilization, presumably due to excluded volume effects that reduce the entropy of the unfolded state. The magnitude of this stabilization, however, is strongly correlated with the charged-residue content of the protein. In particular, we find statistically significant correlations with both the net charge of the protein, with greater negative charge leading to less stabilization by the surface, and with the number of arginines, with more arginines leading to greater stabilization. Such findings refine our understanding of protein-surface interactions, providing in turn a guiding rationale to achieve the functional deposition of proteins on artificial surfaces for implementation in, for example, protein-based biotechnologies.

The ecology of electricity and electroreception

Electricity, the interaction between electrically charged objects, is widely known to be fundamental to the functioning of living systems. However, this appreciation has largely been restricted to the scale of atoms, molecules, and cells. By contrast, the role of electricity at the ecological scale has historically been largely neglected, characterised by punctuated islands of research infrequently connected to one another. Recently, however, an understanding of the ubiquity of electrical forces within the natural environment has begun to grow, along with a realisation of the multitude of ecological interactions that these forces may influence. Herein, we provide the first comprehensive collation and synthesis of research in this emerging field of electric ecology. This includes assessments of the role electricity plays in the natural ecology of predator-prey interactions, pollination, and animal dispersal, among many others, as well as the impact of anthropogenic activity on these systems. A detailed introduction to the ecology and physiology of electroreception - the biological detection of ecologically relevant electric fields - is also provided. Further to this, we suggest avenues for future research that show particular promise, most notably those investigating the recently discovered sense of aerial electroreception.

A rush to explore protein-ligand electrostatic interaction energy with Charger

The mutual penetration of electron densities between two interacting molecules complicates the computation of an accurate electrostatic interaction energy based on a pseudo-atom representation of electron densities. The numerical exact potential and multipole moment (nEP/MM) method is time-consuming since it performs a 3D integration to obtain the electrostatic energy at short interaction distances. Nguyen et al. [(2018), Acta Cryst. A74, 524-536] recently reported a fully analytical computation of the electrostatic interaction energy (aEP/MM). This method performs much faster than nEP/MM (up to two orders of magnitude) and remains highly accurate. A new program library, Charger, contains an implementation of the aEP/MM method. Charger has been incorporated into the MoProViewer software. Benchmark tests on a series of small molecules containing only C, H, N and O atoms show the efficiency of Charger in terms of execution time and accuracy. Charger is also powerful in a study of electrostatic symbiosis between a protein and a ligand. It determines reliable protein-ligand interaction energies even when both contain S atoms. It easily estimates the individual contribution of every residue to the total protein-ligand electrostatic binding energy. Glutathione transferase (GST) in complex with a benzophenone ligand was studied due to the availability of both structural and thermodynamic data. The resulting analysis highlights not only the residues that stabilize the ligand but also those that hinder ligand binding from an electrostatic point of view. This offers new perspectives in the search for mutations to improve the interaction between the two partners. A proposed mutation would improve ligand binding to GST by removing an electrostatic obstacle, rather than by the traditional increase in the number of favourable contacts.

Highly Basic Clusters in the Herpes Simplex Virus 1 Nuclear Egress Complex Drive Membrane Budding by Inducing Lipid Ordering

During replication of herpesviruses, capsids escape from the nucleus into the cytoplasm by budding at the inner nuclear membrane. This unusual process is mediated by the viral nuclear egress complex (NEC) that deforms the membrane around the capsid by oligomerizing into a hexagonal, membrane-bound scaffold. Here, we found that highly basic membrane-proximal regions (MPRs) of the NEC alter lipid order by inserting into the lipid headgroups and promote negative Gaussian curvature. We also find that the electrostatic interactions between the MPRs and the membranes are essential for membrane deformation. One of the MPRs is phosphorylated by a viral kinase during infection, and the corresponding phosphomimicking mutations block capsid nuclear egress. We show that the same phosphomimicking mutations disrupt the NEC-membrane interactions and inhibit NEC-mediated budding in vitro, providing a biophysical explanation for the in vivo phenomenon. Our data suggest that the NEC generates negative membrane curvature by both lipid ordering and protein scaffolding and that phosphorylation acts as an off switch that inhibits the membrane-budding activity of the NEC to prevent capsid-less budding. IMPORTANCE Herpesviruses are large viruses that infect nearly all vertebrates and some invertebrates and cause lifelong infections in most of the world's population. During replication, herpesviruses export their capsids from the nucleus into the cytoplasm by an unusual mechanism in which the viral nuclear egress complex (NEC) deforms the nuclear membrane around the capsid. However, how membrane deformation is achieved is unclear. Here, we show that the NEC from herpes simplex virus 1, a prototypical herpesvirus, uses clusters of positive charges to bind membranes and order membrane lipids. Reducing the positive charge or introducing negative charges weakens the membrane deforming ability of the NEC. We propose that the virus employs electrostatics to deform nuclear membrane around the capsid and can control this process by changing the NEC charge through phosphorylation. Blocking NEC-membrane interactions could be exploited as a therapeutic strategy.

Comparison of the MSMS and NanoShaper molecular surface triangulation codes in the TABI Poisson-Boltzmann solver

The Poisson-Boltzmann (PB) implicit solvent model is a popular framework for studying the electrostatics of solvated biomolecules. In this model the dielectric interface between the biomolecule and solvent is often taken to be the molecular surface or solvent-excluded surface (SES), and the quality of the SES triangulation is critical in boundary element simulations of the model. This work compares the performance of the MSMS and NanoShaper surface triangulation codes for a set of 38 biomolecules. While MSMS produces triangles of exceedingly small area and large aspect ratio, the two codes yield comparable values for the SES surface area and electrostatic solvation energy, where the latter calculations were performed using the treecode-accelerated boundary integral (TABI) PB solver. However we found that NanoShaper is computationally more efficient and reliable than MSMS, especially when parameters are set to produce highly resolved triangulations.