Electrostatics in the active site of an alpha-amylase.

The importance of electrostatics in catalysis has been emphasized in the literature for a large number of enzymes. We examined this hypothesis for the Bacillus licheniformis alpha-amylase by constructing site-directed mutants that were predicted to change the pKa values of the catalytic residues and thus change the pH-activity profile of the enzyme. To change the pKa of the catalytic residues in the active site, we constructed mutations that altered the hydrogen bonding network, mutations that changed the solvent accessibility, and mutations that altered the net charge of the molecule. The results show that changing the hydrogen bonding network near an active site residue or changing the solvent accessibility of an active site residue will very likely result in an enzyme with drastically reduced activity. The differences in the pH-activity profiles for these mutants were modest. pH-activity profiles of mutants which change the net charge on the molecule were significantly different from the wild-type pH-activity profile. The differences were, however, difficult to correlate with the electrostatic field changes calculated. In several cases we observed that pH-activity profiles shifted in the opposite direction compared to the shift predicted from electrostatic calculations. This strongly suggests that electrostatic effects cannot be solely responsible for the pH-activity profile of the B. licheniformis alpha-amylase.     

Protein effects in non-heme iron enzyme catalysis: insights from multiscale models

Many non-heme iron enzymes have similar sets of ligands but still catalyze widely different reactions. A key question is, therefore, the role of the protein in controlling reactivity and selectivity. Examples from multiscale simulations, primarily QM/MM, of both mono- and binuclear non-heme iron enzymes are used to analyze the stability of these models and what they reveal about the protein effects. Consistent results from QM/MM modeling are the importance of the hydrogen bond network to control reactivity and electrostatic stabilization of electron transfer from second-sphere residues. The long-range electrostatic effects on reaction barriers are small for many systems. In the systems where large electrostatic effects have been reported, these lead to higher barriers. There is thus no evidence of any significant long-range electrostatic effects contributing to the catalytic efficiency of non-heme iron enzymes. However, the correct evaluation of electrostatic contributions is challenging, and the correlation between calculated residue contributions and the effects of mutation experiments is not very strong. The largest benefits of QM/MM models are thus the improved active-site geometries, rather than the calculation of accurate energies. Reported differences in mechanistic predictions between QM and QM/MM models can be explained by differences in hydrogen bonding patterns in and around the active site. Correctly constructed cluster models can give results with similar accuracy as those from multiscale models, but the latter reduces the risk of drawing the wrong mechanistic conclusions based on incorrect geometries and are preferable for all types of modeling, even when using very large QM parts.     

Discovering protein−ligand chalcogen bonding in the protein data bank using endocyclic sulfur-containing heterocycles as ligand search subsets

The chalcogen bond, the noncovalent, electrostatic attraction between covalently bonded atoms in group 16 and Lewis bases, is present in protein?ligand interactions based on X-ray structures deposited in the Protein Data Bank (PDB). Discovering protein?ligand chalcogen bonding in the PDB employed a strategy that focused on searching the database for protein complexes of five-membered, heterocyclic ligands containing endocyclic sulfur with endo electron-withdrawing groups (isothiazoles; thiazoles; 1,2,3-, 1,2.4-, 1,2,5-, 1,3,4-thiadiazoles) and thiophenes with exo electron-withdrawing groups, e.g., 2-chloro, 2-bromo, 2-amino, 2-alkylthio. Out of 930 ligands investigated, 33 or 3.5% have protein?ligand S---O interactions of which 31 are chalcogen bonds and two appear to be S---HO hydrogen bonds. The bond angles for some of the chalcogen bonds found in the PDB are less than 90°, and an electrostatic model is proposed to explain this phenomenon.     

Influence of hydrogen bonds on edge-to-face interactions between pyridine molecules

Edge-to-face interactions between two pyridine molecules and the influence of simultaneous hydrogen bonding of one or both of the pyridines to water on those interactions were studied by analyzing data from ab initio calculations. The results show that the edge-to-face interactions of pyridine dimers that are hydrogen bonded to water are generally stronger than those of non-H-bonded pyridine dimers, especially when the donor pyridine forms a hydrogen bond. The binding energy of the most stable edge-to-face interacting H-bonded pyridine dimer is ?5.05 kcal/mol, while that for the most stable edge-to-face interacting non-H-bonded pyridine dimer is ?3.64 kcal/mol. The interaction energy data obtained in this study cannot be explained solely by the differences in electrostatic potential between pyridine and the pyridine–water dimer. However, the calculated cooperative effect can be predicted using electrostatic potential maps.     

