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

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 piecewise 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 as in 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 toward 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.

Modeling of the Electrostatic Interaction and Catalytic Activity of [NiFe] Hydrogenases on a Planar Electrode

Hydrogenases are a group of enzymes that have caught the interest of researchers in renewable energies, due to their ability to catalyze the redox reaction of hydrogen. The exploitation of hydrogenases in electrochemical devices requires their immobilization on the surface of suitable electrodes, such as graphite. The orientation of the enzyme on the electrode is important to ensure a good flux of electrons to the catalytic center, through an array of iron-sulfur clusters. Here we present a computational approach to determine the possible orientations of a [NiFe] hydrogenase (PDB 1e3d) on a planar electrode, as a function of pH, salinity, and electrode potential. The calculations are based on the solution of the linearized Poisson-Boltzmann equation, using the PyGBe software. The results reveal that electrostatic interactions do not truly immobilize the enzyme on the surface of the electrode, but there is instead a dynamic equilibrium between different orientations. Nonetheless, after averaging over all thermally accessible orientations, we find significant differences related to the solution's salinity and pH, while the effect of the electrode potential is relatively weak. We also combine models for the protein adsoption-desorption equilibria and for the electron transfer between the proteins and the electrode to arrive at a prediction of the electrode's activity as a function of the enzyme concentration.

Quantitative electrostatic force tomography for virus capsids in interaction with an approaching nanoscale probe

Electrostatic interactions are crucial for the assembly, disassembly and stability of proteinaceous viral capsids. Moreover, at the molecular scale, elucidating the organization and structure of the capsid proteins in response to an approaching nanoprobe is a major challenge in biomacromolecular research. Here, we report on a generalized electrostatic model, based on the Poisson-Boltzmann equation, that quantifies the subnanometric electrostatic interactions between an AFM tip and a proteinaceous capsid from molecular snapshots. This allows us to describe the contributions of specific amino acids and atoms to the interaction force. We show validation results in terms of total electrostatic forces with previous semi-empirical generalized models at available length scales (d > 1 nm). Then, we studied the interaction of the Zika capsid with conical and spherical AFM tips in a tomography-type analysis to identify the most important residues and atoms, showing the localized nature of the interaction. This method can be employed for the interpretation of force microscopy experiments in fundamental virological characterization and in diverse nanomedicine applications, where specific regions of the protein cages are aimed to electrostatically interact with molecular sized functionalized inhibitors, or tailoring protein-cage functional properties for nucleic acid delivery.

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 a spontaneous process that leads to a rather compact layer of randomly oriented molecules. However, controlling such orientation is critically important for the development of catalytic surfaces. In this regard, the use of electric fields is one of the most promising alternatives. Our work is motivated by experimental observations that show important differences in catalytic activity of a trypsin-covered surface, which depended on the applied potential during the adsorption. Even though adsorption results from the combination of several processes, we were able to determine that (under the selected conditions) mean-field electrostatics play a dominant role, determining the orientation and yielding a difference in catalytic activity. We simulated the electrostatic potential numerically, using an implicit-solvent model based on the linearized Poisson-Boltzmann equation. This was implemented in an extension of the code PyGBe that included an external electric field, and rendered the electrostatic component of the solvation free energy. Our model (extensions available at the Github repository) allowed estimating the overall affinity of the protein with the surface, and their most likely orientation as a function of the potential applied. Our results show that the active sites of trypsin are, on average, more exposed when the electric field is negative, which agrees with the experimental results of catalytic activity, and confirm the premise that electrostatic interactions can be used to control the orientation of adsorbed proteins.

Dielectric Spectroscopy can Predict the Effect of External AC Fields on the Dynamic Adsorption of Lysozyme

This report describes the application of dielectric spectroscopy as a simple and fast way to guide protein adsorption experiments. Specifically, the polarization behavior of a layer of adsorbed lysozyme was investigated using a triangular‐wave signal with frequencies varying from 0.5 to 2 Hz. The basic experiment, which can be performed in less than 5 min and with a single sample, not only allowed confirming the susceptibility of the selected protein towards the electric signal but also identified that this protein would respond more efficiently to signals with lower frequencies. To verify the validity of these observations, the adsorption behavior of lysozyme onto optically transparent carbon electrodes was also investigated under the influence of an applied alternating potential. In these experiments, the applied signal was defined by a sinusoidal wave with an amplitude of 100 mV and superimposed to +800 mV (applied as a working potential) and varying the frequency in the 0.1–10000 Hz range. The experimental data showed that the greatest adsorbed amounts of lysozyme were obtained at the lowest tested frequencies (0.1–1.0 Hz), results that are in line with the corresponding dielectric features of the protein.

Design, Fabrication, and Characterisation of a Label-Free Nanosensor for Bioapplications

In this paper, we present a hybrid semiconductor structure for biosensing applications that features the co-integration of nanoelectromechanical systems with the well-known metal oxide semiconductor technology. The proposed structure features an MOSFET as a readout element, and a doubly clamped beam that is isolated from the substrate by a thin air gap, as well as by a tunnel oxide layer. The beam structure is functionalised by a thin layer of biotargets, and the main aim is to detect a particular set of biomolecules, such as enzymes, bacteria, viruses, and DNA/RNA chains, among others. In here, a three-dimensional finite element analysis is performed in order to study the behaviour of the functionalised, doubly clamped beam. Preliminary results for the fabrication and characterisation processes show good agreement between the simulated and measured characteristics.

