Co-reporter:Liping Xu ; Margaret J. Hilton ; Xinhao Zhang ; Per-Ola Norrby ; Yun-Dong Wu ; Matthew S. Sigman ;Olaf Wiest
Journal of the American Chemical Society 2014 Volume 136(Issue 5) pp:1960-1967
Publication Date(Web):January 11, 2014
DOI:10.1021/ja4109616
The enantioselective Pd-catalyzed redox-relay Heck arylation of acyclic alkenyl alcohols allows access to various useful chiral building blocks from simple olefinic substrates. Mechanistically, after the initial migratory insertion, a succession of β-hydride elimination and migratory insertion steps yields a saturated carbonyl product instead of the more general Heck product, an unsaturated alcohol. Here, we investigate the reaction mechanism, including the relay function, yielding the final carbonyl group transformation. M06 calculations predict a ΔΔG⧧ of 1 kcal/mol for the site selectivity and 2.5 kcal/mol for the enantioselectivity, in quantitative agreement with experimental results. The site selectivity is controlled by a remote electronic effect, where the developing polarization of the alkene in the migratory insertion transition state is stabilized by the C–O dipole of the alcohol moiety. The enantioselectivity is controlled by steric repulsion between the oxazoline substituent and the alcohol-bearing alkene substituent. The relay efficiency is due to an unusually smooth potential energy surface without high barriers, where the hydroxyalkyl-palladium species acts as a thermodynamic sink, driving the reaction toward the carbonyl product. Computational predictions of the relative reactivity and selectivity of the double bond isomers are validated experimentally.
Co-reporter:Yun-Fang Yang ; Gui-Juan Cheng ; Peng Liu ; Dasheng Leow ; Tian-Yu Sun ; Ping Chen ; Xinhao Zhang ; Jin-Quan Yu ; Yun-Dong Wu ;K. N. Houk
Journal of the American Chemical Society 2013 Volume 136(Issue 1) pp:344-355
Publication Date(Web):December 8, 2013
DOI:10.1021/ja410485g
Density functional theory investigations have elucidated the mechanism and origins of meta-regioselectivity of Pd(II)-catalyzed C–H olefinations of toluene derivatives that employ a nitrile-containing template. The reaction proceeds through four major steps: C–H activation, alkene insertion, β-hydride elimination, and reductive elimination. The C–H activation step, which proceeds via a concerted metalation–deprotonation (CMD) pathway, is found to be the rate- and regioselectivity-determining step. For the crucial C–H activation, four possible active catalytic species—monomeric Pd(OAc)2, dimeric Pd2(OAc)4, heterodimeric PdAg(OAc)3, and trimeric Pd3(OAc)6—have been investigated. The computations indicated that the C–H activation with the nitrile-containing template occurs via a Pd–Ag heterodimeric transition state. The nitrile directing group coordinates with Ag while the Pd is placed adjacent to the meta-C–H bond in the transition state, leading to the observed high meta-selectivity. The Pd2(OAc)4 dimeric mechanism also leads to the meta-C–H activation product but with higher activation energies than the Pd–Ag heterodimeric mechanism. The Pd monomeric and trimeric mechanisms require much higher activation free energies and are predicted to give ortho products. Structural and distortion energy analysis of the transition states revealed significant effects of distortions of the template on mechanism and regioselectivity, which provided hints for further developments of new templates.
Co-reporter:Fan Jiang, Wei Han and Yun-Dong Wu
Physical Chemistry Chemical Physics 2013 vol. 15(Issue 10) pp:3413-3428
Publication Date(Web):12 Dec 2012
DOI:10.1039/C2CP43633G
The local conformational (ϕ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the PII conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The ϕ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ1 rotamer preference and χ1-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the ϕ, ψ plot of alanine with the influence of side-chain χ1 rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental 3J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix–coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.
Co-reporter:Ping Chen, Brenda A. Dougan, Xinhao Zhang, Yun-Dong Wu, Zi-Ling Xue
Polyhedron 2013 58() pp: 30-38
Publication Date(Web):
DOI:10.1016/j.poly.2012.07.042
Co-reporter:Liangchun Li;Hongmei Zhan;Pengfei Duan;Jian Liao;Junming Quan;Yu Hu;Zhongzhu Chen;Jin Zhu;Minghua Liu;Jingen Deng
Advanced Functional Materials 2012 Volume 22( Issue 14) pp:3051-3056
Publication Date(Web):
DOI:10.1002/adfm.201200488
Abstract
Self-assembling cyclic peptide nanotubes (SPNs) have been extensively studied due to their potential applications in biology and material sciences. Cyclic γ-peptides, which have a larger conformational space, have received less attention than the cyclic α- and β-peptides. The self-assembly of cyclic homo-γ-tetrapeptide based on cis-3-aminocyclohexanecarboxylic acid (γ-Ach) residues, which can be easily synthesized by a one-pot process is investigated. Fourier transform infrared (FTIR) and NMR analysis along with density functional theory (DFT) calculations indicate that the cyclic homo-γ-tetrapeptide, with a non-planar conformation, can self-assemble into nanotubes through hydrogen-bond-mediated parallel stacking. Atomic force microscopy (AFM) and transmission electron microscopy (TEM) experiments reveal the formation of bundles of nanotubes in CH2Cl2/hexane, but individual nanotubes and bundles of only two nanotubes are obtained in water. The integration of TEG (triethylene glycol) monomethyl ether chains and cyclopeptide backbones may allow the control of width of single nanotubes.
Co-reporter:Kun Zhao, Gui-Juan Cheng, Hongzhi Yang, Hai Shang, Xinhao Zhang, Yun-Dong Wu, and Yefeng Tang
Organic Letters 2012 Volume 14(Issue 18) pp:4878-4881
Publication Date(Web):September 4, 2012
DOI:10.1021/ol302205w
A concise biomimetic total synthesis of incarvilleatone and incarviditone is achieved in one pot via the highly stereoselective hetero- and homodimerization of (±)-rengyolone, respectively. The structure of incarviditone is revised on the basis of spectroscopic and computational evidence.
Co-reporter:Cheuk-Kin Wan, Wei Han, and Yun-Dong Wu
Journal of Chemical Theory and Computation 2012 Volume 8(Issue 1) pp:300-313
Publication Date(Web):November 8, 2011
DOI:10.1021/ct2004275
The recently developed PACE force field was further parametrized so that it can be applied to the studies of membrane systems. Parameters for the interactions between united-atom protein particles and lipid hydrophobic tails were developed by reproducing the solvation free energies of small organic molecules in hexadecane. Interactions between protein particles and lipid heads were parametrized by fitting the potential of mean force of the corresponding all-atom simulation. The force field was applied to the study of five helical peptides in membrane environments. The calculated tilt angles of WALP and GWALP and their mutations are in good agreement with experimental data. The association of two glycophorin A (GpA) helices was simulated for 6 μs. Root-mean-square-deviation of the simulated dimer from the nuclear magnetic resonance structure was found to be 0.272 nm, better than all results obtained so far. These findings demonstrate the high accuracy and applicability of the PACE force field in studying membrane proteins.
Co-reporter:Yun-Fang Yang;Gui-Juan Cheng;Dr. Jun Zhu;Dr. Xinhao Zhang;Dr. Shigeyoshi Inoue;Dr. Yun-Dong Wu
Chemistry - A European Journal 2012 Volume 18( Issue 24) pp:7516-7524
Publication Date(Web):
DOI:10.1002/chem.201103443
Abstract
Density functional theory calculations (B3LYP) have been carried out to investigate the 4π-electron systems of 2,4-disila-1,3-diphosphacyclobutadiene (compound 1) and the tetrasilacyclobutadiene dication (compound 2). The calculated nucleus-independent chemical shift (NICS) values for these two compounds are negative, which indicates that the core rings of compounds 1 and 2 have a certain amount of aromaticity. However, deep electronic analysis reveals that neither of these two formal 4π-electron four-membered ring systems is aromatic. Compound 1 has very weak, almost negligible antiaromaticity, and the amidinate ligands attached to the Si atoms play an important role in stabilizing this conjugated 4π-electron system. The monoanionic bidentate ligand interacts with the conjugated π system to cause π-orbital splitting. This ligand-induced π-orbital splitting effect provides an opportunity to manipulate the gap between occupied and unoccupied π orbitals in conjugated systems. Conversely, compound 2 is nonaromatic because its core ring does not have a conjugated π ring system and does not fulfill the requirements of a Hückel system.
Co-reporter:Qian Peng, Hong Yan, Xinhao Zhang, and Yun-Dong Wu
The Journal of Organic Chemistry 2012 Volume 77(Issue 17) pp:7487-7496
Publication Date(Web):August 9, 2012
DOI:10.1021/jo301319j
Density functional theory studies have been carried out to investigate the mechanism of the Pd(II)(bpy)- and Rh(I)(bpy)-catalyzed conjugate additions and their competitive Heck reactions involving α,β-unsaturated carbonyl compounds. The critical steps of the mechanism are insertion and termination. The insertion step favors 1,2-addition of the vinyl-coordinated species to generate a stable C-bound enolate intermediate, which then may isomerize to either an oxa-π-allyl species or an O-bound enolate. The termination step involves a competition between β-hydride elimination, leading to a Heck reaction product, and protonolysis reaction that gives a conjugate addition product. These two pathways are competitive in the Pd(II)-catalyzed reaction, while a preference for protonolysis has been found in the Rh(I)-catalyzed reaction. The calculations are in good agreement with the experimental observations. The potential energy surface and the rate-determining step of the β-hydride elimination are similar for both Pd(II)- and Rh(I)-catalyzed processes. The rate-determining steps of the Pd(II)- and Rh(I)-catalyzed protonolysis are different. Introduction of an N- or P-ligand significantly stabilizes the protonolysis transition state via the O-bound enolate or oxa-π-allyl complex intermediate, resulting in a reduced free energy of activation. However, the barrier of the β-hydride elimination is less sensitive to ligands. For the Rh(I)-catalyzed reaction, protonolysis is calculated to be more favorable than the β-hydride elimination for all investigated N and P ligands due to the significant ligand stabilization to the protonolysis transition state. For the Pd(II)-catalyzed reaction, the complex with monodentate pyridine ligands prefers the Heck-type product through β-hydride elimination, while the complex with bidentate N and P ligands favors the protonolysis. The theoretical finding suggests the possibility to control the selectivity between the conjugate addition and the Heck reaction by using proper ligands.
Co-reporter:Jian-Ping Chen ; Qian Peng ; Bai-Lin Lei ; Xue-Long Hou
Journal of the American Chemical Society 2011 Volume 133(Issue 36) pp:14180-14183
Publication Date(Web):August 9, 2011
DOI:10.1021/ja2039503
α-Carbanions of cyclic and acyclic imines have been successfully applied as nucleophiles in the Pd-catalyzed allylic alkylation reaction. Tuning of chemo- and regioselectivity has been realized by using t-BuOK/THF and LDA/toluene to give branched and linear products, respectively, with high regio- and diastereoselectivities. A plausible mechanism is proposed on the basis of the experimental results and DFT calculations.
Co-reporter:Bo Chen ; Xue-Long Hou ; Yu-Xue Li
Journal of the American Chemical Society 2011 Volume 133(Issue 20) pp:7668-7671
Publication Date(Web):April 27, 2011
DOI:10.1021/ja201425e
DFT calculations suggest that the unexpected meta product in the copper-catalyzed arylation of anilide is formed via a Heck-like four-membered-ring transition state involving a CuIII–Ph species. A competitive electrophilic substitution mechanism delivers the ortho product when a methoxy group is present at the meta position of pivanilide. A series of experiments including kinetic studies support the involvement of a CuI catalyst.
Co-reporter:Yun-Fang Yang, Ting Shi, Xin-Hao Zhang, Zong-Xun Tang, Zhen-Yi Wen, Jun-Min Quan and Yun-Dong Wu
Organic & Biomolecular Chemistry 2011 vol. 9(Issue 16) pp:5845-5855
Publication Date(Web):23 May 2011
DOI:10.1039/C1OB05501A
Density functional theory calculations (B3LYP) have been carried out to understand the mechanism and stereochemistry of an asymmetric reductive aldol reaction of benzaldehyde and tert-butyl acrylate with hydrosilanes catalyzed by Rh(Phebox-ip)(OAc)2(OH2). According to the calculations, the reaction proceeds via five steps: (1) oxidative addition of hydrosilane, (2) hydride migration to carbon–carbon double bond of tert-butyl acrylate, which determines the chirality at C2, (3) tautomerization from rhodium bound C-enolate to rhodium bound O-enolate, (4) intramolecular aldol reaction, which determines the chirality at C3 and consequently the anti/syn-selectivity, and (5) reductive elimination to release aldol product. The hydride migration is the rate-determining step with a calculated activation energy of 23.3 kcal mol−1. In good agreement with experimental results, the formation of anti-aldolates is found to be the most favorable pathway. The observed Si-facial selectivity in both hydride migration and aldol reaction are well-rationalized by analyzing crucial transition structures. The Re-facial attack transition state is disfavored because of steric hindrance between the isopropyl group of the catalyst and the tert-butyl acrylate.
Co-reporter:Chun-Shan Zuo;Olaf Wiest
Journal of Physical Organic Chemistry 2011 Volume 24( Issue 12) pp:1157-1165
Publication Date(Web):
DOI:10.1002/poc.1840
Abstract
The structures and conformational energies of twelve MeN-, O-, S-, MeP-, CO-bridged homocalix[4]arenes and two kinds of O-bridged alternate hybrid-calix[4]arenes have been calculated at the B3LYP/6-31G* level of theory. The 1,3-alternate or twisted-1,3-alternate conformations are found to be the lowest energy structures in all cases except for MeP-bridged calix[4]pyridine and calix[4]benzene, for which the twist-pinched cone and partial cone are the most stable conformations, respectively. The conformational energy differences calculated between the lowest energy and the next conformation are on the order of 2.0–3.0 kcal/mol, but smaller for the S- and MeP-bridged compounds. Analysis of the structures and relative energies shows that the phenyl hydrogens have electrostatic attractions with the lone pairs of the heteroatoms and steric repulsions with the methyl groups on bridgehead groups. Conversely, the lone electron pair in the pyridyl compounds engage in a repulsive electrostatic interaction. The interactions of 1,3-aromatic rings play a minor but important role in the relative stability sequence. This detailed understanding of the factors governing the conformational space of hetero calixarences can be used to design conformationally biased analogs of these interesting compounds. Copyright © 2011 John Wiley & Sons, Ltd.
Co-reporter:Wei Han, Cheuk-Kin Wan, and Yun-Dong Wu
Journal of Chemical Theory and Computation 2010 Volume 6(Issue 11) pp:3390-3402
Publication Date(Web):October 25, 2010
DOI:10.1021/ct100313a
We present the application of our recently developed PACE force field to the folding of peptides. These peptides include α-helical (AK17 and Fs), β-sheet (GB1m2 and Trpzip2), and mixed helical/coil (Trp-cage) peptides. With replica exchange molecular dynamics (REMD), our force field can fold the five peptides into their native structures while maintaining their stabilities reasonably well. Our force field is also able to capture important thermodynamic features of the five peptides that have been observed in previous experimental and computational studies, such as different preferences for a helix−turn−helix topology for AK17 and Fs, the relative contribution of four hydrophobic side chains of GB1p to the stability of β-hairpin, and the distinct role of a hydrogen bond involving Trp-Hε and a D9/R16 salt bridge in stabilizing the Trp-cage native structure. Furthermore, multiple folding and unfolding events are observed in our microsecond-long normal MD simulations of AK17, Trpzip2, and Trp-cage. These simulations provide mechanistic information such as a “zip-out” pathway of the folding mechanism of Trpzip2 and the folding times of AK17 and Trp-cage, which are estimated to be about 51 ± 43 ns and 270 ± 110 ns, respectively. A 600 ns simulation of the peptides can be completed within one day. These features of our force field are potentially applicable to the study of thermodynamics and kinetics of real protein systems.
Co-reporter:Wei Han, Cheuk-Kin Wan, Fan Jiang, and Yun-Dong Wu
Journal of Chemical Theory and Computation 2010 Volume 6(Issue 11) pp:3373-3389
Publication Date(Web):October 25, 2010
DOI:10.1021/ct1003127
A further parametrization of a united-atom protein model coupled with coarse-grained water has been carried out to cover all amino acids (AAs). The local conformational features of each AA have been fitted on the basis of restricted coil-library statistics of high-resolution X-ray crystal structures of proteins. Potential functions were developed on the basis of combined backbone and side chain rotamer conformational preferences, or rotamer Ramachandran plots (ϕ, Ψ, χ1). Side chain−side chain and side chain−backbone interaction potentials were parametrized to fit the potential mean forces of corresponding all-atom simulations. The force field has been applied in molecular dynamics simulations of several proteins of 56−108 AA residues whose X-ray crystal and/or NMR structures are available. Starting from the crystal structures, each protein was simulated for about 100 ns. The Cα RMSDs of the calculated structures are 2.4−4.2 Å with respect to the crystal and/or NMR structures, which are still larger than but close to those of all-atom simulations (1.1−3.6 Å). Starting from the PDB structure of malate synthase G of 723 AA residues, the wall-clock time of a 30 ns simulation is about three days on a 2.65 GHz dual-core CPU. The RMSD to the experimental structure is about 4.3 Å. These results implicate the applicability of the force field in the study of protein structures.
Co-reporter:Qiang Li;Kittiya Wongkhan;XianCai Luo;Andrei S. Batsanov
Science Bulletin 2010 Volume 55( Issue 25) pp:2794-2798
Publication Date(Web):2010 September
DOI:10.1007/s11434-010-3270-9
While stable in CH2Cl2, hexane or THF, in the presence of MeOH, self-promoted dimerization of the triarylphosphine-alkene 1, a ligand for Pd-catalyzed reactions, produced an unusual racemic bis(phosphine) 2 in high yield. The reaction of 2 with Pd(dba)2, followed by oxidative addition of p-IC6H4NO2, yielded a trans-chelated Pd(II) aryl iodide complex.
Co-reporter:Fan Jiang, Wei Han and Yun-Dong Wu
The Journal of Physical Chemistry B 2010 Volume 114(Issue 17) pp:5840-5850
Publication Date(Web):April 14, 2010
DOI:10.1021/jp909088e
Statistical analysis of coil regions in protein structures has been used to obtain the local backbone ϕ, ψ preferences of amino acids, which agree well with the NMR experiments of unfolded peptides and proteins. We analyzed the conformational features of amino acid residues in a restricted coil library of 4220 high-resolution protein crystal structures. In addition to Gly, Ala, and Pro, the ϕ, ψ distribution (Ramachandran plot) of each amino acid is analyzed with respect to three side chain conformers: g+ (χ1 ∼ −60°), g− (χ1 ∼ 60°), and t (χ1 ∼ 180°). The statistical study indicates that the effect of side chain conformations on ϕ, ψ distributions is even greater than the effect of amino acid types. On the basis of the χ1, ϕ, ψ conformational preferences, the amino acids in addition to Gly, Pro, and Ala can be divided into five types: (1) ordinary amino acids, (2) Ser, (3) Asp and Asn, (4) Val and Ile, and (5) Thr, each with distinguished χ1 rotamers. The α-helix, β-sheet, and type-I β-turn preferences of the different rotamers of various amino acid types can be captured by their intrinsic ϕ, ψ preferences from our coil library. Molecular dynamics simulations of dipeptide Ac−X−NHMe and tetrapeptide Ac−A−X−A−NHMe models give nearly the same side chain rotamer distributions. However, for many amino acids, both OPLS−AA/L and AMBER−FF03 force fields give very different χ1 rotamer distributions from the coil library. This may partially explain why dipeptide models sometimes cannot reproduce those of protein structures well. The current coil library analysis may be valuable in improving the force field for protein simulations.
Co-reporter:Lin-Tai Da, Jun-Min Quan and Yun-Dong Wu
The Journal of Physical Chemistry B 2009 Volume 113(Issue 43) pp:14536-14543
Publication Date(Web):October 8, 2009
DOI:10.1021/jp9081239
As the initial step of the entry of HIV-1 into cells, the interaction of CD4 with gp120 is a central area of concern in HIV-1 biology and intervention studies. CD4 binding induces large conformational changes to gp120, such as the formation of the bridging sheet between the V1/V2 stem and β20/β21. Understanding the dynamic process and the mechanism that leads to the formation of the bridging sheet is important. We have studied the formation of the bridging sheet via extensive molecular dynamics simulations on a modeled intermediate state. The intermediate state is derived from the crystal structure of the gp120/CD4 complex with rotation of the α1 helix and separation of the V1/V2 stem from β20/β21. The molecular dynamics simulations reveal that CD4-bound gp120 leads to the refolding of the bridging sheet but the CD4-free gp120 leads to structures similar to unliganded structures of SIV gp120. The bridging sheet refolds with the S375W mutant, but it does not refold with the W112A and S375W/T257S mutants. Our simulation results are in agreement with experimental observations, and they reveal the importance of the formation of the hydrophobic core to the conformational change of CD4-induced gp120. Other experimental observations, including variations in the binding entropy in different mutants, are also qualitatively reproduced.
Co-reporter:Yu Lan, Lujiang Deng, Jing Liu, Can Wang, Olaf Wiest, Zhen Yang and Yun-Dong Wu
The Journal of Organic Chemistry 2009 Volume 74(Issue 14) pp:5049-5058
Publication Date(Web):May 26, 2009
DOI:10.1021/jo900919v
Density functional theory calculations and experimental studies have been carried out on the intramolecular Pauson−Khand-Type reaction mediated by a PdCl2-thiourea catalyst, which proceeds under mild reaction conditions and provides a useful alternative to traditional Pauson−Khand reactions. The classical mechanism of the Pauson−Khand reaction involving the alkyne/alkene C−C bond formation as the key step has been found to be energetically unfavorable and is not in line with the experimental observations. A novel reaction mechanism has been proposed for the reaction. The first step involves the cis-halometalation of the alkyne, followed by sequential alkene and carbonyl insertion. The rate-determining fourth step is an intramolecular C−Cl oxidative addition, leading to a PdIV intermediate. A C−C bond formation by reductive elimination completes the reaction. The mechanism is in agreement with the key experimental observations including (1) the need of a chloride for catalytic activity and the absence of catalysis with Pd(OAc)2 alone; (2) the rate acceleration by the addition of LiCl; both with PdCl2 and Pd(OAc)2 catalysts; and (3) the preferred formation of the trans diastereomer in substituted cases. The cis halometalation and the formation and stability of the PdIV intermediate is studied in detail and provides general insights into these novel steps.
Co-reporter:Chun-Shan Zuo, Olaf Wiest and Yun-Dong Wu
The Journal of Physical Chemistry A 2009 Volume 113(Issue 43) pp:12028-12034
Publication Date(Web):August 31, 2009
DOI:10.1021/jp905865g
The use of correct ion radii is essential for the calculation of free energies of hydration using continuum models. A simple method for the fitting of the ion radii for ions in aqueous solution, which is a consistently difficult problem for implicit solvent models, is described. A new set of ionic radii based on experimental ionic hydration free energies for use in the integral equation formalism of polarizable continuum model (IEFPCM) is derived using B3LYP calculations with a 6-311++g** basis set for Li, Na, K, Be, Mg, and Ca and a SDD basis set for all other metals. The new radii reproduce the experimental stability constants of metal ions and their pyridine, 2,2′-bipyridine, and 1,10-phenanthroline complexes in aqueous solution significantly better than the results obtained using the default UAHF ion radii. The standard deviation (SD) of binding free energies between the calculations and experiments for the metal−ligand complexes in aqueous solution is 3.7 kcal/mol, while the mean unsigned error (MUE) is 3.1 kcal/mol. These results improve on the standard UFF radii for metal atoms, in which the MUE and the SD are 30.4 and 16.9 kcal/mol, respectively. The new ionic radii greatly improve the computational tools for the study of a variety of metals with ligands such as pyridines or calixarenes, which have found significant interest in materials science and for the removal of toxic metals.
Co-reporter:Xiang Li, Yun-Dong Wu and Dan Yang
Accounts of Chemical Research 2008 Volume 41(Issue 10) pp:1428
Publication Date(Web):September 12, 2008
DOI:10.1021/ar8001393
Naturally occurring peptides serve important functions as enzyme inhibitors, hormones, neurotransmitters, and immunomodulators in many physiological processes including metabolism, digestion, pain sensitivity, and the immune response. However, due to their conformational flexibility and poor bioavailability, such peptides are not generally viewed as useful therapeutic agents in clinical applications. In an effort to solve these problems, chemists have developed peptidomimetic foldamers, unnatural oligomeric molecules that fold into rigid and well-defined secondary structures mimicking the structures and biological functions of these natural peptides. We have designed peptidomimetic foldamers that give predictable, backbone-controlled secondary structures irrespective of the nature of the side chains. This Account presents our efforts to develop a novel class of peptidomimetic foldamers comprising α-aminoxy acids and explore their applications in the simulation of ion recognition and transport processes in living systems. Peptides constructed from α-aminoxy acids fold according to the following rules: (1) A strong intramolecular eight-membered-ring hydrogen bond forms between adjacent α-aminoxy acid residues (the α N−O turn). The chirality of the α-carbon, not the nature of the side chains, determines the conformation of this chiral N−O turn. (2) While homochiral oligomers of α-aminoxy acids form an extended helical structure (1.88 helix), heterochiral ones adopt a bent reverse turn structure. (3) In peptides of alternating α-amino acids and α-aminoxy acids, the seven-membered-ring intramolecular hydrogen bond, that is, the γ-turn, is initiated by a succeeding α N−O turn. Thus, this type of peptide adopts a novel 7/8 helical structure. In investigating the potential applications of α-aminoxy acids, we have found that the amide NH units of α-aminoxy acids are more acidic than are regular amide NH groups, which makes them better hydrogen bond donors when interacting with anions. This property makes α-aminoxy acids ideal building blocks for the construction of anion receptors. Indeed, we have constructued both cyclic and acyclic anion receptors that have strong affinities and good (enantio-)selectivities toward chloride (Cl−) and chiral carboxylate ions. Taking advantage of these systems’ preference for Cl− ions, we have also employed α-aminoxy acid units to construct a synthetic Cl− channel that can mediate the passage of Cl− ions across cell membranes. Continued studies of these peptidomimetic systems built from α-aminoxy acids should lead to a broad range of applications in chemistry, biology, medicine, and materials science.
Co-reporter:Yun-Dong Wu, Wei Han, De-Ping Wang, Yi Gao and Yi-Lei Zhao
Accounts of Chemical Research 2008 Volume 41(Issue 10) pp:1418
Publication Date(Web):October 2, 2008
DOI:10.1021/ar800070b
Unlike α-amino acids, peptides formed from β-amino acids (β-peptides) display stability toward enzymatic degradation and may form turns and helices with as few as four residues. Because both the Cα and Cβ of the β-amino acid may bear substituents, a large number of β-amino acids can be synthesized. β-Peptides form various well-defined secondary structures, including 14-helix, 12-helix, 10/12-helix, 10-helix, 8-helix, turn structures, sheets, and hairpins. For all of these reasons, β-amino acids have been increasingly used as building blocks for molecular design and pharmaceutical applications. To explain the conformational features of β-peptides, several quantum mechanics and molecular dynamics studies that rationalize the observed conformational features have been reported. However, a systematic account that unifies various factors critical to the conformational features is still lacking. In this Account, we present a detailed analysis of the conformational features of various β-peptides. We start by studying the basic local conformational features of β-peptides using di- and tripeptide models. Then, various secondary structures of unsubstituted β-peptides with differing numbers of residues are investigated using a repeating unit approach to derive the intrinsic backbone conformational features. We find that the 10/12-helix is intrinsically most stable for the β-peptide backbone. The 14-helix, 12-helix, and 10-helix structures have similar stabilities for β-peptide backbones of four to six residues. The substituent effects on the stabilities of β-peptide secondary structures are then analyzed. Combined with the substituent effect and the intrinsic backbone preferences, all experimental observations of secondary structure formation can be understood. For example, the 10/12-helix is favored for like-β2/β3-peptides, unlike-β3/β3-peptides, and β3/β-hGly-peptides because these substitution patterns do not cause steric problems for the 10/12-helix. β3-Peptides, β2-peptides, and β2,3-peptides favor the 14-helix because the substituents in these peptides benefit the 14-helix the most but significantly destabilize the 10/12-helix. Because the 10/12-helix is intrinsically favored and has two favorable positions in each residue for substituents, many more hybrid β-peptides are predicted to exist in this secondary structure, which suggests the need for further experiments. These results are valuable for determining the best use of these building blocks in the design of well-structured molecules with desirable chemical functions.
Co-reporter:Wei Han, Cheuk-Kin Wan and Yun-Dong Wu
Journal of Chemical Theory and Computation 2008 Volume 4(Issue 11) pp:1891-1901
Publication Date(Web):October 15, 2008
DOI:10.1021/ct800184c
Recently, we reported that molecular dynamics (MD) simulations using a coarse-grained (CG) peptide model coupled with a CG water model are able to reproduce many of the structural and thermodynamic features of short peptides with nonpolar side chains at 103 times the normal speed (JCTC, 2007, 3, 2146−2161). To further develop a CG protein model for MD simulations, we systematically parametrized the side chains of all 20 naturally occurring amino acids. We developed the parameters by fitting the dihedral potentials of 13 small molecules, the densities and self-solvation free energies of liquids of eight organic molecules, and the hydration free energies of 35 small organic molecules. In a set of 11 classes of compounds (105 in total) including alkanes, alcohols, ethers, ketones/aldehydes, amines, amides, aromatics, carboxylic acids, sulfides/thiols, alkyl ammoniums, and carboxylate ions, the average error in the calculated hydration free energies compared with experimental results is about 1.4 kJ/mol. The average error in the calculated transfer free energies of the 19 side-chain analogues of amino acids from cyclohexane to water is about 2.2 kJ/mol. These results are comparable to the results of all-atom models.
Co-reporter:Can Wang and Yun-Dong Wu
Organometallics 2008 Volume 27(Issue 23) pp:6152-6162
Publication Date(Web):November 6, 2008
DOI:10.1021/om8004178
Density functional theory (DFT) calculations have been carried out to understand the mechanism of the Ru3(CO)12-catalyzed Pauson−Khand-type [2+2+1] reaction and related [2+2+1+1] cycloadditions. The geometries were optimized using the BP86/6-31G*(SDD for Ru) method, and the energies were evaluated with the 6-311+G*(SDD) basis set. We found that these reactions are initiated by a CO−alkyne coupling, forming a ruthenacyclobutenone intermediate, and the widely accepted alkene−alkyne coupling pathway has a much higher activation energy. In the intermolecular reaction between alkene and alkyne, the formation of quinones and hydroquinones through [2+2+1+1] cycloadditions is more favorable than the Pauson−Khand-type reaction, while the intramolecular reaction with 1,6-enyne leads to a favorable Pauson−Khand-type reaction. These results are in agreement with experimental observations. For the [2+2+1+1] cycloadditions we found that the formation of quinones is favored over the formation of hydroquinones due to the preferred insertion of alkynes, which can be attributed to the preferred orbital interaction between the π orbital of the alkyne moiety and the d orbital of the metal center.
Co-reporter:Xin-Hao Zhang, Shu-Jian Chen, Hu Cai, Hee-Jung Im, Tianniu Chen, Xianghua Yu, Xuetai Chen, Zhenyang Lin, Yun-Dong Wu and Zi-Ling Xue
Organometallics 2008 Volume 27(Issue 6) pp:1338-1341
Publication Date(Web):February 19, 2008
DOI:10.1021/om701173m
Reactions of MCl5 (M = Nb, Ta) with LiNMe2 have been found to yield M(NMe2)4(η2-MeNCH2NMe2) (M = Nb, 2a; Ta, 2b) containing a chelating ligand (dimethylaminomethylene)methylamide, as confirmed by NMR spectroscopy, DFT calculations, and their reactivity studies.
Co-reporter:Hai-Chuan Liu, Xin-Hao Zhang, Changsheng Wang, Yun-Dong Wu and Shihe Yang
Physical Chemistry Chemical Physics 2007 vol. 9(Issue 5) pp:607-615
Publication Date(Web):11 Dec 2006
DOI:10.1039/B614445D
Gas phase complexes Mg˙+(2,6-difluoropyridine) (1) and Mg˙+(pentafluoropyridine) (2) have been subjected to photodissociation in the spectral range of ∼230–440 nm. Except for the evaporative photofragment Mg˙+, the primary photoproduct for 1 is C5H3N˙+, which is associated with the rupture of two C–F bonds by the photoexcited Mg˙+, forming very stable MgF2. In contrast, the direct loss of MgF+ is more favorable for 2 due to fluorine substitution. Given enough energy, C5H3N˙+ can undergo decomposition to form C4H2˙+ and HCN. These results are very different from those for Mg˙+(2-fluoropyridine), highlighting the significance of the additional F at C6 of 1 and 2. Density functional theory (DFT) calculations have been employed to examine the geometries and energetics of the complexes as well as relevant reaction mechanisms. All of the complexes feature the direct attachment of Mg˙+ to the N atom. The key intermediate is found to be FMg+(C5HxF4−xN) (x = 3 or 0), which can lead to the formation of MgF+ directly or MgF2 through activation of another C–F bond adjacent to N, producing the pyridyne radical cations. However, hydrogen-transfer prior to the rupture of the second C–F bond followed by ring-opening of C5H3N˙+ may result in the formation of chain forms of C5H3N˙+. The influence of the fluorine substitution on the competition of the two routes have been demonstrated.
Co-reporter:Yun-Dong Wu, Zhi-Hui Peng
Inorganica Chimica Acta 2003 Volume 345() pp:241-254
Publication Date(Web):10 March 2003
DOI:10.1016/S0020-1693(02)01353-1
The mechanism of ring-opening metathesis reaction of molybdenum alkylidene complexes with norbornadiene has been studied by quantum mechanics method. Both cycloaddition and ring-opening reactions have a twisted trigonal bipyramidal transition structure. The cycloaddition is the rate-determining step while the ring-opening step has very low activation energy. When the alkoxyl ligands are small (OCH3, OCF3), there is only a small preference for norbornadiene attacking syn to the imido group. When the size of the alkoxyl ligands increases, the syn preference for norbornadiene addition becomes more favorable. Our calculations strongly support Schrock's proposal for the stereoselectivity of ring-opening metathesis polymerization. That is, the stereochemistry is determined by the relative rates of alkylidene syn/anti interconversion and polymerization reaction.The mechanism of ring-opening metathesis polymerization of norbornadiene mediated by molybdenum alkylidene complexes has been studied by quantum mechanics method. Both cycloaddition and ring-opening reactions have a twisted trigonal bipyramidal transition structure. Our study supports Schrock's proposal for the stereoselectivity of ROMP. That is, the stereochemistry is determined by the relative rates of alkylidene syn/anti interconversion and polymerization.
Co-reporter:Hai-Chuan Liu, Xin-Hao Zhang, Changsheng Wang, Yun-Dong Wu and Shihe Yang
Physical Chemistry Chemical Physics 2007 - vol. 9(Issue 5) pp:NaN615-615
Publication Date(Web):2006/12/11
DOI:10.1039/B614445D
Gas phase complexes Mg˙+(2,6-difluoropyridine) (1) and Mg˙+(pentafluoropyridine) (2) have been subjected to photodissociation in the spectral range of ∼230–440 nm. Except for the evaporative photofragment Mg˙+, the primary photoproduct for 1 is C5H3N˙+, which is associated with the rupture of two C–F bonds by the photoexcited Mg˙+, forming very stable MgF2. In contrast, the direct loss of MgF+ is more favorable for 2 due to fluorine substitution. Given enough energy, C5H3N˙+ can undergo decomposition to form C4H2˙+ and HCN. These results are very different from those for Mg˙+(2-fluoropyridine), highlighting the significance of the additional F at C6 of 1 and 2. Density functional theory (DFT) calculations have been employed to examine the geometries and energetics of the complexes as well as relevant reaction mechanisms. All of the complexes feature the direct attachment of Mg˙+ to the N atom. The key intermediate is found to be FMg+(C5HxF4−xN) (x = 3 or 0), which can lead to the formation of MgF+ directly or MgF2 through activation of another C–F bond adjacent to N, producing the pyridyne radical cations. However, hydrogen-transfer prior to the rupture of the second C–F bond followed by ring-opening of C5H3N˙+ may result in the formation of chain forms of C5H3N˙+. The influence of the fluorine substitution on the competition of the two routes have been demonstrated.
Co-reporter:Yun-Fang Yang, Ting Shi, Xin-Hao Zhang, Zong-Xun Tang, Zhen-Yi Wen, Jun-Min Quan and Yun-Dong Wu
Organic & Biomolecular Chemistry 2011 - vol. 9(Issue 16) pp:NaN5855-5855
Publication Date(Web):2011/05/23
DOI:10.1039/C1OB05501A
Density functional theory calculations (B3LYP) have been carried out to understand the mechanism and stereochemistry of an asymmetric reductive aldol reaction of benzaldehyde and tert-butyl acrylate with hydrosilanes catalyzed by Rh(Phebox-ip)(OAc)2(OH2). According to the calculations, the reaction proceeds via five steps: (1) oxidative addition of hydrosilane, (2) hydride migration to carbon–carbon double bond of tert-butyl acrylate, which determines the chirality at C2, (3) tautomerization from rhodium bound C-enolate to rhodium bound O-enolate, (4) intramolecular aldol reaction, which determines the chirality at C3 and consequently the anti/syn-selectivity, and (5) reductive elimination to release aldol product. The hydride migration is the rate-determining step with a calculated activation energy of 23.3 kcal mol−1. In good agreement with experimental results, the formation of anti-aldolates is found to be the most favorable pathway. The observed Si-facial selectivity in both hydride migration and aldol reaction are well-rationalized by analyzing crucial transition structures. The Re-facial attack transition state is disfavored because of steric hindrance between the isopropyl group of the catalyst and the tert-butyl acrylate.
Co-reporter:Fan Jiang, Wei Han and Yun-Dong Wu
Physical Chemistry Chemical Physics 2013 - vol. 15(Issue 10) pp:NaN3428-3428
Publication Date(Web):2012/12/12
DOI:10.1039/C2CP43633G
The local conformational (ϕ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the PII conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The ϕ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ1 rotamer preference and χ1-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the ϕ, ψ plot of alanine with the influence of side-chain χ1 rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental 3J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix–coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.