Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

NMR evidence for the participation of a low-barrier hydrogen bond in the mechanism of delta 5-3-ketosteroid isomerase.

Delta 5-3-Ketosteroid isomerase (EC 5.3.3.1) promotes an allylic rearrangement involving intramolecular proton transfer via a dienolic intermediate. This enzyme enhances the catalytic rate by a factor of 10(10). Two residues, Tyr-14, the general acid that polarizes the steroid 3-carbonyl group and facilitates enolization, and Asp-38 the general base that abstracts and transfers the 4 beta-proton to the 6 beta-position, contribute 10(4.7) and 10(5.6) to the rate increase, respectively. A major mechanistic enigma is the huge disparity between the pKa values of the catalytic groups and their targets. Upon binding of an analog of the dienolate intermediate to isomerase, proton NMR detects a highly deshielded resonance at 18.15 ppm in proximity to aromatic protons, and with a 3-fold preference for protium over deuterium (fractionation factor, phi = 0.34), consistent with formation of a short, strong (low-barrier) hydrogen bond to Tyr-14. The strength of this hydrogen bond is estimated to be at least 7.1 kcal/mol. This bond is relatively inaccessible to bulk solvent and is pH insensitive. Low-barrier hydrogen bonding of Tyr-14 to the intermediate, in conjunction with the previously demonstrated tunneling contribution to the proton transfer by Asp-38, provide a plausible and quantitative explanation for the high catalytic power of this isomerase.

Kinetics of hydrogen bond breakage in the process of unfolding of ribonuclease A measured by pulsed hydrogen exchange.

A sensitive test for kinetic unfolding intermediates in ribonuclease A (EC 3.1.27.5) is performed under conditions where the enzyme unfolds slowly (10 degrees C, pH 8.0, 4.5 M guanidinium chloride). Exchange of peptide NH protons (2H-1H) is used to monitor structural opening of individual hydrogen bonds during unfolding, and kinetic models are developed for hydrogen exchange during the process of protein unfolding. The analysis indicates that the kinetic process of unfolding can be monitored by EX1 exchange (limited by the rate of opening) for ribonuclease A in these conditions. Of the 49 protons whose unfolding/exchange kinetics was measured, 47 have known hydrogen bond acceptor groups. To test whether exchange during unfolding follows the EX2 (base-catalyzed) or the EX1 (uncatalyzed) mechanism, unfolding/exchange was measured both at pH 8.0 and at pH 9.0. A few faster-exchanging protons were found that undergo exchange by both EX1 and EX2 processes, but the 43 slower-exchanging protons at pH 8 undergo exchange only by the EX1 mechanism, and they have closely similar rates. Thus, it is likely that all 49 protons undergo EX1 exchange at the same rate. The results indicate that a single rate-limiting step in unfolding breaks the entire network of peptide hydrogen bonds and causes the overall unfolding of ribonuclease A. The additional exchange observed for some protons that follows the EX2 mechanism probably results from equilibrium unfolding intermediates and will be discussed elsewhere.

2-(3-Hydroxy-propyl)isoindoline-1,3-dione:Competition among hydrogen-bond acceptors

The title compound, C11H11NO3, has two mol-ecules in the asymmetric unit, which differ in the orientation of their side-chain OH groups, allowing them to form inter-molecular O - H⋯O hydrogen bonds to different acceptors. In one case, the acceptor is the OH group of the other mol-ecule, and in the other case it is an imide O=C group. This is the first example in the N-substituted phthalimide series in which independent mol-ecules have different types of acceptor. Mol-ecular-orbital calculations place the greatest negative charge on the OH group. © 2008 International Union of Crystallography.

Experimental and theoretical studies on six-membered ring hydrogen bond systems in crystalline drug substances

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

Observed and predicted hydrogen bond motifs in crystal structures of hydantoins, dihydrouracils and uracils

