Understanding the Boron–Nitrogen Interaction and ... - ACS Publications

Oct 16, 2015 - Jianwei Sun,. ∥ ... Institute for Computational Molecular Science, Temple University, 1900 North 12th Street, Philadelphia, Pennsylva...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV LAVAL

Article

Understanding the Boron - Nitrogen Interaction and Its Possible Implications in Drug Design Hao Dong, Wei Li, Jianwei Sun, Shuhua Li, and Michael L. Klein J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.5b07783 • Publication Date (Web): 16 Oct 2015 Downloaded from http://pubs.acs.org on October 17, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Understanding the Boron - Nitrogen Interaction and Its Possible Implications in Drug Design Hao Dong,*,†,‡ Wei Li,‖ Jianwei Sun,§ Shuhua Li



and Michael L. Klein





Kuang Yaming Honors School, Nanjing University, P.R. China;



Institute for Computational Molecular Science, Temple University, 1900 North 12th

Street, Philadelphia, Pennsylvania 19122-6078, United States; ‖

School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic

Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, P.R. China; §

Department of Physics, Temple University, 1900 North 12th Street, Philadelphia,

Pennsylvania 19122-6078, United States

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT 2-Aminoethoxydiphenylborate (2-APB) is a broad-spectrum modulator of various membrane proteins. Specifically, it exhibits concentration dependent modulation of calcium signaling through store-operated calcium (SOC) channels: low micro-molar concentration of 2-APB stimulates SOC entry while higher concentration induces complete inhibition. ab initio quantum chemical calculations show the relative stability of the two major isomers of 2-APB (cyclic and extended), is about 8 kcal/mol. The dual functionality of 2-APB for SOC channels is thus likely associated with its ability to switch among isomeric forms, suited to different binding sites in the SOC channels with distinct binding affinities. Importantly, the moderate relative stability of different isomers results from a delicate balance between the intramolecular boron-nitrogen coordinate bond with strength about -45 kcal/mol, and ring strain engendered by cyclic oligomerization. The synergistic effect of these two factors likely makes 2-APB an ideal dual effect drug.

2

ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

INTRODUCTION 2-Aminoethoxydiphenylborate (2-APB), shown in Figure 1, is a membrane permeable reagent that has extensive interactions with membrane proteins, including InsP3 receptors,1 transient receptor potential (TRP) channels,2 specific gap junctions (composed of connexin26 and/or connexin32),3 SERCA Ca2+ pumps,4 voltage-gated potassium channels,5 and store-operated calcium (SOC) channels.6 Especially for the study of calcium signaling through SOC channels, 2-APB has been shown to exhibit concentration dependent modulation behavior: low micromolar concentration(1~20 μM) of 2-APB can stimulate the store-operated Ca2+ entry, while higher concentration (25~100 μM) induces a complete inhibition.6-8 This bimodal effect is a characteristic of 2-APB function on SOC entry (SOCE), though its regulation mechanism is still unclear.6 Boron and nitrogen atoms in 2-APB are typical components of coordinate bonding. The boron atom, which has an open shell and is electron deficient, acts as a Lewis acid and provides a vacant coordination site, while the nitrogen contributes an electron pair as a Lewis base.9 Replacing the terminal NH2 by a NHCH3 (monomethyl-APB) or a N(CH3)2 (dimethyl-APB) does not hurt its inhibitory activity.10 However, further substituting the boron by a carbon on dimethyl-APB result in the loss of inhibition on SOCE.11 The substitute of boron-oxygen core with carbon-phosphorus core induced the loss of potentiation capacity.12 In contrast, methyl-diethylborinate (MDEB), a 2-APB analogue containing the boron-oxygen core alone is capable of potentiating the SOCE by increasing the Ca2+ influx amplitude.13 Interestingly, trimethylborate (TMB), an analog of MDEB, containing the boron-oxygen core but with the two ethyl groups being replaced by methoxyl groups, was devoid of any regulation effect on SOCE.13 Even though it was proposed that neither the boron center is the prerequisite for SOCE inhibitory activity,10 nor does the aminoethyl group being critical for SOCE inhibition,12 and some compounds are more potent than 2-APB in inhibiting SOCE, the analogues based on the mother compound 2-APB having both functional groups -- the boron center and amine group -- were reported to exhibit this dual effect.12-14 We also want to mention that, as an 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

exception, another 2-APB analog, diphenylborinic anhydride (DPBA), without the aminoethyl group, has this bimodal effect as well.12 Seemingly, the structure-function relationship analysis of these modulators leads an elusive conclusion. Yet, 2-APB was usually used as a prototype compound for designing new modulators, and a number of 2-APB analogues have been synthesized and characterized, most of which containing both B and N centers.10, 14-16 Therefore, we focused our attention to the B…N interaction in the present work. The pKa of 2-APB was determined to be 9.6,17 suggesting that the protonated form is dominant in aqua phase at neutral pH condition. However, large down shifts in the pKa value of the amine group (up to 4 pH units) in a protein interior has been observed,18-19 mainly because of the low polarizability and polarity of the low dielectric constant environment relative to the aqua phase. Therefore, the amine group, the potential Lewis base, may adopt the deprotonated form as well. For example, 2-APB in dimethyl sulfoxide (DMSO), an aprotic solvent, was detected as a cyclic form where the amine is deprotonated to serve as the Lewis base.10 In addition, the formation of dimeric 2-APB, in which the deprotonated amine group on one 2-APB interacts with the boron center on a neighboring 2-APB and thus forming a 10-member ring, was already proposed based on the mass spectral analysis.20 Presumably, the population of the deprotonated 2-APB will increase in the protein interior or in lipid bilayers. More importantly, it is likely that the different functional abilities of 2-APB for activating or blocking the SOCE is associated with its ability to switch among distinct forms, suited to different binding sites in the channel with different binding affinities (see Figure 1), in which the dynamic bonding between boron and nitrogen plays a critical role. Here we report quantum mechanics calculations designed to explore the bonding between B and N in 2-APB. The purpose of the present work is two-fold: firstly, we want to understand the nature of B…N bonding and the flexibility of the molecule, which may be helpful to better design drugs that are capable of modulating membrane proteins. Secondly, for 2-APB, we want to identify a density functional that is both reliable and more computational efficient than the accurate 4

ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

