Contribution of intramolecular C=O...H-N hydrogen bonding to the solvent-induced reentrant phase separation of poly(N-isopropylacrylamide)

Changes in the local environment around amide groups of poly(N-isopropylacrylamide) (PNiPA) during a solvent-induced reentrant phase separation have been investigated by infrared spectroscopy combined with quantum chemical calculations. The addition of methanol or tetrahydrofuran as a cosolvent to an aqueous solution of PNiPA causes spectral changes in the amide I regions. By preparing a dimer model compound for PNiPA, we can establish the assignment of the amide I bands for the polymer in solutions. Hydrogen-deuterium exchange experiments of the amide protons of PNiPA and its dimer models have revealed that the amide groups of PNiPA form an intramolecular C=O...H-N hydrogen bond even in a good solvent. The result has suggested that the change in the amide I envelope of PNiPA observed during the solvent-induced phase transition reflects the modification of the intramolecular C=O...H-N hydrogen bond of PNiPA as well as the variation in solvation state of the amide groups. On the basis of the assignment, we have discussed contributions of the intramolecular C=O...H-N hydrogen bond to the phase behavior of PNiPA.     

An application of coupled reference interaction site model/molecular dynamics to the conformational analysis of the alanine dipeptide

We present an application of our recently proposed coupled reference interaction site model (RISM) molecular dynamics (MD) solvation free energy methodology [Freedman and Truong, Chem. Phys. Lett. 381, 362 (2003); J. Chem. Phys. 121, 2187 (2004)] to study the conformational stability of alanine dipeptide in aqueous solution. In this methodology, radial distribution functions obtained from a single MD simulation are substituted into a RISM expression for solvation free energy. Consequently, iterative solution of the RISM equation is not needed. The relative solvation free energies of seven different conformations of the alanine dipeptide in aqueous solution are calculated. Results from the coupled RISM/MD methodology are in good agreement with those from earlier simulations using the accurate free energy perturbation approach, showing that the alphaR conformation is most stabilized by solution. This study establishes a framework for applying this coupled RISM/MD method to larger biological systems.     

Solvation properties of N-acetyl-β-glucosamine: molecular dynamics study incorporating electrostatic polarization

N-Acetyl-β-glucosamine (NAG) is an important moiety of glycoproteins and is involved in many biological functions. However, conformational and dynamical properties of NAG molecules in aqueous solution, the most common biological environment, remain ambiguous due to limitations of experimental methods. Increasing efforts are made to probe structural properties of NAG and NAG-containing macromolecules, like peptidoglycans and polymeric chitin, at the atomic level using molecular dynamics simulations. In this work, we develop a polarizable carbohydrate force field for NAG and contrast simulation results of various properties using this novel force field and an analogous nonpolarizable (fixed charge) model. Aqueous solutions of NAG and its oligomers are investigated; we explore conformational properties (rotatable bond geometry), electrostatic properties (dipole moment distribution), dynamical properties (self-diffusion coefficient), hydrogen bonding (water bridge structure and dynamics), and free energy of hydration. The fixed-charge carbohydrate force field exhibits deviations from the gas phase relative rotation energy of exocyclic hydroxymethyl side chain and of chair/boat ring distortion. The polarizable force field predicts conformational properties in agreement with corresponding first-principles results. NAG-water hydrogen bonding pattern is studied through radial distribution functions (RDFs) and correlation functions. Intermolecular hydrogen bonding between solute and solvent is found to stabilize NAG solution structures while intramolecular hydrogen bonds define glycosidic linkage geometry of NAG oligomers. The electrostatic component of hydration free energy is highly dependent on force field atomic partial charges, influencing a more favorable free energy of hydration in the fixed-charge model compared to the polarizable model.     

Carbohydrate intramolecular hydrogen bonding cooperativity and its effect on water structure

