Theoretical studies on the conformation of saccharides. VIII. Solvent effect on the stability of β-cellobiose conformers

The conformational equilibria of five β-cellobiose conformers have been studied theoretically in 10 solvents. The stability of the conformers in dilute solution has been compared by using the method that has already been tested for 2-methoxytetrahydropyran, β-maltose, and D -glucose. The solvation energy consists of electrostatic, dispersion, and cavity terms which have been determined from the properties of the solute calculated by the PCILO quantum-chemical method and physicochemical properties of the solvents. The calculated abundance of conformers depends on the solvent (e.g., in dioxane C1:C2:C3:C4:C5 = 60.0:34.1:2.9:2.0:1.0; in dimethylsulfoxide, 75.5:22.1:1.8:0.5:0.2; and in water, 82.2:16.2:1.3:0.2:0.1). The results obtained indicate that the preponderant conformer in the aqueous solution is similar to the one adopted by β-cellobiose in the crystalline form. The role of individual contributions to the solvation energy have been analyzed. Based on the determined abundance of conformers, averaged residual optical activity and nmr parameters have been calculated and compared with observable properties. The marked differences observed between solvent-induced conformational changes for β-cellobiose and β-maltose have been discussed from the viewpoint of the solubility of the cellulose.     

Solvent effect on the stability of mannobiose conformers

The conformational equilibria of seven methyl β-D -mannobioside conformers have been studied theoretically in five solvents. The structure of each individual conformer has been refined by the PCILO quantum-chemical method from the seven distinct low-energy regions determined from (Φ, Ψ) maps calculated by a potential function method. The stability of the conformers in dilute solution has been evaluated by using a method that consists of electrostatic, dispersion, and cavity terms. The calculated abundance of conformers depends on the solvent and results indicate that the preponderant conformer in the solution may not be the one adopted by mannobiose in the crystalline form. Based on the determined abundance of conformers, thermodynamically averaged nmr parameters, dipole moment, and linkage rotation have been calculated. The solvation behavior of methyl β-D -mannobioside is compared to those previously estimated for cellobiose and maltose.     

Theoretical studies on the conformation of saccharides. XIV. Structure and conformational properties of the glycosylamines

The 2-methylaminotetrahydropyran was used as a model to study conformational properties of the N-glycosidic linkage in glycosylamines. Relaxed two-dimensional conformational (phi, psi) maps in 20 solvents were calculated by a method in which the total energy is divided into the energy of the isolated molecule and the solvation energy. Molecular geometry optimization has been carried out for each conformer using the quantum chemical method PCILO. The calculated variations of the geometry are consistent with the results obtained by the statistical analysis of available experimental data retrieved from the Cambridge Structural Database. The calculated abundances of conformers show that the polarity of the solvent has little effect on the anomeric ratio, and the form having the methylamino group equatorial is favored in all considered solvents.     

Effects of electric field on alamethicin bound at the lipid-water interface: a molecular mechanics study.

A systematic molecular mechanics study of the alamethicin molecule was made to determine a set of low-energy conformers in vacuo and in aqueous environment. The behavior of these conformers was investigated at the phase boundary which was modeled as a plane dividing two compartments with solvation properties of water and octanol with a constant electric field applied normal to the boundary. The calculations were performed with a molecular mechanics program for calculation of stable conformations at the phase boundary utilizing the Empiric Conformational Energy Program for Peptides force field and the Hopfinger-Scheraga solvation model. 371 minimum energy conformers of alamethicin, determined in vacuo with the build-up procedure, were used as starting conformations for energy minimization in aqueous environment and at the phase boundary. Only 49 interphase-bound structures were within 12 kcal/mol of the minima which was found. No helical structures having values close to the canonical parameters for an alpha- or 3(10)-helix were found despite the presence of eight alpha-methylalanine residues which favor the formation of these helices; four helix-like structures were found, having all negative phi, psi values. All the helical conformers have very high energies in water (approximately 14 kcal/mol), but are quite stable at the phase boundary (3.7-6.8 kcal/mol above the lowest minima found). The implications of these results for proposed mechanisms for membrane-binding and voltage-dependent gating are considered.     

Effect of partial atomic charges on the calculated free energy of solvation of poly(vinyl alcohol) in selected solvents