but expensive post-Hartree-Fock (post-HF) methods. The former in turn, could then be used to parameterize a series of compounds containing the B…N bond suitable to investigate drug-channel interactions with molecular dynamics simulations.

COMPUTATIONAL DETAILS Based on the following considerations, the medium size basis set 6-311++G(d,p) was employed: 1) one aim of the present work is to identify an appropriate and efficient density functional theory (DFT) method suitable to describe B…N bonding, and the parameterization of a force field using a large training set, 2)the energy decomposition analysis (based on localized molecular orbitals)21 we used to explore the B…N bonding is basis set insensitive, and 3)the computational cost will increase significantly with the oligomeric state (for example, the number of basis functions is 1944 for the tetrameric structure at the 6-311++G(d,p) level). We also tested the basis set effect, and found that increasing of basis set does not change the general conclusions. Benchmark calculations were carried out with generalized energy-based fragmentation (GEBF)22 method at the level of explicitly correlated coupled-cluster singles and doubles with noniterative triples corrections [CCSD(T)-F12a],23 e.g. GEBF-CCSD(T)-F12a.24 To be specific, the structures were optimized at either the M06-2X25 or the second order Møller-Plesset perturbation (MP2) theory levels with the 6-311++G(d,p) basis set. The energy difference between two representative states (A1 and A2, Figure 1) was calculated at the GEBF-CCSD(T)-F12a level with the aug-cc-pVDZ basis set, and was taken as the reference for the subsequent calculations. In the GEBF calculations, the canonical HF energies were combined with the GEBF correlation energies. Due to the fast basis set convergence, the GEBF-CCSD(T)-F12a results with the aug-cc-pVDZ basis set are comparable to the corresponding high-level CCSD(T) results with aug-cc-pVQZ basis set.24 Then, for the sake of compromised accuracy and efficiency, several density functionals were benchmarked against this standard, including the pure functional M06L,26 the hybrid functionals B3LYP,27-28 X3LYP,29 M05,30-31 M05-2X,31 M06,25 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

M06-2X,25 M06-HF,32 the D2 version of Grimme’s dispersion corrected33 functionals B97D,34 ωB97XD,35 mPW2PLYPD,36 B2PLYPD,37 and the D3 version of dispersion corrected38 functionals TPSS-D3,39 revTPSS-D3,40 PBE-D3,41 PBE0-D3,42 BLYP-D3,27-28 B3LYP-D3,27-28 as well as the newly developed “meta-GGA made simple” and its hybrid (MGGA-MS2 and MGGA-MS2h)43-44 and the hybrid “meta-GGA made very simple” (MGGA-MVSh)45 functionals. By using each functional, the two states were optimized individually to calculate the energy difference. With the advantage of yielding quantitative agreement with high-level benchmark calculations, M06-2X is used in all of the following calculations (further described below). Due to the inherent difficulties in estimating intramolecular bonding, we built some model systems containing intermolecular B…N bonding (Figure 2) to mimic the intramolecular one in 2-APB. Some of them (for example, B1 and B2) have already been extensively studied by computation, however the results vary from different levels of theory.46-49 We further tested the performance of M06-2X by calculating the binding energy of model complexes. The interaction was characterized as follows:50  =  + 

(1)

where  is the interaction energy,  represents the distortion energy for each fragment to be promoted from isolation to the complex, and  is the instantaneous interaction energy between the two fragments within the complex. In other words,  and  are the energies of the fragment in either fully relaxed (equation 2) or fixed in the complex structure (equation 3). In the fragmentation scheme, the molecule was separated into two neutral fragments by cutting a C-C bond (homolytic bond cleavage), so the spin multiplicity of each fragment is 2.  is used for benchmark in this section, while  is calculated to estimate the ring steric effect, as it is compatible with the energy decomposition analysis (EDA) we used. ( ) = ( ) − ( )  − ( )

(2)

( ) = ( ) − ( ) − ( )

(3)

With each fragment relaxed (equation 2), the calculated binding enthalpies of B1 is -28.24 kcal/mol, which is quite close to the experimental measurement of -27.5 6

ACS Paragon Plus Environment

Page 7 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

kcal/mol,51 and comparable with those computationally expensive high-level results (for

example,

-27.5

kcal/mol

obtained

by

QCISD(T)/aug-cc-pV5Z//MP2/

6-311++G(3df,2p) calculation),52 but slightly higher than Haaland’s estimation of -31.1±1.0 kcal/mol.53 Similar result can be found for B2, where the binding enthalpies is -17.92 kcal/mol (calculation, present work) vs -17.6 kcal/mol in literature (experiment).51 Therefore, the M06-2X scheme is able to describe well the systems containing either the intramolecular or intermolecular B…N interaction. However, we stress that, even though M06-2X exhibited good performance in the present study, carefully calibration of this functional is necessary, as its construction depends on a relative large number of empirical parameters, which is likely to suffer from numerical instabilities.54 Moreover, the M06-2X parameterization is not good for transition metals,55 where the increased occurrence of metal-to-boron coordinate bonding complexes have been observed.56 For different oligomeric states, the cyclic form of the monomer was taken as the reference, and the relative stability of others were calculated by using equation 4, where n represents the oligomeric state: ()   =

() 

− (1)

(4)

Most of the present calculations were carried out with the Gaussian 09 program.57 For those functionals with empirical dispersion corrections based on Grimme’s scheme,58 GAMESS59 was used. The GEBF results were obtained with the LSQC program,60 in which the HF calculations were performed with Gaussian 0957 and the CCSD(T)-F12a calculations were carried out with the MOLPRO 2012 package.61-62 The electron localization function was generated with Multiwfn.63-64

