Skip to main content

OPINION article

Front. Mol. Biosci., 02 July 2019
Sec. Biological Modeling and Simulation
This article is part of the Research Topic Multiscale Modeling from Macromolecules to Cell: Opportunities and Challenges of Biomolecular Simulations View all 27 articles

Building Minimalist Models for Functionalized Metal Nanoparticles

  • 1Istituto Nanoscienze–CNR-NANO S3, Modena, Italy
  • 2Istituto Nanoscienze–CNR and NEST-Scuola Normale Superiore, Pisa, Italy

The Landscape of Coarse Grained NP Models

Metal nanoparticles (NPs) have been recently proposed for an increasing number of applications in nano-medicine (Vlamidis and Voliani, 2018) and nanotechnology (Chen et al., 2015). For instance, gold NPs (Alex and Tiwari, 2015) allow covalent versatile functionalization via thiol chemistry (Hakkinen, 2012) with different biomolecules or functional groups to selectively favor interactions with proteins or other specific components of the cell milieu. In particular, thiol-protected gold NPs functionalized with phenyl groups, Au25L18- (L = S(CH2)2Ph) were considered capable of interfering with protein aggregation, and therefore viewed as possible therapeutic agents against degenerative diseases due to amyloid fibrils accumulation (Brancolini et al., 2014, 2018; Marcinko et al., 2017; Torsten et al., 2018). The optimization of the size and decoration of the NP for therapy can benefit from computer simulations exploring aggregation in different environmental conditions (relative concentration, temperature, ionic strength). However, such extremely large time and size scale simulations call for the use of super-atomistic representations (low resolution or coarse grained—CG—models) (Brancolini and Tozzini, 2019).

A number of CG models for proteins are available (Seo et al., 2012), even minimalist ones, i.e. with single-bead per amino-acid resolution and implicit solvent (Di Fenza et al., 2009; Tozzini, 2010; Trovato and Tozzini, 2012; Trovato et al., 2013). Conversely, for the NPs, available CG models are rather sparse and diverse. The presence of the gold core suggests treating it at the meso-scale as a single spheroidal object (Vàcha et al., 2014), but the roughness of the surface (Radic et al., 2015), and the specificity of the chemical decoration (Tavanti et al., 2015a; Cantarutti et al., 2017) have fundamental roles in the interaction with proteins and must be treated at a higher resolution (Brancolini et al., 2015; Tavanti et al., 2015b; Charchar et al., 2016; Cardellini et al., 2019). Particular attention must be paid to the representation of hydrophobic character of the chemical groups and to the presence of possible net charges, whose medium- and long-range character, respectively, is the determinant of the macroscopic aggregation properties of the system. Implicit solvent requires the use of accurate screened potentials to account for the ionic strength. Finally, for the NP model to be compatible with the protein counterpart, both resolution and parameterization of the force field (FF) should be well matched.

While these prescriptions are followed inprevious literature in given models (Radic et al., 2015; Charchar et al., 2016), here we outline a general strategy to build models for NPs including all of them. In our view (Brancolini et al., 2018) these should contain the following ingredients: (1) Minimalism, i.e., including the minimum possible amount of degrees of freedom (DoF), and implicit solvent (2) Compatibility with the protein models (3) Transferability to different sizes and chemical decorations. Clearly, each of these characteristics involves one or more among the following actions: (i) choice of the model structure/topology, (ii) choice of the functional forms for the interactions, (iii) optimization of parameterization. (ii) and (iii) are complex tasks which have been addressed using a large number of different methodologies (Bauer et al., 2017; Lin et al., 2018; Brancolini et al., submitted). Particularly effective are usually combinations of bottom up and top-down strategy (Leonarski et al., 2013; Mereghetti et al., 2016) including both atomistic simulations and experimental data (Trovato and Tozzini, 2014) from different sources (e.g., structural, or thermodynamic). Here we focus on a general strategy to address (i) (Brancolini et al., 2018).

Rational Building of a Minimalist NP Model

