Researcher profile

Christopher D. Cooper

Christopher D. Cooper contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

4 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.

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.

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.

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/.