RESULTS AND DISCUSSION Benchmark for functionals. We aim to identify a reliable and computationally efficiency functional, to accurately characterize B…N bonding in 2-APB. To do so, we choose two representative conformers of deprotonated 2-APB: the fully extended state (A1), in which the two atoms are far away from each other; and the cyclic state

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(A2), where B and N form close contact. Potential energy surface scan data shows that they are both local minimum (Figure S1 in Supporting Information). As shown in Figure 3, the reference data is in the range 7.65 to 8.08 kcal/mol, which is insensitive to the subtle difference of the structures optimized with different methods (MP2 or M06-2X). The corresponding data calculated with M06-2X is 8.01 kcal/mol, which is quite close to the reference. This is consistent with the results reported in literature, where by using experimental binding enthalpies as the target, the M06-2X was found to show good performance in describing substituted B…N complexes.47 Besides the M06-2X, some D3-corrected functionals including revTPSS-D3, PBE0-D3, TPSS-D3, PBE-D3 as well as M05-2X, MGGA-MS2h and MGGA-MS2, show reliable performance in the present case as well. With the advantage of yielding quantitative agreement with high-level benchmark calculations, M06-2X is used in all of the following calculations. It is worth pointing out that, the MP2 and some DFT functionals are incapable of describing well the energy difference in the present system. The stabilization energy predicted with MP2 (10.57 kcal/mol) is 2.49 kcal/mol higher than the CCSD(T) value; an overestimation error of 31%. One of the explanations is that the MP2 method overestimates the dispersion energy in some cases.65-66 Interestingly, the good performance of B3LYP in other systems does not carry over to the present system, mainly because of its difficulties to characterize pericyclic reactions.67 Even with the D3 correction, B3LYP-D3 still does not lead to a reliable prediction. In contrast, TPSS-D3 was found to have good performance to describe cyclophanes with notable strain components.68

The strength of the B…N bonding in 2-APB. The calculated B-N distance in the cyclic form of 2-APB is 1.73 Å, which is slightly longer than the x-ray structure of 1.65 Å.69-70 The packing in the crystal environment may partially account for this discrepancy. It should be noted that this bond length is significantly longer than the B-N bond of 1.57 Å71 in cubic or 1.45 Å72 in hexagonal boron nitride, but much shorter than the sum of the van der Waals radii of the two atoms (2.90 Å), indicating the 8

ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

characteristics of a coordinative bond. The calculated instantaneous interaction energy of B1, ( ) , with each fragment remaining at the geometry of the complex (equation 3), is -44.08 kcal/mol at the M06-2X/6-311++G(d,p) level, which is consistent with the high level CCSD(T)/aug-cc-pVQZ data of -44.16 kcal/mol reported in literature.21 Being the simplest model system containing B…N bonding, its interaction energy mainly comes from the formation of the Lewis pair. We further realized that all of the model complexes (B1~B5) have similar instantaneous interaction energies, regardless of the chemical environment of the B…N group (Figure 2). How about the B…N bonding in 2-APB? By comparing the physical descriptors of A2 and model complexes with respect to the electron localization function (ELF), molecular orbital (MO) and electrostatic potential surface (EPS), we will demonstrate that the intramolecular B…N bonding in 2-APB resembles the intermolecular one in model complexes, as further described below (Figure 4). ELF: ELF is a useful tool in interpreting bonding.73-74 According to the definition, ELF = 1 corresponds to a fully localized state and ELF = ½ represents the electron gas.73 As shown in Figure 4a, the x-axis is the coordinate of the B…N bonding with the center of the boron atom being at the origin. The location of the single attractor V(B,N) corresponding to bonding is at ~1 Å position away from the boron atom, which confirms that the bonding pair is more closely associated with the electron donor N atom. The ELF value there is ~0.943, almost the same in B1~B5 and in A2, indicating the well resembled B…N bonding in both intra- and inter-molecular systems. MO: The boron atom is approximating the idealized configuration for sp3 hybridization, that is (2s)0.51(2p)1.54 valence configuration in A2 and (2s)0.51(2p)1.55 in B5. The valance orbital hybridization of boron is consistent with its approximately tetrahedral geometry. The frontier orbital picture shows that HOMO and LUMO orbitals in B5 (Figures 4b and 4b') are quite similar to those in A2 (Figures 4c and 4c'). Specifically, HOMO is delocalized over the two phenyl rings as well as the central boron atom, the Lewis acid part. The partial overlap between the vacant p orbital of 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the boron center and the π orbitals of the two phenyl groups results in a pronounced Lewis acidity. When the two phenyl groups in 2-APB are substituted with two methyl group, the relative stability between the cyclic- and extended-forms drops significantly from 8.01 kcal/mol in 2-APB to 3.67 kcal/mol, indicating a weakened B…N interaction. This may explain the observations that the diphenyl moiety is important for its activity by comparing 2-APB with its structural analogs.11 Meantime, nitrogen on the amine head, the Lewis base part, is also sp3 hybridization, with the configuration of (2s)1.37(2p)4.38 in A2 and (2s)1.36(2p)4.37 in B5. The HOMO-LUMO gaps are similar in A2 and B5, which are -6.16 and -6.34 kcal/mol, respectively. Notably, the delocalization of molecular orbitals on the Lewis acid and Lewis base centers, respectively, is an indication of the non-covalent feature of the B…N bond in 2-APB. Intramolecular electron transfer (~0.1 e) takes place from N on amine head to B, confirming the formation of the adduct, and suggesting that polarization (coming from the electron transfer and charge accumulation between nuclei) could be important for the formation of the Lewis pair. EPS: The molecular electrostatic potential is a rigorously defined property, which can be characterized with quantum mechanics. The calculated EPSs of B5 (Figure 4d) and A2 (Figure 4e) resemble each other, suggesting the similar B…N bonding within the two constructs: the amine group is predominantly positive, and the boron center is almost neutral, while the neighboring oxygen atom shows strong negative distribution of electrostatic potential. Therefore, the boron center is not electrostatically complementary to the amine group even though the B…N bond is polar. However, with the presence of oxygen, the electrostatic component may still be important in such bonding, though it may not be a dominant factor. This may explain the experimental observation that replacement of boron-oxygen core with carbon-phosphorus core leads the loss of potentiation capacity.12 All of the above three independent descriptors indicate that the intramolecular B…N bonding in 2-APB resembles the intermolecular one in model complexes, among which, B5 is the best model to mimic the environment in A2. Therefore, though we do not have a direct estimation on the B…N bond strength, it is a rational speculation 10

ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

that the instantaneous interaction energy of B…N in A2 is -46.75 kcal/mol as well, which is much stronger than those non-bonded interactions (no stronger than -10 kcal/mol in most cases), but notable weaker than the covalent B-N single bond (-88 kcal/mol, 1.47 Å).53

