Diastereoselective [3 + 2] Cycloaddition between Tertiary Amine N-Oxides and Substituted Alkenes to Access 7-Azanorbornanes.

We have developed a diastereoselective synthesis of 43 novel 7-azanorbornanes using tertiary amine N-oxides and substituted alkenes. Our method uses an efficient [3 + 2] cycloaddition, starting from either commercially available or easily accessible precursors to generate yields up to 97% and diastereomeric ratios up to >20:1. Density functional theory (DFT) calculations were performed, suggesting that the observed diastereoselectivity is likely due to steric considerations.

Delta SARS-CoV-2 s2m Structure, Dynamics, and Entropy: Consequences of the G15U Mutation.

Bioinformatic analysis of the Delta SARS-CoV-2 genome reveals a single nucleotide mutation (G15U) in the stem-loop II motif (s2m) relative to ancestral SARS-CoV-2. Despite sequence similarity, unexpected differences between SARS-CoV-2 and Delta SARS-CoV-2 s2m homodimerization experiments require the discovery of unknown structural and thermodynamic changes necessary to rationalize the data. Using our reported SARS-CoV-2 s2m model, we induced the G15U substitution and performed 3.5 microseconds of unbiased molecular dynamics simulation at 283 and 310 K. The resultant Delta s2m adopted a secondary structure consistent with our reported NMR data, resulting in significant deviations in the tertiary structure and dynamics from our SARS-CoV-2 s2m model. First, we find differences in the overall three-dimensional structure, where the characteristic 90° L-shaped kink of the SARS-CoV-2 s2m did not form in the Delta s2m resulting in a "linear" hairpin with limited bending dynamics. Delta s2m helical parameters are calculated to align closely with A-form RNA, effectively eliminating a hinge point to form the L-shape kink by correcting an upper stem defect in SARS-CoV-2 induced by a noncanonical and dynamic G:A base pair. Ultimately, the shape difference rationalizes the migration differences in reported electrophoresis experiments. Second, increased fluctuation of the Delta s2m palindromic sequence, within the terminal loop, compared to SARS-CoV-2 s2m results in an estimated increase of entropy of 6.8 kcal/mol at 310 K relative to the SARS-CoV-2 s2m. The entropic difference offers a unique perspective on why the Delta s2m homodimerizes less spontaneously, forming fewer kissing dimers and extended duplexes compared to SARS-CoV-2. In this work, both the L-shape reduction and palindromic entropic penalty provides an explanation of our reported in vitro electrophoresis homodimerization results. Ultimately, the structural, dynamical, and entropic differences between the SARS-CoV-2 s2m and Delta s2m serve to establish a foundation for future studies of the s2m function in the viral lifecycle.

Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p.

The s2m, a highly conserved 41-nt hairpin structure in the SARS-CoV-2 genome, serves as an attractive therapeutic target that may have important roles in the virus life cycle or interactions with the host. However, the conserved s2m in Delta SARS-CoV-2, a previously dominant variant characterized by high infectivity and disease severity, has received relatively less attention than that of the original SARS-CoV-2 virus. The focus of this work is to identify and define the s2m changes between Delta and SARS-CoV-2 and the subsequent impact of those changes upon the s2m dimerization and interactions with the host microRNA miR-1307-3p. Bioinformatics analysis of the GISAID database targeting the s2m element reveals a >99% correlation of a single nucleotide mutation at the 15th position (G15U) in Delta SARS-CoV-2. Based on 1H NMR spectroscopy assignments comparing the imino proton resonance region of s2m and the s2m G15U at 19°C, we show that the U15-A29 base pair closes, resulting in a stabilization of the upper stem without overall secondary structure deviation. Increased stability of the upper stem did not affect the chaperone activity of the viral N protein, as it was still able to convert the kissing dimers formed by s2m G15U into a stable duplex conformation, consistent with the s2m reference. However, we show that the s2m G15U mutation drastically impacts the binding of host miR-1307-3p. These findings demonstrate that the observed G15U mutation alters the secondary structure of s2m with subsequent impact on viral binding of host miR-1307-3p, with potential consequences on immune responses.

