ABSTRACT
Experimentally determined mean pK(a) values of carboxyl residues located at the N-termini of alpha-helices are lower than their overall mean values. Here, we perform three types of analyses to account for this phenomenon. We estimate the magnitude of the helix macrodipole to determine its potential role in lowering carboxyl pK(a) values at the N-termini. No correlation between the magnitude of the macrodipole and the pK(a) values is observed. Using the pK(a) program propKa we compare the molecular surroundings of 18 N-termini carboxyl residues versus 233 protein carboxyl groups from a previously studied database. Although pK(a) lowering interactions at the N-termini are similar in nature to those encountered in other protein regions, pK(a) lowering backbone and side-chain hydrogen bonds appear in greater number at the N-termini. For both Asp and Glu, there are about 0.5 more hydrogen bonds per residue at the N-termini than in other protein regions, which can be used to explain their lower than average pK(a) values. Using a QM-based pK(a) prediction model, we investigate the chemical environment of the two lowest Asp and the two lowest Glu pK(a) values at the N-termini so as to quantify the effect of various pK(a) determinants. We show that local interactions suffice to account for the acidity of carboxyl residues at the N-termini. The effect of the helix dipole on carboxyl pK(a) values, if any, is marginal. Backbone amide hydrogen bonds constitute the single biggest contributor to the lowest carboxyl pK(a) values at the N-termini. Their estimated pK(a) lowering effects range from about 1.0 to 1.9 pK(a) units.
Subject(s)
Carbon Dioxide/chemistry , Models, Molecular , Peptide Fragments/chemistry , Protein Structure, Secondary , Databases, Protein , Free Radicals/chemistry , Hydrogen Bonding , Hydrogen-Ion Concentration , Protein Structure, TertiaryABSTRACT
The presence of a low-barrier hydrogen bond (LBHB) in aspartyl proteases and its implications in drug design have been the subject of intense study. Here, we present a combined quantum mechanical/molecular mechanical (QM/MM)-Numerov procedure and use it to characterize the O(δ)(1)-H-O(δ)(1) hydrogen bond (HB) in unbound HIV protease. The QM/MM scheme fully traces the shape of the HB's potential energy curve. The potential is used to obtain numerical solutions to the wave functions and vibrational energies of hydrogen, deuterium, and tritium. The vibrational eigenfunctions are used to compute expectation values for interatomic distances and vibrationally and thermally averaged spectroscopic properties of the O(δ)(1)-H-O(δ)(1) HB. Our work corroborates previous results by Piana and Carloni who found a LBHB via an ab initio molecular dynamics simulation (Piana, S.; Carloni, P. Proteins 2000, 39, 26-36). Our predictions of isotope effects on the chemical shift of unbound HIV protease are consistent with experimental measurements in similar HBs. These results support the predictive power of this method and its potential use in screening inhibitors of aspartyl proteases.
ABSTRACT
We describe the development and application of a computational method for the prediction and rationalization of pKa values of ionizable residues in proteins, based on ab initio quantum mechanics (QM) and the effective fragment potential (EFPs) method (a hybrid QM/MM method). The theoretical developments include (1) a covalent boundary method based on frozen localized orbitals, (2) divide-and-conquer methods for the ab initio computation of protein EFPs consisting of multipoles up to octupoles and dipole polarizability tensors, (3) a method for computing vibrational free energies for a localized molecular region, and (4) solutions of the polarized continuum model of bulk solvation equations for protein-sized systems. The QM-based pKa prediction method is one of the most accurate methods currently available and can be used in cases where other pKa prediction methods fail. Preliminary analysis of the computed results indicate that many pKa values (1) are primarily determined by hydrogen bonds rather than long-range charge-charge interactions and (2) are relatively insensitive to large-scale dynamical fluctuations of the protein structure.
Subject(s)
Algorithms , Models, Chemical , Proteins/chemistry , Proteins/metabolism , Quantum Theory , Hydrogen Bonding , Ions , Protein Conformation , ThermodynamicsABSTRACT
Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of N(zeta) of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pK(a) of Lys55, and compared to other computational methodologies in which force field charges are employed.
ABSTRACT
Se realizó una recopilación de información técnica y científica relativa a Lepidium meyenii Chacón (Maca), una planta de la familia Brassicaceae que crece en Los Andes Centrales del Perú, entre los 4.000 y 4.500 msnm. Maca ha sido tradicionalmente empleada en la región altoandina debido a supuestas propiedades afrodisíacas y como estimulante de la actividad reproductiva en el hombre y sus animales. Comparativamente con otros vegetales, el análisis químico de la raíz demuestra un alto contenido de carbohidratos, proteínas, lípidos, vitaminas y minerales. Adicionalmente, aminoácidos esenciales y esteroles de alta especificidad que forman parte de neurotransmisores y hormonas, que a su vez cumplen un importante rol funcional en la actividad reproductiva de los mamíferos, se encuentran en alta concentración. Se concluye, a partir de los datos experimentales logrados por variados autores, que el valor nutricional y el efecto sinérgico de los constituyentes naturales de Lepidium meyenii (Maca) deben ser considerados para explicar su acción preventiva de los efectos deletéreos de la altitud sobre la espermatogénesis y estimulante de la fertilidad en sus variados componentes.