Adsorption of lysozyme into a charged confining pore

Several applications arise from the confinement of proteins on surfaces because their stability and biological activity are enhanced. It is also known that the way in which a protein adsorbs on the surface is important for its biological function since its active sites should not be obstructed. In this study, the adsorption properties of hen egg-white lysozyme, HEWL, into a negatively charged silica pore is examined by employing a coarse-grained model and constant-pH Monte Carlo simulations. The role of electrostatic interactions is taken into account via including the Debye-Hückel potentials into the Cα structure-based model. We evaluate the effects of pH, salt concentration, and pore radius on the protein preferential orientation and spatial distribution of its residues regarding the pore surface. By mapping the residues that stay closer to the pore surface, we find that the increase of pH leads to orientational changes of the adsorbed protein when the solution pH gets closer to the HEWL isoelectric point. Under these conditions, the pK_a shift of these important residues caused by the adsorption into the charged confining surface results in a HEWL charge distribution that stabilizes the adsorption in the observed protein orientation. We compare our observations to the results of the pK_a shift for HEWL available in the literature and to some experimental data.

Tethering of the IgG1 Antibody to Amorphous Silica for Immunosensor Development: A Molecular Dynamics Study

A key factor for improving the sensitivity and performance of immunosensors based on mechanical-plasmonic methods is the orientation of the antibody proteins immobilized on the inorganic surface. Although experimental techniques fail to determine surface phenomena at the molecular level, modern simulations open the possibility for improving our understanding of protein-surface interactions. In this work, replica exchange molecular dynamics (REMD) simulations have been used to model the IgG1 protein tethered onto the amorphous silica surface by considering a united-atom model and a relatively large system (2500 nm² surface). Additional molecular dynamics (MD) simulations have been conducted to derive an atomistic model for the amorphous silica surface using the cristobalite crystal structure as a starting point and to examine the structure of the free IgG1 antibody in the solution for comparison when immobilized. Analyses of the trajectories obtained for the tethered IgG1, which was sampled considering 32 different temperatures, have been used to define the geometry of the protein with respect to the inorganic surface. The tilt angle of the protein with respect to the surface plane increases with temperature, the most populated values being 24, 66, and 87° at the lowest (250 K), room (298 K), and the highest (380 K) temperatures. This variation indicates that the importance of protein-surface interactions decreases with increasing temperature. The influence of the surface on the structure of the antibody is very significant in the constant region, which is directly involved in the tethering process, while it is relatively unimportant for the antigen-binding fragments, which are farthest from the surface. These results are expected to contribute to the development of improved mechanical-plasmonic sensor microarrays in the near future.

Computational nanoplasmonics in the quasistatic limit for biosensing applications

The phenomenon of localized surface plasmon resonance (LSPR) provides high sensitivity in detecting biomolecules through shifts in resonance frequency when a target is present. Computational studies in this field have used the full Maxwell equations with simplified models of a sensor-analyte system, or they neglected the analyte altogether. In the long-wavelength limit, one can simplify the theory via an electrostatics approximation while adding geometrical detail in the sensor and analytes (at moderate computational cost). This work uses the latter approach, 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(NlogN) computational complexity. These features allow PyGBe to handle problems with half a million boundary elements or more. In this work, we demonstrate the suitability of PyGBe, extended to compute LSPR response in the electrostatic limit, for biosensing applications. 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 redshift 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 recreate the figures. PyGBe is open source under a permissive license and openly developed. Documentation is available at http://pygbe.github.io/pygbe/docs/.

Orienting proteins by nanostructured surfaces: evidence of a curvature-driven geometrical resonance

Experimental and theoretical reports have shown that nanostructured surfaces have a dramatic effect on the amount of protein adsorbed and the conformational state and, in turn, on the performances of the related devices in tissue engineering strategies. Here we report an innovative method to prepare silica-based nanostructured surfaces with a reproducible, well-defined local curvature, consisting of ordered hexagonally packed arrays of curved hemispheres, from nanoparticles of different diameters (respectively 147 nm, 235 nm and 403 nm). The nanostructured surfaces have been made chemically homogeneous by partially embedding silica nanoparticles in poly(hydroxymethylsiloxane) films, further modified by means of UV-O3 treatments. This paper has been focused on the experimental and theoretical study of laminin, taken as a model protein, to study the nanocurvature effects on the protein configuration at nanostructured surfaces. A simple model, based on the interplay of electrostatic interactions between the charged terminal domains of laminin and the nanocurved charged surfaces, closely reproduces the experimental findings. In particular, the model suggests that nanocurvature drives the orientation of rigid proteins by means of a "geometrical resonance" effect, involving the matching of dimensions, charge distribution and spatial arrangement of both adsorbed molecules and adsorbent nanostructures. Overall, the results pave the way to unravel the nanostructured surface effects on the intra- and inter-molecular organization processes of proteins.

Site-selective orientated immobilization of antibodies and conjugates for immunodiagnostics development

Immobilized antibody systems are the key to develop efficient diagnostics and separations tools. In the last decade, developments in the field of biomolecular engineering and crosslinker chemistry have greatly influenced the development of this field. With all these new approaches at our disposal, several new immobilization methods have been created to address the main challenges associated with immobilized antibodies. Few of these challenges that we have discussed in this review are mainly associated to the site-specific immobilization, appropriate orientation, and activity retention. We have discussed the effect of antibody immobilization approaches on the parameters on the performance of an immunoassay.