PubMed Health⌕ Search

Biomedical subjects

Weitao Yang

Publications and source records attributed to Weitao Yang.

At least 19 recordsLinked to original sources
PubMed2007

Computational study on the relative acidity of acetic acid by the QM/MM method combined with the theory of energy representation.

We have applied the quantum mechanical/molecular mechanical (QM/MM) method combined with the theory of energy representation (ER) to study the acidity of acetic acid in aqueous solution. We have focused our attention on the relative acidity DeltapK(a) of the molecule with respect to water solvent to circumvent the ambiguity of the solvation free energies of the molecular species referred to as proton. The value of DeltapK(a) for the acetic acid has been computed as -11.5 when we adopt the free energy change in the gas phase obtained by the B3LYP functional, which is in excellent agreement with the experimental value of -11.0. It has been demonstrated that the QM/MM-ER approach recently developed gives an adequate description for the solvation free energies related to the acidity/basicity calculations of organic molecules.

Acetic Acid↗
PubMed2007

Designing molecules with optimal properties using the linear combination of atomic potentials approach in an AM1 semiempirical framework.

The linear combination of atomic potentials (LCAP) approach is implemented in the AM1 semiempirical framework and is used to design molecular structures with optimized properties. The optimization procedure uses property derivative information to search molecular space and thus avoid direct enumeration and evaluation of each molecule in a library. Two tests are described: the optimization of first hyperpolarizabilities of substituted aromatics and the optimization of a figure of merit for n-type organic semiconductors.

Journal Article↗
PubMed2006

Automatic integration of the reaction path using diagonally implicit Runge-Kutta methods.

The diagonally implicit Runge-Kutta framework is shown to be a general form for constructing stable, efficient steepest descent reaction path integrators, of any order. With this framework tolerance driven, adaptive step-size methods can be constructed by embedding methods to obtain error estimates of each step without additional computational cost. There are many embedded and nonembedded, diagonally implicit Runge-Kutta methods available from the numerical analysis literature and these are reviewed for orders two, three, and four. New embedded methods are also developed which are tailored to the application of reaction path following. All integrators are summarized and compared for three systems: the Muller-Brown [Theor. Chem. Acta 53, 75 (1979)] potential and two gas phase chemical reactions. The results show that many of the methods are capable of integrating efficiently while reliably keeping the error bound within the desired tolerance. This allows the reaction path to be determined through automatic integration by only specifying the desired accuracy and transition state.

Algorithms↗
PubMed2006

Many-electron self-interaction error in approximate density functionals.

One of the most important challenges in density functional theory (DFT) is the proper description of fractional charge systems relating to the self-interaction error (SIE). Traditionally, the SIE has been formulated as a one-electron problem, which has been addressed in several recent functionals. However, these recent one-electron SIE-free functionals, while greatly improving the description of thermochemistry and reaction barriers in general, still exhibit many of the difficulties associated with SIE. Thus we emphasize the need to surpass this limit and shed light on the many-electron SIE. After identifying the sufficient condition for functionals to be free from SIE, we focus on the symptoms and investigate the performance of most popular functionals. We show that these functionals suffer from many-electron SIE. Finally, we give a SIE classification of density functionals.

Journal Article↗
PubMed2006

Density functional theory investigation of the polarizability and second hyperpolarizability of polydiacetylene and polybutatriene chains: treatment of exact exchange and role of correlation.

The static polarizability and second hyperpolarizability of increasingly large polydiacetylene and polybutatriene (PBT) chains have been evaluated using the optimized effective potential for exact exchange (OEP-EXX) method developed by Yang and Wu [Phys. Rev. Lett. 89, 143002 (2002)], where the unknown part of the effective potential is expressed as a linear combination of Gaussian functions. Various conventional atomic orbital basis sets were employed for the exchange potential (X basis) as well as for the Kohn-Sham orbitals [molecular orbital (MO) basis]. Our results were compared to coupled-perturbed Hartree-Fock (CPHF) calculations and to ab initio correlated values obtained at various levels of approximation. It turns out that (a) small conventional basis sets are, in general, unsatisfactory for the X basis; (b) the performance of a given X basis depends on the MO basis and is generally improved when using a larger MO basis; (c) these effects are exaggerated for the second hyperpolarizability compared to the polarizability; (d) except for the second hyperpolarizability of PBT chains, using 6-311++G** for the X basis gives reasonable agreement with the CPHF results for all MO basis sets; (e) our results suggest that in the limit of a complete X basis the OEP-EXX values may approach the CPHF data; and (f) in general, the quality of a given conventional X basis degrades with the length of the oligomer, which correlates with the fact that the number of X basis functions becomes a smaller fraction of the number required to reproduce exactly the finite-basis-set Hartree-Fock energies. Linear and especially nonlinear electric field responses constitute a very stringent test for assessing the quality of functionals and potentials; appropriately tailored basis sets are needed to describe the latter. Finally, this study further highlights the importance of electron correlation effects on linear and nonlinear responses, for which correlated functionals with OEP are required.