DFT calculations of the anomeric and exo-anomeric effect of the hydroperoxy and peroxy groups

DFT calculations were carried out on axially and equatorially oriented 2-hydroperoxy and 2-peroxy tetrahydropyran, cyclohexyl hydroperoxide, hydroperoxides of 2,3-unsaturated hexapyranoses, and hydroperoxides of OMe and OBn substituted derivatives of 2-deoxy-glucopyranose and 2-deoxy-galactopyranose to investigate the anomeric and exo-anomeric effects of these groups. The structure and energy of the conformers were calculated at the B3LYP/6-311++G** level. Calculations showed that the peroxy anion group exhibits a strong anomeric effect, comparable in magnitude to the methoxy group, and that the anomeric effect of the hydroperoxy group is similar to the hydroxyl group. These results revealed that hydroperoxy and peroxy anion groups display an exo-anomeric effect, but the orientation around the C1-O1 bond is also affected by hydrogen bonding and electrostatic interactions.     

Metal-activated histidine carbon donor hydrogen bonds contribute to metalloprotein folding and function

Carbon donor hydrogen bonds are typically weak interactions that contribute less than 2 kcal/mol, and provide only modest stabilization in proteins. One exception is the class of hydrogen bonds donated by heterocyclic side chain carbons. Histidine is capable of particularly strong interactions through the Cε¹ and Cδ² carbons when the imidazole is protonated or bound to metal. Given the frequent occurrence of metal-bound histidines in metalloproteins, we characterized the energies of these interactions through DFT calculations on model compounds. Imidazole-water hydrogen bonding could vary from −11.0 to −17.0 kcal/mol, depending on the metal identity and oxidation state. A geometric search of metalloprotein structures in the PDB identified a number of candidate His C-H···O hydrogen bonds which may be important for folding or function. DFT calculations on model complexes of superoxide reductase show a carbon donor hydrogen bond positioning a water molecule above the active site.     

Steric and electrostatic effects in DNA synthesis by the SOS-induced DNA polymerases II and IV of Escherichia coli

The SOS-induced DNA polymerases II and IV (pol II and pol IV, respectively) of Escherichia coli play important roles in processing lesions that occur in genomic DNA. Here we study how electrostatic and steric effects play different roles in influencing the efficiency and fidelity of DNA synthesis by these two enzymes. These effects were probed by the use of nonpolar shape analogues of thymidine, in which substituted toluenes replace the polar thymine base. We compared thymine with nonpolar analogues to evaluate the importance of hydrogen bonding in the polymerase active sites, while we used comparisons among a set of variably sized thymine analogues to measure the role of steric effects in the two enzymes. Steady-state kinetics measurements were carried out to evaluate activities for nucleotide insertion and extension. The results showed that both enzymes inserted nucleotides opposite nonpolar template bases with moderate to low efficiency, suggesting that both polymerases benefit from hydrogen bonding or other electrostatic effects involving the template base. Surprisingly, however, pol II inserted nonpolar nucleotide (dNTP) analogues into a primer strand with high (wild-type) efficiency, while pol IV handled them with an extremely low efficiency. Base pair extension studies showed that both enzymes bypass non-hydrogen-bonding template bases with moderately low efficiency, suggesting a possible beneficial role of minor groove hydrogen bonding interactions at the N-1 position. Measurement of the two polymerases' sensitivity to steric size changes showed that both enzymes were relatively flexible, yielding only small kinetic differences with increases or decreases in nucleotide size. Comparisons are made to recent data for DNA pol I (Klenow fragment), the archaeal polymerase Dpo4, and human pol kappa.     

Hydration of amino acid side chains: dependence on secondary structure.

