Nano-Detoxification of Organophosphate Agents by PAMAM Derivatives

PAMAM G5-Folate derivative (G5-FA) Briefly, 104 mg (0.23 mmol) of FA reacted with 150 mg (0.96 mmol) of EDC·HCl and 150 of HOBt (1.11 mmol) in a mixture of 2.7 mL of dry DMF and 1.0 mL of dry DMSO under a nitrogen atmosphere for 1 h. The reaction mixture was slowly added dropwise to a solution of 40 mg (1.38 × 10−3 mmol) of G5 in 3.0 mL of deionized (DI) water. The reaction mixture was vigorously stirred for 72 h. The final functionalized PAMAM was purified through dialysis (using water) membranes with cut-offs of 500 Da to take off the excess of FA. After lyophilization, the weight of G5-FA was 43 mg, affording 90% of yield.


Introduction
Worldwide, the most employed insecticides for indoor use and agricultural purposes belong to carbamates, pyrethroids or organophosphates (OP). 1 As pesticides, OP have been used since the end of the Second World War.But OP have also been used as lethal nerve and chemical warfare agents (CWAs), such as VX, soman and sarin. 2 OP pesticides and CWAs have similar chemical structures (Table 1) and therefore have similar mechanisms of toxicity. 3This is based on its capacity to bind the active site of acetylcholinesterase (AChE), resulting on an accumulation of acetylcholine at the neuron synapses. 4,5The neurotransmitter accumulation leads to symptoms related to the autonomous nervous system (abdominal cramps, nausea, diarrhea and salivation) and the central nervous system (dizziness, tremor, anxiety and confusion). 6,7Due to this, one of the current greatest threats in modern world is the possibility of OPs use as chemical weapons by terrorist groups. 8he World Health Organization estimates that around 3 million people per year are poisoned by pesticides, most of them by OP pesticides, resulting in around 200,000 deaths. 9oisoning occurs by accidental spillage, suicide attempts and occupational hazards involving workers and farmers. 3anotechnology has become an area of great interest in nanomedicine and life science applications as drug delivery, gene delivery, or nanodiagnosis systems, [10][11][12] but there are few studies that use dendrimers as nano-detoxification agents of harmful compounds in the body.For this reason, OPs' capture, through the use of nanocompounds, is a growing and interesting field of development, due to its capacity for selective capture of toxic compounds in vivo.The aim of this study is to evaluate the ability of dendrimers to capture OPs in a selective manner.Polyamidoamine (PAMAM) dendrimers are a class of artificial macromolecules that have shown excellent potential in biomedical applications due to their unique physical and chemical properties.Among these properties, PAMAM show a precise control on molecular size, shape and density, high solubility, a large number of surface functional groups and hydrophobic cavities that can capture molecules of interest. 13PAMAM are biocompatible, non-immunogenic, water-soluble, and possess terminal modifiable amine functional groups for binding various targeting or guest molecules. 14PAMAM are large, highly branched macromolecules that are synthesized from an ethylenediamine core by a successive addition of methyl acrylate and ethylenediamine.The number of terminal amino groups increase as the generation of PAMAM increases. 15][18] Its easy acquisition and handling can pose a great threat if used as CW, and the potential attack employing these agents gives a strong reason for continuous research on the development of more effective antidotes against them.
Structure-based molecular characterization using semi-empirical quantum mechanical studies 19,20 can be used to accurately understand the properties that govern the interactions between molecular systems.However, the application of these computational chemistry methods to analyze nanostructured systems requires its suitable adaptation and implementation.][24] In this work, a nanoinformatics methodology was used to calculate the interaction energy between molecule(1)molecule(2) pairs. 25,26The first molecule represents the monomeric unit of a specific dendrimer attached to a functional group, and the second molecule represents the OP compound.To obtain an accurate prediction of the interaction energies, the nanoinformatics methodology implements a Metropolis Monte Carlo algorithm 27 that performs a random sampling through Euler angles 25 to generate thousands of different molecular conformations (for a specific pair), then for each new conformation, its interaction energy is calculated using a semi-empirical quantum mechanical method. 19,20The long list of interaction energy values can be used for in silico analysis of the affinity between OP agents and the functionalized PAMAM.Thus, in silico techniques are perfect tools for this task since they can quickly provide structures of several new lethal compounds for further experimental work. 28he aim of the present work is to incorporate active biomolecules such as folic acid, coumarin, lysine, arginine and asparagine as new end-groups in the structures of

