pH Dependence of Charge Multipole Moments in Proteins

Correlation between positron annihilation lifetime and photoluminescence measurements for calcined Hydroxyapatite

Hydroxyapatite (HAp) Ca10(PO4)6(OH)2 is a compound that has stable chemical properties, composition, and an affinity for human bone. As a result, it can be used in odontology, cancer treatment, and orthopedic grafts to repair damaged bone. To produce calcined HAp at 600 °C with different pH values, a wet chemical precipitation method was employed. All synthesized HAp samples were characterized by X-ray diffraction (XRD), Raman spectroscopy, scanning electron microscopy (SEM), transmission electron microscopy (TEM), Fourier-transform infrared spectroscopy (FTIR), photoluminescence (PL), Zeta potential, and positron annihilation lifetime spectroscopy (PALS). The XRD results revealed that all calcined HAp samples were formed in a hexagonal structure with a preferred (002) orientation at different pH values. The crystal size of the samples was determined using the Scherrer equation, which ranged from 16 to 25 nm. The SEM and TEM results showed that the morphology of the samples varied from nanorods to nanospheres and rice-like structures depending on the pH value of the sample. The PL measurements indicated that the blue and green emission peaks of HAp were due to defects (bulk, surface, and interface) in the samples, which created additional energy levels within the band gap. According to Zeta potential measurements, the charge carrier changed from a positive to negative value, ranging from 3.94 mV to - 2.95 mV. PALS was used to understand the relationship between the defects and the photoluminescence (PL) properties of HAp. Our results suggest that HAp nanoparticles have excellent potential for developing non-toxic biomedical and optical devices for phototherapy.

Influence of different pH milieu on the structure and function of human Aurora kinase B protein (AURK-B): Amalgamation of both spectroscopic and computational approach

Aurora kinase B (AURK-B) is a serine/threonine kinase protein that plays an essential role in chromosomal separation during the cell cycle event. AURK-B is highly expressed in various types of cancer such as human seminoma, thyroid carcinoma, non-small cell lung carcinoma (NSCLC), oral carcinoma, and gastric cancer. Hence, it is a potential therapeutic target in the treatment of various cancers. The structure of AURK-B in complex with one of its substrate inner centromeric protein (INCENP) is present, but the structural and functional characterization of native AURK-B at different pH environment is still unexplored.This study determines the effect of different pH milieu on the structure and function of AURK-B protein wherein the influence of pH on the protein conformation was probed using Circular dichroism (CD) and fluorescence spectroscopy. The structural studies were further combined with functional activity assay to observe the change in kinase activity at various pH milieu (2.0-11.0). CD and fluorescence spectroscopy experiments dictate that at high acidic conditions (pH 2.0 - 5.0), the secondary and tertiary structures of AURK-B become distorted, leading to diminished activity. The protein, however, was observed to stabilize towards pH 7.0 - 8.0 with minimal structure alteration over the basic pH range (pH 9.0 -11.0). The measured spectroscopic structural features were found to be in-line with obtained experimental kinase activity assays. Further, in-vitro experiments indicate that the enzyme is maximally active at pH 8.0. More ordered conformation and compact structure was observed at this pH (pH 8.0) as compared to other pH values through molecular dynamics simulation studies (MDS). As AURK-B localizes itself in the intracellular compartment, this study may provide a clue about the role of different pH environments in enhancing cancer growth, proliferation, and invasion.

Changes in total charge on spike protein of SARS-CoV-2 in emerging lineages

Motivation

Charged amino acid residues on the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been shown to influence its binding to different cell surface receptors, its non-specific electrostatic interactions with the environment, and its structural stability and conformation. It is therefore important to obtain a good understanding of amino acid mutations that affect the total charge on the spike protein which have arisen across different SARS-CoV-2 lineages during the course of the virus' evolution.

Results