Energy calculations have been used to study the hydration sites around the polar groups of serine, threonine and tyrosine side chains. These hydration sites depend not only on the hybridization of the polar group but also on the local secondary structure, the chi 1 side chain torsion angle and the position of the hydroxyl hydrogen atom. For tyrosine side chains, two solvent sites are found approximately in the plane of the ring. Even for serine and threonine side chains only two minimum energy sites are found in general of which one is in an expected position within hydrogen bonding of the hydroxyl hydrogen atom (unless this is blocked from interaction with solvent molecules by, for example, Oi-4 or Oi-3. The position of the second of these sites depends not only on the position of the hydroxyl oxygen but also on neighbouring main chain atoms to which it can also hydrogen bond. There is good agreement with the solvent distributions obtained from crystallographic data.     

Improved side-chain prediction accuracy using an ab initio potential energy function and a very large rotamer library

Accurate prediction of the placement and comformations of protein side chains given only the backbone trace has a wide range of uses in protein design, structure prediction, and functional analysis. Prediction has most often relied on discrete rotamer libraries so that rapid fitness of side-chain rotamers can be assessed against some scoring function. Scoring functions are generally based on experimental parameters from small-molecule studies or empirical parameters based on determined protein structures. Here, we describe the NCN algorithm for predicting the placement of side chains. A predominantly first-principles approach was taken to develop the potential energy function incorporating van der Waals and electrostatics based on the OPLS parameters, and a hydrogen bonding term. The only empirical knowledge used is the frequency of rotameric states from the PDB. The rotamer library includes nearly 50,000 rotamers, and is the most extensive discrete library used to date. Although the computational time tends to be longer than most other algorithms, the overall accuracy exceeds all algorithms in the literature when placing rotamers on an accurate backbone trace. Considering only the most buried residues, 80% of the total residues tested, the placement accuracy reaches 92% for chi(1), and 83% for chi(1 + 2), and an overall RMS deviation of 1 A. Additionally, we show that if information is available to restrict chi(1) to one rotamer well, then this algorithm can generate structures with an average RMS deviation of 1.0 A for all heavy side-chains atoms and a corresponding overall chi(1 + 2) accuracy of 85.0%.     

The effects of acetate anion on cellulose dissolution and reaction in imidazolium ionic liquids

Quantum mechanical calculations were carried out to determine the mechanisms for the superiority of the imidazolium acetate-based ionic liquids to the corresponding chloride-based ionic liquids. Our results indicate that the imidazolium cation can react with the acetate anion to generate a carbene, a highly reactive intermediate. The carbene produced then reacts with cellulose to facilitate its dissolution in the ionic liquid solvents in addition to the stronger hydrogen bonds formed between the acetate anion and the hydroxyl groups on cellulose. The mechanisms for the imidazolium cation and acetate anion reactions involve the initial ion pairing of the cation and anion via hydrogen bonding and electrostatic interactions. The hydrogen bond formed between the C2–H on the imidazolium cation and COO⁻ of the anion facilitates the transfer of the H⁺ to the anion to form a carbene intermediate.     

Rational design of carbonitrile-carboxaldehyde cation receptor models: probing the nature of the heteroatom–metal interaction

Amino acids are constituents of proteins and enzymes which take part almost in all metabolic reactions. Glutamic acid, with an ability to form a negatively charged side chain, plays a major role in intra and intermolecular interactions of proteins, peptides, and enzymes. An exhaustive conformational analysis has been performed for all eight possible forms at B3LYP/cc-pVTZ level. All possible neutral, zwitterionic, protonated, and deprotonated forms of glutamic acid structures have been investigated in solution by using polarizable continuum model mimicking water as the solvent. Nine families based on the dihedral angles have been classified for eight glutamic acid forms. The electrostatic effects included in the solvent model usually stabilize the charged forms more. However, the stability of the zwitterionic form has been underestimated due to the lack of hydrogen bonding between the solute and solvent; therefore, it is observed that compact neutral glutamic acid structures are more stable in solution than they are in vacuum. Our calculations have shown that among all eight possible forms, some are not stable in solution and are immediately converted to other more stable forms. Comparison of isoelectronic glutamic acid forms indicated that one of the structures among possible zwitterionic and anionic forms may dominate over the other possible forms. Additional investigations using explicit solvent models are necessary to determine the stability of charged forms of glutamic acid in solution as our results clearly indicate that hydrogen bonding and its type have a major role in the structure and energy of conformers.     

Origin of sequence-specific recognition of DNA by non-intercalating anti-tumor antibiotics

Partitioning of energy in the interaction of non-intercalating antibiotics (netropsin, netropsin without its cationic ends and two analogs of distamycin A) with different base sequences of B-DNA is studied here by the atom-atom potential technique and geometry optimization procedures. The results show that electrostatic forces contribute substantially to the stabilization energy as well as to the sequence specificity. The hydrogen-bonding term is also sequence specific and is significant in properly orienting the drug molecule. Relative roles of the hydrogen bonding and electrostatic interactions depend on the dielectric property of the medium.     

Common Hydrogen Bond Interactions in Diverse Phosphoryl Transfer Active Sites

Phosphoryl transfer reactions figure prominently in energy metabolism, signaling, transport and motility. Prior detailed studies of selected systems have highlighted mechanistic features that distinguish different phosphoryl transfer enzymes. Here, a top-down approach is developed for comparing statistically the active site configurations between populations of diverse structures in the Protein Data Bank, and it reveals patterns of hydrogen bonding that transcend enzyme families. Through analysis of large samples of structures, insights are drawn at a level of detail exceeding the experimental precision of an individual structure. In phosphagen kinases, for example, hydrogen bonds with the O_3β of the nucleotide substrate are revealed as analogous to those in unrelated G proteins. In G proteins and other enzymes, interactions with O_3β have been understood in terms of electrostatic favoring of the transition state. Ground state quantum mechanical calculations on model compounds show that the active site interactions highlighted in our database analysis can affect substrate phosphate charge and bond length, in ways that are consistent with prior experimental observations, by modulating hyperconjugative orbital interactions that weaken the scissile bond. Testing experimentally the inference about the importance of O_3β interactions in phosphagen kinases, mutation of arginine kinase Arg₂₈₀ decreases k_cat, as predicted, with little impact upon K_M.     

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).     