Journal Article↗
PubMed2006

Direct detection of the formation of V-amylose helix by single molecule force spectroscopy.

An important polysaccharide, amylose crystallizes as a regular single left-handed helix from a propanol, butanol, or iodine solution. However, its solution structure remains elusive because amylose does not form molecular solutions in these solvents, and standard spectroscopic techniques cannot be exploited to determine its structure. Using AFM, we forced individual amylose chains adsorbed to a surface to enter these poor solvents and carried out stretch-release measurements on them in solution. In this manner, we directly captured the formation of individual amylose helices induced by butanol and iodine. With an accuracy approaching that of X-ray diffraction on amylose crystals, we determined that the pitch of the helix in solution is 1.3 angstroms/ring. We also directly measured the force driving the formation of the helix in solution to be 50 pN. SMD simulations in explicit butanol reproduced the AFM-measured force-extension curves and revealed that the long plateau feature is caused by the rupture of O(2)n-O(6)(n+6) and O(3)n-O(6)(n+6) hydrogen bonds and by the unwinding of the helix. We also found that amylose helices formed in iodine solution are more compliant and hysteretic as compared to helices in butanol, which extend/relax reversibly. In iodine solution, the formation of the helix is inhibited by force and limited by the slow kinetics of the amylose-iodine complex. By forcing individual molecules into poor solvents and performing force spectroscopy measurements in solution, our AFM approach uniquely supplements X-ray diffraction and NMR methods for investigating solution conformations of insoluble biopolymers.

Amylose↗
PubMed2006

Simulating force-induced conformational transitions in polysaccharides with the SMD replica exchange method.

Conventional steered molecular dynamics (SMD) simulations do not readily reproduce equilibrium conditions of atomic force microscopy (AFM) stretch and release measurements of polysaccharides undergoing force-induced conformational transitions because of the gap between the timescales of computer simulations ( approximately 1 mus) and AFM measurements ( approximately 1 s). To circumvent this limitation, we propose using the replica exchange method (REM) to enhance sampling during SMD simulations. By applying REM SMD to a small polysaccharide system and comparing the results with those from AFM stretching experiments, we demonstrate that REM SMD reproduces the experimental results not only qualitatively but quantitatively, approaching near equilibrium conditions of AFM measurements. As tested in this work, hysteresis and computational time of REM SMD simulations of short polysaccharide chains are significantly reduced as compared to regular SMD simulations, making REM SMD an attractive tool for studying forced-induced conformational transitions of small biopolymer systems.

Computer Simulation↗
PubMed2006

Legendre-transform functionals for spin-density-functional theory.

We provide a rigorous proof that the Hohenberg-Kohn theorem holds for spin densities by extending Lieb's Legendre-transform formulation to spin densities. The resulting spin-density-functional theory resolves several troublesome issues. Most importantly, the present paper provides an explicit construction for the spin potentials at any point along the adiabatic connection curve, thus providing a formal basis for the use of exchange-correlation functionals of the spin density in the Kohn-Sham density-functional theory (DFT). The practical implications of this result for unrestricted Kohn-Sham DFT calculations is considered, and the existence of holes below the Fermi level is discussed. We argue that an orbital's energy tends to increase as its occupation number increases, which provides the basis for a computational algorithm for determining the occupation numbers in Kohn-Sham DFT and helps explain the origin of Hund's rules and holes below the Fermi level.

Journal Article↗
PubMed2006

A combined explicit-implicit method for high accuracy reaction path integration.