We analyse the change in the number of ionizable amino acids and the corresponding total charge on the spike proteins of almost 2200 SARS-CoV-2 lineages that have emerged over the span of the pandemic. Our results show that the previously observed trend toward an increase in the positive charge on the spike protein of SARS-CoV-2 variants of concern has essentially stopped with the emergence of the early omicron variants. Furthermore, recently emerged lineages show a greater diversity in terms of their composition of ionizable amino acids. We also demonstrate that the patterns of change in the number of ionizable amino acids on the spike protein are characteristic of related lineages within the broader clade division of the SARS-CoV-2 phylogenetic tree. Due to the ubiquity of electrostatic interactions in the biological environment, our findings are relevant for a broad range of studies dealing with the structural stability of SARS-CoV-2 and its interactions with the environment.

Availability and implementation

The data underlying the article are available in the Supplementary material.

Interaction of Polyanionic and Polycationic Brushes with Globular Proteins and Protein-like Nanocolloids

A large number of experimental studies have demonstrated that globular proteins can be absorbed from the solution by both polycationic and polyanionic brushes when the net charge of protein globules is of the same or of the opposite sign with respect to that of brush-forming polyelectrolyte chains. Here, we overview the results of experimental studies on interactions between globular proteins and polycationic or polyanionic brushes, and present a self-consistent field theoretical model that allows us to account for the asymmetry of interactions of protein-like nanocolloid particles comprising weak (pH-sensitive) cationic and anionic groups with a positively or negatively charged polyelectrolyte brush. The position-dependent insertion free energy and the net charge of the particle are calculated. The theoretical model predicts that if the numbers of cationic and anionic ionizable groups of the protein are approximately equal, then the interaction patterns for both cationic and anionic brushes at equal offset on the "wrong side" from the isoelectric point (IEP), i.e., when the particle and the brush charge are of the same sign, are similar. An essential asymmetry in interactions of particles with polycationic and polyanionic brushes is predicted when fractions of cationic and anionic groups differ significantly. That is, at a pH above IEP, the anionic brush better absorbs negatively charged particles with a larger fraction of ionizable cationic groups and vice versa.

Conformational and Structural Changes in Chickpea Proteins Caused by Simulated Salivary Alterations in the Elderly

The impact of salivary alterations on chickpea protein structure in the elderly has not been well documented. This study aimed to understand the role of simulated salivary alterations in the conformational properties and secondary structure of the chickpea protein isolate (CPI). Whey protein isolate (WPI) was used as the reference. Protein dispersions (10%) were subjected to in vitro oral processing under simulated salivary conditions in both the elderly and adult subjects. Proteins and their oral counterparts were characterized in terms of their composition, charge, size, solubility, water absorption, molecular weight (MW), and secondary structure (Circular Dichroism and Raman spectroscopy). Under condition of simulated oral digestion in the elderly population, the ordered secondary protein structure was significantly affected, decreasing α-helix by ~36% and ~29% in CPI and WPI compared to the control (adult) population, respectively. An increase in the unordered random coil state was observed. These results could be attributed to an increase in electrolytes in the salivary composition. The structure of CPI is more stable than that of WPI because of its higher MW, more rigid structure, less charged surface, and different amino acid compositions. This study is meaningful in understanding how alterations in the elderly oral system affect protein conformation and is expected to improve the understanding of plant-based protein digestibility.

Integration of Cypoviruses into polyhedrin matrix

Unlike in other viruses, in Cypoviruses the genome is doubly protected since their icosahedral capsids are embedded into a perfect polyhedrin crystal. Current experimental methods cannot resolve the resulting interface structure and we propose a symmetry-based approach to predict it. We reveal a remarkable match between the surfaces of Cypovirus and the outer polyhedrin matrix. The match arises due to the preservation of the common tetragonal symmetry, allowing perfect contacts of polyhedrin trimers with VP1 and VP5 capsid proteins. We highlight a crucial role of the VP5 proteins in embedding the Cypovirus into the polyhedrin matrix and discuss the relationship between the nucleoside triphosphatase activity of the proteins and their role in the superstructure formation. Additionally, we propose an electrostatic mechanism that drives the viral superstructure disassembly occurring in the alkaline environment of the insect intestines. Our study may underpin novel strategies for engineering proteinaceous nanocontainers in diverse biotechnological and chemical applications.

