Atomic structure of a short alpha-helix stabilized by a main chain hydrogen-bond surrogate

Herein we report the 1.15 A X-ray crystal structure of a short alpha-helix in which the N-terminal i and i + 4 main chain hydrogen bond is replaced with a carbon-carbon bond. The structure shows that the hydrogen-bond surrogate (HBS) derived alpha-helix truly resembles the structure of canonical alpha-helices and provides unequivocal support for our helix nucleation strategy.     

Nature of the Ga-Ga bonding in Na2[Arx*GaGaArx*]（Arx*=C6H3 -2,6-（C6H5）2）:Electron localization function analysis

The nature of the Ga-Ga bonding in Na2[Arx*GaGaArx*](Arx* = C6H3-2,6-(C6H5)2) has been investigated and compared with that of in H2[Arx*GaGaArx*] using electron localization function(ELF) and orbital analysis.The calculation results show that in Na2[Arx*GaGaArx*],the Ga-Ga interaction is a non-classical triple bond,the heart of Na2[Arx*GaGaArx*] is the Ga2Na2 cluster rather than a simple Ga-Ga bond,and the contribution of the sodium atoms to the short Ga-Ga bond length is considerable.As the two sodium atoms are substituted by two hydrogen atoms,the Ga-Ga bond is replaced by two 3-center,2-electron(3c-2e) Ga-H-Ga covalent bridged bonding.     

FALDI‐based decomposition of an atomic interaction line leads to 3D representation of the multicenter nature of interactions

Atomic interaction lines (AILs) and the QTAIM's molecular graphs provide a predominantly two‐center viewpoint of interatomic interactions. While such a bicentric interpretation is sufficient for most covalent bonds, it fails to adequately describe both formal multicenter bonds as well as many non‐covalent interactions with some multicenter character. We present an extension to our Fragment, Atomic, Localized, Delocalized and Interatomic (FALDI) electron density (ED) decomposition scheme, with which we can measure how any atom‐pair's delocalized density concentrates, depletes or reduces the electron density in the vicinity of a bond critical point. We apply our method on five classical bonds/interactions, ranging from formal either two‐ or three‐center bonds, a non‐covalent interaction (an intramolecular hydrogen bond) to organometallic bonds with partial multicenter character. By use of 3D representation of specific atom‐pairs contributions to the delocalized density we (i) fully recover previous notion of multicenter bonding in diborane and predominant bicentric character of a single covalent C C bond, (ii) reveal a multicenter character of an intramolecular H‐bond and (iii) illustrate, relative to a Schrock carbene, a larger degree of multicenter M C interaction in a Fischer carbene (due to a presence of a heteroatom), whilst revealing the holistic nature of AILs from multicenter ED decomposition. © 2018 Wiley Periodicals, Inc.     

The initiation mechanisms for surface hydrosilylation with 1-alkenes

Hydrosilylation provides a route to form substituted silanes in solution. A similar reaction has been observed in the formation of covalent organic monolayers on a hydrogen-terminated silicon surface and is called thermal hydrosilylation. In solution, the mechanism requires a catalyst to add the basal silicon and saturating hydrogen to the C=C double bond. On the silicon surface, however, the reaction proceeds efficiently at 200 °C, initiated by visible light, and more slowly at room temperature in the dark. Such low activation energy barriers for the reactions on a surface relative to that required for solution hydrosilylation are remarkable, and although many explanations have been suggested, controversy still exists. In this work using a constrained molecular dynamics approach within the density functional theory framework, we show that the free energy activation barrier for abstraction of a hydrogen from silicon by an alkene molecule can be overcome by visible light or thermal excitation. Furthermore, we show that by concerted transfer of a hydrogen from the α-carbon to the β-carbon, a 1-alkene can insert its α-carbon into a surface Si-H bond to accomplish hydrosilylation.     

Calculation of carbon-carbon bond deformation curves for fluorinated ethylenes and the corresponding carbocations

