Computation of the permanent dipole moment of α-chymotrypsin from linear-scaling semiempirical quantum mechanical methods

Signal-Transducing Proteins for Nanoelectronics

This aim of this article is to provide novel paradigms for 21st century nanoelectronics by taking inspiration from the biology of signal transduction events where Nature has solved many complex problems, particularly those concerned with signal integration and amplification.

RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I

Twenty years ago, the landmark AM1 was introduced, and has since had an increasingly wide following among chemists due to its consistently good results and time-tested reliability--being presently available in countless computational quantum chemistry programs. However, semiempirical molecular orbital models still are of limited accuracy and need to be improved if the full potential of new linear scaling techniques, such as MOZYME and LocalSCF, is to be realized. Accordingly, in this article we present RM1 (Recife Model 1): a reparameterization of AM1. As before, the properties used in the parameterization procedure were: heats of formation, dipole moments, ionization potentials and geometric variables (bond lengths and angles). Considering that the vast majority of molecules of importance to life can be assembled by using only six elements: C, H, N, O, P, and S, and that by adding the halogens we can now build most molecules of importance to pharmaceutical research, our training set consisted of 1736 molecules, representative of organic and biochemistry, containing C, H, N, O, P, S, F, Cl, Br, and I atoms. Unlike AM1, and similar to PM3, all RM1 parameters have been optimized. For enthalpies of formation, dipole moments, ionization potentials, and interatomic distances, the average errors in RM1, for the 1736 molecules, are less than those for AM1, PM3, and PM5. Indeed, the average errors in kcal x mol(-1) of the enthalpies of formation for AM1, PM3, and PM5 are 11.15, 7.98, and 6.03, whereas for RM1 this value is 5.77. The errors, in Debye, of the dipole moments for AM1, PM3, PM5, and RM1 are, respectively, 0.37, 0.38, 0.50, and 0.34. Likewise, the respective errors for the ionization potentials, in eV, are 0.60, 0.55, 0.48, and 0.45, and the respective errors, in angstroms, for the interatomic distances are 0.036, 0.029, 0.037, and 0.027. The RM1 average error in bond angles of 6.82 degrees is only slightly higher than the AM1 figure of 5.88 degrees, and both are much smaller than the PM3 and PM5 figures of 6.98 degrees and 9.83 degrees, respectively. Moreover, a known error in PM3 nitrogen charges is corrected in RM1. Therefore, RM1 represents an improvement over AM1 and its similar successor PM3, and is probably very competitive with PM5, which is a somewhat different model, and not fully disclosed. RM1 possesses the same analytical construct and the same number of parameters for each atom as AM1, and, therefore, can be easily implemented in any software that already has AM1, not requiring any change in any line of code, with the sole exception of the values of the parameters themselves.

The electronic structure of human erythropoietin as an aid in the design of oxidation-resistant therapeutic proteins

The electronic structure of human erythropoietin (HuEPO) has been investigated with the aid of quantum mechanical calculations. The results indicate that the protein is highly polarized and its permanent dipole moment has a magnitude of 471 D. The HOMO of HuEPO is localized on Trp51, which stays in close proximity to Met54. Three oxidation-resistant mutants of HuEPO (W51F, M54V, and W51F-M54V) have been modeled and their electronic structures are compared to that of the native protein. Among them, the W51F mutation is predicted to be the most effective in increasing the oxidation potential of the protein.

A quantum mechanical study on phosphotyrosyl peptide binding to the SH2 domain of p56lck tyrosine kinase with insights into the biochemistry of intracellular signal transduction events

A study on the interaction between a phosphotyrosyl peptide with the SH2 domain of Lck kinase has been undertaken with the aid of semiempirical linear-scaling quantum mechanical methods. The structure of this complex has been solved at atomic resolution and, hence, it represents the ideal candidate for studying the charge deformation effects induced by the phosphopeptide on the binding site. Substantial changes in the charge of amino acid residues located in the binding pocket of the protein are observed upon ligand binding. More specifically, our quantum chemical calculations indicate that H-bonds involving charged side-chains are subject to consistent charge deformation effects whereas those forming salt bridges are unaffected by ligand binding. Furthermore, ligand binding has the effect of changing both the magnitude and direction of the protein's macrodipole, which rotates approximately 150 degrees with respect that of the unliganded protein. This suggests that a change in the polarization state of the protein might acts as a switch during the transmission of intracellular signals. The binding energy calculated with the aid of the COSMO solvation model corresponds to about -200 kcal/mol, most of which is attributed to the interaction of the phosphotyrosine head with the amino acid chains located in the binding site of the SH2 domain.