We present the use of an optimal combined explicit-implicit method for following the reaction path to high accuracy. This is in contrast to most purely implicit reaction path integration algorithms, which are only efficient on stiff ordinary differential equations. The defining equation for the reaction path is considered to be stiff, however, we show here that the reaction path is not uniformly stiff and instead is only stiff near stationary points. The optimal algorithm developed in this work is a combination of explicit and implicit methods with a simple criterion to switch between the two. Using three different chemical reactions, we combine and compare three different integration methods: the implicit trapezoidal method, an explicit stabilized third order algorithm implemented in the code DUMKA3 and the traditional explicit fourth order Runge-Kutta method written in the code RKSUITE. The results for high accuracy show that when the implicit trapezoidal method is combined with either explicit method the number of energy and gradient calculations can potentially be reduced by almost a half compared with integrating either method alone. Finally, to explain the improvements of the combined method we expand on the concepts of stability and stiffness and relate them to the efficiency of integration methods.

Acetaldehyde↗
PubMed2006

Nanotube-metal junctions: 2- and 3-terminal electrical transport.

We address the quality of electrical contact between carbon nanotubes and metallic electrodes by performing first-principles calculations for the electron transmission through ideal 2- and 3-terminal junctions, thus revealing the physical limit of tube-metal conduction. The structural model constructed involves surrounding the tube by the metal atoms of the electrode as in most experiments; we consider metallic (5,5) and n-doped semiconducting (10,0) tubes surrounded by Au or Pd. In the case of metallic tubes, the contact conductance is shown to approach the ideal 4e2/h in the limit of large contact area. For three-terminals, the division of flux among the different transmission channels depends strongly on the metal material. A Pd electrode has nearly perfect tube-electrode transmission and therefore turns off the straight transport along the tube. Our results are in good agreement with some recent experimental reports and clarify a fundamental discrepancy between theory and experiment.

Journal Article↗
PubMed2006

Nuclear quantum effects on an enzyme-catalyzed reaction with reaction path potential: proton transfer in triosephosphate isomerase.

Nuclear quantum mechanical effects have been examined for the proton transfer reaction catalyzed by triosephosphate isomerase, with the normal mode centroid path integral molecular dynamics based on the potential energy surface from the recently developed reaction path potential method. In the simulation, the primary and secondary hydrogens and the C and O atoms involving bond forming and bond breaking were treated quantum mechanically, while all other atoms were dealt classical mechanically. The quantum mechanical activation free energy and the primary kinetic isotope effects were examined. Because of the quantum mechanical effects in the proton transfer, the activation free energy was reduced by 2.3 kcal/mol in comparison with the classical one, which accelerates the rate of proton transfer by a factor of 47.5. The primary kinetic isotope effects of kH/kD and kH/kT were estimated to be 4.65 and 9.97, respectively, which are in agreement with the experimental value of 4+/-0.3 and 9. The corresponding Swain-Schadd exponent was predicted to be 3.01, less than the semiclassical limit value of 3.34, indicating that the quantum mechanical effects mainly arise from quantum vibrational motion rather than tunneling. The reaction path potential, in conjunction with the normal mode centroid molecular dynamics, is shown to be an efficient computational tool for investigating the quantum effects on enzymatic reactions involving proton transfer.

Carbon↗
PubMed2006

Designing molecules by optimizing potentials.

The astronomical number of accessible discrete chemical structures makes rational molecular design extremely challenging. We formulate the design of molecules with specific tailored properties as performing a continuous optimization in the space of electron-nuclear attraction potentials. The optimization is facilitated by using a linear combination of atomic potentials (LCAP), a general framework that creates a continuous property landscape from an otherwise unlinked set of discrete molecular-property values. A demonstration of this approach is given for the optimization of molecular electronic polarizability and hyperpolarizability. We show that the optimal structures can be determined without enumerating and separately evaluating the characteristics of the combinatorial number of possible structures, a process that would be much slower. The LCAP approach may be used with quantum or classical Hamiltonians, suggesting possible applications to drug design and new materials discovery.

Journal Article↗
PubMed2006

Self-interaction-free exchange-correlation functional for thermochemistry and kinetics.

We develop a self-interaction-free exchange-correlation functional which is very accurate for thermochemistry and kinetics. This is achieved by theoretical construction of the functional form and nonlinear fitting. We define a simple interpolation of the adiabatic connection that uses exact exchange, generalized gradient approximation (GGA) and meta-GGA functionals. The performance is optimized by fitting a small number of empirical parameters. Overall the new functional improves significantly upon hybrids and meta-GGAs while correctly describing one-electron systems. The mean absolute error on a large set of reaction barriers is reduced to 1.99 kcal/mol.

