BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christopher D. Cooper

Christopher D. Cooper appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phphysics.comp-phBiological Physicscond-mat.soft

Catalog footprint

What is connected

7works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2023arXiv

Accurate boundary-integral formulations for the calculation of electrostatic forces with an implicit-solvent model

An accurate force calculation with the Poisson-Boltzmann equation is challenging, as it requires the electric field on the molecular surface. Here, we present a calculation of the electric field on the solute-solvent interface that is exact for piece-wise linear variations of the potential and analyze four different alternatives to compute the force using a boundary element method. We performed a verification exercise for two cases: the isolated and two interacting molecules. Our results suggest that the boundary element method outperforms the finite difference method, as the latter needs a much finer mesh than in solvation energy calculations to get acceptable accuracy in the force, whereas the same surface mesh than a standard energy calculation is appropriate for the boundary element method. Among the four evaluated alternatives of force calculation, we saw that the most accurate one is based on the Maxwell stress tensor. However, for a realistic application, like the barnase-barstar complex, the approach based on variations of the energy functional, which is less accurate, gives equivalent results. This analysis is useful towards using the Poisson-Boltzmann equation for force calculations in applications where high accuracy is key, for example, to feed molecular dynamics models or to enable the study of the interaction between large molecular structures, like viruses adsorbed onto substrates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predicting the orientation of adsorbed proteins steered with electric fields using a simple electrostatic model

Under the most common experimental conditions, the adsorption of proteins to solid surfaces is an spontaneous process that leads to a rather compact layer of randomly oriented molecules. Due to the importance of this process for the development of catalytic surfaces, a number of existing computational and experimental approaches try to predict and control the orientation of such molecules. However, and despite their own advantages, these tend to be either too expensive computationally, or oversimplified, undermining their ability to predict the most appropriate experimental conditions to maximize the catalytic activity of adsorbed proteins. To address this current need, we present an efficient computational approach to model the behavior of proteins near surfaces in the presence of an external electric field, based on continuum electrostatics. Our model can not only estimate the overall affinity of the protein with the surface, but also their most likely orientation as a function of the potential applied. In this way, a rational selection of the potential can be performed to maximize the accessibility of the protein's active site to the solvent. The model relies on the Poisson-Boltzmann equation and was implemented in an extension of the code PyGBe that includes an external electric field, and renders the electrostatic component of the solvation free energy. To demonstrate the feasibility of this technique, we investigate the adsorption of trypsin onto a carbon electrode under potentiostatic conditions both numerically and experimentally. We found that even though the adsorption process is largely dominated by hydrophobic effects, the orientation of trypsin can be controlled through an external potential, influencing the position of the active sites, and resulting in an important change in the catalytic activity of the surface.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Simple Electrostatic Model for the Hard-Sphere Solute Component of Nonpolar Solvation

We propose a new model for estimating the free energy of forming a molecular cavity in a solvent, by assuming this energy is dominated by the electrostatic energy associated with creating the static (interface) potential inside the cavity. The new model approximates the cavity-formation energy as that of a shell capacitor: the inner, solute-shaped conductor is held at the static potential, and the outer conductor (at the first solvation shell) is held at zero potential. Compared to cavity energies computed using free-energy pertubation with explicit-solvent molecular dynamics, the new model exhibits surprising accuracy (Mobley test set, RMSE 0.45 kcal/mol). Combined with a modified continuum model for solute-solvent van der Waals interactions, the total nonpolar model has RMSE of 0.55 kcal/mol on this test set, which is remarkable because the two terms largely cancel. The overall nonpolar model has a small number of physically meaningful parameters and compares favorably to other published models of nonpolar solvation. Finally, when the proposed nonpolar model is combined with our solvation-layer interface condition (SLIC) continuum electrostatic model, which includes asymmetric solvation-shell response, we predict solvation free energies with an RMS error of 1.35 kcal/mol relative to experiment, comparable to the RMS error of explicit-solvent FEP (1.26 kcal/mol). Moreover, all parameters in our model have a clear physical meaning, and employing reasonable temperature dependencies yields remarkable correlation with solvation entropies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computational nanoplasmonics in the quasistatic limit for biosensing applications

This work uses the long-wavelength limit to compute LSPR response of biosensors, expanding the open-source PyGBe code to compute the extinction cross-section of metallic nanoparticles in the presence of any target for sensing. The target molecule is represented by a surface mesh, based on its crystal structure. PyGBe is research software for continuum electrostatics, written in Python with computationally expensive parts accelerated on GPU hardware, via PyCUDA. It is also accelerated algorithmically via a treecode that offers O(N log N) computational complexity. These features allow PyGBe to handle problems with half a million boundary elements or more. Using a model problem consisting of an isolated silver nanosphere in an electric field, our results show grid convergence as 1/N, and accurate computation of the extinction cross-section as a function of wavelength (compared with an analytical solution). For a model of a sensor-analyte system, consisting of a spherical silver nanoparticle and a set of bovine serum albumin (BSA) proteins, our results again obtain grid convergence as 1/N (with respect to the Richardson extrapolated value). Computing the LSPR response as a function of wavelength in the presence of BSA proteins captures a red-shift of 0.5 nm in the resonance frequency due to the presence of the analytes at 1-nm distance. The final result is a sensitivity study of the biosensor model, obtaining the shift in resonance frequency for various distances between the proteins and the nanoparticle. All results in this paper are fully reproducible, and we have deposited in archival data repositories all the materials needed to run the computations again and re-create the figures. PyGBe is open source under a permissive license and openly developed. Documentation is available at http://barbagroup.github.io/pygbe/docs/.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Poisson-Boltzmann model for protein-surface electrostatic interactions and grid-convergence study using the PyGBe code