Carbon-carbon bond deformation curves for fluorinated ethylene molecules and the corresponding carbocations were calculated by the ab initio self-consistent field method in the 5-31G basis set. The maximum force required for bond cleavage was taken as a criterion for bond strength. It has been found that for ethylene, replacement of hydrogen by fluorine insignificantly strengthens the C=C bonds in symmetric molecules. However, in molecules with an asymmetric arrangement of fluorine atoms, the bond is slightly weakened due to different charges on the carbon atoms. The configuration of the corresponding carbocations also depends on the bond polarity: an assymmetric distribution of electron density in the C=C bond region leads to the formation of σ-complexes, while a symmetric distribution of electron density (pure covalent bonding) gives π-complexes. Since the carbon-carbon bond in the σ-complexes is essentially weaker, one should expect significant weakening of the bond in high-acidity media if the bond exhibits any kind of asymmetry (chain branching, defects, etc.). For the considered molecules, an antibatic correlation has been established between the strength criterion F_max (unlike the dissociation energy) and the bond length. Institute of Physical Chemistry, Russian Academy of Sciences, Moscow. Translated fromZhurnal Strukturnoi Khimii, Vol. 36, No. 1, pp. 34–41, January–February, 1995. Translated by I. Izvekova     

Interfacially Bridging Covalent Network Yields Hyperstable and Ultralong Virus-Based Fibers for Engineering Functional Materials

We present a strategy of interfacially bridging covalent network within tobacco mosaic virus (TMV) virus-like particles (VLPs). We arranged T103C cysteine to laterally conjugate adjacent subunits. In the axis direction, we set A74C mutation and systematically investigated candidate from E50C to P54C as the other thiol function site, for forming longitudinal disulfide bond chains. Significantly, the T103C-TMV-E50C-A74C shows the highest robustness in assembly capability and structural stability with the largest length, for TMV VLP to date. The fibers with lengths from several to a dozen of micrometers even survive under pH 13. The robust nature of this TMV VLP allows for reducer-free synthesis of excellent electrocatalysts for application in harshly alkaline hydrogen evolution.     

The charge density distribution in a model compound of the catalytic triad in serine proteases

Combined low temperature (28(1) K) X-ray and neutron diffraction measurements were carried out on the co-crystallised complex of betaine, imidazole, and picric acid (1). The experimental charge density was determined and compared with ab initio theoretical calculations at the B3LYP/6-311G(d,p) level of theory. The complex serves as a model for the active site in, for example, the serine protease class of enzymes, the so-called catalytic triad. The crystal contains three short strong N-H...O hydrogen bonds (HBs) with dN...O < 2.7 A. The three HBs have energies above 13 kcalmol(-1), although the hydrogen atoms are firmly localized in the "nitrogen wells". This suggests that low-barrier hydrogen bonding in catalytic enzyme reactions may be a sufficient, but not a necessary, condition for obtaining transition-state stabilization. Structural analysis (e.g., covalent N-H bond lengthening) indicates that the hydrogen bond between H3A and 08 of imidazole and betaine respectively (HB2) is slightly stronger than the bond between H1A and O1A of imidazole and picric acid (HB1), although HB1 is shorter than HB2: (dN...O(HB1)= 2.614(1) A, dN...O(HB2) = 2.684(1) A, dH...O(HB1) = 1.630(1) A, dH...O(HB2)= 1.635(1) A, dN-H(HB1) = 1.046(1) A, dN-H(HB2) = 1.057(1) A). Furthermore, the charge density analysis reveals that HB2 has a larger covalent character than HB1, with considerable polarization of the density towards the acceptor atom. The Gatti and Bader source function (S) is introduced to the analysis of strong HBs. The source function is found to be a sensitive measure for the nature of a hydrogen bond, and comparison with low-barrier and single-well hydrogen bonding systems (e.g., benzoylacetone and nitromalonamide) shows that the low-barrier hydrogen bond (LBHB) state is characterized by an enormously increased hydrogen atom source contribution to the bond critical point in the HB. In this context, HB2 can be characterized as intermediate between localized HBs and delocalized LBHBs.     

取代基对N—H…O＝C氢键三聚体中氢键强度的影响

使用MP2方法研究了氢键三聚体中N-H…O=C氢键强度,探讨了氢键受体分子中不同取代基对N-H…O=C氢键强度的影响.研究表明,不同取代基对氢键三聚体中N-H…O=C氢键强度的影响是不同的:取代基为供电子基团,氢键键长r(H…O)缩短,氢键强度增强;取代基为吸电子基团,氢键键长r(H…O)伸长,氢键强度减弱.自然键轨道(NBO)分析表明,N-H…O=C氢键强度越强,氢键中氢原子的正电荷越多,氧原子的负电荷越多,质子供体和受体分子间的电荷转移越多.供电子基团使N-H…O=C氢键中氧原子的孤对电子n(O)对N-H的反键轨道σ~*(N-H)的二阶相互作用稳定化能增加,吸电子基团使这种二阶相互作用稳定化能减小.取代基对与其相近的N-H…O=C氢键影响更大.     