Uridylation of the histone mRNA stem-loop weakens binding interactions with SLBP while maintaining interactions with 3'hExo.

Histone mRNA degradation is controlled by the unique 3' stem-loop of histone mRNA and the stem-loop binding protein (SLBP). As part of this process, the 3' stem-loop is trimmed by the histone-specific 3' exonuclease (3'hExo) and uridylated by the terminal uridylyl transferase 7 (TUT7), creating partially degraded intermediates with short uridylations. The role of these uridylations in degradation is not fully understood. Our work examines changes in the stability of the ternary complex created by trimming and uridylation of the stem-loop to better understand the role of this process in the histone mRNA life cycle. In this study, we used fluorescence polarization and electrophoretic mobility shift assays to demonstrate that both SLBP and 3'hExo can bind to uridylated and partially degraded stem-loop intermediates, although with lower affinity. We further characterized this complex by performing 1-µs molecular dynamics simulations using the AMBER force field and Nanoscale Molecular Dynamics (NAMD). These simulations show that while uridylation helps maintain the overall shape of the stem-loop, the combination of uridylation and dephosphorylation of the TPNK motif in SLBP disrupts key RNA-protein interactions. They also demonstrate that uridylation allows 3'hExo to maintain contact with the stem-loop after partial degradation and plays a role in disrupting key base pairs in partially degraded histone mRNA intermediates. Together, these experiments and simulations suggest that trimming by 3'hExo, uridylation, and SLBP dephosphorylation weakens both RNA-protein interactions and the stem-loop itself. Our results further elucidate the role of uridylation and SLBP dephosphorylation in the early stages of histone mRNA degradation.

Parameterization of the miniPEG-Modified γPNA Backbone: Toward Induced γPNA Duplex Dissociation.

γ-Modified peptide nucleic acids (γPNAs) serve as potential therapeutic agents against genetic diseases. Miniature poly(ethylene glycol) (miniPEG) has been reported to increase solubility and binding affinity toward genetic targets, yet details of γPNA structure and dynamics are not understood. Within our work, we parameterized missing torsional and electrostatic terms for the miniPEG substituent on the γ-carbon atom of the γPNA backbone in the CHARMM force field. Microsecond timescale molecular dynamics simulations were carried out on six miniPEG-modified γPNA duplexes from NMR structures (PDB ID: 2KVJ). Three NMR models for the γPNA duplex (PDB ID: 2KVJ) were simulated as a reference for structural and dynamic changes captured for the miniPEG-modified γPNA duplex. Principal component analysis performed on the γPNA backbone atoms identified a single isotropic conformational substate (CS) for the NMR simulations, whereas four anisotropic CSs were identified for the ensemble of miniPEG-modified γPNA simulations. The NMR structures were found to have a 23° helical bend toward the major groove, consistent with our simulated CS structure of 19.0°. However, a significant difference between simulated methyl- and miniPEG-modified γPNAs involved the opportunistic invasion of miniPEG through the minor and major groves. Specifically, hydrogen bond fractional analysis showed that the invasion was particularly prone to affect the second G-C base pair, reducing the Watson-Crick base pair hydrogen bond by 60% over the six simulations, whereas the A-T base pairs decreased by only 20%. Ultimately, the invasion led to base stack reshuffling, where the well-ordered base stacking was reduced to segmented nucleobase stacking interactions. Our 6 µs timescale simulations indicate that duplex dissociation suggests the onset toward γPNA single strands, consistent with the experimental observation of decreased aggregation. To complement the insight of miniPEG-modified γPNA structure and dynamics, the new miniPEG force field parameters allow for further exploration of such modified γPNA single strands as potential therapeutic agents against genetic diseases.

Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p.