The nature of the B…N bonding in 2-APB It is interesting to understand the nature of B…N bonding in 2-APB, which has strength intermediate between covalent and non-covalent bonds. Among different protocols, EDA was developed as a powerful tool for the quantitative interpretation of chemical bonds.75 Here we used the localized molecular orbital based EDA protocol21 implemented in GAMESS59 to explore the B…N bonding. Besides model complexes B1~B5, HO-CH2-CH2-NH2 (B6, aiming at the C-C bond) and H2O…H2O (B7, aiming at the hydrogen bond) were used as the representatives for covalent and non-covalent bonding models. The individual contributions from the total interaction energies, including electrostatic, exchange, repulsion, and polarization terms, were normalized for direct comparison. As seen from Figure 5, the contributions of each component towards attractions are quite different among the complexes B6, B7 and B1~B5. Polarization term, accounting for the additional electrostatic stabilization due to the induced multipoles as well as the charge accumulation between nuclei, is dominant in the covalent bond (B6), which counts for ~33% of the attractive energy, but only contributes for ~10% for the non-covalent bond (B7). B…N bonding in B1~B5 has a moderate polarization component (25%) in between of the covalent and non-covalent bond. Electrostatic term is dominant in non-covalent bond (B7), and is decreasing in those B1~B5 complexes, and is further decreasing in the covalent bond (B6). Interestingly, contributions from dispersion are more important in the non-covalent bond than others. All of the three components, including polarization, electrostatic, and dispersion, vary monotonically from B6 to B1~B5, and to B7, indicating that the B…N bonding falls in between of that covalent and non-covalent bond. By calculating the localized resonance structures and the corresponding weighting factors based on the 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 28

natural resonance theory, it was proposed that the B…N bond in B1 complex is 65% ionic.52

Ring strain Given the significant contribution from the instantaneous intramolecular B…N interaction (-46.75 kcal/mol), why is the cyclic form of 2-APB only 8.08 kcal/mol more stable than the extended form? When compared to the extended form, it is very likely that strain is involved in ring formation in the cyclic form and that this mainly accounts for the “extra” instability. So we propose the following energy partition protocol (protocol-I) to roughly estimate the energy of ring strain. By breaking the C-C covalent bond to divide the 2-APB into two fragments with either B or N atoms (two radical monomers), as shown in Figure 6, we get the instantaneous interaction energy between two fragments of -129.96 and -104.47 kcal/mol, respectively, for the cyclic- and extended-forms. We further propose equations (5) and (6) to partition the interaction energy into different components, including the steric of the ring, where the superscript "CS" and "ES" mean the cyclicand extended-state, respectively. Eint, EB…N, EC-C and Ering are the instantaneous interaction energy between two fragments, the B…N interaction energy, the C-C bond strength, and steric energy of the ring, respectively. The Eerror term is introduced to account for the error in the energy partition protocol we proposed here.      =  …" + # +  + 

(5)

    =  …" + # +  +  

(6)

   By using (5)-(6), and assuming  =  , and both 

and  …" equal

to 0 as there is no ring steric energy and B…N bonding in the fully extended structure, we get:       = ( −  ) −  …" − (# − # )

(7)

Eint could be calculated with EDA, and we already have an estimation on  …" . So

12

ACS Paragon Plus Environment

Page 13 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

once the # is known, we are able to evaluate the ring steric. To estimate # , another model containing this covalent C-C bond is constructed, as shown in Figure 6.  Taking # as an example, we built A1CC by keeping the fragment (O-CH2-CH2-NH2)

as in the original structure of A1. A hydrogen atom is added to saturate the broken bond at oxygen atom, and its position is freely optimized. We then estimate the C-C  in the same way. The bond strength by breaking the bond and calculate the #   calculated # and # are -107.15 kcal/mol in A1CC and -102.87 kcal/mol in

A2CC, respectively. By using this energy partition protocol, the ring strain in the cyclic form of 2-APB (A2) is estimated to be 16.98 kcal/mol. Similarly, we construct analogs (C1 and C2) to 2-APB (A1 and A2, respectively) with one more CH2 group inserted (Figure 6) to make a longer chain. After structure optimization, the C1CC and C2CC are also constructed to estimate the # in the analogs (Figure 6). The ring strain in C2 is then estimated as 12.25 kcal/mol. Therefore, increasing from 5-members A2 (2-APB) to 6-members C2 (2-APB analog) ring is estimated to gain 4.73 kcal/mol energy. To further validate the above protocol-I to estimate ring strain, we calculate the relative stability of the cyclic- and extended-forms (energy difference protocol, protocol-II). The energy difference (∆E5-members) in 2-APB is 8.01 kcal/mol, and the cyclic form of the aforementioned analog is 12.57 kcal/mol (∆E6-members) lower than its extended form, mainly because of the released ring strain. Therefore, the energy difference of ring strain (∆∆E = ∆E6-members - ∆E5-members) is 4.56 kcal/mol. The consistency of data obtained from both protocol-I and II justifies the effect of ring strain we proposed above. However, by using protocol-II, we can only estimate the relative steric energy between different ring sizes (for example, 4.56 kcal/mol between 5- and 6-members rings). In contrast, protocol-I provides an opportunity to quantitatively evaluate the ring steric effect in each specific system (for example, 16.98 kcal/mol in the 5-members ring of 2-APB, and 12.25 kcal/mol in its 6-members analog). It is clear that the ring steric effect partially counterbalance the strong attraction from the B...N bonding, and makes it more dynamic. The reversible and dissociative 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

characteristics of such bonding was already observed using fs kinetic-energy resolved time-of-light mass spectrometry.76 To quantitatively show this offset, we manually change the B…N distance in B1, B5 (model compounds) and A2 (cyclic 2-APB), and calculate the relative stability by using the energy minimum of each as the reference. The former two show quite similar disassociation curves (Figure 7). The stability of the cyclic 2-APB, on the other hand, decreases slower than B1 and B5. Clearly, the release of ring strain in the cyclic 2-APB partially compensates the energy lost from B…N separation. Therefore, there is a delicate balance between ring strain and B…N bonding in the cyclic form of 2-APB, and the reversible forming and breaking of the B…N bond is possible, as disclosed by the only 8.08 kcal/mol energy difference between the cyclic- and the extended-forms. Consequently, 2-APB is free to switch among different functional forms, and may have different binding affinities in the target channels.

