Hydrogen bonding in complex of serine with histidine: computational and spectroscopic study of model compounds

Experimental and DFT dimer modeling studies of the H-bond induced-vibration modes of l-β-Homoserine

The vibrational spectra for l-β-Homoserine have been measured (IR absorption: 4000-400cm^-1/Raman spectra: 4000-200cm^-1). Characteristic vibrational modes of ammonium (-NH₃⁺), carboxylate (-CO₂^-) and hydroxyl (-OH) groups across the 3700-1400cm^-1 are all identified to have originated in inter-molecular hydrogen bonding involving these functional groups. DFT calculations at B3LYP/6-311++G(d, p) level have yielded a single neutral monomer in the gas phase. Since as a member of the amino acids which are known to possess zwitterionic structure in condensed phase, the neutral monomer of l-β-Homoserine is optimized to a zwitterionic structure in a water medium. Consideration of two dimer structures, one dimer with -NH‧‧‧O bond and another -OH‧‧‧O bond, has given rise to vibrational modes that satisfactorily fit to all the observed absorption and Raman bands. It is found that the dimer with -OH‧‧‧O bond (binding energy, 8.896kcal/mol) is more tightly bound than the dimer with -NH‧‧‧O bond (8.363kcal/mol).

Alcohol dimers – how much diagonal OH anharmonicity?

The massive infrared attenuation and increased anharmonicity of OH stretching overtones in alcohol dimers is experimentally quantified.

Biochemical and biophysical analysis of a chiral PqsD inhibitor revealing tight-binding behavior and enantiomers with contrary thermodynamic signatures

Antivirulence strategies addressing bacterial pathogenicity without exhibiting growth inhibition effects represent a novel approach to overcome today's crisis in antibiotic development. In recent studies, we examined various inhibitors of PqsD, an enzyme involved in formation of Pseudomonas aeruginosa cell-to-cell signaling molecules, and observed desired cellular effects for 2-nitrophenyl derivatives. Herein, we investigated the binding characteristics of this interesting compound class using several biochemical and biophysical methods. The inhibitors showed time-dependent activity, tight-binding behavior, and interactions with the catalytic center. Furthermore, isothermal titration calorimetry (ITC) experiments with separated enantiomers revealed contrary thermodynamic signatures showing either enthalpy- or entropy-driven affinity. A combination of site-directed mutagenesis and thermodynamic profiling was used to identify key residues involved in inhibitor binding. This information allowed the proposal of experimentally confirmed docking poses. Although originally designed as transition state analogs, our results suggest an altered position for both enantiomers. Interestingly, the main difference between stereoisomers was found in the orientation of the hydroxyl group at the stereogenic center. The predicted binding modes are in accordance with experimental data and, thus, allow future structure-guided optimization.

Computational chemical analysis of unconjugated bilirubin anions and insights into pKa values clarification

The pKa, the negative logarithm of the acid dissociation equilibrium constant, of the carboxylic acid groups of unconjugated bilirubin in water is a discussed issue because there are quite different experimental values reported. Using quantum mechanical calculations we have studied the conformational behavior of unconjugated bilirubin species (in gas phase and in solution modeled implicitly and explicitly) to provide evidence that may clarify pKa values because of its pathophysiological relevance. Our results show that rotation of carboxylate group, which is not restricted, settles it in a suitable place to establish stronger interactions that stabilizes the monoanion and the dianion to be properly solvated, demonstrating that the rationalization used to justify the high pKa values of unconjugated bilirubin is inappropriate. Furthermore, low unconjugated bilirubin (UCB) pKa values were estimated from a linear regression analysis.

The structure of geminal imidazolium bis(trifluoromethylsulfonyl)amide ionic liquids: a theoretical study of the gas phase ionic complexes

In this work we report molecular mechanics and ab initio calculations on the geminal di-imidazolium bis(trifluoromethylsulfonyl)amide ionic liquid in the gas phase. We report the likely energetically preferred geometries of the ionic complex and its main features in terms of charge distribution, electronic density, structure, and energetics. We find that the gas phase structure of the ionic complex is quite compact and that the alkyl chain connecting the two imidazolium charged rings is strongly bent in order to maximize their electrostatic interactions with the two anions.

