Gas-Phase and Microsolvated Glycine Interacting with Boron Nitride Nanotubes . A B 3 LYP-D 2 * Periodic Study

The adsorption of glycine (Gly) both in gas-phase conditions and in a microsolvated state on a series of zig-zag (n,0) single-walled boron nitride nanotubes (BNNTs, n = 4, 6, 9 and 15) has been studied by means of B3LYP-D2* periodic calculations. Gas-phase Gly is found to be chemisorbed on the (4,0), (6,0) and (9,0) BNNTs by means of a dative interaction between the NH2 group of Gly and a B atom of the BNNTs, whose computed adsorption energies are gradually decreased by increasing the tube radius. On the (15,0) BNNT, Gly is found to be physisorbed with an adsorption driving force mainly dictated by π-stacking dispersion interactions. Gly adsorption in a microsolvated environment has been studied in the presence of seven water molecules by progressively microsolvating the dry Gly/BNNT interface. The most stable structures on the (6,0), (9,0) and (15,0) BNNTs present the Gly/BNNT interface fully bridged by the water solvent molecules; i.e., no direct contact between Gly and the BNNTs takes place, whereas on the (4,0) BNNT the most stable structure presents a unique direct interaction between the COO Gly group and a B atom of the nanotube. Further energetic analyses indicate that the (6,0), (9,0) and (15,0) BNNTs exhibit a low water affinity, which favors the Gly/water interactions upon BNNT coadsorption. In contrast, the (4,0) BNNT has been found to show a large water affinity, bringing the replacement of adsorbed water by a microsolvated glycine molecule as an unfavorable process.