The stem loop 2 motif (s2m), a highly conserved 41-nucleotide hairpin structure in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome, serves as an attractive therapeutic target that may have important roles in the virus life cycle or interactions with the host. However, the conserved s2m in Delta SARS-CoV-2, a previously dominant variant characterized by high infectivity and disease severity, has received relatively less attention than that of the original SARS-CoV-2 virus. The focus of this work is to identify and define the s2m changes between Delta and SARS-CoV-2 and subsequent impact of those changes upon the s2m dimerization and interactions with the host microRNA miR-1307-3p. Bioinformatics analysis of the GISAID database targeting the s2m element reveals a greater than 99% correlation of a single nucleotide mutation at the 15 th position (G15U) in Delta SARS-CoV-2. Based on 1 H NMR assignments comparing the imino proton resonance region of s2m and the G15U at 19°C, we find that the U15-A29 base pair closes resulting in a stabilization of the upper stem without overall secondary structure deviation. Increased stability of the upper stem did not affect the chaperone activity of the viral N protein, as it was still able to convert the kissing dimers formed by s2m G15U into a stable duplex conformation, consistent with the s2m reference. However, we find that the s2m G15U mutation drastically reduces the binding affinity of the host miR-1307-3p. These findings demonstrate that the observed G15U mutation alters the secondary structure of s2m with subsequent impact on viral binding of host miR-1307-3p, with potential consequences on the immune response.

Structural, Dynamical, and Entropic Differences between SARS-CoV and SARS-CoV-2 s2m Elements Using Molecular Dynamics Simulations.

The functional role of the highly conserved stem-loop II motif (s2m) in SARS-CoV and SARS-CoV-2 in the viral lifecycle remains enigmatic and an intense area of research. Structure and dynamics of the s2m are key to establishing a structure-function connection, yet a full set of atomistic resolution coordinates is not available for SARS-CoV-2 s2m. Our work constructs three-dimensional coordinates consistent with NMR solution phase data for SARS-CoV-2 s2m and provides a comparative analysis with its counterpart SARS-CoV s2m. We employed initial coordinates based on PDB ID 1XJR for SARS-CoV s2m and two models for SARS-CoV-2 s2m: one based on 1XJR in which we introduced the mutations present in SARS-CoV-2 s2m and the second based on the available SARS-CoV-2 NMR NOE data supplemented with knowledge-based methods. For each of the three systems, 3.5 µs molecular dynamics simulations were used to sample the structure and dynamics, and principal component analysis (PCA) reduced the ensembles to hierarchal conformational substates for detailed analysis. Dilute solution simulations of SARS-CoV s2m demonstrate that the GNRA-like terminal pentaloop is rigidly defined by base stacking uniquely positioned for possible kissing dimer formation. However, the SARS-CoV-2 s2m simulation did not retain the reported crystallographic SARS-CoV motifs and the terminal loop expands to a highly dynamic "nonaloop." Increased flexibility and structural disorganization are observed for the larger terminal loop, where an entropic penalty is computed to explain the experimentally observed reduction in kissing complex formation. Overall, both SARS-CoV and SARS-CoV-2 s2m elements have a similarly pronounced L-shape due to different motif interactions. Our study establishes the atomistic three-dimensional structure and uncovers dynamic differences that arise from s2m sequence changes, which sets the stage for the interrogation of different mechanistic pathways of suspected biological function.

The Scientific Method as a Scaffold to Enhance Communication Skills in Chemistry.

Scientific success in the field of chemistry depends upon the mastery of a wide range of soft skills, most notably scientific writing and speaking. However, training for scientific communication is typically limited at the undergraduate level, where students struggle to express themselves in a clear and logical manner. The underlying issue is deeper than basic technical skills; rather, it is a problem of students' unawareness of a fundamental and strategic framework for writing and speaking with a purpose. The methodology has been implemented for individual mentorship and in our regional summer research program to deliver a blueprint of thought and reasoning that endows students with the confidence and skills to become more effective communicators. Our didactic process intertwines undergraduate research with the scientific method and is partitioned into six steps, referred to as "phases", to allow for focused and deep thinking on the essential components of the scientific method. The phases are designed to challenge the student in their zone of proximal development so they learn to extract and ultimately comprehend the elements of the scientific method through focused written and oral assignments. Students then compile their newly acquired knowledge to create a compelling and logical story, using their persuasive written and oral presentations to complete a research proposal, final report, and formal 20 min presentation. We find that such an approach delivers the necessary guidance to promote the logical framework that improves writing and speaking skills. Over the past decade, we have witnessed both qualitative and quantitative gains in the students' confidence in their abilities and skills (developed by this process), preparing them for future careers as young scientists.