Oligomerization Given the ~45 kcal/mol instantaneous strength of the B…N bonding in 2-APB, it is very likely such favorable interaction can be formed either intra- or inter-molecularly. One of the potential advantages of forming a molecular complex that consists of a few monomer units is to form a more stable cluster with higher molecular volume, which may have distinct modulating behavior on the target membrane protein. Formation of higher oligomeric structures gradually releases the ring strain and therefore makes the B...N bonding more favorable, as indicated by the bond length. It is 1.73 Å in the monomer (A2), and the averaged values are 1.67 Å in the dimer (A3), 1.65 Å in the trimer (A4) and 1.63 Å in the tetramer (A5), respectively. Consequently, aggregation of the monomers increase the stability of the complex (tetramer > trimer > dimer > monomer, Figure 8). For example, the dimer is 8.8 kcal/mol per monomer lower than the monomer. However, with the increased oligomerization status, the favorable entropy changes may partly offset by the unfavorable enthalpy changes, and the steric repulsion from the bulky phenyl groups between different units is also introduced. Therefore, the relative stability of 14

ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

oligomers with increasing number of monomers converges rapidly.

Possible 2-APB modulation mechanism for CRAC channel gating The detailed information about 2-APB regulation is still unclear. Presumably, the stimulatory effect of 2-APB on SOCE comes from its induced pore dilation. Specially, the interaction between the extended form of 2-APB (so that the B and/or N centers are free to bind to the protein) and Orai (at the N- and/or C-terminals) facilities the opening of the central pore.77 Another plausible scenario is that the presence of 2-APB (or its analog) helps the formation of the Orai-STIM complex.16, 78-79 With regard to inhibition, an intuitive picture is that the volume-occupied cyclic 2-APB or its oligomer with larger size (which is more likely, because of the higher stability as shown in Figure 8) physically blocks the permeation pathway of ions through the channel, while other cases, including the intracellular acidification mechanism,80 are also possible. All of the above speculations are subject to further validations from both molecular modeling and experimental measurements.

CONCLUSIONS In this work, we propose that the B…N interaction, a reversible coordinate bonding in 2-APB, is critical for its known dual pharmaceutical effect. Importantly, the character of the intramolecular B…N bonding, with respective to the electron localization function, molecular orbitals and electrostatic potential surface, resembles that found in complexes containing intermolecular B…N bonding. The instantaneous strength of the B…N bond is 45 kcal/mol, which lies between a covalent and a non-covalent bond. The delicate balance between B…N bonding and ring steric strain makes 2-APB in principle an ideal drug to present distinct forms and/or binding affinities at different binding sites. Seemingly facile formation and easy dissociation of oligomers is the key factor enabling 2-APB to exhibit broad-spectrum activity on membrane proteins. Understanding the subtlety of oligomer formation and stability obtained herein should provide new insights for generating novel, improved clinical candidates for the CRAC channel. 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 28

Supporting Information 2D potential energy surfaces of the cyclic- and extended-forms of 2-APB (A1 and A2). This material is available free of charge via the Internet at http://pubs.acs.org.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]

Notes The authors declare no competing financial interest.

ACKNOWLEDGEMENTS This work was supported by the National Natural Science Foundation of China (Grant No. 21503107) and the "Fundamental Research Funds for the Central Universities". Part of the calculations were performed using resources on IBM Blade cluster system from the High Performance Computing Center (HPCC) of Nanjing University.

REFERENCES (1) Maruyama, T.; Kanaji, T.; Nakade, S.; Kanno, T.; Mikoshiba, K. 2APB, 2-Aminoethoxydiphenyl Borate, a Membrane-Penetrable Modulator of Ins(1,4,5)P-3-Induced Ca2+ Release. J. Biochem. 1997, 122, 498-505. (2) Hu, H. Z.; Gu, Q. H.; Wang, C. B.; Colton, C. K.; Tang, J. S.; Kinoshita-Kawada, M.; Lee, L. Y.; Wood, J. D.; Zhu, M. X. 2-Aminoethoxydiphenyl Borate Is a Common Activator of TRPV1, TRPV2, and TRPV3. J. Biol. Chem. 2004, 279, 35741-35748. (3) Tao, L.; Harris, A. L. 2-Aminoethoxydiphenyl Borate Directly Inhibits Channels Composed of Connexin26 and/or Connexin32. Mol. Pharmacol. 2007, 71, 570-579. (4) Missiaen, L.; Callewaert, G.; De Smedt, H.; Parys, J. B. 2-Aminoethoxydiphenyl Borate Affects the 2+

2+

Inositol 1,4,5-Trisphosphate Receptor, the Intracellular Ca Pump and the Non-Specific Ca Leak from the Non-Mitochondrial Ca2+ Stores in Permeabilized A7r5 Cells. Cell Calcium 2001, 29, 111-116. (5) Wang, Y.; Deshpande, M.; Payne, R. 2-Aminoethoxydiphenyl Borate Inhibits Phototransduction and Blocks Voltage-Gated Potassium Channels in Limulus Ventral Photoreceptors. Cell Calcium 2002, 32, 209-216. 2+

2+

(6) Prakriya, M.; Lewis, R. S. Potentiation and Inhibition of Ca Release-Activated Ca Channels by 2-Aminoethyldiphenyl Borate (2-APB) Occurs Independently of IP3 Receptors. J. Physiol. (Lond.) 2001,

16

ACS Paragon Plus Environment

Page 17 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

536, 3-19. (7) Ma, H. T.; Patterson, R. L.; van Rossum, D. B.; Birnbaumer, L.; Mikoshiba, K.; Gill, D. L. 2+

