Thermostabilizing mechanisms of canonical single amino acid substitutions at a GH1 ß-glucosidase probed by multiple MD and computational approaches.

ß-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable ß-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 ß-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized ß-glucosidases for 2G-biofuel production for industry, biotechnology, and science.

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.

Unravelling the Interaction of Piperlongumine with the Nucleotide-Binding Domain of HSP70: A Spectroscopic and In Silico Study.

Piperlongumine (PPL) is an alkaloid extracted from several pepper species that exhibits anti-inflammatory and anti-carcinogenic properties. Nevertheless, the molecular mode of action of PPL that confers such powerful pharmacological properties remains unknown. From this perspective, spectroscopic methods aided by computational modeling were employed to characterize the interaction between PPL and nucleotide-binding domain of heat shock protein 70 (NBD/HSP70), which is involved in the pathogenesis of several diseases. Steady-state fluorescence spectroscopy along with time-resolved fluorescence revealed the complex formation based on a static quenching mechanism. Van't Hoff analyses showed that the binding of PPL toward NBD is driven by equivalent contributions of entropic and enthalpic factors. Furthermore, IDF and Scatchard methods applied to fluorescence intensities determined two cooperative binding sites with Kb of (6.3 ± 0.2) × 104 M-1. Circular dichroism determined the thermal stability of the NBD domain and showed that PPL caused minor changes in the protein secondary structure. Computational simulations elucidated the microenvironment of these interactions, showing that the binding sites are composed mainly of polar amino acids and the predominant interaction of PPL with NBD is Van der Waals in nature.

Efficient mesh refinement for the Poisson-Boltzmann equation with boundary elements.

The Poisson-Boltzmann equation is a widely used model to study electrostatics in molecular solvation. Its numerical solution using a boundary integral formulation requires a mesh on the molecular surface only, yielding accurate representations of the solute, which is usually a complicated geometry. Here, we utilize adjoint-based analyses to form two goal-oriented error estimates that allow us to determine the contribution of each discretization element (panel) to the numerical error in the solvation free energy. This information is useful to identify high-error panels to then refine them adaptively to find optimal surface meshes. We present results for spheres and real molecular geometries, and see that elements with large error tend to be in regions where there is a high electrostatic potential. We also find that even though both estimates predict different total errors, they have similar performance as part of an adaptive mesh refinement scheme. Our test cases suggest that the adaptive mesh refinement scheme is very effective, as we are able to reduce the error one order of magnitude by increasing the mesh size less than 20% and come out to be more efficient than uniform refinement when computing error estimations. This result sets the basis toward efficient automatic mesh refinement schemes that produce optimal meshes for solvation energy calculations.

Molecular Dynamics Gives New Insights into the Glucose Tolerance and Inhibition Mechanisms on ß-Glucosidases.

ß-Glucosidases are enzymes with high importance for many industrial processes, catalyzing the last and limiting step of the conversion of lignocellulosic material into fermentable sugars for biofuel production. However, ß-glucosidases are inhibited by high concentrations of the product (glucose), which limits the biofuel production on an industrial scale. For this reason, the structural mechanisms of tolerance to product inhibition have been the target of several studies. In this study, we performed in silico experiments, such as molecular dynamics (MD) simulations, free energy landscape (FEL) estimate, Poisson-Boltzmann surface area (PBSA), and grid inhomogeneous solvation theory (GIST) seeking a better understanding of the glucose tolerance and inhibition mechanisms of a representative GH1 ß-glucosidase and a GH3 one. Our results suggest that the hydrophobic residues Y180, W350, and F349, as well the polar one D238 act in a mechanism for glucose releasing, herein called "slingshot mechanism", dependent also on an allosteric channel (AC). In addition, water activity modulation and the protein loop motions suggest that GH1 ß-Glucosidases present an active site more adapted to glucose withdrawal than GH3, in consonance with the GH1s lower product inhibition. The results presented here provide directions on the understanding of the molecular mechanisms governing inhibition and tolerance to the product in ß-glucosidases and can be useful for the rational design of optimized enzymes for industrial interests.

Synergistic long-range effects of mutations underlie aggregation propensities of amylin analogues.