Introduction
Bioconjugated nanostructured materials resulting from the coupling of biomolecules with inorganic nanomaterials including nanotubes, nanowires, nanoparticles and nanosheets have attracted much attention during the last years as they exhibit unique features derived from combining synergistically the properties of the interacting components.These exclusive physico-chemical properties render these materials as suitable substrates with potential applications in diverse biological- [1,2] and material-related [3] areas such as biocatalysis, [4,5], drug delivery [6][7][8], biosensing [9][10][11][12][13] and medical diagnostics [14,15].The functionalities resulting from these biohybrid materials are largely mediated by the biomolecule/inorganic surface interactions, which in turn are dictated by the structure-specific binding properties of the two partners.Accordingly, precise knowledge on the interactions between the biomolecule and the inorganic components is of fundamental relevance.
Among the different nanostructured materials, boron nitride nanotubes (BNNTs) have been proposed to be suitable candidates to be combined with biomolecules [16].BNNTs are isosteres and structurally similar to carbon nanotubes (CNTs), in which alternating B and N atoms substitute for C atoms.However, these two type of nanotubes exhibit different physico-chemical properties.Whereas CNTs exhibit metallic or semiconducting behavior, which moreover is strongly dependent on the tube diameter, helicity and concentric layers, BNNTs are electrical insulators with a band gap of ca.5.5 eV regardless of the tube geometry features [17].Moreover, at variance with the non-polar C-C bonds in CNTs, the B-N bonds of BNNTs exhibit a certain polar character, the degree of which depends on the curvature of the nanotube.That is, the increase of the tube curvature induces the transformation of the sp 2 hybrid character of the B and N atoms in large diameter BNNTs into a sp 3 one in small diameter BNNTs.As recently shown by us [18], this has important consequences for the nature of interaction between functional molecules and the BNNTs walls; i.e., polar molecules strongly chemisorb on small radius BNNTs, whereas interaction of non-polar molecules are energetically more favourable when physisorbed on large radius BNNTs.Furthermore, unlike CNTs, which present an inherent cytotoxicity [19], BNNTs have been found to be nontoxic [20] due to their high chemical and structural stability and high oxidation resistance, which alongside their uniformity and stability in dispersion in solution [21] make them suitable for biomedical applications.
Different experimental studies have focused on the interaction of peptides and proteins with BNNTs, showing a natural affinity between the two conjugates, which allows a direct immobilization of proteins on the BNNTs [22] as well as the isolation of individual BNNTs through a novel pathway based on peptide wrapping [23].Moreover, biofunctionalized BNNTs via glycine interaction are good reactant substrates to obtain polysaccharide-coated BNNTs under mild conditions, in which the role of glycine is crucial during the interfacial process.The interactions of DNA and RNA with BNNTs have also been addressed and exploited to obtain nematic ordered ensembles of BNNT [24].Other works have been devoted to assess the cytotoxicity of BNNTs when in contact with cells.Chen et al. [20] concluded that pristine BNNTs are inherently non-cytotoxic in view of the non-altered growth of human embryonic kidney (HEK) cells when cultured with BNNTs.Similar results were found by Ciofani and coworkers, in which coated-BNNTs presented a good cytocompatibility with human cells [25][26][27][28].However, Goldberg and coworkers more recently found that BNNTs are actually cytotoxic for cells present in the lung alveoli and for HEK, in which the discrepancies with the other works were discussed and suggested to be due to the different morphology and size distribution of the BNNTs tested and the different assay techniques [29].
Theoretical works, mainly based on density functional theory (DFT) methods, have also studied the interaction of biomolecules with boron nitride nanostructures, most studies being limited to biomolecule building blocks, (amino acids and DNA and RNA nucleobases) due to the demanding computational cost of these calculations.Works on the gas-phase interaction of nucleobases using the local density approximation (LDA) and generalized gradient approximation (GGA) DFT levels of theory showed that this depends on the individual polarizations of the nucleobases [30][31][32].The interaction of BNNTs with glycine (Gly) among other different amines has been studied in the gas-phase revealing an affinity of the BNNT with the NH 2 Gly group [33].Study on the gas-phase interaction of the arginine (Arg), aspartic acid (Asp) and tryptophane (Trp) amino acids, with basic, acidic and aromatic side chain functionalities, respectively, at the LDA DFT level revealed that the binding is accompanied by charge transfer following the trend of Arg > Asp > Trp [34].The binding of different biomolecules inside the cavity of BNNTs has also been studied at the LDA level [35].The calculated weak interactions led the authors to suggest BNNTs to be suitable biological carriers due to the limited delivery kinetic barrier.
All these works focus on the intrinsic adsorption properties; i.e., they are limited to the gas-phase and, accordingly, solvation effects were not accounted for.Moreover, each work addresses the interaction of biomolecules with a particular BNNT.Since it has been shown that the tubular radius can modulate the adsorption properties of BNNTs [18], which is also applicable for biomolecules, the nature of interaction can significantly be different depending on the radius of the BNNT.Moreover, for the particular works addressing the interaction of amino acids, no conformational exploration to find out the most stable amino acid/BNNT adduct was performed (i.e., the initial amino acid conformation guesses were the most stable gas-phase structure), which is an important drawback due to the large conformational mobility of these molecules.In order to provide a more complete atomic-scale description of the interaction of amino acids with BNNTs, the present work reports a systematic periodic B3LYP-D2* study, using a hybrid functional and including dispersion corrections, on the interaction of Gly with different zig-zag (n,0) single-walled BNNTs (n = 4, 6, 9 and 15) rendering nanotubes of different radius.Note that dispersive effects, not included in previous works, are expected to play a role in these systems.Moreover, with the aim to study in a more realistic way the interactions between biological systems and BNNTs, the very same Gly/BNNTs interaction study has been addressed considering a microsolvated environment modeled by the presence of seven water molecules.The effect of water has been analyzed from a structural and energetic point of view, with particular attention paid to the Gly/BNNT interface to determine whether the interaction is direct or bridged by the water molecules.