Requirement of the Inositol Trisphosphate Receptor for Activation of Store-Operated Ca Channels. Science 2000, 287, 1647-1651. (8) Ma, H. T.; Venkatachalam, K.; Parys, J. B.; Gill, D. L. Modification of Store-Operated Channel Coupling and Inositol Trisphosphate Receptor Function by 2-Aminoethoxydiphenyl Borate in Dt40 Lymphocytes. J. Biol. Chem. 2002, 277, 6915-6922. (9) Staubitz, A.; Robertson, A. P.; Sloan, M. E.; Manners, I. Amine- and Phosphine-Borane Adducts: New Interest in Old Molecules. Chem. Rev. 2010, 110, 4023-4078. (10) Dobrydneva, Y.; Abelt, C. J.; Dovel, B.; Thadigiri, C. M.; Williams, R. L.; Blackmore, P. F. 2-Aminoethoxydiphenyl Borate as a Prototype Drug for a Group of Structurally Related Calcium Channel Blockers in Human Platelets. Mol. Pharmacol. 2006, 69, 247-256. (11) Dobrydneva, Y.; Blackmore, P. 2-Aminoethoxydiphenyl Borate Directly Inhibits Store-Operated Calcium Entry Channels in Human Platelets. Mol. Pharmacol. 2001, 60, 541-552. (12) Dellis, O.; Mercier, P.; Chomienne, C. The Boron-Oxygen Core of Borinate Esters Is Responsible for the Store-Operated Calcium Entry Potentiation Ability. BMC Pharmacol. 2011, 11, 1. (13) Djillani, A.; Doignon, I.; Luyten, T.; Lamkhioued, B.; Gangloff, S. C.; Parys, J. B.; Nusse, O.; Chomienne, C.; Dellis, O. Potentiation of the Store-Operated Calcium Entry (SOCE) Induces Phytohemagglutinin-Activated Jurkat T Cell Apoptosis. Cell Calcium 2015, 58, 171-185. (14) Goto, J.; Suzuki, A. Z.; Ozaki, S.; Matsumoto, N.; Nakamura, T.; Ebisui, E.; Fleig, A.; Penner, R.; Mikoshiba, K. Two Novel 2-Aminoethyl Diphenylborinate (2-APB) Analogues Differentially Activate and Inhibit Store-Operated Ca(2+) Entry via STIM Proteins. Cell Calcium 2010, 47, 1-10. (15) Ozaki, S.; Suzuki, A. Z.; Bauer, P. O.; Ebisui, E.; Mikoshiba, K. 2-Aminoethyl Diphenylborinate (2-Apb) Analogues: Regulation of Ca2+ Signaling. Biochem. Biophys. Res. Commun. 2013, 441, 286-290. (16) Hendron, E.; Wang, X. Z.; Zhou, Y. D.; Cai, X. Y.; Goto, J. I.; Mikoshiba, K.; Baba, Y.; Kurosaki, T.; Wang, Y. J.; Gill, D. L. Potent Functional Uncoupling between STIM1 and Orai1 by Dimeric 2-Aminodiphenyl Borinate Analogs. Cell Calcium 2014, 56, 482-492. 2+

(17) Bilmen, J. G.; Wootton, L. L.; Godfrey, R. E.; Smart, O. S.; Michelangeli, F. Inhibition of Serca Ca

Pumps by 2-Aminoethoxydiphenyl Borate (2-APB) - 2-APB Reduces Both Ca2+ Binding and Phosphoryl Transfer from Atp, by Interfering with the Pathway Leading to the Ca2+-Binding Sites. Eur. J. Biochem. 2002, 269, 3678-3687. (18) Isom, D. G.; Cannon, B. R.; Castaneda, C. A.; Robinson, A.; Bertrand, G. M. E. High Tolerance for Ionizable Residues in the Hydrophobic Interior of Proteins. Proc. Natl. Acad. Sci. USA 2008, 105, 17784-17788. (19) Isom, D. G.; Castaneda, C. A.; Cannon, B. R.; Garcia-Moreno, B. Large Shifts in Pka Values of Lysine Residues Buried inside a Protein. Proc. Natl. Acad. Sci. USA 2011, 108, 5260-5265. (20) van Rossum, D. B.; Patterson, R. L.; Ma, H. T.; Gill, D. L. Ca2+ Entry Mediated by Store Depletion, S-Nitrosylation, and TRP3 Channels. Comparison of Coupling and Function. J. Biol. Chem. 2000, 275, 28562-28568. (21) Su, P.; Li, H. Energy Decomposition Analysis of Covalent Bonds and Intermolecular Interactions. J. Chem. Phys. 2009, 131, 014102. (22) Li, S. H.; Li, W.; Ma, J. Generalized Energy-Based Fragmentation Approach and Its Applications to Macromolecules and Molecular Aggregates. Acc. Chem. Res. 2014, 47, 2712-2720. (23) Werner, H. J.; Schutz, M. An Efficient Local Coupled Cluster Method for Accurate

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 28

Thermochemistry of Large Systems. J. Chem. Phys. 2011, 135, 144116 (24) Wang, K. D.; Li, W.; Li, S. H. Generalized Energy-Based Fragmentation CCSD(T)-F12a Method and Application to the Relative Energies of Water Clusters (H2O)(20). J. Chem. Theory Comput. 2014, 10, 1546-1553. (25) Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215-241. (26) Zhao, Y.; Truhlar, D. G. A New Local Density Functional for Main-Group Thermochemistry, Transition Metal Bonding, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2006, 125, 194101. (27) Becke,

A.

D.

Density-Functional

Exchange-Energy

Approximation

with

Correct

Asymptotic-Behavior. Phys. Rev. A 1988, 38, 3098-3100. (28) Lee, C. T.; Yang, W. T.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron-Density. Phys. Rev. B 1988, 37, 785-789. (29) Xu, X.; Goddard, W. A. The X3lyp Extended Density Functional for Accurate Descriptions of Nonbond Interactions, Spin States, and Thermochemical Properties. Proc. Natl. Acad. Sci. USA 2004, 101, 2673-2677. (30) Zhao, Y.; Schultz, N. E.; Truhlar, D. G. Exchange-Correlation Functional with Broad Accuracy for Metallic and Nonmetallic Compounds, Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2005, 123, 161103. (31) Zhao, Y.; Schultz, N. E.; Truhlar, D. G. Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Theory Comput. 2006, 2, 364-382. (32) Zhao, Y.; Truhlar, D. G. Density Functional for Spectroscopy: No Long-Range Self-Interaction Error, Good Performance for Rydberg and Charge-Transfer States, and Better Performance on Average Than B3LYP for Ground States. J. Phys. Chem. A 2006, 110, 13126-13130. (33) Grimme, S. Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787-1799. (34) Becke,

A.

D.

Density-Functional

Thermochemistry

.5.

Systematic

Optimization

of