The USFDA has approved pramlintide, commercially named Symlin (sIAPP), as adjunctive therapy for type 2 diabetes (T2D). This analogue of the human amylin peptide (hIAPP) has triple proline substitutions typical of the rat isoform (rIAPP). Recently, it was proposed that pramlintide solubility and aggregation resistance might be improved by incorporating further mutations, as S20R, screened from the wild-type porcine isoform (pIAPP), which leads to the variant named sIAPP+. To better elucidate how such properties might be systematically induced in rationally designed analogues, we performed comparative assessments of rIAPP, sIAPP, and sIAPP+ using replica-exchange molecular dynamics (REMD) with an accurate combination of force field Charmm22* and explicit aqueous solvation TIP4P/Ew. Our thermo-structural analyses show that sIAPP exhibits a thermal conversion channel of helices[Formula: see text]-sheets resembling hIAPP. This channel is depleted in rIAPP and is absent in sIAPP+. As a consequence, sIAPP+ presents an overall decrease of ß-like secondary structures and an overstabilization of α-helices. Additionally, we observed in rIAPP and sIAPP+ an increase in the backbone RMSF of molecular terminals and the exposed area of key residues. These structural features of sIAPP+ suggest a nonamyloidogenic character, which is corroborated by our judicious estimate of the electrostatic component of the solvation free energy using a generalized Born model, and so it may constitute an alternative strategy to sIAPP as a peptide analogue of hIAPP. Furthermore, our findings confirm that different aggregation propensities of amylin and its analogues are synergistically modulated by long-range effects of key mutations. Graphical Abstract S20R-Pramlintide.

Comparative Analysis of Electrostatic Models for Ligand Docking.

The precise modeling of molecular interactions remains an important goal among molecular modeling techniques. Some of the challenges in the field include the precise definition of a Hamiltonian for biomolecular systems, together with precise parameters derived from Molecular Mechanics Force Fields, for example. The problem is even more challenging when interaction energies from different species are computed, such as the interaction energy involving a ligand and a protein, given that small differences must be computed from large energies. Here we evaluated the effects of the electrostatic model for ligand binding energy evaluation in the context of ligand docking. For this purpose, a classical Coulomb potential with distance-dependent dielectrics was compared with a Poisson-Boltzmann (PB) model for electrostatic potential computation, based on DelPhi calculations. We found that, although the electrostatic energies were highly correlated for the Coulomb and PB models, the ligand pose and the enrichment of actual ligands against decoy compounds, were improved when binding energies were computed using PB as compared to the Coulomb model. We observed that the electrostatic energies computed with the Coulomb model were, on average, ten times larger than the energies computed with the PB model, suggesting a strong overestimation of the polar interactions in the Coulomb model. We also found that a slightly smoothed Lennard-Jones potential combined with the PB model resulted in a good compromise between ligand sampling and energetic scoring.

A Boundary-Integral Approach for the Poisson-Boltzmann Equation with Polarizable Force Fields.

Implicit-solvent models are widely used to study the electrostatics in dissolved biomolecules, which are parameterized using force fields. Standard force fields treat the charge distribution with point charges; however, other force fields have emerged which offer a more realistic description by considering polarizability. In this work, we present the implementation of the polarizable and multipolar force field atomic multipole optimized energetics for biomolecular applications (AMOEBA), in the boundary integral Poisson-Boltzmann solver PyGBe. Previous work from other researchers coupled AMOEBA with the finite-difference solver APBS, and found difficulties to effectively transfer the multipolar charge description to the mesh. A boundary integral formulation treats the charge distribution analytically, overlooking such limitations. This becomes particularly important in simulations that need high accuracy, for example, when the quantity of interest is the difference between solvation energies obtained from separate calculations, like happens for binding energy. We present verification and validation results of our software, compare it with the implementation on APBS, and assess the efficiency of AMOEBA and classical point-charge force fields in a Poisson-Boltzmann solver. We found that a boundary integral approach performs similarly to a volumetric method on CPU. Also, we present a GPU implementation of our solver. Moreover, with a boundary element method, the mesh density to correctly resolve the electrostatic potential is the same for standard point-charge and multipolar force fields. Finally, we saw that for binding energy calculations, a boundary integral approach presents more consistent results than a finite difference approximation for multipolar force fields. © 2019 Wiley Periodicals, Inc.

Cooperative hydrogen bonds and mobility of the non-aromatic ring as selectivity determinants for human acetylcholinesterase to similar anti-Alzheimer's galantaminics: a computational study.