Langmuir monolayer of lysozyme at variable subphase pH conditions: a comprehensive study on structure, morphology and hysteresis behaviour

Formation of a pure Langmuir monolayer of lysozyme at the air-water interface and its investigation by means of a surface pressure (π)-mean molecular area (A) isotherm has been accomplished under different subphase pH conditions. A normalized area-time curve confirms the stable nature of the lysozyme monolayer whose compressibility variation with an increased surface pressure at specific subphase pH has also been studied from π-A isotherms. The monolayers exhibit irreversible hysteresis behaviour irrespective of subphase pH conditions, as evidenced from successive compression-expansion π-A isotherm cycles. Comparison of surface thermodynamics under hysteresis with subphase pH variation confirms that the monolayer at subphase pH ≈ 4.0 involves a greater amount of energy to attain and retain the ordered and compact monolayer than the other two pH conditions (pH ≈ 7.0 and 9.5). In situ visualization of lysozyme monolayers by Brewster angle microscopy suggests the homogeneous and stripe-like pattern formation at lower and higher surface pressure respectively. Further investigations of lysozyme films at solid surfaces have been carried out with atomic force microscopy and X-ray reflectivity (XRR) analysis. Structural reversibility of lysozyme molecules under compression-expansion-compression of the monolayer is revealed from the comparison of height profiles of AFM images and electron density profiles as extracted from XRR analysis of the films deposited during both first and second compressions of the monolayer. The mechanism of the structural rearrangement of lysozyme molecules with surface pressure variation at different subphase pH is explored, correlating macroscopic and microscopic information.

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

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

Unusual Aspects of Charge Regulation in Flexible Weak Polyelectrolytes

This article reviews the state of the art of the studies on charge regulation (CR) effects in flexible weak polyelectrolytes (FWPE). The characteristic of FWPE is the strong coupling of ionization and conformational degrees of freedom. After introducing the necessary fundamental concepts, some unconventional aspects of the the physical chemistry of FWPE are discussed. These aspects are: (i) the extension of statistical mechanics techniques to include ionization equilibria and, in particular, the use of the recently proposed Site Binding-Rotational Isomeric State (SBRIS) model, which allows the calculation of ionization and conformational properties on the same foot; (ii) the recent progresses in the inclusion of proton equilibria in computer simulations; (iii) the possibility of mechanically induced CR in the stretching of FWPE; (iv) the non-trivial adsorption of FWPE on ionized surfaces with the same charge sign as the PE (the so-called "wrong side" of the isoelectric point); (v) the influence of macromolecular crowding on CR.

Unfolding and Aggregation Pathways of Variable Domains from Immunoglobulin Light Chains

Light chain amyloidosis is the most common form of systemic amyloidosis. This disease is caused by the formation and deposition of amyloid fibers made from immunoglobulin light chains. Environmental conditions such as pH and temperature can affect protein structure and induce the development of these fibers. Several studies have shed light on the native state, stability, dynamics, and final amyloid state of these proteins; however, the initiation process and the fibril formation pathway remain poorly understood structurally and kinetically. To study this, we analyzed the unfolding and aggregation process of the 6aJL2 protein under acidic conditions, with temperature changes, and upon mutation, using biophysical and computational techniques. Our results suggest that the differences in amyloidogenicity displayed by 6aJL2 under these conditions are caused by traversing different aggregation pathways, including unfolded intermediates and the formation of oligomers.

The significance of bioelectricity on all levels of organization of an organism. Part 1: From the subcellular level to cells

Bioelectricity plays an essential role in the structural and functional organization of biological organisms. In this first article of our three-part series, we summarize the importance of bioelectricity for the basic structural level of biological organization, i.e. from the subcellular level (charges, ion channels, molecules and cell organelles) to cells.