Identification of hydrogen bonds using quantum electrodynamics

A method for the identification of hydrogen bonds was investigated from the viewpoint of the stress tensor density proposed by Tachibana and following other works in this field. Hydrogen bonds are known to exhibit common features with ionic and covalent bonds. In quantum electrodynamics, the covalent bond has been demonstrated to display a spindle structure of the stress tensor density. Importantly, this spindle structure is also seen in the hydrogen bond, although the covalency is considerably weaker than in a typical covalent bond. Distinguishing it from the ionic bond is most imperative for the identification of the hydrogen bond. In the present study, the directionality of the hydrogen bond is investigated as the ionic bond is nearly isotropic, while the hydrogen bond exhibits the directionality. It was demonstrated that the hydrogen bond can be distinguished from the ionic bond using the angle dependence of the largest eigenvalue of the stress tensor density.     

Interfacially Bridging Covalent Network Yields Hyperstable and Ultralong Virus‐Based Fibers for Engineering Functional Materials

We present a strategy of interfacially bridging covalent network within tobacco mosaic virus (TMV) virus‐like particles (VLPs). We arranged T103C cysteine to laterally conjugate adjacent subunits. In the axis direction, we set A74C mutation and systematically investigated candidate from E50C to P54C as the other thiol function site, for forming longitudinal disulfide bond chains. Significantly, the T103C‐TMV‐E50C‐A74C shows the highest robustness in assembly capability and structural stability with the largest length, for TMV VLP to date. The fibers with lengths from several to a dozen of micrometers even survive under pH 13. The robust nature of this TMV VLP allows for reducer‐free synthesis of excellent electrocatalysts for application in harshly alkaline hydrogen evolution.     

Electrostatic potentials and covalent radii

We begin with a brief overview of the electrostatic potential V(r) as a fundamental determinant of the properties of systems of electrons and nuclei. The minimum of V(r) along the internuclear axis between two bonded atoms is a natural and physically meaningful boundary point, at which the electrostatic forces of the two nuclei upon an element of charge exactly cancel. We propose that the distances from nuclei to V(r) bond minima provide the basis for a well-defined set of covalent radii. Density functional calculations at the B3PW91/6-311+G** level were carried out for 59 molecules to locate the V(r) minima in 95 bonds and use these as the basis for determining single- and multiple-bond covalent radii for eight first- and second-row atoms plus hydrogen. It was found to be unrealistic to assign a single covalent radius to each atom; different values are needed for bonds to first- and second-row atoms, as well as to hydrogen. Using these results, we are able to predict the bond lengths of 33 single and multiple bonds with average errors of less than 0.04 A relative to experimental data.     

取代基对N-H…O=C氢键三聚体中氢键强度的影响

使用MP2方法研究了氢键三聚体中N—H…O=C氢键强度, 探讨了氢键受体分子中不同取代基对N—H…O=C氢键强度的影响. 研究表明, 不同取代基对氢键三聚体中N—H…O=C氢键强度的影响是不同的: 取代基为供电子基团, 氢键键长r(H…O)缩短, 氢键强度增强; 取代基为吸电子基团, 氢键键长r(H…O)伸长, 氢键强度减弱. 自然键轨道(NBO)分析表明, N—H…O=C氢键强度越强, 氢键中氢原子的正电荷越多, 氧原子的负电荷越多, 质子供体和受体分子间的电荷转移越多. 供电子基团使N—H…O=C氢键中氧原子的孤对电子n(O)对N—H的反键轨道滓*(N—H)的二阶相互作用稳定化能增加, 吸电子基团使这种二阶相互作用稳定化能减小. 取代基对与其相近的N—H…O=C氢键影响更大.     

Molecular dynamics extended for fluctuating networks: Application to water

Molecular simulation models are increasingly important tools in efforts to understand the role that water plays in biochemical processes. However, existing models of water have limited capacity to deal with the characteristics of hydrogen bond networks. This article proposes a new fluctuating network (FN) algorithm as an extension of the standard molecular dynamics algorithm. The new algorithm allows for the simulation of a molecular system based on an underlying network, such as the hydrogen bond network in water. This algorithm distinguishes strong from weak network connections, applying a potential that best describes the specific connection behavior. We model liquid water with this new technique using a single‐site, isotropic, short‐range potential. We successfully reproduce liquid water's signature molecular spacing (as represented by the radial distribution function) and characterize its dynamic properties including the exponential hydrogen bond lifetime distribution, diffusion rate, and average hydrogen bonds per molecule. The FN algorithm allows exploration of the behavior of networked systems where explicit coordination limits are required. As such it could also be used to model covalent interactions, reaction dynamics, and applied to simulation of cellular networks. © 2012 Wiley Periodicals, Inc.     