Galantamine (Gnt) is a natural alkaloid inhibitor of acetylcholinesterase and is presently one of the most used drugs in the treatment against Alzheimer's disease during both the initial and intermediate stages. Among several natural Gnt derivatives, sanguinine (Sng) and lycoramine (Lyc) attract attention because of the way their subtle chemical differences from Gnt lead to drastic and opposite distinctions in inhibitory effects. However, to date, there is no solved structure for these natural derivatives. In the present study, we applied computational modeling and free energy calculation methods to better elucidate the molecular basis of the subtle distinctions between these derivatives and Gnt. The results showed that differences in the mobility of the non-aromatic ring carried by the Lyc-like sp2-sp3 modification display drastic conformational, vibrational, and entropic penalties at binding compared to Gnt. Additionally, the establishment of a stronger hydrogen bond network added enthalpic advantages for the linkage of the Sng-like methoxy-hydroxy substituted ligands. These results, which suggest an affinity ranking in agreement with that found in the literature, provided insights that are helpful for future planning and development of new anti-Alzheimer's disease drugs.

A modified Poisson-Boltzmann equation applied to protein adsorption.

Ion-exchange chromatography has been widely used as a standard process in purification and analysis of protein, based on the electrostatic interaction between the protein and the stationary phase. Through the years, several approaches are used to improve the thermodynamic description of colloidal particle-surface interaction systems, however there are still a lot of gaps specifically when describing the behavior of protein adsorption. Here, we present an improved methodology for predicting the adsorption equilibrium constant by solving the modified Poisson-Boltzmann (PB) equation in bispherical coordinates. By including dispersion interactions between ions and protein, and between ions and surface, the modified PB equation used can describe the Hofmeister effects. We solve the modified Poisson-Boltzmann equation to calculate the protein-surface potential of mean force, treated as spherical colloid-plate system, as a function of process variables. From the potential of mean force, the Henry constants of adsorption, for different proteins and surfaces, are calculated as a function of pH, salt concentration, salt type, and temperature. The obtained Henry constants are compared with experimental data for several isotherms showing excellent agreement. We have also performed a sensitivity analysis to verify the behavior of different kind of salts and the Hofmeister effects.

The electrostatic behavior of the bacterial cell wall using a smoothing function to describe the charge-regulated volume charge density profile.

The Donnan potential can be observed in many biological systems due to the presence of polyelectrolytes as proteins and nucleic acids. The aim of this work was to present a useful tool to describe the fixed and charge-regulated volume charge density profile through the use of a smoothing function and to obtain the electrostatic potential profile as well as the Donnan potential of this system by solving Poisson-Boltzmann (PB) equation. When we use the smoothing function, the Donnan potential arises automatically from the solution of only one Poisson-Boltzmann equation and it is not necessary to impose this potential for treating charged system in the presence of a membrane. The electrostatic behavior across the Bacillus brevis wall considering the dependence on the ionization of the cell wall functional groups as a function of the solution pH was analyzed. An important issue was to show that potentiometric titration data could be used together with the Poisson-Boltzmann equation to predict the electrostatic behavior (e.g., zeta potential) of the bacterial cell surface.

Modeling the electrostatic potential of asymmetric lipopolysaccharide membranes: the MEMPOT algorithm implemented in DelPhi.

Four chemotypes of the rough lipopolysaccharides (LPS) membrane from Pseudomonas aeruginosa were investigated by a combined approach of explicit water molecular dynamics (MD) simulations and Poisson-Boltzmann continuum electrostatics with the goal to deliver the distribution of the electrostatic potential across the membrane. For the purpose of this investigation, a new tool for modeling the electrostatic potential profile along the axis normal to the membrane, MEMbrane POTential (MEMPOT), was developed and implemented in DelPhi. Applying MEMPOT on the snapshots obtained by MD simulations, two observations were made: (a) the average electrostatic potential has a complex profile but is mostly positive inside the membrane due to the presence of Ca(2+) ions, which overcompensate for the negative potential created by lipid phosphate groups; and (b) correct modeling of the electrostatic potential profile across the membrane requires taking into account the water phase, while neglecting it (vacuum calculations) results in dramatic changes including a reversal of the sign of the potential inside the membrane. Furthermore, using DelPhi to assign different dielectric constants for different regions of the LPS membranes, it was investigated whether a single frame structure before MD simulations with appropriate dielectric constants for the lipid tails, inner, and the external leaflet regions, can deliver the same average electrostatic potential distribution as obtained from the MD-generated ensemble of structures. Indeed, this can be attained by using smaller dielectric constant for the tail and inner leaflet regions (mostly hydrophobic) than for the external leaflet region (hydrophilic) and the optimal dielectric constant values are chemotype-specific.