Highly conserved s2m element of SARS-CoV-2 dimerizes via a kissing complex and interacts with host miRNA-1307-3p.

The ongoing COVID-19 pandemic highlights the necessity for a more fundamental understanding of the coronavirus life cycle. The causative agent of the disease, SARS-CoV-2, is being studied extensively from a structural standpoint in order to gain insight into key molecular mechanisms required for its survival. Contained within the untranslated regions of the SARS-CoV-2 genome are various conserved stem-loop elements that are believed to function in RNA replication, viral protein translation, and discontinuous transcription. While the majority of these regions are variable in sequence, a 41-nucleotide s2m element within the genome 3' untranslated region is highly conserved among coronaviruses and three other viral families. In this study, we demonstrate that the SARS-CoV-2 s2m element dimerizes by forming an intermediate homodimeric kissing complex structure that is subsequently converted to a thermodynamically stable duplex conformation. This process is aided by the viral nucleocapsid protein, potentially indicating a role in mediating genome dimerization. Furthermore, we demonstrate that the s2m element interacts with multiple copies of host cellular microRNA (miRNA) 1307-3p. Taken together, our results highlight the potential significance of the dimer structures formed by the s2m element in key biological processes and implicate the motif as a possible therapeutic drug target for COVID-19 and other coronavirus-related diseases.

Multi-Ion Bridged Pathway of N-Oxides to 1,3-Dipole Dilithium Oxide Complexes.

Roussi's landmark work on the generation of 1,3-dipoles from tertiary amine N-oxides has not reached its full potential since its underlying mechanism is neither well explored nor understood. Two competing mechanisms were previously proposed to explain the transformation involving either an iminium ion or a diradical intermediate. Our investigation has revealed an alternative mechanistic pathway that explains experimental results and provides significant insights to guide the creation of new N-oxide reagents beyond tertiary alkylamines for direct synthetic transformations. Truhlar's M06-2x functional and Møller-Plesset second-order perturbation theory with Dunning's [jul,aug]-cc-pv[D,T]z basis sets and discrete-continuum solvation models were employed to determine activation enthalpies and structures. During these mechanistic explorations, we discovered a unique multi-ion bridged pathway resulting from the rate-determining step, which was energetically more favorable than other alternate mechanisms. This newly proposed mechanism contains no electrophilic intermediates, strengthening the reaction potential by broadening the reagent scope and limiting the possible side reactions. This thoroughly defined general mechanism supports a more direct route for improving the use of N-oxides in generating azomethine ylide-dilithium oxide complexes with expanded functional group tolerance and breadth of chemistry.

Metalated nitriles: SNi' cyclizations with a propargylic electrophile.

A rare 5-exo-dig SNi' cyclization with magnesiated and lithiated nitriles affords a cis-fused hydrindane bearing an exocyclic allene. The cyclization of the dilithiated nitrile pits a stereoelectronic preference for a trans-hydrindane against a cyclization through a less strained transition structure to the corresponding cis-hydrindane. Computational modeling suggests that the dilithiated nitrile cyclizes to a cis-hydrindane because the preferred transition structure positions the lithium cation in a cone of electron density that bridges the nitrile-bearing carbon, an alkoxide, and an electron-rich alkyne functionality.

Intramolecular charge-assisted hydrogen bond strength in pseudochair carboxyphosphate.