Gas-Phase Interaction
An exhaustive description relative to the modeling of the different BNNTs used in this work in their pristine state is available in a very recent work by us [18].The BNNTs are automatically generated by the CRYSTAL code through the wrapping of a periodic boron nitride monolayer (hexagonal P3m1 layer symmetry group) into cylinders of different radius and fully exploiting the symmetry operators of the nanotubes [36][37][38][39][40].In the interest of the present work, it is worth mentioning that the calculated electrostatic potential maps indicate a prominent positive/negative valued potential region for the (4,0) BNNT, which become progressively less pronounced for BNNT with increasing radius, until obtaining a practically shallow electrostatic potential for the (15,0) BNNT.For the interaction with Gly, the geometry optimizations were carried out as 1D polymers within the P1 space group, in which the unit cell parameters have been enlarged twice to avoid lateral interactions between molecules of adjacent unit cells.The different optimized adducts for the adsorption of Gly on the BNNTs are shown in Figure 1, whereas the calculated adsorption energies alongside the pure electronic and dispersion contributions are shown in Table 1.The adsorption energies (∆E ads ) per unit cell of the probe molecules with the BNNTs are computed as: where E(Gly/BNNT) is the energy of a fully relaxed unitary cell containing the BNNT in interaction with Gly, E(BNNT) is the energy of a fully relaxed unitary cell of the BNNTs alone, and E m (Gly) is the molecular energy of the free Gly.For all the BNNT systems, six initial structural guesses were considered (see Figure 1a): (i) pure interaction between the NH 2 group of Gly and one nanotube B atom (hereafter referred to as BN/NH 2 ); (ii) BN/NH 2 interaction plus H-bonding between the Gly OH group and one nanotube N atoms (hereafter referred as BN/NH 2 -OH); (iii) and (iv) interaction between the Gly CO group and one nanotube B atom plus H-bonding between the Gly OH group and one nanotube N atom, with the difference that the nanotube B and N atoms are chemically bonded or not (hereafter referred to as BN/COOH-1 and BN/COOH-2, respectively); (v) interaction of Gly in its zwitterionic form (hereafter referred to as BN/zwitt), in which the Gly COO − group interacts with one B atom and the Gly NH 3 + is H-bonded to one nanotube N atom; and (vi) interaction between Gly and the BNNTs purely through the π system of the Gly COOH group (hereafter referred to as BN/π).Data reported in Table 1 clearly indicate that the most stable adduct for the (4,0) BNNT is the BN/NH 2 -OH whereas for the (6,0) BNNTs both BN/NH 2 and BN/NH 2 -OH are nearly degenerate.These complexes result from dative covalent interactions between the Gly NH 2 group and the B atoms, which act as Lewis acid sites.These findings are consistent with the data reported for the interaction of these BNNTs with probe molecules, in which the interaction of NH 3 with the (4,0) and (6,0) BNNTs was found to be the strongest one among all tested molecules [18].The fact that the BN/NH 2 -OH adduct becomes the most stable one for the Gly/BNNT(4,0) system is consistent with the large polar character of this nanotube, the N atoms acting as H-bonding acceptor groups.This is not in line with the most stable BN/NH 2 adduct found for the (6,0) BNNT and is due to the weaker H-bonding acceptor character of the N nanotube atom when increasing the nanotube radius, which is reflected by an increase of the H-bond distance in the BN/NH 2 -OH adducts (1.686 and 1.748 Å for the (4,0) and (6,0) BNNTs, see Figure 1a).In the same way, the B-N Gly bond lengths of the dative interactions in the BN/NH 2 and BN/NH 2 -OH complexes also increase with increases in the nanotube radius due to the progressive decrease of the Lewis acceptor character of the B atom.Because of that, the calculated adsorption energy (∆E ads ) is more negative and larger for the most stable BN/NH 2 -OH adduct of Gly/BNNT(4,0) than the most stable BN/NH 2 adduct of Gly/BNNT(6,0) (i.e., −33.2 and −18.9 kcal mol −1 , respectively).The calculated energetic contributions; i.e., purely electronic and dispersion (∆E el and ∆E D values of Table 1, respectively) indicate that Gly adsorption on the (4,0) BNNT is largely dictated by the covalent dative interaction, whereas on the (6,0) BNNT ∆E el decreases in favor of dispersion.An analysis of the Mulliken charges (Q values of Table 2, only limited to the most stable Gly/BNNT complexes) confirms the formation of charge transfer complexes for both the Gly/BNNT(4,0) and Gly/BNNT(6,0) systems, the computed charge transfer values from Gly to the BNNTs being 0.30 e and 0.22 e, respectively.The most stable Gly/BNNT(15,0) adduct has been found to be the BN/π one.This is in perfect agreement with the EPM results, which point out the (15,0) BNNT as a practically non-polar nanomaterial.Data reported in Table 1 clearly indicate that the binding mechanism involved in this adduct is mainly based on dispersive forces (calculated ∆E ads is practically equal to ∆E D ; i.e., −10.2 and −11.7 kcal mol −1 , respectively) dictated by π-stacking interactions between the π systems of the COOH group and the B-N hexagon rings of the (15,0) BNNT.Because of the presence of only non-covalent interactions, ∆E ads is less negative compared to Gly interaction on (4,0) and (6,0) BNNTs.It is worth mentioning that the very same BN/π complexes have also been calculated for the Gly/BNNT (4,0) and (6,0) systems, meaning that, in the former case, the structure collapses onto the BN/NH 2 complex, whereas for the latter case the calculated ∆E ads is found to be 9.0 kcal mol −1 above the most stable one, due to the lower propensity of this BNNT to establish π-stacking interactions.The interaction of Gly with the (9,0) BNNT is a frontier case between small radius (i.e., (4,0) and (6,0)) and large radius (i.e., (15,0)) BNNTs.Although the BN/NH 2 adduct has been found to be the most stable one, the BN/π complex is the second most stable one lying 1.9 kcal mol −1 above.Calculated ∆E ads values, however, can suffer from the basis set superposition error (BSSE).Indeed, upon correction, results indicate that these two complexes are nearly degenerate (BSSE-corrected ∆E ads values being −9.4 and −8.9 kcal mol −1 for the BN/NH 2 and BN/π, respectively).
The interaction of Gly through the COOH group by means of a simultaneous CO-B dative bond and a OH•••N (BNNT) H-bond has also been considered (see BN/COOH-1 and BN/COOH-2 adducts).Although none of the calculated complexes are the most stable ones, important structural and energetic features deserve to be mentioned.For the BN/COOH-1 adduct, Gly adsorption on the (4,0) BNNT results in a spontaneous proton transfer from the Gly COOH group to the N nanotube atom, hence forming a COO-/BNNT-H + ion pair (see Figure 1b, structure of right).Such a proton transfer was already observed for the adsorption of HCOOH on the very same (4,0) BNNT and is attributed to the net charge transfer occurring from Gly to the BNNT, which induces an increase of the COOH acidity and the nanotube basicity up to the point of promoting the proton transfer to a nearby N atom of the nanotube.Moreover, for this adduct the CO-B distance is significantly shorter than for the other BN/COOH-1 adducts (1.487 Å versus 2.514-2.921Å, respectively), which results in a stronger interaction (∆E ads = −31.7 kcal mol −1 and ≈ −9.6-−7.5 kcal mol −1 , respectively).For the BN/COOH-2 complex on the (4,0) BNNT, no proton transfer has been found; although the OH•••N (BNNT) H-bond and the CO-B dative bond are actually shorter than those present in the other BNNTs, in line with what has been described for the BN/COOH-1 cases.Interestingly, the difference between BN/COOH-1 and BN/COOH-2 is that in the former the COOH interaction occurs on B and N atoms that are chemically bonded to each other, whereas in the latter this is not the case.Accordingly, the fact that the spontaneous proton transfer only occurs in the former system seems to indicate that the charge transfer in enhanced by a cooperative effect between the OH•••N (BNNT) and the CO-B interactions when the B and the N atoms are chemically bonded, which is in agreement with the larger and more negative calculated ∆E ads values (−31.7 and −17.6 kcal mol −1 for BN/COOH-1 and BN/COOH-2, respectively).
Finally, it is worth mentioning that the interaction of Gly in its zwitterionic state has also been computed (BN/zwitt).On the (4,0) and (6,0) BNNTs, a spontaneous proton transfer from the NH 3 + group to the N nanotube atom has been found, whereas for the (9,0) and (15,0) BNNTs the zwitterionic form is maintained.Consistently, calculated ∆E ads values are negative for the two former adducts (−31.2 and −13.5 kcal mol −1 , respectively), whereas for the two latter ones they have been found to be positive (+1.5 and +7.6 kcal mol −1 , respectively) and, accordingly, are not stable complexes.