Is Asp-His-Ser/Thr-Trp tetrad hydrogen-bond network important to WD40-repeat proteins: a statistical and theoretical study

WD40-repeat proteins are abundant and play important roles in forming protein complexes. The domain usually has seven WD40 repeats, which folds into a seven beta-sheet propeller with each beta-sheet in a four-strand structure. An analysis of 20 available WD40-repeat proteins in Protein Data Bank reveals that each protein has at least one Asp-His-Ser/Thr-Trp (D-H-S/T-W) hydrogen-bonded tetrad, and some proteins have up to six or seven such tetrads. The relative positions of the four residues in the tetrads are also found to be conserved. A sequence alignment analysis of 560 WD40-repeat protein sequences in human reveals very similar features, indicating that such tetrad may be a general feature of WD40-repeat proteins. We carried out density functional theory and found that these tetrads can lead to significant stabilization including hydrogen-bonding cooperativity. The hydrogen bond involving Trp is significant. These results lead us to propose that the tetrads may be critical to the stability and the mechanism of folding of these proteins.

Combined QM/MM mechanistic study of the acylation process in furin complexed with the H5N1 avian influenza virus hemagglutinin's cleavage site

Combined quantum mechanical/molecular mechanical (QM/MM) techniques have been applied to investigate the detailed reaction mechanism of the first step of the acylation process by furin in which the cleavage site of the highly pathogenic avian influenza virus subtype H5N1 (HPH5) acts as its substrate. The energy profile shows a simultaneous mechanism, known as a concerted reaction, of the two subprocesses: the proton transfer from Ser368 to His194 and the nucleophilic attack on the carbonyl carbon of the scissile peptide of the HPH5 cleavage site with a formation of tetrahedral intermediate (INT). The calculated energy barrier for this reaction is 16.2 kcal.mol(-1) at QM/MM B3LYP/6-31+G*//PM3-CHARMM22 level of theory. Once the reaction proceeds, the ordering of the electrostatic stabilization by protein environment is of the enzyme-substrate < transition state < INT complexes. Asp153 was found to play the most important role in the enzymatic reaction by providing the highest degree of intermediate complex stabilization. In addition, the negatively charged carbonyl oxygen of INT is well stabilized by the oxyanion hole constructed by Asn295's carboxamide and Ser368's backbone.

Gas-Phase Structure of Protonated Histidine and Histidine Methyl Ester: Combined Experimental Mass Spectrometry and Theoretical ab Initio Study

Gas-phase H/D exchange experiments with CD3OD and D2O and quantum chemical ab initio G3(MP2) calculations were carried out on protonated histidine and protonated histidine methyl ester in order to elucidate their bonding and structure. The H/D exchange experiments show that both ions have three equivalent fast hydrogens and one appreciably slower exchangeable hydrogen assigned to the protonated amino group participating in a strong intramolecular hydrogen bond (IHB) with the nearest N(sp2) nitrogen of the imidazole fragment and to the distal ring NH-group, respectively. It is taken for granted that the proton exchange in the IHB is much faster than the H/D exchange. Unlike in other protonated amino acids (glycine, proline, phenylalanine, tyrosine, and tryptophan) studied earlier, the exchange rate of the carboxyl group in protonated histidine is slower than that of the amino group. The most stable conformers and the enthalpies of neutral and protonated histidine and its methyl ester are calculated at the G3(MP2) level of theory. It is shown that strong intramolecular hydrogen bonding between the amino group and the imidazole ring nitrogen sites is responsible for the stability and specific properties of the protonated histidine. It is found that the proton fluctuates between the amino and imidazole groups in the protonated form across an almost vanishing barrier. Proton affinity (PA) of histidine calculated by the G3(MP2) method is 233.2 and 232.4 kcal mol(-1) for protonation at the imidazole ring and at the amino group nitrogens, respectively, which is about 3-5 kcal mol(-1) lower than the reported experimental value.