No abstract

A new betacoronavirus named SARS-CoV-2 has emerged as a new threat to global health and economy. A promising target for both diagnosis and therapeutics treatments of the new disease named COVID-19 is the coronavirus (CoV) spike (S) glycoprotein. By constant-pH Monte Carlo simulations and the PROCEEDpKa method, we have mapped the electrostatic epitopes for four monoclonal antibodies and the angiotensin-converting enzyme 2 (ACE2) on both SARS-CoV-1 and the new SARS-CoV-2 S receptor binding domain (RBD) proteins. We also calculated free energy of interactions and shown that the S RBD proteins from both SARS viruses binds to ACE2 with similar affinities. However, the affinity between the S RBD protein from the new SARS-CoV-2 and ACE2 is higher than for any studied antibody previously found complexed with SARS-CoV-1. Based on physical chemical analysis and free energies estimates, we can shed some light on the involved molecular recognition processes, their clinical aspects, the implications for drug developments, and suggest structural modifications on the CR3022 antibody that would improve its binding affinities for SARS-CoV-2 and contribute to address the ongoing international health crisis.

Both natural and synthetic polyelectrolytes form strong complexes with a variety of proteins. One peculiar phenomenon is that association can take place even when the protein and the polyelectrolyte carry the same charge. This has been interpreted as if the ion-dipole interaction can overcome the repulsive ion-ion interaction. On the basis of Monte Carlo simulations and perturbation theory, we propose a different explanation for the association, namely, charge regulation. We have investigated three different protein-polymer complexes and found that the induced ionization of amino acid residues due to the polyelectrolyte leads to a surprisingly strong attractive interaction between the protein and the polymer. The extra attraction from this charge-induced charge interaction can be several kT and is for the three cases studied here, lysozyme, alpha-lactalbumin, and beta-lactoglobulin, of the same magnitude or stronger than the ion-dipole interaction. The magnitude of the induced charge is governed by a response function, the protein charge capacitance Z2-Z2. This fluctuation term can easily be calculated in a simulation or measured in a titration experiment.

Viruses are enthusiastically studied due to the great impact that these organisms can have on human health. Computational approaches can contribute offering tools that can shed light on important molecular mechanisms that help to design new diagnostic procedures. Several cellular processes between the immune-host system and the pathogenic organism are dependent on specific intermolecular interactions. In this study, we evaluated theoretical approaches to understand some properties of the antigen−antibody interactions considering the titratable properties of all ionizable residues of the nonstructural viral protein 1 (NS1) of the West Nile virus (WNV) and the Zika virus (ZIKV). Constant-pH Monte Carlo simulations were performed to estimate electrostatic properties such as the pK a shifts (ΔpK a ). We proposed an alternative criterion for the discrimination of antigenic residues based on ΔpK a s. Our outcomes were analyzed by an evaluation of the sensitivity and specificity through a receiver operating characteristic (ROC). As a starting point, we used the known crystallographic structure for the complex of NS1 WNV(176−352) and the specific antibody 22NS1 (PDB ID 4OII) to differentiate the residues belonging to that interface. With an optimal threshold for the absolute value of the pK a shifts, we found that is possible to predict antigenic epitopes reproducing the interfaces as defined by the X-ray structure. After this validation, we evaluated theoretical predictions based on protein− protein (PP) complexation simulations. From them, we observe amino acids with an antigenic potential and defined the optimum threshold that was applied for two strains of ZIKV (i.e., Uganda and Brazil). Several ionizable residues with antigenic capacity were identified. This is favorably related to some studies that show the high immunogenicity of secreted NS1. This approach opens up an important discussion about what are termed here "electrostatic epitopes" and how they work as an important reference in the paratope−epitope interaction for viral systems.

Proton exchange between titratable amino acid residues and the surrounding solution gives rise to exciting electric processes in proteins. We present a proton titration scheme for studying acid-base equilibria in Metropolis Monte Carlo simulations where salt is treated at the Debye-Hückel level. The method, rooted in the Kirkwood model of impenetrable spheres, is applied on the three milk proteins α-lactalbumin, β-lactoglobulin, and lactoferrin, for which we investigate the net-charge, molecular dipole moment, and charge capacitance. Over a wide range of pH and salt conditions, excellent agreement is found with more elaborate simulations where salt is explicitly included. The implicit salt scheme is orders of magnitude faster than the explicit analog and allows for transparent interpretation of physical mechanisms. It is shown how the method can be expanded to multiscale modeling of aqueous salt solutions of many biomolecules with nonstatic charge distributions. Important examples are protein-protein aggregation, protein-polyelectrolyte complexation, and protein-membrane association.