Molecular dynamics (MD) simulations combined with water-water H-bond angle analysis and calculation of solvent accessible surface area and approximate free energy of solvation were used to determine the influence of hydroxyl orientation on solute hydration and surrounding water structure for a group of chemically identical solutes-the aldohexopyranose sugars. Intramolecular hydrogen bond cooperativity was closely associated with changes in water structure surrounding the aldohexopyranose stereoisomers. The OH-4 group played a pivotal role in hydration as it was able to participate in a number of hydrogen bond networks utilizing the OH-6 group. Networks that terminated within the molecule (OH-4 --> OH-6 --> O-5) had relatively more nonpolar-like hydration than those that ended in a free hydroxyl group (OH-6 --> OH-4 --> OH-3). The OH-2 group modulated the strength of OH-4 networks through syndiaxial OH-2/4 intramolecular hydrogen bonding, which stabilized and induced directionality in the network. Other syndiaxial interactions, such as the one between OH-1 and OH-3, only indirectly affected water structure. Water structure surrounding hydrogen bond networks is discussed in terms of water-water hydrogen bond populations. The impact of syndiaxial versus vicinal hydrogen bonds is also reviewed. The results suggest that biological events such as protein-carbohydrate recognition and cryoprotection by carbohydrates may be driven by intramolecular hydrogen bond cooperativity.     

Solvation free energy of amino acids and side-chain analogues

The solvation free energies of amino acids and their side-chain analogues in water and cyclohexane are calculated by using Monte Carlo simulation. The molecular interactions are described by the OPLS-AA force field for the amino acids and the TIP4P model for water, and the free energies are determined by using the Bennett acceptance method. Results for the side-chain analogues in cyclohexane and in water are used to evaluate the performance of the force field for the van der Waals and the electrostatic interactions, respectively. Comparison of the calculated hydration free energies for the amino acid analogues and the full amino acids allows assessment of the additivity of the side chain contributions on the number of hydrating water molecules. The hydration free energies of neutral amino acids can be reasonably approximated by adding the contributions of their side chains to that of the hydration of glycine. However, significant nonadditivity in the free energy is found for the zwitterionic form of amino acids with polar side chains. In serine and threonine, intramolecular hydrogen bonds are formed between the polar side chains and backbone groups, leading to weaker solvation than for glycine. In contrast, such nonadditivity is not observed in tyrosine, in which the hydroxyl group is farther separated from, and therefore cannot form an intramolecular hydrogen bond with, the backbone. For histidine we find that a water molecule can form a bridge when the intramolecular hydrogen bond between the polar group and the backbone is broken.     

Molecular mechanics of organic halides. V. Conformational equilibria in solution

With a view of using data on solutions and liquids for parameter fitting in molecular mechanical force fields, Abraham's theory of solvation is incorporated in the force field procedure. Geometries and bond moments are estimated internally, partial account being taken of bond–bond induction, and used to calculate the intramolecular electrostatic energy, dipole moment, and the dipole and quadrupole terms in the solvation energy. Three dielectric constants are used, one for the solute in the vapor, one for the solution, and one for the intramolecular space through which dipole–dipole interactions take place. Examples are given, including such where computation differs with measurement, to illustrate the performance of the scheme.     

Theoretical studies of the tautomeric equilibria for five-member N-heterocycles in the gas phase and in solution

Tautomeric equilibria have been studied for five-member N-heterocycles and their methyl derivatives in the gas phase and in different solvents with dielectric constants of epsilon = 4.7-78.4. The free energy changes differently for tautomers upon solvation as compared to the gas phase, resulting in a shift of the equilibrium constant in solution. Solvents with increasing dielectric constant produce more negative solute-solvent interaction energies and increasing internal energies. The methyl-substituted imidazole and pyrrazole form delicate equilibria between two tautomeric forms. Depending on the solvent, the methyl-substituted triazoles and tetrazole have one or two major tautomers in solution. When estimating the relative solvation free energies by means of an explicit solvent model and using the FEP/MC method, one observes that the preferred tautomers differ in several cases from those predicted by the continuum solvent model. The 1,2-prototropic shift, as an intramolecular tautomerization path, requires about 50 kcal/mol activation energy for imidazole in the gas phase, and this route is also disfavored in a solution. The calculated activation free energy along the intramolecular path is 48-50 kcal/mol in chloroform and water as compared to a literature value of 13.6 kcal/mol for pyrrazole in DMSO. A molecular dynamics computer experiment favors the formation of an imidazole chain in chloroform, making the 1,3-tautomerization feasible along an intermolecular path in nonprotic solvents. In aqueous solution, one strong N-H...Ow hydrogen bond is formed for each species, whereas all other nitrogens in the ring form weaker, N...HwOw type hydrogen bonds. The tetrahydrofuran solvent acts as a hydrogen bond acceptor and forms N-H...Oether bonds. Molecules of the dichloromethane solvent are in favorable dipole-dipole interactions with the solute. The results obtained are useful in the design of N-heterocyclic ligands forming specified hydrogen bonds with protein side chains.     