Exchange-Correlation Functionals. J. Chem. Phys. 1997, 107, 8554-8560. (35) Chai, J. D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615-6620. (36) Schwabe, T.; Grimme, S. Towards Chemical Accuracy for the Thermodynamics of Large Molecules: New Hybrid Density Functionals Including Non-Local Correlation Effects. Phys. Chem. Chem. Phys. 2006, 8, 4398-4401. (37) Schwabe, T.; Grimme, S. Double-Hybrid Density Functionals with Long-Range Dispersion Corrections: Higher Accuracy and Extended Applicability. Phys. Chem. Chem. Phys. 2007, 9, 3397-3406. (38) Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456-1465. (39) Tao, J.; Perdew, J. P.; Staroverov, V. N.; Scuseria, G. E. Climbing the Density Functional Ladder: Nonempirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids. Phys. Rev. Lett. 2003, 91, 146401. (40) Perdew, J. P.; Ruzsinszky, A.; Csonka, G. I.; Constantin, L. A.; Sun, J. Workhorse Semilocal Density

18

ACS Paragon Plus Environment

Page 19 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Functional for Condensed Matter Physics and Quantum Chemistry. Phys. Rev. Lett. 2009, 103, 026403. (41) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (42) Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158-6170. (43) Sun, J. W.; Xiao, B.; Ruzsinszky, A. Communication: Effect of the Orbital-Overlap Dependence in the Meta Generalized Gradient Approximation. J. Chem. Phys. 2012, 137, 051101. (44) Sun, J. W.; Haunschild, R.; Xiao, B.; Bulik, I. W.; Scuseria, G. E.; Perdew, J. P. Semilocal and Hybrid Meta-Generalized Gradient Approximations Based on the Understanding of the Kinetic-Energy-Density Dependence. J. Chem. Phys. 2013, 138, 044113. (45) Sun, J. W.; Perdew, J. P.; Ruzsinszky, A. Semilocal Density Functional Obeying a Strongly Tightened Bound for Exchange. Proc. Natl. Acad. Sci. USA 2015, 112, 685-689. (46) LeTourneau, H. A.; Birsch, R. E.; Korbeck, G.; Radkiewicz-Poutsma, J. L. Study of the Dative Bond in 2-Aminoethoxydiphenyl Borate at Various Levels of Theory: Another Poor Performance of the B3LYP Method for B-N Dative Bounds. J. Phys. Chem. A 2005, 109, 12014-12019. (47) Plumley, J. A.; Evanseck, J. D. Hybrid Meta-Generalized Gradient Functional Modeling of Boron-Nitrogen Coordinate Covalent Bonds. J. Chem. Theory Comput. 2008, 4, 1249-1253. (48) Huang, F.; Jiang, J. L.; Wen, M. W.; Wang, Z. X. Assessing the Performance of Commonly Used DFT Functionals in Studying the Chemistry of Frustrated Lewis Pairs. J. Theor. Comput. Chem. 2014, 13, 1350074. (49) Mitoraj, M. P. Bonding in Ammonia Borane: An Analysis Based on the Natural Orbitals for Chemical Valence and the Extended Transition State Method (ETS-NOCV). J. Phys. Chem. A 2011, 115, 14708-14716. (50) Bessac, F.; Frenking, G. Chemical Bonding in Phosphane and Amine Complexes of Main Group Elements and Transition Metals. Inorg. Chem. 2006, 45, 6956-6964. (51) Brown, H. C.; Bartholomay, H.; Taylor, M. D. Acid-Base Studies in Gaseous Systems Ii the Anomalous Base Strength of the Methylamines, a New Manifestation of Steric Strain. J. Am. Chem. Soc. 1944, 66, 435-442. (52) Plumley, J. A.; Evanseck, J. D. Covalent and Ionic Nature of the Dative Bond and Account of Accurate Ammonia Borane Binding Enthalpies. J. Phys. Chem. A 2007, 111, 13472-13483. (53) Haaland, A. Covalent Versus Dative Bonds to Main Group Metals, a Useful Distinction. Angew. Chem. Int. Ed. 1989, 28, 992-1007. (54) Wheeler, S. E.; Houk, K. N. Integration Grid Errors for Meta-GGA-Predicted Reaction Energies: Origin of Grid Errors for the M06 Suite of Functionals. J. Chem. Theory Comput. 2010, 6, 395-404. (55) Zhao, Y.; Truhlar, D. G. Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157-167. (56) Braunschweig, H.; Dewhurst, R. D. Transition Metals as Lewis Bases: "Z-Type" Boron Ligands and Metal-to-Boron Dative Bonding. Dalton Trans. 2011, 40, 549-558. (57) Frisch, M. J.; Trucks, G. W. S., H. B. ; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H., et al. Gaussian 09, Revision C.01; Gaussian, Inc.; Wallingford, CT, 2010. (58) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(59) Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.; Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S. J., et al. General Atomic and Molecular Electronic-Structure System. J. Comput. Chem. 1993, 14, 1347-1363. (60) Li, W.; Chen, C. H.; Zhao, D. B.; Li, S. H. Lsqc: Low Scaling Quantum Chemistry Program. Int. J. Quantum. Chem. 2015, 115, 641-646. (61) Werner, H. J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schutz, M. Molpro: A General-Purpose Quantum Chemistry Program Package. Wiley Interdiscipl. Rev. - Comput. Mol. Sci. 2012, 2, 242-253. (62) Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G., et al. MOLPRO, version 2012.1, a package of ab initio programs, 2012. (63) Lu, T.; Chen, F. W. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem. 2012, 33, 580-592. (64) Lu, T.; Chen, F. W. Quantitative Analysis of Molecular Surface Based on Improved Marching Tetrahedra Algorithm. J. Mol. Graphics Model. 2012, 38, 314-323. (65) Hobza, P.; Selzle, H. L.; Schlag, E. W. Potential Energy Surface for the Benzene Dimer. Results of ab Initio CCSD(T) Calculations Show Two Nearly Isoenergetic Structures: T-Shaped and Parallel-Displaced. J. Phys. Chem. 1996, 100, 18790-18794. (66) Sinnokrot, M. O.; Sherrill, C. D. Highly Accurate Coupled Cluster Potential Energy Curves for the Benzene Dimer: Sandwich, T-Shaped, and Parallel-Displaced Configurations. J. Phys. Chem. A 2004, 108, 10200-10207. (67) Hao, P.; Sun, J. W.; Xiao, B.; Ruzsinszky, A.; Csonka, G. I.; Tao, J. M.; Glindmeyer, S.; Perdew, J. P. Performance of Meta-GGA Functionals on General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions. J. Chem. Theory Comput. 2013, 9, 355-363. (68) Grimme, S.; Muck-Lichtenfeld, C. Accurate Computation of Structures and Strain Energies of Cyclophanes with Modern DFT Methods. Isr. J. Chem. 2012, 52, 180-192. (69) Rettig, S. J.; Trotter, J. Crystal and Molecular-Structure of B,B-Diphenylboroxazolidine (2-Aminoethyl Diphenylborinate), Ph2BO(CH2)2NH2. Can. J. Chem. 1973, 51, 1288-1294. (70) Rettig, S. J.; Trotter, J. Crystal and Molecular-Structure of B,B-Bis(P-Tolyl)Boroxazolidine and Orthorhombic Form of B,B-Diphenylboroxazolidine. Can. J. Chem. 1976, 54, 3130-3141. (71) Fermi, E. A Statistical Method for Determining Some Properties of the Atoms and Its Application to the Theory of the Periodic Table of Elements. Zeitschrift Fur Physik 1928, 48, 73-79. (72) Wirtz, L.; Rubio, A.; de la Concha, R. A.; Loiseau, A. ab Initio Calculations of the Lattice Dynamics of Boron Nitride Nanotubes. Phys. Rev. B 2003, 68, 045425. (73) Becke, A. D.; Edgecombe, K. E. A Simple Measure of Electron Localization in Atomic and Molecular-Systems. J. Chem. Phys. 1990, 92, 5397-5403. (74) Silvi, B.; Savin, A. Classification of Chemical-Bonds Based on Topological Analysis of Electron Localization Functions. Nature 1994, 371, 683-686. (75) von Hopffgarten, M.; Frenking, G. Energy Decomposition Analysis. Wiley Interdiscipl. Rev. Comput. Mol. Sci. 2012, 2, 43-62. (76) Zhong, D. P.; Zewail, A. H. Femtosecond Dynamics of Dative Bonding: Concepts of Reversible and Dissociative Electron Transfer Reactions. Proc. Natl. Acad. Sci. USA 1999, 96, 2602-2607. (77) Amcheslavsky, A.; Safrina, O.; Cahalan, M. D. State-Dependent Block of Orai3 TM1 and TM3 Cysteine Mutants: Insights into 2-APB Activation. J. Gen. Physiol. 2014, 143, 621-631. (78) Wang, X.; Wang, Y.; Zhou, Y.; Hendron, E.; Mancarella, S.; Andrake, M. D.; Rothberg, B. S.; Soboloff, J.; Gill, D. L. Distinct Orai-Coupling Domains in STIM1 and STIM2 Define the Orai-Activating Site. Nat.