Attraction between two similar particles in an electrolyte: effects of Stern layer absorption

When Debye length is comparable or larger than the distance between two identical particles, the overlapping among the particles double-layers can play an important role in their interactions. This paper presents a theoretical analysis of the interaction among two identical particles with overlapped double-layers. We particularly focus on the effect of a Stern electro static condition from linearization of the adsorption isotherm near the isoelectric (neutrality) point in order to capture how polyvalent ion condensation affect sand reverses the surface charge. The stationary potential problem is solved within the framework of an asymptotic lubrication approach for a mean-field Poisson-Boltzmann model. Both spherical and cylindrical particles are analyzed. The results are finally discussed in the context of Debye-Hückel (D-H) limit and beyond it.

Quando o comprimento de Debye é comparável ou maior do que as distâncias entre duas partículas idênticas, a interseção entre as duplas camadas pode desempenhar papel importante na interação entre elas. Este artigo apresenta uma análise teórica da interação entre duas partículas idênticas as quais apresentam interseção entre as camadas duplas. Nós particularmente focamos a análise sobre o efeito da condição eletrostática de Stern a partir da linearização da isoterma de adsorção perto do ponto isoelétrico para capturar como a condensação do íon polivalente afeta e reverte a densidade de carga superficial. O problema que governa o potencial estacionário é resolvido no contexto de uma teoria assintótica de lubrificação para o modelo de Poisson-Boltzmann. O modelo é analisado para partículas cilíndricas e esféricas. Os resultados são finalmente discutidos no contexto do limite de Debye-Hückel e além dele.

Electro-osmosis in kaolinite with pH-dependent surface charge modelling by homogenization

A new three-scale model to describe the coupling between pH-dependent flows and transient ion transport, including adsorption phenomena in kaolinite clays, is proposed. The kaolinite is characterized by three separate nano/micro and macroscopic length scales. The pore (micro)-scale is characterized by micro-pores saturated by an aqueous solution containing four monovalent ions and charged solid particles surrounded by thin electrical double layers. The movement of the ions is governed by the Nernst-Planck equations, and the influence of the double layers upon the flow is dictated by the Helmholtz-Smoluchowski slip boundary condition on the tangential velocity. In addition, an adsorption interface condition for the Na+ transportis postulated to capture its retention in the electrical double layer. Thetwo-scalenano/micro model including salt adsorption and slip boundary condition is homogenized to the Darcy scale and leads to the derivation of macroscopic governing equations. One of the notable features of the three-scale model is there construction of the constitutive law of effective partition coefficient that governs the sodium adsorption in the double layer. To illustrate the feasibility of the three-scale model in simulating soil decontamination by electrokinetics, the macroscopic model is discretized by the finite volume method and the desalination of a kaolinite sample by electrokinetics is simulated.

Neste artigo propomos um modelo em três escalas para descrever o acoplamento entre o fluxo eletroosmótico e o transporte de íons incluindo fenômenos de adsorção em uma caulinita. A argila é caracterizada por três escalas nano/micro e macroscópica. A escala microscópica é constituída por micro-poros saturados por uma solução aquosa contendo quatro íons monovalentes e partículas sólidas carregadas eletricamente circundadas por uma dupla camada elétrica fina. O movimento dos íons é governado pelas equações de Nernst-Planck e a influência da dupla camada sobre o fluxo aquoso é modelada por uma condição de contorno de deslizamento da componente tangencial do campo de velocidade (condição de Helmholtz-Smoluchowski). Além disso, uma condição de adsorção na interface fluido-sólido para os íons Na+ é postulada capturando a retenção do sódio na dupla camada elétrica. O modelo em duas escalas nano/micro incluindo a adsorção do sal e a condição de deslizamento da velocidade é homogeneizado levando a derivação das equações macroscópicas na escala de Darcy. Um dos aspectos inovadores do modelo em três escalas é a reconstrução da lei constitutiva para o coeficiente de partição que governa a adsorção do Na+ na dupla camada elétrica. Para ilustrar as potencialidades do modelo em três escalas na simulação da eletroremediação de solos argilosos, o modelo macroscópico é discretizado utilizando o método de volumes finitos no intuito de simular a dessalinização de uma amostra de caulinita por técnica de eletrocinética.