The starting point is an atomistic structure of the functionalized NP (Figure 1A). The minimalism requirement suggests using a single large interacting center (“bead”) for the gold core, which is, in fact, a common feature to most of the NP models (Charchar et al., 2016; Shao and Hall, 2017). The chemical decoration is accounted for in several models by covering the central bead with smaller beads (Radic et al., 2015). The compatibility criterion can be satisfied choosing in specific ways the number and location of the decoration beads. For instance, when the functionalizing groups resemble in size and shape the side chains of amino-acids, this choice is rather straightforward: each of the functional group can be represented using the same representation of the protein amino-acids, i.e., 2–4 beads in MARTINI-like models (Seo et al., 2012), or a single bead for the minimalist models (Figure 1B). Remarkably, the model will include a number of DoF (Degrees of Freedom) proportional to the number of functional groups, i.e., will scale proportionally the surface of the NP, rather than to the volume.

FIGURE 1
www.frontiersin.org

Figure 1. Illustration of the phases of the model building, taking as a test case the Au25(S(CH2)2Ph)18. (A) Atomistic resolution, (B) coarse graining of functionalizing beads (C) analysis of the statistical structural data, clustering (D) evaluation of centroids of the space distribution, (E) topology and FF building (F) transfer of topology/structure to larger models. (G) Different symmetry of ligands are highlighted in the atomistic structure: 12 inner ligands (red circle) with sulfur bond to one core Au atom and one staple Au atom and 6 outer ligands (green circle) bound to two staple Au atoms. In principle, inner and outer beads can have different parameters, if needed to better balance the structure and electrostatics. (H) The electrostatic potential generated by the RESP atomistic charges using a representative configuration (UHBD isosurfaces drawn at +0.5 (red) and −0.5 (blue) kcal/e electrostatic potential). (I) The different types of CG beads reflecting the different symmetries of ligands are colored in different shades of green and assigned different charges [B2 = −0.230e (green), B3 = −0.128e (light green), the gold bead has charge B1 = +2.525e; the total net charge −1 is reproduced]. (J) CG electrostatic potential generated by the CG charges [UHBD isosurfaces drawn at +0.5 (red) and −0.5 (blue) kcal/e electrostatic potential].

An important point is how to choose the relative location of the decorating beads. Clearly, the thermal fluctuations of the group that they represent will determine the space distribution of the bead locations, which can be evaluated by means of atomistic simulations of small NPs (Maccari et al., 2014) (Figure 1C). The volume map build using this space distribution will form lobes, whose centroid and dispersion can be determined by clustering procedures (Arkhipov et al., 2006) (Figure 1D). This information can be used to build the starting location and topology of the model, and to parameterize the force field (FF) describing its internal dynamics (Figure 1E). Those parameters will then be transferred to larger NPs, once an average position of the functionalizing groups is determined, either from an atomistic model or from structural data (Figure 1F).

Distributing masses and effective charges among the beads is a non-trivial point. Considering masses, for instance, an obvious way would be to assign to each bead the sum of masses of their constituting elements. This, however, might not preserve the rotational inertia of the NP: in fact, being the total mass of the metal core concentrated in the center, it does not contribute, resulting in too small total rotational inertia. The problem can be solved by attributing larger masses to the peripheral beads. The proper balance of masses can be found by imposing that the total mass and the total rotational inertia correspond to that of the atomistic NP (Bauer et al., 2017).

The problem of charges is analogous: in this case an accurate charge distribution might be adjusted to reproduce the electrostatic potential, besides the net charge. The reference electrostatic potential can be generated from the RESP derived atomistic charges (Heaven et al., 2008), based on ab initio calculations (Figures 1G,H). Deriving the CG charges based on the atomistic components (Baker et al., 2001; Terakawa and Takada, 2014; McCullagh et al., 2016) results in effective charges depending on the bead type (gold or ligand) and symmetry (Figure 1I). The electrostatic potential generated by these can be compared with its atomistic counterpart, showing that the general shape of the iso-surfaces is preserved (Figure 1J): although of course the atomistic detail is lost, the CG model reproduce the global net prevalence of negative character (in blue), which however uncovers some positive areas (in red) for given directions, as in the atomistic case.