Packing and trimer-to-dimer protein reconstruction in icosahedral viral shells with a single type of symmetrical structural unit

Understanding the principles of protein packing and the mechanisms driving morphological transformations in virus shells (capsids) during their maturation can be pivotal for the development of new antiviral strategies. Here, we study how these principles and mechanisms manifest themselves in icosahedral viral capsids assembled from identical symmetric structural units (capsomeres). To rationalize such shells, we model capsomers as symmetrical groups of identical particles interacting with a short-range potential typical of the classic Tammes problem. The capsomere particles are assumed to retain their relative positions on the vertices of planar polygons placed on the spherical shell and to interact only with the particles from other capsomeres. Minimization of the interaction energy enforces equal distances between the nearest particles belonging to neighboring capsomeres and minimizes the number of different local environments. Thus, our model implements the Caspar and Klug quasi-equivalence principle and leads to packings strikingly similar to real capsids. We then study a reconstruction of protein trimers into dimers in a Flavivirus shell during its maturation, connecting the relevant structural changes with the modifications of the electrostatic charges of proteins, wrought by the oxidative switch in the bathing solution that is essential for the process. We highlight the key role of pr peptides in the shell reconstruction and show that the highly ordered arrangement of these subunits in the dimeric state is energetically favored at a low pH level. We also discuss the electrostatic mechanisms controlling the release of pr peptides in the last irreversible step of the maturation process.

Negatively Charged Composite Nanofibrous Hydrogel Membranes for High-Performance Protein Adsorption

Nanofibrous materials are considered as promising candidates for fabricating high-efficiency chromatography media, which are urgently needed in protein pharmaceuticals purification and biological research, yet still face several bottlenecks. Herein, novel negatively charged composite nanofibrous hydrogel membranes (NHMs) are obtained by a facile combination of electrospinning and surface coating modification. The resulting NHMs exhibit controllable morphologies and chemical structures. Benefitting from the combined effect of the stable framework of silicon dioxide (SiO₂) nanofiber and the function layer of negatively charged hydrogel, as well as good pore connectivity among nanofibers, NHMs exhibit a high protein adsorption capacity of around 1000 mg g^-1, and are superior to the commercial cellulose fibrous adsorbent (Sartobind^®) and the reported nanofibrous membranous adsorbents. Moreover, due to their relatively stable physicochemical and mechanical properties, NHMs possess comprehensive adsorption performance, favorable resistance to acid and solvents, good selectivity, and excellent regenerability. The designed NHMs composite adsorbents are expected to supply a new protein chromatography platform for effective protein purification in biopharmaceuticals and biochemical reagents.

Adsorption of flexible proteins in the 'wrong side' of the isoelectric point: Casein macropeptide as a model system

We analyze the conditions of the adsorption of a flexible peptide onto a charged substrate in the 'wrong side' of the isoelectric point (WSIP), i.e. when surface and peptide charges have the same sign. As a model system, we focus on the casein macropeptide (CMP), both in the aglycosylated (aCMP) and fully glycosydated (gCMP) forms. We model the substrate as a uniformly charged plane while CMP is treated as a bead-and-spring model including electrostatic interactions, excluded volume effects and acid/base equilibria. Adsorption coverage, aminoacid charges and concentration profiles are computed by means of Monte Carlo simulations at fixed pH and salt concentration. We conclude that for different reasons the CMP can be adsorbed to both positively and negatively charged surfaces in the WSIP. For negatively charged surfaces, WSIP adsorption is due to the patchy distribution of charges: the peptide is attached to the surface by the positively charged end of the chain, while the repulsion of the surface for the negatively charged tail is screened by the small ions of the added salt. This effect increases with salt concentration. Conversely, a positively charged substrate induces strong charge regulation of the peptide: the acidic groups are deprotonated, and the peptide becomes negatively charged. This effect is stronger at low salt concentrations and it is more intense for gCMP than for aCMP, due to the presence of the additional sialic groups in gCMP.