The interplay between the biocolloidal characteristics (especially size and charge), pH, salt concentration and the thermal energy results in a unique collection of mesoscopic forces of importance to the molecular organization and function in biological systems. By means of Monte Carlo simulations and semi-quantitative analysis in terms of perturbation theory, we describe a general electrostatic mechanism that gives attraction at low electrolyte concentrations. This charge regulation mechanism due to titrating amino acid residues is discussed in a purely electrostatic framework. The complexation data reported here for interaction between a polyelectrolyte chain and the proteins albumin, goat and bovine a-lactalbumin, b-lactoglobulin, insulin, k-casein, lysozyme and pectin methylesterase illustrate the importance of the charge regulation mechanism. Special attention is given to pH y pI where iondipole and charge regulation interactions could overcome the repulsive ion-ion interaction. By means of protein mutations, we confirm the importance of the charge regulation mechanism, and quantify when the complexation is dominated either by charge regulation or by the ion-dipole term.

pH is a key parameter for technological and biological processes, intimately related to biomolecular charge. As such, it controls biomolecular conformation and intermolecular interactions, for example, protein/RNA stability and folding, enzyme activity, regulation through conformational switches, protein-polyelectrolyte association, and protein-RNA interactions. pH also plays an important role in technological systems in food, brewing, pharma, bioseparations, and biomaterials in general. Predicting the structure of large proteins and complexes remains a great challenge experimentally, industrially, and theoretically, despite the variety of numerical schemes available ranging from Poisson-Boltzmann approaches to explicit solvent based methods. In this work we benchmark a fast proton titration scheme against experiment and several theoretical methods on the following set of representative proteins: [HP36, BBL, HEWL (triclinic and orthorhombic), RNase, SNASE (V66K/WT, V66K/PHS, V66K/Δ+PHS, L38D/Δ+PHS, L38E/Δ+PHS, L38K/Δ+PHS), ALAC, and OMTKY3]; routinely used in similar tests due to the diversity of their structural features. Our scheme is rooted in the classical Tanford-Kirkwood model of impenetrable spheres, where salt is treated at the Debye-Hückel level. Treating salt implicitly dramatically reduces the computation time, thereby circumventing sampling difficulties faced by other numerical schemes. In comparison with experimental measurements, our calculated pK values have the average, maximum absolute, and root-mean-square deviations of 0.4-0.9, 1.0-5.2, and 0.5-1.2 pH units, respectively. These values are within the ranges commonly observed in theoretical models. They are also in the large majority of the cases studied here more accurate than the NULL model. For BBL, ALAC, and OMTKY3, the predicted pK are closer to experimental results than other analyzed theoretical data. Despite the intrinsic approximations of the fast titration scheme, its robustness and ability to properly describe the main system physics is confirmed.

Septins are GTP binding proteins considered to be novel components of the cytoskeleton. They polymerize into filaments based on hexameric or octameric core particles in which two copies of either three or four different septins, respectively, assemble into a specific sequence. Viable combinations of the 13 human septins are believed to obey substitution rules in which the different septins involved must come from distinct subgroups. The hexameric assembly, for example, has been reported to be SEPT7–SEPT6–SEPT2–SEPT2–SEPT6–SEPT7. Here, we have replaced SEPT2 by SEPT5 according to the substitution rules and used transmission electron microscopy to demonstrate that the resulting recombinant complex assembles into hexameric particles which are inverted with respect that predicted previously. MBP‐SEPT5 constructs and immunostaining show that SEPT5 occupies the terminal positions of the hexamer. We further show that this is also true for the assembly including SEPT2, in direct contradiction with that reported previously. Consequently, both complexes expose an NC interface, as reported for yeast, which we show to be more susceptible to high salt concentrations. The correct assembly for the canonical combination of septins 2‐6‐7 is therefore established to be SEPT2–SEPT6–SEPT7–SEPT7–SEPT6–SEPT2, implying the need for revision of the mechanisms involved in filament assembly.