Kinetics↗
PubMed2006

Quadratic string method for determining the minimum-energy path based on multiobjective optimization.

Based on a multiobjective optimization framework, we develop a new quadratic string method for finding the minimum-energy path. In the method, each point on the minimum-energy path is minimized by integration in the descent direction perpendicular to path. Each local integration is done on a quadratic surface approximated by a damped Broyden-Fletcher-Goldfarb-Shanno updated Hessian, allowing the algorithm to take many steps between energy and gradient calls. The integration is performed with an adaptive step-size solver, which is restricted in length to the trust radius of the approximate Hessian. The full algorithm is shown to be capable of practical superlinear convergence, in contrast to the linear convergence of other methods. The method also eliminates the need for predetermining such parameters as step size and spring constants, and is applicable to reactions with multiple barriers. The effectiveness of this method is demonstrated for the Muller-Brown potential, a seven-atom Lennard-Jones cluster, and the enolation of acetaldehyde to vinyl alcohol.

Acetaldehyde↗
PubMed2006

Theoretical and experimental determination on two substrates turned over by 4-oxalocrotonate tautomerase.

Quantum mechanical/molecular mechanical (QM/MM) calculations and experimental kinetic studies have been performed on 4-oxalocrotonate tautomerase (4OT) for two different substrates, 2-hydroxymuconate (2HM) and 2-oxo-4-hexenedioate (2o4hex). Potential (deltaE) and free energy (deltaG) paths for both steps of the reaction using both substrates were calculated to determine the free energy barriers and compared to the experimental values obtained from the kinetic studies via the transition state theory. In the first step, a proton from the hydroxyl oxygen on the second carbon of 2HM, or from the third carbon of 2o4hex, is abstracted by Pro-1. In the second step, the proton is transferred to the fifth carbon of the substrate to form the product, 2-oxo-3-hexenedioate (2o3hex). For both substrates we obtain a calculated deltaG of approximately 13 kcal/mol, in agreement with experimental determinations. The calculated free energy barrier difference deltaG2o4hex - deltaG2HM (deltadeltaG) is 0.87 kcal/mol. We obtained an experimental deltadeltaG of 0.85 kcal/mol. These results suggest that 2HM is turned over faster than 2o4hex by 4OT. However, these energy differences are so small that both 2HM and 2o4hex need to be taken into account in considering the mechanism of catalysis of 4OT.

Isomerases↗
PubMed2006

Organometallic molecular rectification.

We study the rectification of current through a single molecule with an intrinsic spatial asymmetry. The molecule contains a cobaltocene moiety in order to take advantage of its relatively localized and high-energy d states. A rectifier with large voltage range, high current, and low threshold can be realized. The evolution of molecular orbitals under both forward and reverse biases is captured in a self-consistent nonequilibrium Green function plus density functional theory description. Our calculations demonstrate the plausibility of making excellent molecular diodes by using metallocenes, pointing to a fruitful class of molecules.

Journal Article↗
PubMed2005

Experimental validation of the docking orientation of Cdc25 with its Cdk2-CycA protein substrate.

Cdc25 phosphatases are key activators of the eukaryotic cell cycle and compelling anticancer targets because their overexpression has been associated with numerous cancers. However, drug discovery targeting these phosphatases has been hampered by the lack of structural information about how Cdc25s interact with their native protein substrates, the cyclin-dependent kinases. Herein, we predict a docked orientation for Cdc25B with its Cdk2-pTpY-CycA protein substrate by a rigid-body docking method and refine the docked models with full-scale molecular dynamics simulations and minimization. We validate the stable ensemble structure experimentally by a variety of in vitro and in vivo techniques. Specifically, we compare our model with a crystal structure of the substrate-trapping mutant of Cdc25B. We identify and validate in vivo a novel hot-spot residue on Cdc25B (Arg492) that plays a central role in protein substrate recognition. We identify a hot-spot residue on the substrate Cdk2 (Asp206) and confirm its interaction with hot-spot residues on Cdc25 using hot-spot swapping and double mutant cycles to derive interaction energies. Our experimentally validated model is consistent with previous studies of Cdk2 and its interaction partners and initiates the opportunity for drug discovery of inhibitors that target the remote binding sites of this protein-protein interaction.

Binding Sites↗