Nanopore Current Enhancements Lack Protein Charge Dependence and Elucidate Maximum Unfolding at Protein's Isoelectric Point

Protein sequencing, as well as protein fingerprinting, has gained tremendous attention in the electrical sensing realm of solid-state nanopores and is challenging due to fast translocations and the use of high molar electrolytes. Despite providing an appreciable signal-to-noise ratio, high electrolyte concentrations can have adverse effects on the native protein structure. Herein, we present a thorough investigation of low electrolyte sensing conditions across a broad pH and voltage range generating conductive pulses (CPs) irrespective of protein net charge. We used Cas9 as the model protein and demonstrated that unfolding is noncooperative, represented by the gradual elongation or stretching of the protein, and sensitive to both the applied voltage and pH (i.e., charge state). The magnitude of unfolding and the isoelectric point (pI) of Cas9 was found to be correlated and a critical factor in our experiments. Electroosmotic flow (EOF) was always aligned with the transit direction, whereas electrophoretic force (EPF) was either reinforcing (pH < pI) or opposing (pH > pI) the protein's movement, which led to slower translocations at higher pH values. Further exploration of higher pH values led to slowing down of protein with > 30% of the population being slower than 0.5 ms. Our results would be critical for protein sensing at very low electrolytes and to retard their translocation speed without resorting to high-bandwidth equipment.

Site Correlations, Capacitance, and Polarizability From Protein Protonation Fluctuations

We generalize the Kirkwood-Shumaker theory of protonation fluctuation for an anisotropic distribution of dissociable charges on a globular protein. The fluctuations of the total charge and the total dipole moment, in contrast to their average values, depend on the same proton occupancy correlator, thus exhibiting a similar dependence also on the solution pH. This has important consequences for the Kirkwood-Shumaker interaction and its dependence on the bathing solution conditions.

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

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

Simulation and machine learning modelling based comparative study of InAlGaN and AlGaN high electron mobility transistors for the detection of HER-2

The detection of the cancer biomarker human epidermal growth factor receptor 2 (HER-2) has always been challenging at the early stages of cancer due to its very small presence. A systematic study of biosensors to achieve optimum sensitivity is of paramount significance. Thus, in this paper, we report a simulation study and machine learning (ML) based model for the comparative analysis of indium aluminum gallium nitride (InAlGaN) and aluminum gallium nitride (AlGaN) based high electron mobility transistors (HEMTs) for the detection of HER-2. The sensing performance of the InAlGaN based HEMT exhibits 1.8 times higher sensitivity as compared to that of the AlGaN based HEMT. The presented work also provides insights into the importance of the pH of the medium of HER-2. The results produced by the developed ML-based model are in good agreement with the simulation results. The model is not only capable of predicting within the trained range but also it can predict reasonably well even beyond the range of the training data. The introduction of a ML-based model significantly reduces the computational cost, time to perform similar type of simulations and, unlike the physics-based modelling, it also eliminates the need for empirical fitting of the model parameters.

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.

Microchip-based structure determination of low-molecular weight proteins using cryo-electron microscopy

Interest in cryo-Electron Microscopy (EM) imaging has skyrocketed in recent years due to its pristine views of macromolecules and materials. As advances in instrumentation and computing algorithms spurred this progress, there is renewed focus to address specimen-related challenges. Here we contribute a microchip-based toolkit to perform complementary structural and biochemical analysis on low-molecular weight proteins. As a model system, we used the SARS-CoV-2 nucleocapsid (N) protein (48 kDa) due to its stability and important role in therapeutic development. Cryo-EM structures of the N protein monomer revealed a flexible N-terminal "top hat" motif and a helical-rich C-terminal domain. To complement our structural findings, we engineered microchip-based immunoprecipitation assays that led to the discovery of the first antibody binding site on the N protein. The data also facilitated molecular modeling of a variety of pandemic and common cold-related coronavirus proteins. Such insights may guide future pandemic-preparedness protocols through immuno-engineering strategies to mitigate viral outbreaks.