It is well-known that properties of poly(vinyl alcohol) (PVA) in the pure and solution states depend largely on the hydrogen bonding networks formed. In the context of molecular simulation, such networks are handled through the Coulombic interactions. Therefore, a good set of partial atom charges (PACs) for simulations involving PVA is highly desirable. In this work, we calculated the PACs for PVA using a few commonly used population analysis schemes with a hope to identify an accurate set of PACs for PVA monomers. To evaluate the quality of the calculated parameters, we have benchmarked their predictions for free energy of solvation (FES) in selected solvents by molecular dynamics simulations against the ab initio calculated values. Selected solvents were water, ethanol and benzene as they covered a range of size and polarity. Also, PVA with different tacticities were used to capture their effect on the calculated FESs. Based on our results, neither PACs nor FESs are affected by the chain tacticity. While PACs predicted by the Merz-Singh-Kollman scheme were close to original values in the OPLS-AA force field in way that no significant difference in properties of pure PVA was observed, free energy of solvation calculated using such PACs showed greater agreement with ab initio calculated values than those calculated by OPLS-AA (and all other schemes used in this work) in all three solvents considered.     

Thermochemistry of dissolving glycine, glycyl-glycine, and diglycyl-glycine in a mixed water-dimethylsulfoxide solvent at 298.15 K