Carboxyphosphate, a suspected intermediate in ATP-dependent carboxylases, has not been isolated nor observed directly by experiment. Consequently, little is known concerning its structure, stability, and ionization state. Recently, carboxyphosphate as either a monoanion or dianion has been shown computationally to adopt a novel pseudochair conformation featuring an intramolecular charge-assisted hydrogen bond (CAHB). In this work, additive and subtractive correction schemes to the commonly employed open-closed method are used to estimate the strength of the CAHB. Truhlar's Minnesota M06-2X functional with Dunning's aug-cc-pVTZ basis set has been used for geometry optimization, energy evaluation, and frequency analysis. The CHARMM force field has been used to approximate the Pauli repulsive terms in the closed and open forms of carboxyphosphate. From our additive correction scheme, differential Pauli repulsion contributions between the pseudochair (closed) and open conformations of carboxyphosphate are found to be significant in determining the CAHB strength. The additive correction modifies the CAHB prediction (ΔEclosed-open) of -14 kcal/mol for the monoanion and -12 kcal/mol for the dianion to -22.9 and -18.4 kcal/mol, respectively. Results from the subtractive technique reinforce those from our additive procedure, where the predicted CAHB strength ranges from -17.8 to -25.4 kcal/mol for the monoanion and from -15.7 to -20.9 kcal/mol for the dianion. Ultimately, we find that the CAHB in carboxyphosphate meets the criteria for short-strong hydrogen bonds. However, carboxyphosphate has a unique energy profile that does not result in the symmetric double-well behavior of low-barrier hydrogen bonds. These findings provide deeper insight into the pseudochair conformation of carboxyphosphate, and lead to an improved mechanistic understanding of this intermediate in ATP-dependent carboxylases.

Common hydrogen bond interactions in diverse phosphoryl transfer active sites.

Phosphoryl transfer reactions figure prominently in energy metabolism, signaling, transport and motility. Prior detailed studies of selected systems have highlighted mechanistic features that distinguish different phosphoryl transfer enzymes. Here, a top-down approach is developed for comparing statistically the active site configurations between populations of diverse structures in the Protein Data Bank, and it reveals patterns of hydrogen bonding that transcend enzyme families. Through analysis of large samples of structures, insights are drawn at a level of detail exceeding the experimental precision of an individual structure. In phosphagen kinases, for example, hydrogen bonds with the O3ß of the nucleotide substrate are revealed as analogous to those in unrelated G proteins. In G proteins and other enzymes, interactions with O3ß have been understood in terms of electrostatic favoring of the transition state. Ground state quantum mechanical calculations on model compounds show that the active site interactions highlighted in our database analysis can affect substrate phosphate charge and bond length, in ways that are consistent with prior experimental observations, by modulating hyperconjugative orbital interactions that weaken the scissile bond. Testing experimentally the inference about the importance of O3ß interactions in phosphagen kinases, mutation of arginine kinase Arg280 decreases kcat, as predicted, with little impact upon KM.

Hyperconjugation-mediated solvent effects in phosphoanhydride bonds.

Density functional theory and natural bond orbital analysis are used to explore the impact of solvent on hyperconjugation in methyl triphosphate, a model for "energy rich" phosphoanhydride bonds, such as found in ATP. As expected, dihedral rotation of a hydroxyl group vicinal to the phosphoanhydride bond reveals that the conformational dependence of the anomeric effect involves modulation of the orbital overlap between the donor and acceptor orbitals. However, a conformational independence was observed in the rotation of a solvent hydrogen bond. As one lone pair orbital rotates away from an optimal antiperiplanar orientation, the overall magnitude of the anomeric effect is compensated approximately by the other lone pair as it becomes more antiperiplanar. Furthermore, solvent modulation of the anomeric effect is not restricted to the antiperiplanar lone pair; hydrogen bonds involving gauche lone pairs also affect the anomeric interaction and the strength of the phosphoanhydride bond. Both gauche and anti solvent hydrogen bonds lengthen nonbridging O-P bonds, increasing the distance between donor and acceptor orbitals and decreasing orbital overlap, which leads to a reduction of the anomeric effect. Solvent effects are additive with greater reduction in the anomeric effect upon increasing water coordination. By controlling the coordination environment of substrates in an active site, kinases, phosphatases, and other enzymes important in metabolism and signaling may have the potential to modulate the stability of individual phosphoanhydride bonds through stereoelectronic effects.

Periodic trends and index of boron LEwis acidity.