Charge regulation of colloidal particles in aqueous solutions

We study the charge regulation of colloidal particles inside aqueous electrolyte solutions. To stabilize a colloidal suspension against precipitation, colloidal particles are synthesized with either acidic or basic groups on their surface. On contact with water, these surface groups undergo proton transfer reactions, resulting in colloidal surface charge. The charge is determined by the condition of local chemical equilibrium between hydronium ions inside the solution and at the colloidal surface. We use a model of Baxter sticky spheres to explicitly calculate the equilibrium dissociation constants and to construct a theory which is able to quantitatively predict the effective charge of colloidal particles with either acidic or basic surface groups. The predictions of the theory for the model are found to be in excellent agreement with the results of Monte Carlo simulations. This theory is further extended to treat colloidal particles with a mixture of both acidic and basic surface groups.

Dual promoter strategy enhances co-expression of α-L-rhamnosidase and enhanced fluorescent protein for whole-cell catalysis and bioresource valorization

Developing circular economy is the only way to improve the efficiency of resource utilization. Whole-cell catalysis is an effective method to recycle enzymes, improve catalytic efficiency, and reduce production costs. The enzyme, α-L-rhamnosidase has considerable application prospects in the field of biocatalysis as it can hydrolyze a variety of α-L rhamnoses. In the present study, the genes for α-L-rhamnosidase (rhaB1) and enhanced fluorescent protein (EGFP) were co-expressed using a bi-promoter expression vector pRSFDuet1 and their enzymatic properties were evaluated. To our knowledge, this study has established an effective rhamnosidase-fluorescent indicator and whole-cell catalytic system for the first time. Moreover, we analyzed the change in the activity of the crude rhaB1-EGFP as well as its whole-cell during the biocatalysis process using fluorescence intensity. Recombinant rhaB1-EGFP as a product which contains rhaB1 and EGFP showed higher thermal stability, pH stability, and conversion efficiency than rhaB1, and its optimum temperature for rutin catalysis was ideal for industrial applications. Moreover, under the optimal conditions of a rutin concentration of 0.05 g/L, pH of 6.0, temperature of 40 °C, a yield of 92.5% was obtained. Furthermore, we demonstrated the relationship between the fluorescence intensity and enzyme activity. This study established a highly efficient whole-cell catalytic system whose activity can be evaluated by fluorescence intensity, providing a reference for enzyme recycling.

Multiscale Coarse-Grained Approach to Investigate Self-Association of Antibodies

Self-association of therapeutic monoclonal antibodies (mabs) are thought to modulate the undesirably high viscosity observed in their concentrated solutions. Computational prediction of such a self-association behavior is advantageous early during mab drug candidate selection when material availability is limited. Here, we present a coarse-grained (CG) simulation method that enables microsecond molecular dynamics simulations of full-length antibodies at high concentrations. The proposed approach differs from others in two ways: first, charges are assigned to CG beads in an effort to reproduce molecular multipole moments and charge asymmetry of full-length antibodies instead of only localized charges. This leads to great improvements in the agreement between CG and all-atom electrostatic fields. Second, the distinctive hydrophobic character of each antibody is incorporated through empirical adjustments to the short-range van der Waals terms dictated by cosolvent all-atom molecular dynamics simulations of antibody variable regions. CG simulations performed on a set of 15 different mabs reveal that diffusion coefficients in crowded environments are markedly impacted by intermolecular interactions. Diffusion coefficients computed from the simulations are in correlation with experimentally measured observables, including viscosities at a high concentration. Further, we show that the evaluation of electrostatic and hydrophobic characters of the mabs is useful in predicting the nonuniform effect of salt on the viscosity of mab solutions. This CG modeling approach is particularly applicable as a material-free screening tool for selecting antibody candidates with desirable viscosity properties.

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.