A survey of crystal structures containing hydantoin, dihydrouracil and uracil derivatives in the Cambridge Structural Database revealed four main types of hydrogen bond motifs when derivatives with extra substituents able to interfere with the main motif are excluded. All these molecules contain two hydrogen bond donors and two hydrogen bond acceptors in the sequence of NH, C = O, NH, and C=O groups within a 5-membered ring (hydantoin) and two 6-membered rings (dihydrouracil and uracil). In all cases, both ring NH groups act as donors in the main hydrogen bond motif but there is an excess of hydrogen bond acceptors (two C=O able to accept twice each) and so two possibilities are found: (i) each carbonyl O atom may accept one hydrogen bond or (ii) one carbonyl O atom may accept two hydrogen bonds while the other does not participate in the hydrogen bonding. We observed different preferences in the type and symmetry of the motifs adopted by the different derivatives, and a good agreement is found between motifs observed experimentally and those predicted using computational methods. We identified certain molecular factors such as chirality, substituent size and the possibility of C-H⋯O interactions as important factors influencing the motif observation. © 2012 The Royal Society of Chemistry and the Centre National de la Recherche Scientifique.

Stereochemical studies on cyclic peptides. XIII. Energy minimization studies on cyclic hexapeptides having hydrogen bonds

The basic cyclic hexapeptide conformations which accommodate hydrogen bonded β and γ turns in the backbone have been worked out using stereochemical criteria and energy minimization procedures. It was found that cyclic hexapeptides can be made up of all possible combinations of 4 ± 1 hydrogen bonded types I, I', II and II' β turns, giving rise to symmetric conformations having twofold and inversion symmetries as well as nonsymmetric structures. Conformations having exclusive features of 3 ± 1 hydrogen bonded γ turns were found to be possible in threefold and S6 symmetric cyclic hexapeptides. The results show that the cyclic hexapeptides formed by the linking of two β turn tripeptide fragments differ mainly in (a) the hydrogen bonding scheme present in the β turn tripeptides and (b) the conformation at the α-carbon atoms where the two tripeptide fragments link. The different hydrogen bonding schemes found in the component β turns are: 1) a β turn with only a 4 ± 1 hydrogen bond, 2) a type I or I' β turn with 4 ± 1 and 3 ± 1 hydrogen bonds occurring in a bifurcated form and 3) a type II or II' β turn having both the 4 ± 1 and the 3 ± 1 hydrogen bonds with the same acceptor oxygen atom. The conformation at the linking α-carbon atoms was found to lie either in the extended region or in the 3 ± 1 hydrogen bonded γ turn or inverse γ turn regions. Further, the threefold and the S6 symmetric conformations have three γ turns interleaved by three extended regions or three inverse γ turns, respectively. The feasibility of accommodating alanyl residues of both isomeric forms in the CHP minima has been explored. Finally, the available experimental data are reviewed in the light of the present results.

Structure of the hydrogen-bonded water molecule in crystals

A correlation of the structural data on IS hydrates obtained by x-ray diffraction, neutron diffraction, and proton magnetic resonance reveals that when a water molecule is hydrogen bonded into a crystal structure and the angle subtended at the donor water oxygen by the acceptor atoms deviates from the vapor H-O-H angle, bent hydrogen bonds are formed in preference to distortion of the H-O-H angle. Theoretical justification for this result is obtained from energy considerations by calculating the energy of formation of bent hydrogen bonds on the basis of the Lippincott-Schroeder potential function model for the hydrogen bond and the energy of deformation of the H-O-H angle from spectroscopic force constants.

The Mechanism Of The Hydrogen Abstraction Process By Photo-Excited Ketones

Qualitative potential energy surfaces for hydrogen abstraction from alkanes containing primary, secondary and tertiary C-H bonds by a photo-excited ketone have been reported, The results suggest that the activation barriers for these processes decrease in the order primary > secondary > tertiary in agreement with the observed trend in the rate constants. The analysis of the electronic structure of the transition-state reveal that electron-transfer from hydrocarbon to ketone and formation of a new bond are almost synchronous in the hydrogen transfer process. The tunneling of hydrogen is not important in the normal temperature region even though the barriers are small.

The Nature of Bond Critical Points in Dinuclear Copper(I) Complexes