20

ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Commun. 2014, 5, 3183. (79) Wang, Y.; Deng, X.; Zhou, Y.; Hendron, E.; Mancarella, S.; Ritchie, M. F.; Tang, X. D.; Baba, Y.; Kurosaki, T.; Mori, Y., et al. STIM Protein Coupling in the Activation of Orai Channels. Proc. Natl. Acad. Sci. USA 2009, 106, 7391-7396. (80) Chokshi, R.; Fruasaha, P.; Kozak, J. A. 2-Aminoethyl Diphenyl Borinate (2-APB) Inhibits TRPM7 Channels through an Intracellular Acidification Mechanism. Channels (Austin) 2012, 6, 362-369.

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Different oligomeric states of 2-APB (left) and its potential binding sites (right) in the SOC channel.

Figure 2. Shown are selected model complexes exhibiting intermolecular B…N bonding. The calculated interaction energies between pairs of fragments are quite similar for each of these dimeric systems. Energy is shown in kcal/mol.

22

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. Shown are the estimated energy differences between the cyclic- and extended- forms of 2-APB. Using GEBF-CCSD(T)-F12a/aug-cc-pVDZ data as the benchmark, different functional as well as MP2 are tested at the 6-311++G(d,p) level. M06-2X gives the most accurate result in this case, and some other functionals highlighted in the red frame also show good performance. Ref[MP2] and Ref[M06-2X] are benchmark calculations based on the MP2 and M06-2X optimized structures, respectively. Energy is shown in kcal/mol.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Physical descriptors show high similarity between the inter- and intra-molecular B…N bonding, indicating the medium strength of the B…N interaction in 2-APB. (a) 1-D electron localization function along the B…N bond in model complexes B1~B5 (black, yellow, green, blue and red lines) and 2-APB (pink line). The HOMO (b) and LUMO (b') orbitals in B5, which mimics the 2-APB (A2). The molecular orbitals (c for HOMO and c' for LUMO) of A2 are quite similar to those of B5. Electrostatic potential surfaces of B5 (d) and A2 (e) resemble each other as well.

24

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. Shown is the energy decomposition analysis of B…N bonding. Data from both covalent and non-covalent bonds is presented for comparison. The contributions to attraction from each component are normalized for comparison. Polarization, electrostatic, and dispersion terms vary monotonically from B6 (covalent bond) to B1~B5, A2 and to B7 (non-covalent bond), indicating that the B…N bonding in A2 (as well as in B1~B5) falls in between a covalent and non-covalent bond.

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Shown is the energy partition protocol to estimate the ring steric effect and the instantaneous energy of the C-C bond. C1 and C2 are analogs to A1 and A2, respectively, with one more CH2 group inserted. In the energy partition protocol, each molecule was separated into two neutral pieces by cutting the C-C bond (as indicated by the red arrowhead). The Eint between the two fragments was then calculated with the energy decomposition analysis. Model compound XCC was constructed to estimate EC-C, where its geometry is kept in the original position as in the primary compound X. With known Eint and EC-C, the effect of ring strain (Ering) could be estimated using equation 7.

26

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 7. Shown are the relative stabilities of the compounds (B1, B5 and A2) at different B…N distances. For A2 (blue curve), the release of ring strain partially offset the energy lost when the B…N distance is increasing, showing the delicate balance between B…N bonding and ring strain in the cyclic form of 2-APB.

Figure 8. The computed stability of different oligomers indicates that aggregation increases the stability of the complex.

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table of Contents Image

28

ACS Paragon Plus Environment

Page 28 of 28