High-throughput computational pipeline for 3-D structure preparation and in silico protein surface property screening: A case study on HBcAg dimer structures

Knowledge-based experimental design can aid biopharmaceutical high-throughput screening (HTS) experiments needed to identify critical manufacturability parameters. Prior knowledge can be obtained via computational methods such as protein property extraction from 3-D protein structures. This study presents a high-throughput 3-D structure preparation and refinement pipeline that supports structure screenings with an automated and data-dependent workflow. As a case study, three chimeric virus-like particle (VLP) building blocks, hepatitis B core antigen (HBcAg) dimers, were constructed. Molecular dynamics (MD) refinement quality, speed, stability, and correlation to zeta potential data was evaluated using different MD simulation settings. Settings included 2 force fields (YASARA2 and AMBER03) and 2 pKa computation methods (YASARA and H++). MD simulations contained a data-dependent termination via identification of a 2 ns Window of Stability, which was also used for robust descriptor extraction. MD simulation with YASARA2, independent of pKa computation method, was found to be most stable and computationally efficient. These settings resulted in a fast refinement (6.6-37.5 h), a good structure quality (-1.17--1.13) and a strong linear dependence between dimer surface charge and complete chimeric HBcAg VLP zeta potential. These results indicate the computational pipeline's applicability for early-stage candidate assessment and design optimization of HTS manufacturability or formulability experiments.

pH-induced morphological changes of proteinaceous viral shells

Changes in environmental pH can induce morphological changes in empty proteinaceous shells of bacteriophages in vitro that are very similar to changes occurring in viral capsids in vivo after encapsidation of DNA. These changes in capsid shape and size cannot be explained with a simple elastic model alone. We propose a new theoretical framework that combines the elasticity of thin icosahedral shells with the pH dependence of capsid charge distribution. Minimization of the sum of elastic and electrostatic free energies leads to equilibrium shapes of viral shells that depend on a single elastic parameter and the detailed configuration of the imbedded protein charges. Based on the in vitro shell reconstructions of bacteriophage HK97 we elucidate the details of how the reversible transition between Prohead II and Expansion Intermediate II states of the HK97 procapsid is induced by pH changes, as well as some other features of the bacteriophage maturation.

Hydrodynamic and Electrophoretic Properties of Trastuzumab/HER2 Extracellular Domain Complexes as Revealed by Experimental Techniques and Computational Simulations

The combination of hydrodynamic and electrophoretic experiments and computer simulations is a powerful approach to study the interaction between proteins. In this work, we present hydrodynamic and electrophoretic experiments in an aqueous solution along with molecular dynamics and hydrodynamic modeling to monitor and compute biophysical properties of the interactions between the extracellular domain of the HER2 protein (eHER2) and the monoclonal antibody trastuzumab (TZM). The importance of this system relies on the fact that the overexpression of HER2 protein is related with the poor prognosis breast cancers (HER2++ positives), while the TZM is a monoclonal antibody for the treatment of this cancer. We have found and characterized two different complexes between the TZM and eHER2 proteins (1:1 and 1:2 TZM:eHER2 complexes). The conformational features of these complexes regulate their hydrodynamic and electrostatic properties. Thus, the results indicate a high degree of molecular flexibility in the systems that ultimately leads to higher values of the intrinsic viscosity, as well as lower values of diffusion coefficient than those expected for simple globular proteins. A highly asymmetric charge distribution is detected for the monovalent complex (1:1 complex), which has strong implications in correlations between the experimental electrophoretic mobility and the modeled net charge. In order to understand the dynamics of these systems and the role of the specific domains involved, it is essential to find biophysical correlations between dynamics, macroscopic transport and electrostatic properties. The results should be of general interest for researchers working in this area.

Insights into the ZIKV NS1 Virology from Different Strains through a Fine Analysis of Physicochemical Properties