Summary and Perspectives: The Next Steps

In our opinion, the presented strategy includes all the crucial elements of an optimal low resolution model: the choice of the minimal possible resolution, compatibility between different levels of resolution, a parameterization including the specific coating present on the NP by means of superficial higher resolution interacting sites. The effective charges could be further optimized by directly adopting a RESP procedure for their fitting. This task and the model validation at different concentrations and ionic strengths vs. the aggregation tendency are currently in due course (Brancolini et al., submitted). The following steps will be the use of the model in combination with proteins models at the same CG level (minimalist), to verify their effective capability of preventing the amyloids aggregation. Furthermore, the strategy here outlined is extensible to larger NPs and different functionalization, which opens the possibility of in silico optimization of the NPs size and chemistry for therapeutic use.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

The authors acknowledge SEED project granted by CNR-Istituto Nanoscienze, Italy (GAE PUSEED04), titled LOPE-DeveLopment of a Coarse Grained MOdel forNanoparticle-Protein IntEractions for financial support.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to thank Dr. Hender Lopez for useful discussions. Oak Ridge National Laboratory is acknowledged for the supercomputing project CNMS2018-338, through the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. Department of Energy. Facilities of the National Energy Research Scientific Computing Center (NERSC), which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, are acknowledged.

References

Alex, S., and Tiwari, A. (2015). Functionalized gold nanoparticles: synthesis, properties and applications—a review. J. Nanosci. Nanotechnol. 15, 1869–1894. doi: 10.1166/jnn.2015.9718

PubMed Abstract | CrossRef Full Text | Google Scholar

Arkhipov, A., Freddolino, P. L., and Schulten, K. (2006). Stability and dynamics of virus capsids described by coarse-grained modeling. Structure 14, 1767–1777. doi: 10.1016/j.str.2006.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Baker, N. A., Sept, D., Joseph, S., Holst, M. J., and McCammon, J. A. (2001). Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl. Acad. Sci. U.S.A. 98, 10037–10041. doi: 10.1073/pnas.181342398

PubMed Abstract | CrossRef Full Text | Google Scholar

Bauer, G., Gribova, N., Lange, A., Holm, C., and Gross, J. (2017). Three-body effects in triplets of capped gold nanocrystals. Mol. Phys. 115, 1031–1040. doi: 10.1080/00268976.2016.1213909

CrossRef Full Text | Google Scholar

Brancolini, G., Corazza, A., Vuano, M., Fogolari, F., Mimmi, M. C., Bellotti, V., et al. (2015). Probing the influence of citrate-capped gold nanoparticles on an amyloidogenic protein. ACS Nano 9, 2600–2613. doi: 10.1021/nn506161j

PubMed Abstract | CrossRef Full Text | Google Scholar

Brancolini, G., Maschio, M. C., Cantarutti, C., Corazza, A., Fogolari, F., Corni, S., et al. (2018). Citrate stabilized gold nanoparticles interfere with amyloid fibril formation: D76N and ΔN6 variants. Nanoscale 10, 4793–4806. doi: 10.1039/c7nr06808e

PubMed Abstract | CrossRef Full Text | Google Scholar

Brancolini, G., Toroz, D., and Corni, S. (2014). Can small hydrophobic gold nanoparticles inhibit β2-microglobulin fibrillation? Nanoscale 6, 7903–7911. doi: 10.1039/C4NR01514B

CrossRef Full Text | Google Scholar

Brancolini, G., and Tozzini, V. (2019). Multiscale modeling of proteins interaction with functionalized nanoparticles. COCIS 41, 66–73. doi: 10.1016/j.cocis.2018.12.001

CrossRef Full Text

Cantarutti, C., Raimondi, S., Brancolini, G., Corazza, A., Giorgetti, S., Ballico, M., et al. (2017). Citrate-stabilized gold nanoparticles hinder fibrillogenesis of a pathological variant of β2-microglobulin. Nanoscale 9, 3941–3951. doi: 10.1039/c6nr09362k

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardellini, A., Alberghini, M., Rajan, A. G., Misra, R. P., Blankschtein, D., and Asinari, P. (2019). Multi-scale approach for modeling stability, aggregation, and network formation of nanoparticles suspended in aqueous solutions. Nanoscale 11, 3979–3992. doi: 10.1039/C8NR08782B