Theoretical analysis on changes in thermodynamic quantities upon protein folding: essential role of hydration

The free energy change associated with the coil-to-native structural transition of protein G in aqueous solution is calculated by using the molecular theory of solvation, also known as the three-dimensional reference interaction site model theory, to uncover the molecular mechanism of protein folding. The free energy is decomposed into the protein intramolecular energy, the hydration energy, and the hydration entropy. The folding is accompanied with a large gain in the protein intramolecular energy. However, it is almost canceled by the correspondingly large loss in the hydration energy due to the dehydration, resulting in the total energy gain about an order of magnitude smaller than might occur in vacuum. The hydration entropy gain is found to be a substantial driving force in protein folding. It is comparable with or even larger than the total energy gain. The total energy gain coupled with the hydration entropy gain is capable of suppressing the conformational entropy loss in the folding. Based on careful analysis of the theoretical results, the authors present a challenging physical picture of protein folding where the overall folding process is driven by the water entropy effect.     

Study of intramolecular interactions in aspirin

A detailed quantum chemical study of the solvent effects in the intramolecular hydrogen bonding, conformational stability, and reactivity of aspirin has been performed using density functional theory (DFT) at the B3LYP/6‐31G(d) theory level. Seven conformational isomers, three of them presenting intramolecular hydrogen bonds, have been located. Thermochemical functions have been computed, and relative energies and free enthalpies have been determined in gas and aqueous phases. Several molecular properties have been calculated to predict the ability of aspirin to acylate cyclooxygenase (COX) enzymes. A six‐membered‐ring hydrogen‐bonded conformer was found to be the most reactive species. The solvation in aqueous phase increases the reactivity and strengthens intramolecular hydrogen bonding.     

MD simulation of the Na+-phenylalanine complex in water: competition between cation-pi interaction and aqueous solvation

The competition between cation-pi interaction and aqueous solvation for the Na+ ion has been investigated by molecular dynamics simulations, using the phenylalanine amino acid as the test pi system. Starting from one of the best standard force fields, we have developed new parameters that significantly improve the agreement with experimental and high quality quantum mechanical results for the complexes of Na+ with phenylalanine, benzene, and water. The modified force field performs very well in forecasting energy and geometry of cation coordination for the complexes. Next, analysis of MD trajectories and steered MD simulations indicate that the Na+-phenylalanine complex survives for a significant time in aqueous solution and that the free energy barrier opposing dissociation of the complex is sizable. Finally, we analyze the role of different intermolecular interactions in determining the preference for cation-pi bonding with respect to aqueous solvation. We thus confirm that the Na+-phenylalanine stabilization energy may overcome the interactions with water.     

Conformation, Inversion Barrier, and Solvent-Induced Conformational Shift in Exo- and Endo/Exo-Calix[4]arenes

Calixarenes 4a and 4b having hydroxyl groups in endo and exo positions and the ethanediyl-bridged exo-calixarene 5a were synthesized by a stepwise strategy. Single-crystal X-ray structures were obtained for 4a and for the exo-calixarene 3d, showing the molecules to exist in the 1,2-alternate conformation which is also found for 4a,b in solution. The inversion barriers of 4a and 4b (10.3 and 10.8 kcal mol(-1)) are similar to that determined for the endo-dihydroxycalixarene 12, indicating that the additional intramolecular hydrogen bond between the exo OH groups does not decrease the flexibility of the molecule. In CDCl(3) solution exo-calixarene 5a adopts a 1,2-alternate conformation with the methyl group at the bridge located in an axial position, while in DMSO-d(6) the conformation adopted is the partial cone. Similar solvent-induced conformational shifts were found for the exo-calixarenes 3b and 3d. MM3 calculations predict that the cone form is the lowest energy conformation of 4 and the exo-calixarenes 3 and 5. The calculations suggest that the conformational preferences of the methyl group at the bridge for either the axial or equatorial positions are in large part determined by the repulsive steric interactions with the hydroxyl groups. The inversion barrier of 4b is satisfactorily reproduced by calculations, which indicate that the rotation of the exo rings is less energetically demanding than the rotation of the endo rings.     

Involvement of water in carbohydrate-protein binding.