Contribution of multicentered short hydrogen bond arrays to potency of active site-directed serine protease inhibitors

We describe and compare the pH dependencies of the potencies and of the bound structures of two inhibitor isosteres that form multicentered short hydrogen bond arrays at the active sites of trypsin, thrombin, and urokinase type plasminogen activator (urokinase or uPA) over certain ranges of pH. Depending on the pH, short hydrogen bond arrays at the active site are mediated by two waters, one in the oxyanion hole (H(2)O(oxy)) and one on the other (S2) side of the inhibitor (H(2)O(S2)), by one water (H(2)O(oxy)), or by no water. The dramatic variation in the length of the active site hydrogen bonds as a function of pH, of inhibitor, and of enzyme, along with the involvement or absence of ordered water, produces a large structural manifold of active site hydrogen bond motifs. Diverse examples of multicentered and two-centered short hydrogen bond arrays, both at and away from the active site, recently discovered in several protein crystal systems, suggest that short hydrogen bonds in proteins may be more common than has been recognized. The short hydrogen bond arrays resemble one another with respect to ionic nature, highly polar environment, multitude of associated ordinary hydrogen bonds, and disparate pK(a) values of participating groups. Comparison of structures and K(i) values of trypsin complexes at pH values where the multicentered short hydrogen bond arrays mediating inhibitor binding are present or absent indicate that these arrays have a minor effect on inhibitor potency. These features suggest little covalent nature within the short hydrogen bonds, despite their extraordinary shortness (as short as 2.0 A).     

Barriers of hydrogen abstraction vs halogen exchange: an experimental manifestation of charge-shift bonding

This paper shows that the differences between the barriers of the halogen exchange reactions, in the H + XH systems, and the hydrogen abstraction reactions, in the X + HX systems (X = F, Cl, Br), measure the covalent-ionic resonance energies of the corresponding X-H bonds. These processes are investigated using CCSD(T) calculations as well as the breathing-orbital valence bond (BOVB) method. Thus, the VB analysis shows that (i) at the level of covalent structures the barriers are the same for the two series and (ii) the higher barriers for halogen exchange processes originate solely from the less efficient mixing of the ionic structures into the respective covalent structures. The barrier differences, in the HXH vs XHX series, which decrease as X is varied from F to I, can be estimated as one-quarter of the covalent-ionic resonance energy of the H-X bond. The largest difference (22 kcal/mol) is calculated for X = F in accord with the finding that the H-F bond possesses the largest covalent-ionic resonance energy, 87 kcal/mol, which constitutes the major part of the bonding energy. The H-F bond belongs to the class of "charge-shift" bonds (Shaik, S.; Danovich, D.; Silvi, B.; Lauvergnat, D. L.; Hiberty, P. C. Chem. Eur. J. 2005, 21, 6358), which are all typified by dominant covalent-ionic resonance energies. Since the barrier difference between the two series is an experimental measure of the resonance energy quantity, in the particular case of X = F, the unusually high barrier for the fluorine exchange reaction emerges as an experimental manifestation of charge-shift bonding.     

Topology-Templated Synthesis of Crystalline Porous Covalent Organic Frameworks

A strategy is presented for the synthesis of crystalline porous covalent organic frameworks via topology-templated polymerization. The template is based on imine-linked frameworks and their (001) facets seed the C=C bond formation reaction to constitute 2D sp² carbon-conjugated frameworks. This strategy is applicable to templates with different topologies, enables designed synthesis of frameworks that cannot be prepared via direct polymerization, and creates a series of sp² carbon frameworks with tetragonal, hexagonal, and kagome topologies. The sp² carbon frameworks are highly luminescent even in the solid state and exhibit topology-dependent π transmission and exciton migration; these key fundamental π functions are unique to sp² carbon-conjugated frameworks and cannot be accessible by imine-linked frameworks, amorphous analogues, and 1D conjugated polymers. These results demonstrate an unprecedented strategy for structural and functional designs of covalent organic frameworks.     