PubMed Abstract | CrossRef Full Text | Google Scholar

Charchar, P., Christofferson, A. J., Todorova, N., and Yarovsky, I. (2016). Understanding and designing the goldbio interface: insights from simulations. Small 12, 2395–2418. doi: 10.1002/smll.201503585

CrossRef Full Text | Google Scholar

Chen, L. Y., Wang, C. W., Yuan, Z., and Chang, H. T. (2015). Fluorescent gold nanoclusters: recent advances in sensing and imaging. Anal. Chem. 6, 216–229. doi: 10.1021/ac503636j

CrossRef Full Text | Google Scholar

Di Fenza, A., Rocchia, A. W., and Tozzini, V. (2009). Complexes of HIV-1 integrase with HAT proteins: Multiscale models, dynamics, and hypotheses on allosteric sites of inhibition. Proteins Struct. Funct. Bioinform. 76, 946–958. doi: 10.1002/prot.22399

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakkinen, H. (2012). The gold-sulfur interface at the nanoscale. Nat. Chem. 4, 443–455. doi: 10.1038/nchem.1352

PubMed Abstract | CrossRef Full Text

Heaven, M. W., Dass, A., White, P. S., Holt, K. M., and Murray, R. W. (2008). Crystal Structure of the Gold Nanoparticle [N(C8H17)4][Au25(SCH2CH2Ph)18]. J. Am. Chem. Soc. 130, 3754–3755. doi: 10.1021/ja800561b

PubMed Abstract | CrossRef Full Text | Google Scholar

Leonarski, F., Trovato, F., Tozzini, V., Les, A., and Trylska, J. (2013). Evolutionary algorithm in the optimization of a coarse-grained force field. J. Chem. Theor. Comput. 9, 4874–4889. doi: 10.1021/ct4005036

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y., Charchar, P., Christofferson, A. J., Thomas, M. R., Todorova, N., Mazo, M. M., et al. (2018). Surface dynamics and ligand–core interactions of quantum sized photoluminescent gold nanoclusters. J. Am. Chem. Soc. 140, 18217–18226. doi: 10.1021/jacs.8b04436

PubMed Abstract | CrossRef Full Text | Google Scholar

Maccari, G., Spampinato, G. L. B., and Tozzini, V. (2014). SecStAnT: Secondary Structure Analysis Tool for data selection, statistics and models building. Bioinformatics 30, 668–674. doi: 10.1093/bioinformatics/btt586

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcinko, T. M., Dong, J., LeBlanc, R., Daborowski, K. V., and Vachet, R. W. (2017). Small moleculemediatedInhibition of β-2-microglobulin amyloid fibril formation. J. Biol. Chem. 292, 10630–10638. doi: 10.1074/jbc.M116.774083

PubMed Abstract | CrossRef Full Text | Google Scholar

McCullagh, P., Lake, P. T., and Mccullagh, M. (2016). Deriving coarse-grained charges from all-atom systems: an analytic solution. J. Chem. Theor. Comput. 12, 4390–4399. doi: 10.1021/acs.jctc.6b00507

PubMed Abstract | CrossRef Full Text | Google Scholar

Mereghetti, P., Maccari, G., Spampinato, G. L. B., and Tozzini, V. (2016). Optimization of analytical potentials for coarse-grained biopolymers models. J. Phys. Chem. B 120, 8571–8579. doi: 10.1021/acs.jpcb.6b02555

CrossRef Full Text | Google Scholar

Radic, S., Davis, T. P., Keb, P. C., and Ding, F. (2015). Contrasting effects of nanoparticle-protein attraction on amyloid aggregation. RSC Adv. 5, 105489–105498. doi: 10.1039/C5RA20182A

PubMed Abstract | CrossRef Full Text | Google Scholar