Probing multiple effects on 15N, 13C alpha, 13C beta,and 13C' chemical shifts in peptides using density functional theory

We have used density functional calculations on model peptides to study conformational effects on (15)N, (13)C alpha, (13)C beta, and (13)C' chemical shifts, associated with hydrogen bonding, backbone conformation, and side-chain orientation. The results show a significant dependence on the backbone torsion angles of the nearest three residues. Contributions to (15)N chemical shifts from hydrogen bonding (up to 8 ppm), backbone conformation (up to 13 ppm), side-chain orientation and neighborhood residue effects (up to 22 ppm) are significant, and a unified theory will be required to account for their behavior in proteins. In contrast to this, the dependence on sequence and hydrogen bonding is much less for (13)C alpha and (13)C beta chemical shifts (<0.5 ppm), and moderate for carbonyl carbon shifts (<2 ppm). The effects of side-chain orientation are mainly limited to the residue itself for both nitrogen and carbon, but the chi(1) effect is also significant for the nitrogen shift of the following residue and for the (13)C' shift of the preceding residue. The calculated results are used, in conjunction with an additive model of chemical shift contributions, to create an algorithm for prediction of (15)N and (13)C shifts in proteins from their structure; this includes a model to extrapolate results to regions of torsion angle space that have not been explicitly studied by density functional theory (DFT) calculations. Crystal structures of 20 proteins with measured shifts have been used to test the prediction scheme. Root mean square deviations between calculated and experimental shifts 2.71, 1.22, 1.31, and 1.28 ppm for N, C alpha, C beta, and C', respectively. This prediction algorithm should be helpful in NMR assignment, crystal and solution structure comparison, and structure refinement.     

Dynamics of Cleft Closure of the GluA2 Ligand-binding Domain in the Presence of Full and Partial Agonists Revealed by Hydrogen-Deuterium Exchange