The interactions of trimannosides 1 and 2 with Con A were studied to reveal the effects of displacement of well-ordered water molecules on the thermodynamic parameters of protein-ligand complexation. Trisaccharide 2 is a derivative of 1, in which the hydroxyl at C-2 of the central mannose unit is replaced by a hydroxyethyl moiety. Upon binding, this moiety displaces a conserved water molecule present in the Con A binding site. Structural studies by NMR spectroscopy and MD simulations showed that the two compounds have very similar solution conformational properties. MD simulations of the complexes of Con A with 1 and 2 demonstrated that the hydroxyethyl side chain of 2 can establish the same hydrogen bonds in a low energy conformation with the protein binding site as those mediated by the water molecule in the complex of 1 with Con A. Isothermal titration microcalorimetry (ITC) measurements showed that 2 has a more favorable entropy of binding compared to 1. This term, which was expected, arises from the return of the highly ordered water molecule to bulk solution. The favorable entropy term was, however, offset by a relatively large unfavorable enthalpy term. This observation was rationalized by comparing the extent of hydrogen bond and solvation changes during binding. It is proposed that an indirect interaction through a water molecule will provide a larger number of hydrogen bonds in the complex that have higher occupancies than in bulk solution, thereby stabilizing the complex.     

PEARLS: program for energetic analysis of receptor-ligand system

Analysis of the energetics of small molecule ligand-protein, ligand-nucleic acid, and protein-nucleic acid interactions facilitates the quantitative understanding of molecular interactions that regulate the function and conformation of proteins. It has also been extensively used for ranking potential new ligands in virtual drug screening. We developed a Web-based software, PEARLS (Program for Energetic Analysis of Ligand-Receptor Systems), for computing interaction energies of ligand-protein, ligand-nucleic acid, protein-nucleic acid, and ligand-protein-nucleic acid complexes from their 3D structures. AMBER molecular force field, Morse potential, and empirical energy functions are used to compute the van der Waals, electrostatic, hydrogen bond, metal-ligand bonding, and water-mediated hydrogen bond energies between the binding molecules. The change in the solvation free energy of molecular binding is estimated by using an empirical solvation free energy model. Contribution from ligand conformational entropy change is also estimated by a simple model. The computed free energy for a number of PDB ligand-receptor complexes were studied and compared to experimental binding affinity. A substantial degree of correlation between the computed free energy and experimental binding affinity was found, which suggests that PEARLS may be useful in facilitating energetic analysis of ligand-protein, ligand-nucleic acid, and protein-nucleic acid interactions. PEARLS can be accessed at http://ang.cz3.nus.edu.sg/cgi-bin/prog/rune.pl.     

Analyzing the robustness of the MM/PBSA free energy calculation method: application to DNA conformational transitions

The ability to predict and characterize free energy differences associated with conformational equilibria or the binding of biomolecules is vital to understanding the molecular basis of many important biological functions. As biological studies focus on larger molecular complexes and properties of the genome, proteome, and interactome, the development and characterization of efficient methods for calculating free energy becomes increasingly essential. The aim of this study is to examine the robustness of the end-point free energy method termed the molecular mechanics Poisson-Boltzmann solvent accessible surface area (MM/PBSA) method. Specifically, applications of MM/PBSA to the conformational equilibria of nucleic acid (NA) systems are explored. This is achieved by comparing A to B form DNA conformational free energy differences calculated using MM/PBSA with corresponding free energy differences determined with a more rigorous and time-consuming umbrella sampling algorithm. In addition, the robustness of NA MM/PBSA calculations is also evaluated in terms of the sensitivity towards the choice of force field and the choice of solvent model used during conformational sampling. MM/PBSA calculations of the free energy difference between A-form and B-form DNA are shown to be in very close agreement with the PMF result determined using an umbrella sampling approach. Further, it is found that the MM/PBSA conformational free energy differences were also in agreement using either the CHARMM or AMBER force field. The influence of ionic strength on conformational stability was particularly insensitive to the choice of force field. Finally, it is also shown that the use of a generalized Born implicit solvent during conformational sampling results in free energy estimates that deviate slightly from those obtained using explicitly solvated MD simulations in these NA systems.     

SAMPL5: 3D-RISM partition coefficient calculations with partial molar volume corrections and solute conformational sampling