MALDI analysis
To confirm the molecular weight of surface modified PAMAM, mass spectra analyses of the dendrimers were performed on matrix-assisted laser desorption ionizationtime of flight (MALDI-TOF) equipment (Bruker Autoflex) with a pulsed nitrogen laser (337 nm), operating in positive ion reflector mode, using 19 kV acceleration voltage and 2,5-dihydroxybenzoic acid (DHB) as matrix.

Synthesis of PAMAM derivatives
Six different PAMAM formulations varying size and surface charges were synthetized as described below.
PAMAM-G5-Folate derivative (G5-FA): conjugation of G5 with FA was carried out by a condensation between the γ-carboxyl group of FA and the primary amino group of PAMAM, according to the previously reported method 12,29 [available as Supplementary Information (SI)].
PAMAM-Coumarine derivative (G5-Cou): conjugation of G5 with Cou was carried out by a condensation between the carboxyl group of Cou and the primary amino group of PAMAM, according to the previously reported method 12,30 (available as SI).
PAMAM G4-Arginine(Tos)-OH (G4-Arg): conjugation of G4 with Boc-Arg(Tos)-OH was carried out by a condensation between the carboxyl group of arginine and the primary amino group of PAMAM, according to the previously reported method 12,30 (available as SI).
PAMAM G4-Asn(Trt)-OH) (G4-Asn): conjugation of G4 with Boc-Asn(Trt)-OH was carried out by a condensation between the carboxyl group of Boc-asn(Trt)-OH and the primary amino group of PAMAM, according to the previously reported method 12,30 (available as SI).

Affinity assays
For each PAMAM and functionalized PAMAM, experiments were carried out to determine the extent of binding and fractional binding of MMP in methanol solutions.Briefly, functionalized PAMAM adjusted at 5.83 × 10 -3 mmol L -1 was mixed in a 1:1 v/v ratio with a 0.07 mmol L -1 (10 ppm) of MMP to give a final MMP/PAMAM mol L -1 ratio of 12:1.Assays were performed at pH 4.5-5.5.The samples were mixed for 45 min with the PAMAM derivatives at constant room temperature (25 °C), and then centrifuged at 10000 rpm for 10 min.The concentrations of MMP in supernatants separated from precipitated were analyzed by HPLC.The adsorption efficiency of each OP compound by PAMAM were evaluated by determining the percentage decrease in the absorbance at each specific maximum absorbance wavelength using the following equation: where A0 is the initial absorbance at specific wavelength and A is the final absorbance at the same wavelength of each OP. 31

Chromatographic method
An Agilent 1260 Infinity HPLC system (USA) with a photodiode array detector and a C-18 Lichrospher 100 RP-18 (250 mm × 4 mm i.d.× 5 μm) column was used for the analysis of column eluents.
A 40 μL solution of the studied supernatant was injected into the HPLC apparatus.An isocratic elution with methanol-water mixture (70:30, v/v) at a flow rate of 1.5 mL min -1 was used as mobile phase.The eluents were monitored at 230 nm by absorbance detection.For data processing, the chromatographic software ChemStation was used.The chromatographic determination was performed at room temperature.Reproducibility of the method was determined for individual standards.Three replicates per pesticide and concentration level were prepared.

Calibration procedures
AZM and MMP standards were prepared in methanol.The concentrations ranged from 10, 5.0, 2.5 and 1.25 ppm.Linear calibration curves were obtained by plotting the peak areas of the individual chemicals as a function of the concentration using the GraphPad Prism program for Windows (GraphPad Software, Inc., San Diego, CA, U.S.A.).

Statistical analysis
All the experiments were carried out in triplicate and the student t-test was used to calculate the statistical differences between the experimental compound and the negative control.A p-value of < 0.05 was considered significant.