Microsolvated Interaction
Here, results on the interaction of Gly with the BNNTs in the presence of seven water molecules are reported.We have introduced seven water molecules since this is the minimum number of water molecules to have a relatively complete first-solvation shell of Gly upon adsorption; i.e., three water molecules interacting with the NH 3 + group, two water molecules interacting with the COO − group and two more water molecules to complete the solvation shell.For these cases, the unit cell parameters of the BNNTs have been enlarged thrice to avoid lateral interactions between water molecules of adjacent unit cells.It is worth mentioning that a statistical sampling of the hypersurface of these systems can be carried out adopting either the molecular dynamics or the Monte Carlo approaches [41].However, these calculations are extremely expensive at the ab-initio level adopting realistic models for BNNTs like the (9,0) and the (15,0) ones, which contain 129 and 201 atoms, respectively.For the present work, we have followed a different approach consisting of a progressive microsolvation of the dry interface at the Gly/BNNT structures.This microsolvation procedure consists of adding water molecules at the dry Gly/BNNT interface in such a way that Gly progressively loses direct contact with the BNNTs up to a point in which the interaction is fully bridged by water.This procedure has already been performed by some of us in other works for the interaction of Gly with silica [42] and hydroxyapatite [43] surfaces.Since the most stable state of Gly in water is the zwitterionic one, we considered the BN/zwitt adducts as initial guesses for the progressive microsolvation.For the (4,0), (6,0) and (9,0) BNNTs, the resulting structures are shown in Figure 2a.In the BN/CONH adduct, Gly directly interacts with the BNNTs in a similar fashion as in the gas-phase (i.e., COO − •••B BNNT dative bond and NH 3 + •••N BNNT H-bond ) , while the seven microsolvating water molecules are simple spectators interacting with available points of the COO − and NH 3 + groups through H-bonding.It is worth remarking that now for the (4,0) and (6,0) BNNTs no H transfer from Gly to the BNNT occurs (at variance with the gas-phase adsorption, vide supra) due to the screening effect of water.The BN/CO-w/NH and BN/NH-w/CO adducts result from moving one spectator water molecule from the outer shell to the inner shell so that the following water mediated interactions NH 3 occur respectively.Finally, From these two adducts, a second water displacement to remove the remaining direct Gly/BNNT interaction gives the w/CONH adduct, in which water fully mediates the Gly/BNNT contact.It is worth mentioning that each of the seven H 2 O molecules can in principle be displaced from their positions to lead to a water mediated contact between Gly and BNNT.We choose the one exhibiting the weakest interaction energy with the other water molecules by computing the cost to remove one water molecule from the BN/CONH adduct by a single point energy evaluation for each H 2 O.
For the (15,0) BNNT, all the optimization calculations collapsed to structures with no direct contact between Gly and BNNT, the most stable one being presented in Figure 2b.In this structure, water fully solvates the Gly molecule and, at variance with the other BNNT cases, no charge transfer between water and BNNT takes place, due to the highly apolar character of this nanotube.values in parenthesis for the (6,0) one; values in brackets for the (9,0) one; and italic underscored values for the (15,0) one.For this latter case, the distance is that between the C atom and the plane defined by the closest B-N hexagon ring.
The relative stabilities between the different calculated adducts for a given microsolvated complex are shown in Table 3.As one can observe, for the (6,0), (9,0) and (15,0) BNNTs the most stable systems are the w/CONH adducts; i.e., those in which no direct Gly/BNNT contact occurs, whereas for the (4,0) one, the BN/CO-w/NH adduct (direct interaction only through the COO − ) was found to be the most energetically stable one.It is worth mentioning, however, that the energy difference between the BN-CO-w/NH and w/CONH adducts for the (4,0) BNNT (the first and second most stable ones) is relatively small (2.6 kcal mol −1 ) and, accordingly, it might be inverted due to entropic effects associated with water rearrangement, as it is shown for peptide adsorption on hydrophobic and polar surfaces [44].To further analyze this point, finite temperature molecular dynamics simulations would be desirable.Besides these results, three different processes have moreover been considered to study the stability of the structures shown in Figure 2. The first one involves Gly solvated by 7 H 2 O molecules being adsorbed on the clean walls of the BNNTs, whose reaction energy was computed as (reported by ∆E R1 of Table 2) where Gly/7w is glycine solvated by the seven H 2 O molecules and Gly/7w/BNNT represents the microsolvated complexes.The ∆E R1 column shows that for all BNNTs the process is exoenergetic, meaning that the structural rearrangement of the seven H 2 O molecules around Gly is compensated by the interaction with the BNNTs.Remarkably, limited to the most favorable adducts per BNNT, ∆E R1 values are less negative with increases in the nanotube radius, consistent with the less polar behavior of the BNNTs.The second process envisages gas phase Gly adsorbed on the already microsolvated BNNTs by the seven water molecules, whereas the third one envisages the solvated Gly being adsorbed on the seven water hydrated BNNTs giving rise to the shown adducts with expulsion of seven water molecules (here considered as a H-bonded cluster), whose reaction energies were computed as (reported by ∆E R2 and ∆E R3 of Table 2) Gly/7w + 7w/BNNT → Gly/7w/BNNT + 7w (4) where 7w/BNNT is the BNNT solvated by the seven H 2 O molecules and 7w the H-bonded cluster made up by seven water molecules.These last two processes require the BNNTs in the presence of seven water molecules, whose optimized structures are given in Figure 3.The initial guess of these systems were the corresponding w/CONH adducts (no direct interaction between Gly and the BNNTs and accordingly the interaction between water molecules and BNNTs is maximum) in which the Gly molecule was removed.Consistent with the polar character of the BNNTs, several interactions between the water molecules and the (4,0) BNNT via covalent dative and H-bond interactions take place, whereas by increasing the BNNT radius these interactions are progressively missed, up to the point in which for the (15,0) no apparent interaction is observed.The energetics of R2 is a tradeoff between the water affinity of the BNNTs and Gly (R2 is in essence the capture of the adsorbed water by Gly).The trend provided by the calculated ∆E R2 values indicates that R2 is more favorable when increasing the radius.That is, since the (15,0) BNNT does not exhibit water affinity, the calculated ∆E R2 values are large and negative due to the strong interaction of water with Gly.In contrast, both the (4,0) BNNT and Gly exhibit large water affinity and accordingly the calculated ∆E R2 values are the less negative ones of the series.R3 is probably the most physically sound process as it involves the replacement of adsorbed water solvent by an already solvated Gly.Calculated ∆E R3 values indicate that such a water replacement is energetically favorable on the (15,0) and (9,0) BNNTs probably due to the low water affinity of these BNNTs.On the (6,0) the process is still favorable by some amount but less than on the other two nanotubes (∆E R3 = −2.5 kcal mol −1 ), whereas on the (4,0) calculated ∆E R1 was found to be positive, indicating that the overall interactions between Gly, water solvent and the nanotube are not as stable as the interaction between this nanotube and water.

