{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Silver nanoparticle response to Bovine Serum Albumina (BSA) \n", "\n", "In section [03_LSPR_sensor_response_POC](http://nbviewer.jupyter.org/github/barbagroup/pygbe_lspr/blob/master/reports/03_LSPR_sensor_response_POC.ipynb) of this report we present a proof of concept of Localized surface plasmon resonance sensor (LSPR) response. We show how [PyGBe](https://github.com/barbagroup/pygbe) can be used to model the response of LSPR biosensors by modeling the BSA as spheres. \n", "\n", "In this report, we present a more realistic approach where we use the geometry of the actual protein and its charges. We obtain the crystal structure of BSA (4f5s) from the [Protein Data Bank](https://www.rcsb.org/) and generate the charge distibution and mesh using [`pdb2pqr`](http://nbcr-222.ucsd.edu/pdb2pqr_2.1.1/) and [NanoShaper](https://www.electrostaticszone.eu/downloads/25-nanoshaper-0-7-linux64-static) respectively. \n", "\n", "We study the interaction between a realistic model of the BSA protein (non-spherical) and the nanoparticle; more specifically how the relative position between them affects sensitivity.\n", "\n", "\n", "
\n", "#### Fig 1. Setup for response calculation using crystal structure of real protein. \n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Mesh convergence\n", "\n", "Before running the response calculations we perform a mesh convergence analysis to ensure that the results are not going to be affected by changing the size of the mesh. \n", "\n", "Since we compute the extinction cross section over the sphere nanoparticle, we pick a fixed mesh density for the protein and varied the refinement over the sphere. After studying different cases we decide to use a mesh of density 2 triangles per Angstrom square for the protein. \n", "\n", "We run the case presented in Fig. 1 for radius of the sphere of $25 \\, nm$, a wavelength of $380 \\, nm$ and a distance $d\\,=\\,1 \\,nm$ for meshes of 512, 2048, 8192 and 32768 elements.\n", "\n", "The values for the dielectric constants for water and silver were obtained by interpolation of experimental data (2,3):\n", "\n", "* Water : 1.7972 + 8.5048e-09j\n", "* Silver: -3.3877 + 0.1922j\n", "\n", "For the protein (modeled as a sphere) we use Bovine Serum Albumina (BSA) data extracted from the functional relationship provided by Pahn, et al. (1)\n", "\n", "* Protein: 2.7514 + 0.2860j\n", "\n", "The parameters used to perform the convergence analysis are:\n", "\n", "```\n", "Precision double\n", "K 4\n", "Nk 9 \n", "K_fine 37\n", "thresold 0.5\n", "BSZ 128\n", "restart 100\n", "tolerance 1e-5\n", "max_iter 1000\n", "P 15\n", "eps 1e-12\n", "NCRIT 500\n", "theta 0.5 \n", "GPU 1\n", "```\n", "\n", "The error calculation uses the Richardson extrapolation value of the extinction cross-section as a reference, $Cext\\,=\\, 6799.27 \\, nm^2$. We obtain an order of convergence of $0.98$" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " \n", "
\n", "#### Fig 2. Convergence of extinction cross-section of a $r=25\\, nm$ silver sphere with BSA protein located at $1\\,nm$ away along the z-axis, versus the number of elements. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Relaxiation of parameters to optimize run time \n", "\n", "We analyzed different combinations of simulation parameters in order to reduce the time of computation and have a relative error (compared to the Richardson extrapolation) smaller %2.\n", "\n", "The final parameters picked are:\n", "\n", "Density of protein mesh: 1 triangle/Ang^2\n", "\n", "```\n", "Precision double\n", "K 4\n", "Nk 9 \n", "K_fine 19\n", "thresold 0.5\n", "BSZ 128\n", "restart 100\n", "tolerance 1e-3\n", "max_iter 1000\n", "P 6\n", "eps 1e-12\n", "NCRIT 500\n", "theta 0.5 \n", "GPU 1\n", "```\n", "\n", "After this relaxation the percentage error compared to the Richardson extrapolation is %0.63 and the time of computation is approximately 7.5 min in a K40. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "##### References\n", "\n", "(1) Phan, Anh D., et al. \"Surface plasmon resonances of protein-conjugated gold nanoparticles on graphitic substrates.\" Applied Physics Letters 103.16 (2013): 163702.\n", "\n", "(2) Hale, G. M. and Querry, M. R. (1972). Optical constants of water in the 200-nm to 200-μm wavelength region. Appl. Opt., 12(3):555–563.\n", "\n", "(3) Johnson, P. B. and Christy, R. W. (1972). Optical constants of nobble metals. Phys. Rev. B, 12(6):4370–4379." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\n", "\n" ], "text/plain": [ "" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#Ignore this cell, It simply loads a style for the notebook.\n", "\n", "from IPython.core.display import HTML\n", "def css_styling():\n", " try:\n", " styles = open(\"styles/custom.css\", \"r\").read()\n", " return HTML(styles)\n", " except:\n", " pass\n", "css_styling()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.5.3" } }, "nbformat": 4, "nbformat_minor": 2 }