The enthalpies of dissolving glycine, glycyl-glycine and diglycyl-glycine (deltaH(soln)0) in a mixed water-dimethylsulfoxide (DMSO) solvent was determined by the calorimetric method in the range of concentrations of the organic component 0 < X2 < 0.4 m.d. at 298.15 K. The enthalpies of solvation ((deltaH(solv)0) and transfer ((deltaH(tr)0) of these compounds from water to a mixed solvent were calculated. The dependencies deltaH(tr)0 =f(X2) were found to be extreme, indicating complex intermolecular interactions between the solution components. The influence of the structure and the properties of the substances dissolved and the composition of the mixture and the nature of organic solvent on their thermochemical characteristics was studied. The coefficients of enthalpy for pair interactions of glycine and its oligomers with DMSO molecules were calculated. These have positive values and increase in the order: glycyl-glycine < glycine < diglycyl-glycine. The changes in the thermochemical characteristics of dissolving, transfer, and solvation of glycine and its olygomers were shown to be determined by the energy of the mixed solvent formation, the nature of the organic solvents, and the structure of amino acids and peptides.     

Combined theoretical studies on solvation and hydrogen bond interactions in glycine tripeptide

The theoretical conformational analysis of glycine tripeptide (GT) has been carried out by molecular dynamics (MD) method in order to find minimum energy conformations. The MD studies on GT with water have been carried out for over 10 ns with a time step of 2 fs using fixed charge force field (AMBER ff03). By adding the solvation effect using water as a solvent, the GT conformers identified in this study exhibit α-helical conformation. Compared with the earlier reports, this MD study is able to identify the energetically favourable GT conformations. The obtained geometry of the five most stable GT conformations was optimised using the density functional theory method at B3LYP/6-311G** level of theory. Subsequently, the effects of solvation on the conformational characteristics of five most stable GT conformers with four water molecules (the number of water molecules in the first solvation shell of GT obtained from MD study) were investigated using the same method and the same level of theory. The effect of microsolvation on the fifth GT conformer has been studied with a cluster of 11 water molecules as the first hydration shell which generates folded structure. The interaction energies of all the complexes are calculated by correcting the basis set superposition error. The strong hydrogen bond mainly contributes to the interaction energies. The atoms in molecules theory and natural bond orbital analysis were used to study the origin of H-bonds. A good correlation between the structural parameters and the properties of charge density is found. NMR calculations show that the C = O carbons of the amine groups of the first and middle glycine fragments have maximum chemical shifts.     

Solution conformation of maltose from optical rotation: A procedure for evaluating carbohydrate force fields

As one means of evaluating the various force fields that are being developed for the molecular modeling of carbohydrate structures, the statistically averaged optical rotation of β-maltose can be calculated using the energy surface to be tested in combination with the detailed conformation dependence of β-maltose optical rotation presented here. Comparison with the observed rotation will indicate the quality of the force field, especially with respect to the parameterization of solvent effects. © 1992 John Wiley & Sons, Inc.     

Intraglobular electrostatic field of an enzyme. Calculation of the dissociation constant of a protein ionogenic group. Dissociation of Asp-102 chymotrypsin

Factors determining the change of dissociation constants of ionogenic groups upon their transfer from water into a protein globule are considered. The change of the short-range interaction is simulated with the aid of two model solvents: dimethylformamide (DMF) and formamide (FA). The change of Bornian solvation energy is calculated taking into account the interaction of ions situated inside a protein globule with the surrounding electrolyte solution. The change of ion energy due to the intraglobular electric field preexisting in the enzyme molecule is calculated. Each of the factors listed above gives a large contribution into the ion energy in the case of Asp-102 these contributions compensate each other to a great extent.     

Empirical solvation models can be used to differentiate native from near-native conformations of bovine pancreatic trypsin inhibitor.

Several hydration models for peptides and proteins based on solvent accessible surface area have been proposed previously. We have evaluated some of these models as well as four new ones in the context of near-native conformations of a protein. In addition, we propose an empirical site-site distance-dependent correction that can be used in conjunction with any of these models. The set of near-native structures consisted of 39 conformations of bovine pancreatic trypsin inhibitor (BPTI) each of which was a local minimum of an empirical energy function (ECEPP) in the absence of solvent. Root-mean-square (rms) deviations from the crystallographically determined structure were in the following ranges: 1.06-1.94 A for all heavy atoms, 0.77-1.36 A for all backbone heavy atoms, 0.68-1.33 A for all alpha-carbon atoms, and 1.41-2.72 A for all side-chain heavy atoms. We have found that there is considerable variation among the solvent models when evaluated in terms of concordance between the solvation free energy and the rms deviations from the crystallographically determined conformation. The solvation model for which the best concordance (0.939) with the rms deviations of the C alpha atoms was found was derived from NMR coupling constants of peptides in water combined with an exponential site-site distance dependence of the potential of mean force. Our results indicate that solvation free energy parameters derived from nonpeptide free energies of hydration may not be transferrable to peptides. Parameters derived from peptide and protein data may be more applicable to conformational analysis of proteins. A general approach to derive parameters for free energy of hydration from ensemble-averaged properties of peptides in solution is described.     

A solvent model for simulations of peptides in bilayers. I. Membrane-promoting alpha-helix formation

We describe an efficient solvation model for proteins. In this model atomic solvation parameters imitating the hydrocarbon core of a membrane, water, and weak polar solvent (octanol) were developed. An optimal number of solvation parameters was chosen based on analysis of atomic hydrophobicities and fitting experimental free energies of gas-cyclohexane, gas-water, and octanol-water transfer for amino acids. The solvation energy term incorporated into the ECEPP/2 potential energy function was tested in Monte Carlo simulations of a number of small peptides with known energies of bilayer-water and octanol-water transfer. The calculated properties were shown to agree reasonably well with the experimental data. Furthermore, the solvation model was used to assess membrane-promoting alpha-helix formation. To accomplish this, all-atom models of 20-residue homopolypeptides-poly-Leu, poly-Val, poly-Ile, and poly-Gly in initial random coil conformation-were subjected to nonrestrained Monte Carlo conformational search in vacuo and with the solvation terms mimicking the water and hydrophobic parts of the bilayer. All the peptides demonstrated their largest helix-forming tendencies in a nonpolar environment, where the lowest-energy conformers of poly-Leu, Val, Ile revealed 100, 95, and 80% of alpha-helical content, respectively. Energetic and conformational properties of Gly in all environments were shown to be different from those observed for residues with hydrophobic side chains. Applications of the solvation model to simulations of peptides and proteins in the presence of membrane, along with limitations of the approach, are discussed.     

Calculation of the thermodynamic quantities of solvation for lysozyme and β-lactoglobulin a in guanidinium chloride

The thermodynamic quantities of solvation, that is, Gibbs free energy, enthalpy, and entropy for the transfer of lysozyme and β-lactoglobulin A from water to 6M guanidinium chloride solution, were calculated from the respective contributions of constituent groups. Comparison of the calculated values with the experimental ones reveals satisfactory agreement for lysozyme, whereas for β-lactoglobulin A, it is semiquantitative. Thus the method of calculation appears to be fundamentally sound. The reasons for the differences in calculated and experimental value are briefly discussed.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

Solvent effect on the synergic extraction of thulium(III) with acetylacetone and 1,10-phenanthroline

The synergic extraction equilibrium of Tm(III) with acetylacetone (Hacac) and 1,10-phenanthroline (phen) in various organic solvents has been studied. The adduct formation constants, β_s, for Tm(acac)₃^? phen, were determined in heptane, cyclohexane, carbon tetrachloride, benzene and chloroform. The solvent effect on β_s is explained in connection with the activity coefficients of the neutral ligand, the chelate, and the adduct in the organic solvent. The activity coefficients can be calculated from the corresponding solubility parameters on the basis of the regular solution theory, and the solubility parameters of the solutes were estimated from their two-phase partition coefficients. It is demonstrated that β_s in different organic solvents except those having a specific interaction with the solute, such as chloroform, can be calculated by the present approach.     

CD conformational studies on synthetic peptides encompassing the processing domain of the ocytocin/neurophysin precursor

Synthetic peptides of different size, reproducing the proteolytic processing site of proocytocin, were studied by CD under several experimental conditions in order to ascertain the ability of different solvents to stabilize secondary structural motifs, such as α-helix tracts and β-turns. A combination of deconvolution methods and empirical calculations subtracting the contributions due to unordered structures from the spectra suggests that in solution (a) mainly two distinct families of ordered conformers containing structurally different β-turns are present, (b) the relative stability of the different conformers depends from the nature of the solvent, and (c) in the case of the larger peptides, a population containing an α-helical conformation is also present. From the biological point of view the presence of at least two families of ordered conformers could be in line with current theories assuming that the catalytic effect of the receptor microenvironment may be determinant in shifting the equilibrium toward the active conformation. © 1997 John Wiley & Sons, Inc. Biopoly 41 : 461–479, 1997     

Computational investigation on microsolvation of the osmolyte glycine betaine [GB (H2O)1-7]

The preferential interactions of glycine betaine (GB) with solvent components and the effect of solvent on its stability have been examined. In particular, the microsolvation of organic osmolyte and widely important osmoprotectant in nature as glycine betaine has been reported by using M06 method. A number of configurations (b_X (a-z)) of the clusters for one to seven water molecules (×?=?1-7) have been considered for the microsolvation. Structures of stable conformers are obtained and denoted as b1a, b2a, b3a, b4a, b5a, b6a and b7a. It is observed from the interaction energy difference (?E) that only seven water molecules can be accommodated in the first solvation shell to stabilize GB. It is also observed that the calculated relative energy using M06 is in close agreement with calculations at the MP2 level of theory.

Calculation of the free energy of association for protein complexes.

We have developed a method for calculating the association energy of quaternary complexes starting from their atomic coordinates. The association energy is described as the sum of two solvation terms and an energy term to account for the loss of translational and rotational entropy. The calculated solvation energy, using atomic solvation parameters and the solvent accessible surface areas, has a correlation of 96% with experimentally determined values. We have applied this methodology to examine intermediates in viral assembly and to assess the contribution isomerization makes to the association energy of molecular complexes. In addition, we have shown that the calculated association can be used as a predictive tool for analyzing modeled molecular complexes.     

Solution structure and dynamics of cyclic and acyclic cholinergic agonists.

Two classes of nicotinic cholinergic agonists, which vary in flexibility and electronegativity, have been synthesized, and their structural and dynamic properties have been studied with nuclear magnetic resonance (NMR) spectroscopy. Although the compounds are chemically identical except for the presence or absence of one cyclicizing C--C bond, single channel recording and radioligand binding studies have shown that the cyclic compounds are considerably more potent than the acyclic derivatives (McGroddy, K.A., A.A. Carter, M.M. Tubbert, and R.E. Oswald. 1993. Biophys. J. 64:325-338). Using one- and two-dimensional NMR spectroscopy, we have shown that these molecules exist in two distinct stable conformers, which differ in the orientation of the amide bond. The cyclic 1,1-dimethyl-4-trifluoroacetyl-piperazinium iodide and its trifluoromethyl derivative compounds are symmetric, and the two conformers are of equal energy. The acyclic N,N,N,N'-tetramethyl-N'-acetylethylene-diamine iodide (TED) and its trifluoromethyl derivative derivatives, however, populate two energetically unequal solution conformations. Using variable temperature NMR spectroscopy on these molecules and their uncharged precursors, we have characterized the energetics of amide bond isomerization and have distinguished steric and electrostatic contributions to the equilibrium between the two conformers. The more populated TED conformer has the amide methyl group trans to the carbonyl oxygen, and it is stabilized by an electrostatic attraction between the partially negative carbonyl oxygen and the positively charged quaternary amine nitrogen. As discussed in the accompanying paper (McGroddy, K.A., A.A. Carter, M.M. Tubbert, and R.E. Oswald. 1993. Biophys. J. 64:325-338), the differences in the stable solution structures of the TED derivatives and their interconversion kinetics may be of biological significance.