Interactions between surfaces and proteins occur in many vital processes and are crucial in biotechnology: the ability to control specific interactions is essential in fields like biomaterials, biomedical implants and biosensors. In the latter case, biosensor sensitivity hinges on ligand proteins adsorbing on bioactive surfaces with a favorable orientation, exposing reaction sites to target molecules. Protein adsorption, being a free-energy-driven process, is difficult to study experimentally. This paper develops and evaluates a computational model to study electrostatic interactions of proteins and charged nanosurfaces, via the Poisson-Boltzmann equation. We extended the implicit-solvent model used in the open-source code PyGBe to include surfaces of imposed charge or potential. This code solves the boundary integral formulation of the Poisson-Boltzmann equation, discretized with surface elements. PyGBe has at its core a treecode-accelerated Krylov iterative solver, resulting in O(N log N) scaling, with further acceleration on hardware via multi-threaded execution on \gpu s. It computes solvation and surface free energies, providing a framework for studying the effect of electrostatics on adsorption. We then derived an analytical solution for a spherical charged surface interacting with a spherical molecule, then completed a grid-convergence study to build evidence on the correctness of our approach. The study showed the error decaying with the average area of the boundary elements, i.e., the method is O(1/N), which is consistent with our previous verification studies using PyGBe. We also studied grid-convergence using a real molecular geometry (protein GB1D4'), in this case using Richardson extrapolation (in the absence of an analytical solution) and confirmed the O(1/N) scaling in this case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probing protein orientation near charged nanosurfaces for simulation-assisted biosensor design

Protein-surface interactions are ubiquitous in biological processes and bioengineering, yet are not fully understood. In biosensors, a key factor determining the sensitivity and thus the performance of the device is the orientation of the ligand molecules on the bioactive device surface. Adsorption studies thus seek to determine how orientation can be influenced by surface preparation. In this work, protein orientation near charged nanosurfaces is obtained under electrostatic effects using the Poisson-Boltzmann equation, in an implicit-solvent model. Sampling the free energy for protein GB1D4' at a range of tilt and rotation angles with respect to the charged surface, we calculated the probability of the protein orientations and observed a dipolar behavior. This result is consistent with published experimental studies and combined Monte Carlo and molecular dynamics simulations using this small protein, validating our method. More relevant to biosensor technology, antibodies such as immunoglobulin G are still a formidable challenge to molecular simulation, due to their large size. We obtained the probability distribution of orientations for the iso-type IgG2a at varying surface charge and salt concentration. This iso-type was not found to have a preferred orientation in previous studies, unlike the iso-type IgG1 whose larger dipole moment was assumed to make it easier to control. We find that the preferred orientation of IgG2a can be favorable for biosensing with positive surface charge of 0.05C/m$^{2}$ or higher and 37mM salt concentration. The results also show that local interactions dominate over dipole moment for this protein. Improving immunoassay sensitivity may thus be assisted by numerical studies using our method (and open-source code), guiding changes to fabrication protocols or protein engineering of ligand molecules to obtain more favorable orientations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A biomolecular electrostatics solver using Python, GPUs and boundary elements that can handle solvent-filled cavities and Stern layers

The continuum theory applied to bimolecular electrostatics leads to an implicit-solvent model governed by the Poisson-Boltzmann equation. Solvers relying on a boundary integral representation typically do not consider features like solvent-filled cavities or ion-exclusion (Stern) layers, due to the added difficulty of treating multiple boundary surfaces. This has hindered meaningful comparisons with volume-based methods, and the effects on accuracy of including these features has remained unknown. This work presents a solver called PyGBe that uses a boundary-element formulation and can handle multiple interacting surfaces. It was used to study the effects of solvent-filled cavities and Stern layers on the accuracy of calculating solvation energy and binding energy of proteins, using the well-known APBS finite-difference code for comparison. The results suggest that if required accuracy for an application allows errors larger than about 2%, then the simpler, single-surface model can be used. When calculating binding energies, the need for a multi-surface model is problem-dependent, becoming more critical when ligand and receptor are of comparable size. Comparing with the APBS solver, the boundary-element solver is faster when the accuracy requirements are higher. The cross-over point for the PyGBe code is in the order of 1-2% error, when running on one GPU card (NVIDIA Tesla C2075), compared with APBS running on six Intel Xeon CPU cores. PyGBe achieves algorithmic acceleration of the boundary element method using a treecode, and hardware acceleration using GPUs via PyCuda from a user-visible code that is all Python. The code is open-source under MIT license.

0 citations0 reviewsNo code linkedOpen access