Computational Details
All calculations were carried out using the periodic ab-initio code [45].All the SCF calculations and geometry optimizations were performed using the B3LYP-D* density functional method, which includes an empirical a posteriori correction term proposed by Grimme [46] to account for dispersion forces (missed in the pure B3LYP [47,48] method), but whose initial parameterization (D) was modified for extended systems (D*) [49], to provide accurate results for the calculations of cohesive energies of molecular crystals and of adsorption processes within a periodic treatment [49][50][51].The adopted Gaussian functions consisted of an all electron triple-ζ 6-311G* standard basis set for the B and N atoms of the BNNTs and a TZP basis set from Ahlrichs and coworkers [52] for the atoms of Gly.This basis set combination has been proved to exhibit small basis set superposition error interaction energies [18,50].The shrinking factor of the reciprocal space net defining the mesh of k points in the irreducible Brillouin zone was set to 5, which requires diagonalizing the Hamiltonian matrix in 3 k points [53].The accuracy of both Coulomb and exchange series was set to values of overlap integrals of 10 −7 and 10 −16 , respectively, which ensure a very good numerical accuracy.A pruned (75, 974) grid has been used for the Gauss-Legendre and Lebedev quadrature schemes in the evaluation of functionals.The condition to achieve SCF convergence between two subsequent cycles was set to 10 −7 Eh.Full relaxations of both lattice parameters and internal atomic coordinates by means of analytical energy gradients [54][55][56] have been carried out.The geometry optimization was performed by means of a quasi-Newton algorithm in which the quadratic step (BFGS Hessian updating scheme) is combined with a linear one (parabolic fit) [57].