The majority of excitatory neurotransmission in the CNS is mediated by tetrameric AMPA receptors. Channel activation begins with a series of interactions with an agonist that binds to the cleft between the two lobes of the ligand-binding domain of each subunit. Binding leads to a series of conformational transitions, including the closure of the two lobes of the binding domain around the ligand, culminating in ion channel opening. Although a great deal has been learned from crystal structures, determining the molecular details of channel activation, deactivation, and desensitization requires measures of dynamics and stabilities of hydrogen bonds that stabilize cleft closure. The use of hydrogen-deuterium exchange at low pH provides a measure of the variation of stability of specific hydrogen bonds among agonists of different efficacy. Here, we used NMR measurements of hydrogen-deuterium exchange to determine the stability of hydrogen bonds in the GluA2 (AMPA receptor) ligand-binding domain in the presence of several full and partial agonists. The results suggest that the stabilization of hydrogen bonds between the two lobes of the binding domain is weaker for partial than for full agonists, and efficacy is correlated with the stability of these hydrogen bonds. The closure of the lobes around the agonists leads to a destabilization of the hydrogen bonding in another portion of the lobe interface, and removing an electrostatic interaction in Lobe 2 can relieve the strain. These results provide new details of transitions in the binding domain that are associated with channel activation and desensitization.     

Binding of transition metal complexes to guanine and guanine–cytosine: hydrogen bonding and covalent effects

Density functional calculations and Atoms in Molecules analysis are used to investigate the role of covalent and hydrogen bondings in determining the binding of transition metal complexes to guanine, and the subsequent effect on pairing with cytosine. Hydrogen bonding is ubiquitous, and typically contributes ca. 10% to overall binding, a value that varies with the coordination site on guanine, as well as metal and ligands. Early transition metals show a clear preference for the O6 position, while later ones prefer N7, the crossover point coming at the vanadium group. Metallation at N7 causes a redistribution of hydrogen bonding strength between guanine and cytosine, but does not greatly affect the overall pairing energy. In contrast, metallation at O6 strongly reduces the pairing energy, as may be expected given the role of O6 in pairing guanine with cytosine. This effect can be quantified using electron density properties, and seems to be due to both electrostatic repulsion from the positive metal centre and a redistribution of electron density within guanine itself. Qualitative agreement with experimental mass spectroscopic results is obtained.     

Effect of solvent environment on the CO band position in the infrared spectrum of trans-[Fe(CN)4(CO)2]

Density functional theory (DFT) is used to understand the effect of hydrogen bonding solvents on the CO band position in the infrared (IR) spectrum of a mono-iron complex, trans-[Fe^II(CN)₄(CO)₂]²⁻. This mono-iron complex has received much attention recently due its potential relation to the biosynthesis of Fe-only hydrogenase enzymes. Our calculations show that the polar solvent molecules preferentially hydrogen bond to the cyano ligands in this complex. The effect of such hydrogen bonding on the electron density distribution is analyzed in terms of the population in natural bond orbitals (NBO). Our results show that the presence of hydrogen bonding to the cyano ligands decreases the extent of back bonding from the metal to the carbonyl ligand. This results in decreased electron density in the π^∗ orbitals of the carbonyl bond leading to a strengthening of the CO bond and a consequent blue shift in the IR band position of the carbonyl group. We also show that the extent of blue shift correlates with the number of nearest neighbor solvent molecules.     

Structure-based virtual screening of essential Mycobacterium tuberculosis enzymes AspS and KatG for potential inhibitors

Mycobacterium tuberculosis (TB) is a leading global cause of disease-related death. Recent works have studied metabolic pathways of the mycobacterium, highlighting essential enzymes to target via competitive inhibition through computational molecular modeling to suppress the organism''s life cycle. We used the Protein Databank (PDB), the UniProt Knowledgebase and the iDock server in this study. In vitro toxicity screening and pharmacokinetic properties were assessed to determine potential ligand safety and drug properties. Our results have revealed five and nine potential ligands for the enzymes AspS and KatG respectively. The KatG active site has displayed binding affinities of -13.443 to -12.895 kcal/mol, while AspS ligands range from -6.580 to -6.490kcal/mol. The intermolecular forces responsible for the differing binding affinities of each enzyme are primarily Coulombic interactions for AspS, versus Coulombic and extensive hydrogen bonding interactions in KatG.