Lewis acidity is customarily gauged by comparing the relative magnitude of coordinate covalent bonding energies, where the Lewis acid moiety is varied and the Lewis base is kept constant. However, the prediction of Lewis acidity from first principles is sometimes contrary to that suggested by experimental bond energies. Specifically, the order of boron trihalide Lewis acidities predicted from substituent electronegativity arguments is opposite to that inferred by experiment. Contemporary explanations for the divergence between theory, computation, and experiment have led to further consternation. Due to the fundamental importance of understanding the origin of Lewis acidity, we report periodic trends for 21 boron Lewis acids, as well as their coordinate covalent bond strengths with NH(3), utilizing ab initio, density functional theory, and natural bond orbital analysis. Coordinate covalent bond dissociation energy has been determined to be an inadequate index of Lewis acid strength. Instead, acidity is measured in the manner originally intended by Lewis, which is defined by the valence of the acid of interest. Boron Lewis acidity is found to depend upon substituent electronegativity and atomic size, differently than for known Brønsted-Lowry periodic trends. Across the second period, stronger substituent electronegativity correlates (R(2) = 0.94) with increased Lewis acidity. However, across the third period, an equal contribution from substituent electronegativity and atomic radii is correlated (R(2) = 0.98) with Lewis acidity. The data suggest that Lewis acidity depends upon electronegativity solely down group 14, while equal contribution from both substituent electronegativity and atomic size are significant down groups 16 and 17. Originally deduced from Pauling's electronegativities, boron's substituents determine acidity by influencing the population of its valence by withdrawing electron density. However, size effects manifest differently than previously considered, where greater sigma bond (not pi bond) orbital overlap between boron and larger substituents increase the electron density available to boron's valence, thereby decreasing Lewis acidity. The computed electronegativity and size effects of substituents establish unique periodic trends that provide a novel explanation of boron Lewis acidity, consistent with first principle predictions. The findings resolve ambiguities between theory, computation, and experiment and provide a clearer understanding of Lewis acidity.

Strong coordination of tetraphenylborate anion to copper(I) bipyridine and phenanthroline-based complexes and its effect on catalytic activity in the cyclopropanation of styrene.

The synthesis, characterization, and cyclopropanation activity of tetrahedral copper(I) complexes with bipyridine- and phenanthroline-based ligands containing strongly coordinated tetraphenylborate anions are reported. Cu(I)(bpy)(BPh(4)), Cu(I)(phen)(BPh(4)), and Cu(I)(3,4,7,8-Me(4)phen)(BPh(4)) complexes are the first examples in which the BPh(4)(-) counterion chelates a transition metal center in bidentate fashion through eta(2) pi interactions with two of its phenyl rings.

Anomeric effect in "high energy" phosphate bonds. Selective destabilization of the scissile bond and modulation of the exothermicity of hydrolysis.

A natural bonding orbital (NBO) analysis of phosphate bonding and connection to experimental phosphotransfer potential is presented. Density functional calculations with the 6-311++G(d,p) basis set carried out on 10 model phosphoryl compounds verify that the wide variability of experimental standard free energies of hydrolysis (a phosphotransfer potential benchmark) is correlated with the instability of the scissile O-P bond through computed bond lengths. NBO analysis is used to analyze all delocalization interactions contributing to O-P bond weakening. Phosphoryl bond lengths are found to correlate strongest (R = 0.90) with the magnitude of the ground-state n(O) --> sigma*(O-P) anomeric effect. Electron-withdrawing interactions of the substituent upon the sigma(O-P) bonding orbital also correlate strongly with O-P bond lengths (R = 0.88). However, an analysis of sigma*(O-P) and sigma(O-P) populations show that the increase in sigma*(O-P) density is up to 6.5 times greater than the decrease in sigma(O-P) density. Consequently, the anomeric effect is more important than other delocalization interactions in impacting O-P bond lengths. Factors reducing anomeric power by diminishing either lone pair donor ability (solvent) or antibonding acceptor ability (substituent) are shown to result in shorter O-P bond lengths. The trends shown in this work suggest that the generalized anomeric effect provides a simple explanation for relating the sensitivity of the O-P bond to diverse environmental and substituent factors. The anomeric n(O) --> sigma*(O-P) interaction is also shown to correlate strongly with experimentally determined standard free energies of hydrolysis (R = -0.93). A causal mechanism cannot be inferred from correlation. Equally, a P-value of 1.2 x 10(-4) from an F-test indicates that it is unlikely that the ground-state anomeric effect and standard free energies of hydrolysis are coincidentally related. It is found that as the exothermicity of hydrolysis increases, the energy stabilization of the ground-state anomeric effect increases with selective destabilization of the high-energy O-P bond to be broken in hydrolysis. The anomeric effect therefore partially counteracts a larger resonance stabilization of products that makes hydrolysis exothermic and needs to be considered in achieving improved agreement between calculated and empirical energies of hydrolysis. The avenues relating the thermodynamic behavior of phosphates to underlying structural factors via the anomeric effect are discussed.