The flavivirus genus has several organisms responsible for generating various diseases in humans. Recently, especially in tropical regions, Zika virus (ZIKV) has raised great health concerns due to the high number of cases affecting the area during the last years that has been accompanied by a rise in the cases of the Guillain-Barré syndrome and fetal and neonatal microcephaly. Diagnosis is still difficult since the clinical symptoms between ZIKV and other flaviviruses (e.g., dengue and yellow fever) are highly similar. The understanding of their common physicochemical properties that are pH-dependent and biomolecular interaction features and their differences sheds light on the relation strain-virulence and might suggest alternative strategies toward differential serological diagnostics and therapeutic intervention. Due to their immunogenicity, the primary focus of this study was on the ZIKV nonstructural proteins 1 (NS1). By means of computational studies and semiquantitative theoretical analyses, we calculated the main physicochemical properties of this protein from different strains that are directly responsible for the biomolecular interactions and, therefore, can be related to the differential infectivity of the strains. We also mapped the electrostatic differences at both the sequence and structural levels for the strains from Uganda to Brazil, which could suggest possible molecular mechanisms for the increase of the virulence of ZIKV in Brazil. Exploring the interfaces used by NS1 to self-associate in some different oligomeric states and interact with membranes and the antibody, we could map the strategy used by the ZIKV during its evolutionary process. This indicates possible molecular mechanisms that can be correlated with the different immunological responses. By comparing with the known antibody structure available for the West Nile virus, we demonstrated that this antibody would have difficulties to neutralize the NS1 from the Brazilian strain. The present study also opens up perspectives to computationally design high-specificity antibodies.

Assessing Force Field Potential Energy Function Accuracy via a Multipolar Description of Atomic Electrostatic Interactions in RNA

Computational investigations of RNA properties often rely on a molecular mechanical approach to define molecular potential energy. Force fields for RNA typically employ a point charge model of electrostatics, which does not provide a realistic quantum-mechanical picture. In reality, electron distributions around nuclei are not spherically symmetric and are geometry dependent. A multipole expansion method which allows for incorporation of polarizability and anisotropy in a force field is described, and its applicability to modeling the behavior of RNA molecules is investigated. Transferability of the model, critical for force field development, is also investigated.

General theory of charge regulation and surface differential capacitance

A generalization of the mean-field approach will be derived that will take into account the ion-ion as well as ion-surface non-electrostatic effects on an equal footing, being based on the bulk and surface equations of state in the absence of electrostatic interactions. This approach will be applied to the analysis of a single planar surface with dissociable sites with several models of the specific ion-surface non-electrostatic interactions, providing a general thermodynamic insight into the characteristics of the surface differential capacitance. The ion-surface interactions and ion-ion packing considerations at the surface will be shown to be more relevant than the bulk packing constraints for ions vicinal to the surface, as well as to set in prior to the conditions where the bulk packing constraints would become relevant.

Charge regulating macro-ions in salt solutions: screening properties and electrostatic interactions

We revisit the charge-regulation mechanism of macro-ions and apply it to mobile macro-ions in a bathing salt solution. In particular, we examine the effects of correlation between various adsorption/desorption sites and analyze the collective behavior in terms of the solution effective screening properties. We show that such a behavior can be quantified in terms of the charge asymmetry of the macro-ions, defined by their preference for a non-zero effective charge, and their donor/acceptor propensity for exchanging salt ions with the bathing solution. Asymmetric macro-ions tend to increase the screening, while symmetric macro-ions can in some cases decrease it. Macro-ions that are classified as donors display a rather regular behavior, while those that behave as acceptors exhibit an anomalous non-monotonic Debye length. The screening properties, in their turn, engender important modifications to the disjoining pressure between two charged surfaces. Our findings are in particular relevant for solutions of proteins, whose exposed amino acids can undergo charge dissociation/association processes to/from the bathing solution, and can be considered as a solution of charged regulated macro-ions, as analyzed here.

From discrete to continuous description of spherical surface charge distributions

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

Varieties of charge distributions in coat proteins of ssRNA+ viruses

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