Conclusions
Periodic quantum mechanical calculations have been used to simulate the adsorption of glycine (Gly) on different zig-zag (n,0) single-walled boron-nitride nanotubes (BNNTs, n = 4, 6, 9 and 15) both in the gas-phase and in a microsolvated state (i.e., modeled by the presence of seven explicit water molecules) with the aim of determining the adsorption properties and the effect exerted by water as a function of surface curvature.These calculations are based on the B3LYP-D2* method, which includes the B3LYP hybrid functional with a revised version of the empirical a posteriori correction term (D2*) to account for dispersion interactions.
Gas-phase results clearly indicate that the most stable interaction between Gly and the (4,0), (6,0) and (9,0) BNNTs takes place through a covalent dative interaction between the NH 2 group of Gly and the B atom of the BNNTs, which produce charge transfers from Gly to the BNNTs.In contrast, the interaction between Gly and the (15,0) BNNT is mainly governed by non-covalent dispersive forces based on a π-stacking between the π systems of the Gly COOH group and the B-N hexagon rings of the nanotube.Remarkably, the energy difference between these two adducts decreases when increasing the BNNT radius, in line with the polar/apolar character of the considered nanotubes.
The adsorption of Gly on the BNNTs in the presence of seven water molecules has been studied adopting a progressive microsolvation procedure, in which water solvent molecules are added at the dry Gly/BNNT interfaces, hence progressively removing the direct interaction between Gly and the BNNTs.The obtained results indicate that for the (6,0), (9,0) and (15,0) BNNTs, the most stable microsolvated systems were found to exhibit no direct contact between the two partners; that is, the Gly/BNNT interfaces are fully bridged by the water solvent molecules.In contrast, for the (4,0) one of the most stable systems shows direct contact between Gly and the BNNT through an interaction between the Gly COO − group and a nanotube B atom, although entropic effects (not accounted for in this work) might favor water mediated interface.Further energetic results point out that the larger the BNNT radius, the less water affinity.Accordingly, for larger radius BNNTs, the interaction between water and Gly was found to be predominant, in detriment to their interaction with the BNNT.However, it is found that the (4,0) BNNT exhibits a large water affinity, which is reflected by the fact that the replacement of seven adsorbed water molecules by a microsolvated Gly has been found to be an unfavorable process.
The results presented here provide evidence that the adsorption properties of the BNNTs as well as their water affinity can significantly be modulated by controlling the tube diameter, as they are expected to exhibit different physico-chemical features, which may be of interest for the design of bioconjugated systems based on boron-nitride nanostructures and their potential bionanotechnological applications.