Theoretical platform
Two OP compounds were chosen for detailed molecular dynamic (MD) investigations, MMP and AZM.G4-Asn 10 , which presented the highest binding capacity for MMP at pH 7, was chosen for build systems and was parameterized according to the rules for all-atom CHARMM27. 32This force field was not included in conventional force field due to the dendritic structures.Charges, missing bond, angle, dihedral parameters of the PAMAM core, as well as intermediate dendrons and the outer surface were estimated from similar terms within the force field.An atomistic model for G4-Asn 10 was built at pH 7 containing 64 protonated amines in total (54 terminal amines from PAMAM and 10 amines from Asn end groups).The 10 units of Asn functional groups have been homogeneously distributed in 64 end groups.The molecular parameters for MMP and AZM were obtained from SwissParam 33 online server (http://www.swissparam.ch/)by using the Merck Molecular Force Field (MMFF) at pH 7 (net charge = 0 for both molecules).The G4-Asn 10 structure and MMP and AZM molecules were implemented for two systems.The first system contained 50 molecules of MMP and 1 molecule of PAMAM (MMP: PAMAM, 50:1) and the second one contained 50 molecules of AZM and 1 molecule of PAMAM (AZM: PAMAM, 50:1) in order to evaluate the interaction of MMP and AZM molecules over the entire surface of the PAMAM.For these models, G4-Asn 10 was set to the origin of the coordinate system.Subsequently, the PAMAM was solvated in a water periodic box (104 × 107 × 101 Å 3 ) of TIP3 water molecules, 34 then the corresponding number of MMP and AZM molecules were added in a random way along the water box (Figure 1).In order to neutralize the system, 64 chloride ions were added for the first and second system, respectively.Both systems were minimized for 10,000 steps and equilibrated through MD simulation.The MD simulations were performed at 298 K in the isobaric-isothermal ensemble for 20 ns.The temperature was maintained with the Langevin thermostat with a damping coefficient of 1 ps -1 whereas the pressure (1 atm) was kept constant using the Langevin piston method. 35Equations of motion were integrated with a 2 fs time step.The van der Waals (vdW) cut-off was set to 12 and the long-range electrostatic forces were computed using the particle-mesh Ewald approach.The MD simulation was performed using the NAMD2.9 31 program and CHARMM27 force field 36 for standard ions and water molecules.All analyses were done with VMD 37 software.For the analysis of the molecular dynamics, the data collected along the trajectories were used to calculate the radius of gyration, R g , of the PAMAM and different concepts were applied such as radial pair distribution function (RDF) of MMP and AZM molecules.Furthermore, several scripts were implemented in order to get relevant information from de MD, such as capture of MMP and AZM within a distance of 4 Å with respect of any atom of the PAMAM for the MMP/AZM-PAMAM systems, the number of hydrogen bond interactions between the PAMAM and the entrapped MMP and AZM molecules.
Building molecular structures: The molecular structures of MMP and the representative units of six functionalized PAMAM (monomer-FA, monomer-Cou, monomer-Arg, monomer-Lysine, monomer-Asparagine and monomer-Amine) were built using GaussView program. 38The geometry of these molecules was optimized at density functional theory (DFT) level using B3LYP method with 6-31G* as the basis set, which has been implemented in Gaussian 03 package program. 39lculation of interaction energies by nanoinformatic methodologies: Metropolis Monte Carlo simulations 27 were used to calculate the interaction energy between conformational pairs [molecule(1)-molecule(2)], in which the molecule(1) represents the monomeric unit of a specific PAMAM attached to a functional group and the molecule(2) represents the OP compound MMP.
The algorithm that generates the conformational sampling includes excluded-volume constraint method. 25rocedures used to calculate the interaction energy were as follows: (i) the mass centers of molecule(1) and molecule(2) are placed near the origin of the Cartesian coordinates frame, (ii) then, the new orientation of molecule (2) [in relation to molecule (1)] is chosen according to a set of three Euler angles (α, β, γ) obtained randomly, (iii) molecule( 2) is then translated along the random vector "n" until the vdW surfaces of each molecule touch each other, (iv) after translation, single-point energy (1SCF) for this specific molecular conformation [molecule(1)-molecule(2)] is calculated using a semi-empirical quantum mechanical method (PM6-DH+) 19 implemented in MOPAC2009™ packaged program, version 11.038L (LINUX), 29 (v) the heat of formation (ΔH f ) is extracted of the previous result and also from its isolated parts.Then, a supermolecular approach was used to obtain directly the interaction energy (ΔE int ) for molecule-molecule pairs, as the difference between the energy of the complex "molecule(1)-molecule(2)" and the sum of the energies of its isolated parts.According to the approach, the Gibbs energy molecule-molecule, molecule-polymer or polymer-polymer is defined by equation 2.
The entropy in the polymer systems 25 has been shown to have a small contribution in the ΔG, thus the term TΔS is negligible and equation 3 can be used.
The above equation is conveniently used to estimate the interaction energy with a standard quantum chemical procedure.In order to solve the superposition error caused by this approach in equation 3, a highly efficient conformational sampling, according to what is described by MD calculations, the ΔG of mixing polymer-molecule depends on the specific theory or model used.Thus, it was employed PM6-DH+ theory that is based on the Flory-Huggins lattice approximation, which is the simplest and most widely used to calculate the free energy of polymermolecule interactions, ΔG.A general expression through Flory-Huggins theory for the ΔG of a binary system (molecule-molecule, molecule-polymer, or polymerpolymer) can be described as: where ΔG is the free energy of polymer-molecule per mole of lattice site, f is the volume fraction, and x is the chain length with each repeating unit defined as occupying a lattice site.The Flory-Huggins c parameter is defined as: where z is the coordination number of the model lattice, and ΔE int is the energy of formation of an unlike pair (ΔE l,2 ).According to equation 5, the energy of a particular [1,2]  pair is defined as: Finally, assuming that ΔG = ΔE int (equation 3), the ΔE int described in this work was defined as: Thus, steps (i) through (v) are repeated until the obtention of 100,000 different molecular configurations and their corresponding interaction energies.Then, from this interaction energies distribution, the average "<ΔE (molecule1 -molecule2) >" was calculated.
The molecule(1) is replaced for each one of the six representative units of the functionalized PAMAM (monomer-FA, monomer-Cou, monomer-Arg, monomer-Lysine, monomer-Asparagine and monomer-Amine) and the above calculations are performed again to finally obtain six-interaction energy averages.
In order to compare and demonstrate the effectiveness of the above methodology, conformations associated with the best interaction energies for each of the 6 analyzed complexes were selected, and energies calculated at DFT level using B3LYP method with 6-31G (d,p) as basis set, which has been implemented in Gaussian 03 package program37.Finally, the interaction energies were calculated based on the supermolecular approach using the equation 7.

AChE inhibition assay
We evaluated in vitro the AChE inhibitory activity of net MMP and encapsulated MMP with G4-Asn and G4 Arg/Lys.The MMP was measured as follows: PAMAM derivatives with MMP were incubated by 45 min on stirring, (molar ratio of MMP:PAMAM = 12:1), after that the AChE inhibitory activity of MMP/PAMAM was measured.Additionally, the AChE inhibitory activity of PAMAM derivatives alone was also evaluated to disregard any AChE inhibitory activity of PAMAM.The assay for measuring AChE inhibitory activity was carried out according to Ellman et al. 40 and adapted to 96-well microtiter plates by Lopez et al. 41 and described by Gutierrez et al. 42 The enzyme was obtained from electric eel (C-3389, type VI-S, Sigma Chemical Co., St. Louis, MO).In brief, 50 mL of the sample dissolved in phosphate buffer (8 mmol L -1 K 2 HPO 4 , 2.3 mmol L -1 NaH 2 PO 4 , 150 mmol L -1 NaCl, and 0.05% Tween 20 at pH 7.6) and 50 mL of the AChE solution (0.25 unit mL -1 ) in the same phosphate buffer were added to the wells.Then, the plates were incubated for 30 min at room temperature.After the incubation time, 100 mL of the substrate solution [40 mmol L -1 Na 2 HPO 4 , 0.2 mmol L -1 5,5'-dithio-bis (2-nitrobenzoic acid) (DTNB), and 0.24 mmol L -1 acetylthiocholine iodide (ACTI) in distilled water at pH 7.5] was added.The absorbance was read in a Bio-Tek Instrument microplate reader at 405 nm after 3 min.The enzyme activity was calculated as a percentage compared to a control using only the buffer and enzyme solution.The compounds were assayed at 100 mg mL -1 , and the alkaloid galanthamine was used as the reference compound.When the enzyme inhibition was > 50% at 100 mg mL -1 , dilutions were performed to determine the corresponding IC 50 values.The IC 50 values were calculated by means of regression analysis from three individual determinations.

Synthesis of PAMAM derivatives
FA and Cou were attached to the amine terminal groups of PAMAM-G5 by EDC•HCl and HOBt coupling reaction. 12,39An excess amount of reagents (1:140) was used to allow functionalization of the 128-terminal amine groups of the PAMAM-G5 (MW ca.28826) in order to ensure the maximum functionalization.However, only 18 FA (MW of PAMAM G5 + 18 FA ca.36800) and 67 Cou molecules (MW of PAMAM G5 + 67 Cou ca.41600) were found attached to the G5, respectively.It is proposed that the difference in the number of molecules coupled to G5 was due to steric hindrance caused by the size of FA and Cou, respectively (Figures 2a and 2b).
Lys, Arg and Asn were coupled to the terminal amines of the G4 by EDC•HCl and HOBt coupling reaction. 12,39An excess amount of reagent (1:70) was used in all three cases to allow maximum possible functionalization.Nevertheless, only 44 molecules of Lys were conjugated to PAMAM (MW G4 + 44 Lys ca.30500), 15 molecules of arginine were conjugated to G4 (MW G4 + 15 Arg ca.19800), and 10 molecules of Asn were found attached to the G4 (MW G4 + 10 Asn ca.18500).Again, it is proposed that the specific number of molecules attached to G4 is determined by the steric hindrance and physicochemical behavior of Arg, Lys and Asn, respectively (Figures 2c-e).
The surface functionalization of PAMAM-G4 with two different amino acids (Arg and Lys) was also achieved.In this case, a small amount of arginine (1:5) was used to accommodate second amino acid lysine by subsequent reaction.With this simple methodology, coupling of two different amino acids on the surfaces of PAMAM-G4 was obtained.The goal of this approach was to increase the capacity of trapping to the OP (synergy) by PAMAM derivatives through more hydrogen bonds, weak interactions such as vdW, hydrophobic bonds, and electrostatic interactions.Thus, it was possible to conjugate 14 molecules of Arg (MW ca.18600) to G4, and additionally introduce 25 molecules of Lys (MW ca.25200) (Figure 2f).
The adsorption behavior of OPs compounds (specifically MMP by PAMAM derivatives) in model solutions was studied through HPLC analysis, but the desorption behavior was not performed.(Figure S1, SI section).It was observed that the adsorption of MMP was higher when it was exposed to G4-Arg/Lys and G4-Asn achieving 45% and 64% of affinity, respectively (Figure 3).This high affinity may be due to the small size of MMP and multiple hydrogen bonds that occur between the tertiary amine and phosphoryl group from MMP and secondary, tertiary amines (interdendritic) and carbonyl groups of PAMAM derivatives (Figure 4).Molecular dynamics calculations can predict and evaluate the contribution of other energies involved in the dendrimer-OP system.
In this system, the hydrophobic interaction strength is obtained through the vdW energy.The variation of the vdW energy for the MMP-PAMAM and AZM-PAMAM systems, are shown in Figure S2 (SI section).The vdW energy pattern do not change significantly with time, it confirms that vdW interaction is less significant than the hydrogen bonding interaction.Another point to consider is that G4-Asn provides greater amount of carbonyl groups to the system than other polymers studied here, allowing the formation of more hydrogen bonds between the PAMAM and the OP.

Bioinformatic studies
Six, three-dimensional atomic models for monomerend-group complexes were generated using the in silico design.These six functionalized PAMAM representative structures were used to estimate the PAMAM/OPs interactions.
The average interaction energy calculated by the nanoinformatic methodologies is depicted in Table 2.The theoretical data was correlated with the experimental percentage of affinity calculated by HPLC (Figure 5).For the six representative units of the studied functionalized PAMAM, the relative order of affinity by MMP was predicted as: monomer-Asparagine > monomer-Lysine > monomer-Folate > monomer-Coumarin > monomer-Arginine > monomer-Amine.As shown in Table 2 and Figure 5b, interaction energies calculated through DFT-6-31G method showed a better correlation (r 2 = 0.995) as expected, corroborating the data obtained at semiempirical level.Furthermore, comparing DFT and PM6-FH+ correlations (r 2 ), it can be noted that semi-empirical method gave an acceptable and highly predictive  correlation.Moreover, the data obtained by DFT showed to be much more accurate than semi-empirical methods, however, considering the efficiency and speed required for filtering several of thousand conformations, the semiempirical method probed to be commendable.
According to the obtained results, the functionalized PAMAM improved its ability to trap OP molecules; thus, strategic changes of the functional groups on the surface of PAMAM have proven to be relevant when defining their interaction with molecules of interest.
Computational chemistry approaches through nanoinformatics platforms have proved to be essential in the estimation of the PAMAM-molecule interactions, due to the speedy results, the low cost of in silico experiments (this is particularly true for those scenarios where the handling of toxic or otherwise hazardous materials is necessary) and the ability to relate the theoretical values with experimental data with high accuracy (Figure 5).To characterize the G4-Asn 10 models and their interaction with MMP and AZM molecules at molecular level, we performed MD simulations to subsequently make a structural analysis of the interactions present in the system.The first analysis performed was the radius of gyration (R g ) calculated as a function of simulated time, as shown in Figure S3 (see SI section).Furthermore, the mean radius of gyration (R g ) was calculated including only the last 8 ns of the 20 ns equilibration simulations, where appeared to take on a steady value.For G4-Asn 10 with MMP molecules, values of 21.0 ± 0.2 Å were observed, while for G4-Asn 10 with AZM molecules, 21.6 ± 0.3 Å.The second analysis was RDF, which represents the probability to find an atom in a shell (dr) at the distance r, of another atom or molecule.Figure S4 (SI section) shows RDF of G4-Asn 10 /MMP (1:50) and G4-Asn 10 /AZM (1:50) systems.The RDF shows a radius r = 4 Å as the distance from which the highest number of MMP molecules are found, whereas  for AZM molecules, the radius was higher and lower its probability, radius r = 6 Å and r =10 Å for AZM molecules.Therefore, the MMP molecules are positioned at locations closer to PAMAM than AZM molecules, allowing greater interaction with the terminal groups of G4-Asn 10 .Distances were calculated between phosphorus atoms (P) of MMP/ AZM molecules and nitrogen atoms (N) of terminal amines and Asn functional groups.The representation of hydrogen bonds between MMP/AZM and G4-Asn 10 are shown in Figure 6.The AZM molecule can only form one hydrogen bond with the G4-Asn 10 monomers, and the majority of the contributions come from the hydroxyl oxygen atom not associated with the phosphate group.On the other hand, a MMP molecule can achieve two hydrogen bonds with G4-Asn 10 , the hydroxyl oxygen atom from phosphate group and hydrogen atoms from amine group.Thus G4-Asn 10 can produce more hydrogen bonds in the presence of MMP molecules than AZM, for instance, in the Figure S5 (SI section) is depicted the number of hydrogen bonds as functions of time for MMP and AZM molecules interacting with G4-Asn 10 .The higher number of hydrogen bonds is generated in the interaction  between G4-Asn 10 with MMP molecules.The capability of entrapment of G4-Asn 10 was calculated through the number of molecules near the PAMAM at a distance of 4 Å or less (value obtained from the RDF analysis).The number of MMP trapped by G4-Asn 10 was higher than AZM: approximately 18-20 MMP versus 10-12 AZM molecules in the last 4 ns of the 20 ns equilibration simulation, in the Figure S6 (available as SI) is depicted the number of MMP and AZM trapped by functionalized PAMAM.Thus, it was rationalized that the entrapment of these molecules was promoted primarily by hydrogen bonds formation, especially by of the MMP molecules ability to form more hydrogen bonds with G4-Asn 10 .

AChE inhibition assay of PAMAM derivatives
To give further support to our model studies for capturing affinity and molecular simulations data, the best PAMAM derivatives (G4-Asn and G4-Arg/Lys, according to percentage of affinity) were employed to evaluate the inhibition of the AChE activity in vitro.Initially, the inhibition activity of MMP to AChE was measured, obtaining an IC 50 of 0.85 μg mL -1 , confirming the higher toxicity of MMP.Next, MMP was exposed (incubated) to PAMAM derivatives and G4-NH 2 (commercial) for a period of 45 min with constant stirring (molar ratio MMP:PAMAM = 12:1).IC 50 values for commercial PAMAM, MMP and derivatized PAMAM are given in Table 3.In the case of MMP trapped by G4-NH 2 , the data obtained indicated that the IC 50 increased slightly, showing IC 50 values of 4.23 μg mL -1 .Whereas the MMP trapped by PAMAM derivatives reached IC 50 values of 65 μg mL -1 and 132 μg mL -1 for G4-Asn and G4-Arg/Lys, respectively.Therefore, the IC 50 of MMP encapsulated PAMAM derivatives was significantly higher than MMP encapsulated PAMAM-NH 2 and free MMP (see Table 3).Finally, the toxicity of PAMAM derivatives was evaluated (G4-Asn and G4-Arg-Lys).In both cases, the values were higher than 300 μmol L -1 (see Table 3), demonstrating that the toxicity of PAMAM is not significant.These data indicate that PAMAM derivatives have the capacity to capture MMP efficiently, and perhaps, could be used in a near future in preventing the poisoning events.Molecular simulation calculations performed in this study predict that there is a good "hydrogen bond" formation and "radius of gyration" between G4-Asn with MMP.This is consistent with the AChE activity assay, where the IC50 of MMP increased significantly when coupled with PAMAM derivatives.We hypothesize it is due to a higher MMP-PAMAM union, due to the formed hydrogen bonds and the closeness of the atoms in each of the molecules involved in the interaction system.Furthermore, if we analyze the in vitro experimental data obtained in the AchE activity assay, we could conclude that the interaction MMP-PAMAM is strong, limiting the contact between MMP and AchE enzyme, presumably by the encapsulation of MMP inside of PAMAM.

Conclusion
Functionalization of PAMAM derivatives with different amino acids such as lysine, arginine, asparagine and others biomolecules (coumarin and folic acid) was demonstrated.These PAMAM derivatives were employed to trap AZM and MMP, OP compounds broadly used as pesticides in the agricultural sector.The molecules were designed and modeled using novel bioinformatic tools, and chemically synthesized, (specifically G4-Asn and G4-Arg/Lys), to demonstrate significant affinity for MMP.The results obtained by HPLC in model solutions, support theoretical chemistry through nanoinformatic methodologies to estimate the interaction energies between MMP and six functionalized PAMAM.In silico results have proven to be highly predictive to estimate a priori, a tendency of affinity from a group of MPP by specific PAMAM, and allowed observing that the PAMAM-molecule interactions play an important role for the development of chemical encapsulation and attachment systems.Moreover, by comparing theoretical affinity of PAMAM derivatives with AZM and MMP, it was demonstrated that the better trapping of MMP was produced by PAMAM-Asn primarily due to hydrogen bonds formation, especially by the MMP ability to form more hydrogen bonds with G4-Asn 10 .Furthermore, our results of affinity by HPLC and molecular dynamics correlated highly with the in vitro enzyme AChE activity assays.The inhibition of AChE by MMP in solution, highly decreased when G4-Arg/Lys and G4-Asn were used compared with commercial PAMAM G4.This information allows us to conclude that PAMAM derivatives G4-Arg/ Lys and G4-Asn are capable of capturing MMP in solution efficiently, and greatly reduce the MMP action in the enzyme and could be used as possible nano-detoxification agents of OPs in the future.In vivo studies for nanodetoxication of OPs are under current investigation.

Figure 1 .
Figure 1.G4-Asn 10 models on a box containing explicit water molecules and Cl -ions.The dendrimer structures are shown by space-filling representation with carbon, nitrogen and oxygen atoms shown in green, blue and red, respectively.The Asn end groups are shown in purple.The Cl -is shown in brown.The chemical structures of MMP and AZM are shown in detail in the center.The distribution of 50 of these molecules in the water box are shown as spacefilling representation with carbon, nitrogen, oxygen, phosphorous and sulfur atoms shown in gray, blue, red, orange and yellow, respectively.

Figure 3 .
Figure 3. Experimental assays indicating the percentage of PAMAM derivatives affinity to MMP.As control, G4-NH 2 and G5-NH 2 were used.

Figure 4 .
Figure 4. Representation view of the interaction between MMP and G4-Asn indicating the role of hydrogen bond in the capture.

Figure 5 .
Figure 5. Correlation between the experimental affinity values obtained by HPLC studies: (a) the interaction energy averages calculated by using semiempirical method (MOPAC); and (b) the interaction energies calculated by DFT (B3LYP/6-31G).

Table 2 .
Values of the theoretical interaction energies (at semi empirical and DFT level) and experimental (%) affinities between MMP and the monomeric unit of PAMAM attached to six different functional groups

Table 3 .
AChE activity assay with MMP, PAMAM derivatives and their formulations a Mass and concentration displayed is in relation to the MMP amount encapsulated by PAMAM (MMP/PAMAM = 12:1).n = 3, p < 0.05.