Closed-shell contacts between two copper(I) ions are expected to be repulsive. However, such contacts are quite frequent and are well documented. Crystallographic characterization of such contacts in unsupported and bridged multinuclear copper(I) complexes has repeatedly invited debates on the existence of cuprophilicity. Recent developments in the application of Baders theory of atoms-in-molecules (AIM) to systems in which weak hydrogen bonds are involved suggests that the copper(I)copper(I) contacts would benefit from a similar analysis. Thus the nature of electron-density distributions in copper(I) dimers that are unsupported, and those that are bridged, have been examined. A comparison of complexes that are dimers of symmetrical monomers and those that are dimers of two copper(I) monomers with different coordination spheres has also been made. AIM analysis shows that a bond critical point (BCP) between two Cu atoms is present in most cases. The nature of the BCP in terms of the electron density, ?, and its Laplacian is quite similar to the nature of critical points observed in hydrogen bonds in the same systems. The ? is inversely correlated to Cu?Cu distance. It is higher in asymmetrical systems than what is observed in corresponding symmetrical systems. By examining the ratio of the local electron potential-energy density (Vc) to the kinetic energy density (Gc), |Vc|/Gc at the critical point suggests that these interactions are not perfectly ionic but have some shared nature. Thus an analysis of critical points by using AIM theory points to the presence of an attractive metallophilic interaction similar to other well-documented weak interactions like hydrogen bonding.

Intermolecular charge transfer and vibrational analysis of hydrogen bonding in acetazolamide

In the present work the structural and spectral characteristics of acetazolamide have been studied by methods of infrared, Raman spectroscopy and quantum chemistry. Electrostatic potential surface, optimized geometry, harmonic vibrational frequencies, infrared intensities and activities of Raman scattering were calculated by density functional theory (DFT) employing B3LYP with complete relaxation in the potential energy surface using 6-311++G(d,p) basis set. Based on these results, we have discussed the correlation between the vibrational modes and the structure of the dimers of acetazolamide. The calculated vibrational spectra of three dimers of acetazolamide have been compared with observed spectra, and the assignment of observed bands was carried out using potential energy distribution. The observed spectra agree well with the values computed from the OFT. A comparison of observed and calculated vibrational spectra clearly shows the effect of hydrogen bonding. The frequency shifts observed for the different dimers are in accord with the hydrogen bonding in acetazolamide. Natural bond orbital (NBO) analyses reflect the charge transfer interaction in the individual hydrogen bond units and the stability of different dimers of acetazolamide. (C) 2012 Elsevier B.V. All rights reserved.

Fe as Hydrogen/Halogen Bond Acceptor in Square Pyramidal Fe(CO)(5)

Hydrogen bonded complexes formed between the square pyramidal Fe(CO)(5) with HX (X = F, Cl, Br), showing X-H center dot center dot center dot Fe interactions, have been investigated theoretically using density functional theory (DFT) including dispersion correction. Geometry, interaction energy, and large red shift of about 400 cm(-1) in the FIX stretching frequency confirm X-H center dot center dot center dot Fe hydrogen bond formation. In the (CO)(5)Fe center dot center dot center dot HBr complex, following the significant red shift, the HBr stretching mode is coupled with the carbonyl stretching modes. This clearly affects the correlation between frequency shift and binding energy, which is a hallmark of hydrogen bonds. Atoms in Molecule (AIM) theoretical analyses show the presence of a bond critical point between the iron and the hydrogen of FIX and significant mutual penetration. These X-H center dot center dot center dot Fe hydrogen bonds follow most but not all of the eight criteria proposed by Koch and Popelier (J. Phys. Chem. 1995, 99, 9747) based on their investigations on C-H center dot center dot center dot O hydrogen bonds. Natural bond orbital (NBO) analysis indicates charge transfer from the organometallic system to the hydrogen bond donor. However, there is no correlation between the extent of charge transfer and interaction,energy, contrary to what is proposed in the recent IUPAC recommendation (Pure Appl.. Chem. 2011, 83, 1637). The ``hydrogen bond radius'' for iron has been determined to be 1.60 +/- 0.02 angstrom, and not surprisingly it is between the covalent (127 angstrom) and van der Waals (2.0) radii of Fe. DFT and AIM theoretical studies reveal that Fe in square pyramidal Fe(CO)(5) can also form halogen bond with CIF and ClH as ``halogen bond donor''. Both these complexes show mutual penetration as well, though the Fe center dot center dot center dot Cl distance is closer to the sum of van der Waals radii of Fe and Cl in (CO)5Fe center dot center dot center dot ClH, and it is about 1 angstrom less in (CO)(5)Fe center dot center dot center dot ClF.

Quantum chemical study on influence of intermolecular hydrogen bonding on the geometry, the atomic charges and the vibrational dynamics of 2,6-dichlorobenzonitrile