Figure 1 .
Figure 1.(a) B3LYP-D2* optimized structures of the different calculated adducts for glycine interacting with the considered boron nitride nanotubes (BNNTs) in the gas-phase; (b) B3LYP-D2* optimized structures for those complexes in which a spontaneous proton transfer occurs during the geometry optimization (see text).Distances in Å: bare values for the (4,0) BNNT; values in parenthesis for the (6,0) one; values in brackets for the (9,0) one; and italic underlined values for the (15,0) one.

Figure 2 .
Figure 2. B3LYP-D2* optimized structures of the different calculated complexes for glycine interacting with the (4,0), (6,0) and (9,0) BNNTs (a) and with the (15,0) BNNT (b) in the presence of seven water molecules.Distances in Å: bare values for the (4,0) BNNT; values in parenthesis for the (6,0) one; values in brackets for the (9,0) one; and italic underscored values for the (15,0) one.For this latter case, the distance is that between the C atom and the plane defined by the closest B-N hexagon ring.

Figure 3 .
Figure 3. B3LYP-D2* optimized structures of the different calculated complexes for the considered BNNTs in the presence of seven water molecules.Distances in Å: bare values for the (4,0) BNNT; values in parenthesis for the (6,0) one; values in brackets for the (9,0) one; and italic underscored values for the (15,0) one.For this latter case, the distance is that between the closest water O to the plane defined by the B-N hexagon ring.

Table 1 .
Calculated adsorption energies (∆E ads ), including the pure electronic energy contribution (∆E el ) and the contribution of dispersion (∆E D ).The relative electronic energies (∆E rel ) for a given Gly/BNNT system are also included.Values in units of kcal mol −1 .

Table 2 .
Muliken charge (Q) of Gly adsorbed on the BNNTs in the most stable adducts and respective calculated direct band gaps (E g ).