Structural characterization of alpha-helices of implicitly solvated poly-alanine

The structural characteristics of alpha-helices in poly-alanine-based peptides have been investigated via molecular dynamics simulation with the goal of understanding the basic features of peptide simulations within the context of a model system, classical molecular dynamics with generalized Born (GB) solvation, and to shed insight into the formation and stabilization of alpha-helices in short peptides. The effects of peptide length, terminal charges, proline substitution, and temperature on the alpha-helical secondary structure have been studied. The simulations have shown that distinct secondary structure begins to develop in peptides with lengths approaching 10 residues while ambiguous structures occur in shorter peptides. The helical content of peptides with lengths > or =10 amino acids is observed to be nearly constant up to (Ala)(40). Interestingly, terminal charges and proline in the second position from the N-terminus alter the secondary structure locally with little effect on the overall alpha-helical content of the peptide. The free energy profile of helix formation was also investigated. A large increase in free energy accompanying the formation of helices with more than two consecutive hydrogen bonds in the (i, i + 4) pattern was observed while the free energy increases linearly with additional hydrogen bonds. Values for the change in enthalpy and entropy of helix nucleation and propagation are reported. Additionally the results obtained from the GB model are compared to explicit solvent simulations of two synthetic alanine-based peptides.     

Very strong C-H...O, N-H...O, and O-H...O hydrogen bonds involving a cyclic phosphate.

Very short C-H...O, N-H...O, and O-H...O hydrogen bonds have been generated utilizing the cyclic phosphate [CH2(6-t-Bu-4-Me-C6H2O)2]P(O)OH (1). X-ray structures of (i) 1 (unsolvated, two polymorphs), 1...EtOH, and 1...MeOH, (ii) [imidazolium](+)[CH2(6-t-Bu-4-Me-C6H2O)2PO2](-)...MeOH [2], (iii) [HNC5H4-N=N-C5H4NH](2+)[(CH2(6-t-Bu-4-Me-C6H2O)2PO2)2](2-)...4CH3CN...H2O [3], (v) [K, 18-crown-6](+)[(CH2(6-t-Bu-4-Me-C6H2O)2P(O)OH)(CH2(6-t-Bu-4-Me-C6H2O)2PO2)](-)...2THF [4], (vi) 1...cytosine...MeOH [5], (vii) 1...adenine...1/2MeOH [6], and (viii) 1...S-(-)-proline [7] have been determined. The phosphate 1 in both its forms is a hydrogen-bonded dimer with a short O-H...O distance of 2.481(2) [triclinic form] or 2.507(3) A [monoclinic form]. Compound 2 has a helical structure with a very short C-H...O hydrogen bond involving an imidazolyl C-H and methanol in addition to N-H...O hydrogen bonds. A helical motif is also seen in 5. In 3, an extremely short N-H...O hydrogen bond [N...O 2.558(4) A] is observed. Compounds 6 and 7 also exhibit short N-H...O hydrogen bonds. In 1...EtOH, a 12-membered hydrogen-bonded ring motif, with one of the shortest known O-H...O hydrogen bonds [O...O 2.368(4) A], is present. 1...MeOH is a similar dimer with a very short O(-H)...O bond [2.429(3) A]. In 4, the deprotonated phosphate (anion) and the parent acid are held together by a hydrogen bond on one side and a coordinate/covalent bond to potassium on the other; the O-H...O bond is symmetrical and very strong [O...O 2.397(3) A].     

Carbene triel bonds between TrR3 (Tr = B,Al) and N-heterocyclic carbenes

The carbene triel bond is predicted and characterized by theoretical calculations. The C lone pair of N-heterocyclic carbenes (NHCs) is allowed to interact with the central triel atom of TrR₃ (Tr = B and Al; R = H, F, Cl, and Br). The ensuing bond is very strong, with an interaction energy of nearly 90 kcal/mol. Replacement of the C lone pair by that of either N or Si weakens the binding. The bond is strengthened by electron-withdrawing substituents on the triel atom, and the reverse occurs with substitution on the NHC. However, these effects do not strictly follow the typical pattern of F > Cl > Br. The TrR₃ molecule suffers a good deal of geometric deformation, requiring on the order of 30 kcal/mol, in forming the complex. The R(C···Tr) bond is quite short, for example, 1.6 Å for Tr = B, and shows other indications of at least a partially covalent bond, such as a high electron density at the bond critical point and a good deal of intermolecular charge transfer.