Hybrid Meta-Generalized Gradient Functional Modeling of Boron-Nitrogen Coordinate Covalent Bonds.

Truhlar's new generation of hybrid meta-generalized gradient functionals has been evaluated in modeling the binding enthalpies of substituted B-N coordinate covalent bonds. The short-range exchange correlation (XC) energy of coordinate covalent bonding coupled with the medium-range XC energy of noncovalent interactions results in a particularly difficult case for density functional theory (DFT). In this study, M06, M06-2X, M05, M05-2X, MPWB1K, and MPW1B95 with the 6-311++G(3df,2p) basis set have been used to evaluate four methylated ammonia trimethylboranes, (CH3)3B-N(CH3)nH3-n (n = 0 to 3), along with H3B-NH3. The predicted binding enthalpies from the new functionals have been compared to experiment as well as previous DFT (B3LYP, MPW1K) and ab initio (HF, MP2, QCISD, and QCISD(T)) results. Previously, only MP2, QCISD, and QCISD(T) were found to model the experimental energetic trend accurately. The mean absolute deviation (MAD) from experimental binding enthalpies for M06-2X and M05-2X is 0.3 and 1.6 kcal/mol, respectively. M06-2X yields a lower MAD than more expensive ab initio methods (MP2 = 1.9 kcal/mol and QCISD = 2.3 kcal/mol) and a comparable MAD to QCISD(T) (MAD = 0.4 kcal/mol). M06-2X is shown to provide a balanced account of the short- and medium-range XC energies necessary to describe the binding enthalpy of coordinate covalent bonds accurately in sterically congested molecular systems.

Covalent and ionic nature of the dative bond and account of accurate ammonia borane binding enthalpies.

The inherent difficulty in modeling the energetic character of the B-N dative bond has been investigated utilizing density functional theory and ab initio methods. The underlying influence of basis set size and functions, thermal corrections, and basis set superposition error (BSSE) on the predicted binding enthalpy of ammonia borane (H3B-NH3) and four methyl-substituted ammonia trimethylboranes ((CH3)3B-N(CH3)nH3-n; n = 0-3) has been evaluated and compared with experiment. HF, B3LYP, MPW1K, MP2, QCISD, and QCISD(T) have been utilized with a wide range of Pople and correlation-consistent basis sets, totaling 336 levels of theory. MPW1K, B3LYP, and HF result in less BSSE and converge to binding enthalpies with fewer basis functions than post-SCF techniques; however, the methods fail to model experimental binding enthalpies and trends accurately, producing mean absolute deviations (MADs) of 5.1, 10.8, and 16.3 kcal/mol, respectively. Despite slow convergence, MP2, QCISD, and QCISD(T) using the 6-311++G(3df,2p) basis set reproduce the experimental binding enthalpy trend and result in lower MADs of 2.2, 2.6, and 0.5 kcal/mol, respectively, when corrected for BSSE and a residual convergence error of ca. 1.3-1.6 kcal/mol. Accuracy of the predicted binding enthalpy is linked to correct determination of the bond's dative character given by charge-transfer frustration, QCTF = -(Delta QN + Delta QB). Frustration gauges the incompleteness of charge transfer between the donor and the acceptor. The binding enthalpy across ammonia borane and methylated complexes is correlated to its dative character (R2 = 0.91), where a more dative bond (less charge-transfer frustration) results in a weaker binding enthalpy. However, a balance of electronic and steric factors must be considered to explain trends in experimentally reported binding enthalpies. Dative bond descriptors, such as bond ionicity and covalency are important in the accurate characterization of the dative bond. The B-N dative bond in ammonia borane is 65% ionic, moderately strong (-27.5 +/- 0.5 kcal/mol), and structurally flexible on the donor side to relieve steric congestion.