Implicit solvent methods for classical molecular modeling are frequently used to provide fast, physics-based hydration free energies of macromolecules. Less commonly considered is the transferability of these methods to other solvents. The Statistical Assessment of Modeling of Proteins and Ligands 5 (SAMPL5) distribution coefficient dataset and the accompanying explicit solvent partition coefficient reference calculations provide a direct test of solvent model transferability. Here we use the 3D reference interaction site model (3D-RISM) statistical-mechanical solvation theory, with a well tested water model and a new united atom cyclohexane model, to calculate partition coefficients for the SAMPL5 dataset. The cyclohexane model performed well in training and testing (\(R=0.98\) for amino acid neutral side chain analogues) but only if a parameterized solvation free energy correction was used. In contrast, the same protocol, using single solute conformations, performed poorly on the SAMPL5 dataset, obtaining \(R=0.73\) compared to the reference partition coefficients, likely due to the much larger solute sizes. Including solute conformational sampling through molecular dynamics coupled with 3D-RISM (MD/3D-RISM) improved agreement with the reference calculation to \(R=0.93\). Since our initial calculations only considered partition coefficients and not distribution coefficients, solute sampling provided little benefit comparing against experiment, where ionized and tautomer states are more important. Applying a simple \(\hbox {p}K_{\text {a}}\) correction improved agreement with experiment from \(R=0.54\) to \(R=0.66\), despite a small number of outliers. Better agreement is possible by accounting for tautomers and improving the ionization correction.     

On the origins of core-electron chemical shifts of small biomolecules in aqueous solution: insights from photoemission and ab initio calculations of glycine(aq)

The local electronic structure of glycine in neutral, basic, and acidic aqueous solution is studied experimentally by X-ray photoelectron spectroscopy and theoretically by molecular dynamics simulations accompanied by first-principle electronic structure and spectrum calculations. Measured and computed nitrogen and carbon 1s binding energies are assigned to different local atomic environments, which are shown to be sensitive to the protonation/deprotonation of the amino and carboxyl functional groups at different pH values. We report the first accurate computation of core-level chemical shifts of an aqueous solute in various protonation states and explicitly show how the distributions of photoelectron binding energies (core-level peak widths) are related to the details of the hydrogen bond configurations, i.e. the geometries of the water solvation shell and the associated electronic screening. The comparison between the experiments and calculations further enables the separation of protonation-induced (covalent) and solvent-induced (electrostatic) screening contributions to the chemical shifts in the aqueous phase. The present core-level line shape analysis facilitates an accurate interpretation of photoelectron spectra from larger biomolecular solutes than glycine.     

Effects of substituting a OH group by a F atom in D-glucose. Ab initio and DFT analysis.

High-level ab initio and DFT methods up to MP2/6-311++G//B3LYP/6-31G and B3LYP/6-311++G//B3LYP/6-31G levels have been used to assess the relative energies of 17 different structures of D-glucose and 13 different structures of 4-deoxy-4-fluoro-D-glucose. The structures were confirmed to correspond to minima on the potential energy surface at the RHF/6-31G level. Solvation Model 5.4/AM1 was used to calculate the effects of aqueous solution. The substitution of a OH group by a F atom does not much change the shape and electrostatic potential around corresponding conformers, but in the gas phase it destabilizes the cooperative network of intramolecular hydrogen bonds. This destabilization mostly affects structures with a chain of intramolecular hydrogen bonds oriented counterclockwise, as fluorine is unable to donate a hydrogen bond and therefore causes a gap in the chain. In contrast, for clockwise-oriented networks of hydrogen bonds, the fluorine can act as an acceptor at the end of a chain of cooperative hydrogen bonds. A slightly higher energy of anomeric and exo-anomeric stabilization is another effect of substituting the fourth hydroxyl group by a fluorine atom in D-glucose, observed both in the gas phase and in aqueous solution. For this reason, the alpha anomers contribute more to the equilibrium population of structures of 4-deoxy-4-fluoro-D-glucose than D-glucose. In aqueous solution, both D-glucose and its 4-deoxy-4-fluoro analogue are present as a mixture of mainly three corresponding structures. This indicates that 4-deoxy-4-fluoro-D-glucose is a good substitute for D-glucose in terms of its biochemical and biological activity. Moreover, this suggests that, for molecules with limited conformational freedom, the substitution of a OH group by a F atom is very likely to lead to a potential new drug. In contrast, it had already been shown that, for conformationally labile aliphatic compounds, replacement of a hydroxyl by a fluorine increases conformational diversity, so the fluorine-containing aliphatic molecules were not likely to be an example of a successful drug design. On the other hand, this work shows that, among molecules with limited conformational freedom, such as cyclic compounds, one is very likely to find targets for a successful rational drug design.