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, Au25 (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
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.
Statements
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.
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.
Conflict of interest
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.
References
1
AlexS.TiwariA. (2015). Functionalized gold nanoparticles: synthesis, properties and applications—a review. J. Nanosci. Nanotechnol.15, 1869–1894. 10.1166/jnn.2015.9718
2
ArkhipovA.FreddolinoP. L.SchultenK. (2006). Stability and dynamics of virus capsids described by coarse-grained modeling. Structure14, 1767–1777. 10.1016/j.str.2006.10.003
3
BakerN. A.SeptD.JosephS.HolstM. J.McCammonJ. A. (2001). Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl. Acad. Sci. U.S.A. 98, 10037–10041. 10.1073/pnas.181342398
4
BauerG.GribovaN.LangeA.HolmC.GrossJ. (2017). Three-body effects in triplets of capped gold nanocrystals. Mol. Phys.115, 1031–1040. 10.1080/00268976.2016.1213909
5
BrancoliniG.CorazzaA.VuanoM.FogolariF.MimmiM. C.BellottiV.et al. (2015). Probing the influence of citrate-capped gold nanoparticles on an amyloidogenic protein. ACS Nano9, 2600–2613. 10.1021/nn506161j
6
BrancoliniG.MaschioM. C.CantaruttiC.CorazzaA.FogolariF.CorniS.et al. (2018). Citrate stabilized gold nanoparticles interfere with amyloid fibril formation: D76N and ΔN6 variants. Nanoscale10, 4793–4806. 10.1039/c7nr06808e
7
BrancoliniG.TorozD.CorniS. (2014). Can small hydrophobic gold nanoparticles inhibit β2-microglobulin fibrillation?Nanoscale6, 7903–7911. 10.1039/C4NR01514B
8
BrancoliniG.TozziniV. (2019). Multiscale modeling of proteins interaction with functionalized nanoparticles. COCIS41, 66–73. 10.1016/j.cocis.2018.12.001
9
CantaruttiC.RaimondiS.BrancoliniG.CorazzaA.GiorgettiS.BallicoM.et al. (2017). Citrate-stabilized gold nanoparticles hinder fibrillogenesis of a pathological variant of β2-microglobulin. Nanoscale9, 3941–3951. 10.1039/c6nr09362k
10
CardelliniA.AlberghiniM.RajanA. G.MisraR. P.BlankschteinD.AsinariP. (2019). Multi-scale approach for modeling stability, aggregation, and network formation of nanoparticles suspended in aqueous solutions. Nanoscale11, 3979–3992. 10.1039/C8NR08782B
11
CharcharP.ChristoffersonA. J.TodorovaN.YarovskyI. (2016). Understanding and designing the goldbio interface: insights from simulations. Small12, 2395–2418. 10.1002/smll.201503585
12
ChenL. Y.WangC. W.YuanZ.ChangH. T. (2015). Fluorescent gold nanoclusters: recent advances in sensing and imaging. Anal. Chem.6, 216–229. 10.1021/ac503636j
13
Di FenzaA.RocchiaA. W.TozziniV. (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. 10.1002/prot.22399
14
HakkinenH. (2012). The gold-sulfur interface at the nanoscale. Nat. Chem.4, 443–455. 10.1038/nchem.1352
15
HeavenM. W.DassA.WhiteP. S.HoltK. M.MurrayR. W. (2008). Crystal Structure of the Gold Nanoparticle [N(C8H17)4][Au25(SCH2CH2Ph)18]. J. Am. Chem. Soc.130, 3754–3755. 10.1021/ja800561b
16
LeonarskiF.TrovatoF.TozziniV.LesA.TrylskaJ. (2013). Evolutionary algorithm in the optimization of a coarse-grained force field. J. Chem. Theor. Comput.9, 4874–4889. 10.1021/ct4005036
17
LinY.CharcharP.ChristoffersonA. J.ThomasM. R.TodorovaN.MazoM. M.et al. (2018). Surface dynamics and ligand–core interactions of quantum sized photoluminescent gold nanoclusters. J. Am. Chem. Soc.140, 18217–18226. 10.1021/jacs.8b04436
18
MaccariG.SpampinatoG. L. B.TozziniV. (2014). SecStAnT: Secondary Structure Analysis Tool for data selection, statistics and models building. Bioinformatics30, 668–674. 10.1093/bioinformatics/btt586
19
MarcinkoT. M.DongJ.LeBlancR.DaborowskiK. V.VachetR. W. (2017). Small moleculemediatedInhibition of β-2-microglobulin amyloid fibril formation. J. Biol. Chem.292, 10630–10638. 10.1074/jbc.M116.774083
20
McCullaghP.LakeP. T.MccullaghM. (2016). Deriving coarse-grained charges from all-atom systems: an analytic solution. J. Chem. Theor. Comput. 12, 4390–4399. 10.1021/acs.jctc.6b00507
21
MereghettiP.MaccariG.SpampinatoG. L. B.TozziniV. (2016). Optimization of analytical potentials for coarse-grained biopolymers models. J. Phys. Chem. B120, 8571–8579. 10.1021/acs.jpcb.6b02555
22
RadicS.DavisT. P.KebP. C.DingF. (2015). Contrasting effects of nanoparticle-protein attraction on amyloid aggregation. RSC Adv. 5, 105489–105498. 10.1039/C5RA20182A
23
SeoM.RauscherS.PomesR.TielemanD. 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. 10.1021/ct200876v
24
ShaoQ.HallC. K. (2017). Allosteric effects of gold nanoparticles on human serum albumin. Nanoscale9, 380–390. 10.1039/c6nr07665c
25
TavantiF.PedoneA.MenzianiM. C. (2015a). A closer look into the ubiquitin corona on gold nanoparticles by computational studies. New J. Chem.39, 2474–2482. 10.1039/C4NJ01752H
26
TavantiF.PedoneA.MenzianiM. C. (2015b). Competitive binding of proteins to gold nanoparticles disclosed by molecular dynamics simulations. J. Phys. Chem. C119, 22172–22180. 10.1021/acs.jpcc.5b05796
27
TerakawaT.TakadaS. (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. 10.1021/ct4007162
28
TorstenJ.AnikaG.ClemensK.LisandraL. M.HerreJ. R.BerndA. (2018). Impact of nanoparticles on amyloid peptide and protein aggregation: a review with a focus on gold nanoparticles. Nanoscale10, 20894–20913. 10.1039/c8nr04506b
29
TozziniV. (2010). Minimalist models for proteins: a comparative analysis. Q. Rev. Biophys.43, 333–371. 10.1017/S0033583510000132
30
TrovatoF.NifosìRDi FenzaA.TozziniV. (2013). A minimalist model of protein diffusion and interactions: the green fluorescent protein within the cytoplasm. Macromolecules46, 8311–8322. 10.1021/ma401843h
31
TrovatoF.TozziniV. (2012). Minimalist models for biopolymers: Open problems, latest advances and perspectives. AIP Conf. Proc.1456, 187–200. 10.1063/1.4730659
32
TrovatoF.TozziniV. (2014). Diffusion within the cytoplasm: a mesoscale model of interacting macromolecules. Biophys. J.107, 2579–2591. 10.1016/j.bpj.2014.09.043
33
VàchaR.LinseS.LundM. (2014). Surface effects on aggregation kinetics of amyloidogenic peptides. J. Am. Chem. Soc. 136, 11776–11782. 10.1021/ja505502e
34
VlamidisY.VolianiV. (2018). Bringing again noble metal nanoparticles to the forefront of cancer therapy. Front. Bioeng. Biotechnol. 6, 143–148. 10.3389/fbioe.2018.00143
Summary
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
Volume
6 - 2019
Edited by
Edina Rosta, King's College London, United Kingdom
Reviewed by
Bart De Nijs, University of Cambridge, United States
Updates
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
This article was submitted to Biological Modeling and Simulation, a section of the journal Frontiers in Molecular Biosciences
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.