FT-IR (4000-400 cm(-1)) and FT-Raman (4000-200 cm(-1)) spectral measurements on solid 2,6-dichlorobenzonitrile (2,6-DCBN) have been done. The molecular geometry, harmonic vibrational frequencies and bonding features in the ground state have been calculated by density functional theory at the B3LYP/6-311++G (d,p) level. A comparison between the calculated and the experimental results covering the molecular structure has been made. The assignments of the fundamental vibrational modes have been done on the basis of the potential energy distribution (PED). To investigate the influence of intermolecular hydrogen bonding on the geometry, the charge distribution and the vibrational spectrum of 2,6-DCBN; calculations have been done for the monomer as well as the tetramer. The intermolecular interaction energies corrected for basis set superposition error (BSSE) have been calculated using counterpoise method. Based on these results, the correlations between the vibrational modes and the structure of the tetramer have been discussed. Molecular electrostatic potential (MEP) contour map has been plotted in order to predict how different geometries could interact. The Natural Bond Orbital (NBO) analysis has been done for the chemical interpretation of hyperconjugative interactions and electron density transfer between occupied (bonding or lone pair) orbitals to unoccupied (antibonding or Rydberg) orbitals. UV spectrum was measured in methanol solution. The energies and oscillator strengths were calculated by Time Dependent Density Functional Theory (TD-DFT) and matched to the experimental findings. TD-DFT method has also been used for theoretically studying the hydrogen bonding dynamics by monitoring the spectral shifts of some characteristic vibrational modes involved in the formation of hydrogen bonds in the ground and the first excited state. The C-13 nuclear magnetic resonance (NMR) chemical shifts of the molecule were calculated by the Gauge independent atomic orbital (GIAO) method and compared with experimental results. Standard thermodynamic functions have been obtained and changes in thermodynamic properties on going from monomer to tetramer have been presented. (C) 2013 Elsevier B.V. All rights reserved.

Microwave spectroscopic and theoretical investigations of the strongly hydrogen bonded hexafluoroisopropanol...water complex

This paper reports microwave spectroscopic and theoretical investigations on the interaction of water with hexafluoroisopropanol (HFIP). The HFIP monomer can exist in two conformations, antiperiplanar (AP) and synclinical (SC). The former is about 5 kJ mol(-1) more stable than the latter. Theoretical calculations predicted three potential minima for the complex, two having AP and one having SC conformations. Though, the binding energy for the HFIP(SC)...H2O turned out to be larger than that for the other two conformers having HFIP in the AP form, the global minimum for the complex in the potential energy hypersurface had HFIP in the AP form. Experimental rotational constants for four isotopologues measured using a pulsed nozzle Fourier transform microwave spectrometer, correspond to the global minimum in the potential energy hypersurface. The structural parameters and the internal dynamics of the complex could be determined from the rotational spectra of the four isotopologues. The global minimum has the HFIP(AP) as a hydrogen bond donor forming a strong hydrogen bond with H2O. To characterize the strength of the bonding and to probe the other interactions within the complex, atoms in molecules, non-covalent interaction index and natural bond orbital theoretical analyses have been performed.

Hydrogen energy future with formic acid: a renewable chemical hydrogen storage system

Formic acid, the simplest carboxylic acid, is found in nature or can be easily synthesized in the laboratory (major by-product of some second generation biorefinery processes); it is also an important chemical due to its myriad applications in pharmaceuticals and industry. In recent years, formic acid has been used as an important fuel either without reformation (in direct formic acid fuel cells, DFAFCs) or with reformation (as a potential chemical hydrogen storage material). Owing to the better efficiency of DFAFCs compared to several other PEMFCs and reversible hydrogen storage systems, formic acid could serve as one of the better fuels for portable devices, vehicles and other energy-related applications in the future. This perspective is focused on recent developments in the use of formic acid as a reversible source for hydrogen storage. Recent developments in this direction will likely give access to a variety of low-cost and highly efficient rechargeable hydrogen fuel cells within the next few years by the use of suitable homogeneous metal complex/heterogeneous metal nanoparticle-based catalysts under ambient reaction conditions. The production of formic acid from atmospheric CO2 (a greenhouse gas) will decrease the CO2 content and may be helpful in reducing global warming.