Seo, M., Rauscher, S., Pomes, R., and Tieleman, D. P. (2012). Improving internal peptide dynamics in the coarse-grained MARTINI model: toward large-scale simulations of amyloid- and elastin-like peptides. J. Chem. Theor. Comput. 8, 1774–1785. doi: 10.1021/ct200876v

PubMed Abstract | CrossRef Full Text | Google Scholar

Shao, Q., and Hall, C. K. (2017). Allosteric effects of gold nanoparticles on human serum albumin. Nanoscale 9, 380–390. doi: 10.1039/c6nr07665c

PubMed Abstract | CrossRef Full Text | Google Scholar

Tavanti, F., Pedone, A., and Menziani, M. C. (2015a). A closer look into the ubiquitin corona on gold nanoparticles by computational studies. New J. Chem. 39, 2474–2482. doi: 10.1039/C4NJ01752H

CrossRef Full Text | Google Scholar

Tavanti, F., Pedone, A., and Menziani, M. C. (2015b). Competitive binding of proteins to gold nanoparticles disclosed by molecular dynamics simulations. J. Phys. Chem. C 119, 22172–22180. doi: 10.1021/acs.jpcc.5b05796

CrossRef Full Text | Google Scholar

Terakawa, T., and Takada, S. (2014). RESPAC: method to determine partial charges in coarse-grained protein model and its application to DNA-binding proteins. J. Chem. Theor. Comput. 10, 711–721. doi: 10.1021/ct4007162

PubMed Abstract | CrossRef Full Text | Google Scholar

Torsten, J., Anika, G., Clemens, K., Lisandra, L. M., Herre, J. R., and Bernd, A. (2018). Impact of nanoparticles on amyloid peptide and protein aggregation: a review with a focus on gold nanoparticles. Nanoscale 10, 20894–20913. doi: 10.1039/c8nr04506b

CrossRef Full Text | Google Scholar

Tozzini, V. (2010). Minimalist models for proteins: a comparative analysis. Q. Rev. Biophys. 43, 333–371. doi: 10.1017/S0033583510000132

PubMed Abstract | CrossRef Full Text | Google Scholar

Trovato, F., Nifosì, R, Di Fenza, A., and Tozzini, V. (2013). A minimalist model of protein diffusion and interactions: the green fluorescent protein within the cytoplasm. Macromolecules 46, 8311–8322. doi: 10.1021/ma401843h

CrossRef Full Text | Google Scholar

Trovato, F., and Tozzini, V. (2012). Minimalist models for biopolymers: Open problems, latest advances and perspectives. AIP Conf. Proc. 1456, 187–200. doi: 10.1063/1.4730659

CrossRef Full Text | Google Scholar

Trovato, F., and Tozzini, V. (2014). Diffusion within the cytoplasm: a mesoscale model of interacting macromolecules. Biophys. J. 107, 2579–2591. doi: 10.1016/j.bpj.2014.09.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Vàcha, R., Linse, S., and Lund, M. (2014). Surface effects on aggregation kinetics of amyloidogenic peptides. J. Am. Chem. Soc. 136, 11776–11782. doi: 10.1021/ja505502e

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlamidis, Y., and Voliani, V. (2018). Bringing again noble metal nanoparticles to the forefront of cancer therapy. Front. Bioeng. Biotechnol. 6, 143–148. doi: 10.3389/fbioe.2018.00143

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: coarse grained models, molecular dynamics, brownian dynamics, multiscale simulations, gold nanocrystal, macromolecules aggregation

Citation: Brancolini G and Tozzini V (2019) Building Minimalist Models for Functionalized Metal Nanoparticles. Front. Mol. Biosci. 6:50. doi: 10.3389/fmolb.2019.00050

Received: 11 March 2019; Accepted: 17 June 2019;
Published: 02 July 2019.

Edited by:

Edina Rosta, King's College London, United Kingdom

Reviewed by:

Bart De Nijs, University of Cambridge, United States

Copyright © 2019 Brancolini and Tozzini. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Giorgia Brancolini, giorgia.brancolini@nano.cnr.it

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.