Alcohols, ethers, carbohydrates, and related compounds. III. The 1,2-dimethoxyethane system

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Additive and Classical Drude Polarizable Force Fields for Linear and Cyclic Ethers

Empirical force field parameters consistent with the CHARMM additive and classical Drude based polarizable force fields are presented for linear and cyclic ethers. Initiation of the optimization process involved validation of the aliphatic parameters based on linear alkanes and cyclic alkanes. Results showed the transfer to cyclohexane to yield satisfactory agreement with target data; however, in the case of cyclopentane direct transfer of the Lennard-Jones parameters was not sufficient due to ring strain, requiring additional optimization of these parameters for this molecule. Parameters for the ethers were then developed starting with the available aliphatic parameters, with the nonbond parameters for the oxygens optimized to reproduce both gas- and condensed-phase properties. Nonbond parameters for the polarizable model include the use of an anisotropic electrostatic model on the oxygens. Parameter optimization emphasized the development of transferable parameters between the ethers of a given class. The ether models are shown to be in satisfactory agreement with both pure solvent and aqueous solvation properties, and the resulting parameters are transferable to test molecules. The presented force field will allow for simulation studies of ethers in condensed phase and provides a basis for ongoing developments in both additive and polarizable force fields for biological molecules.

Force fields for homology modeling

Accurate all-atom energy functions are crucial for successful high-resolution protein structure prediction. In this chapter, we review both physics-based force fields and knowledge-based potentials used in protein modeling. Because it is important to calculate the energy as accurately as possible given the limitations imposed by sampling convergence, different components of the energy, and force fields representing them to varying degrees of detail and complexity are discussed. Force fields using Cartesian as well as torsion angle representations of protein geometry are covered. Since solvent is important for protein energetics, different aqueous and membrane solvation models for protein simulations are also described. Finally, we summarize recent progress in protein structure refinement using new force fields.

A reoptimized GROMOS force field for hexopyranose-based carbohydrates accounting for the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers

This article presents a reoptimization of the GROMOS 53A6 force field for hexopyranose-based carbohydrates (nearly equivalent to 45A4 for pure carbohydrate systems) into a new version 56A(CARBO) (nearly equivalent to 53A6 for non-carbohydrate systems). This reoptimization was found necessary to repair a number of shortcomings of the 53A6 (45A4) parameter set and to extend the scope of the force field to properties that had not been included previously into the parameterization procedure. The new 56A(CARBO) force field is characterized by: (i) the formulation of systematic build-up rules for the automatic generation of force-field topologies over a large class of compounds including (but not restricted to) unfunctionalized polyhexopyranoses with arbritrary connectivities; (ii) the systematic use of enhanced sampling methods for inclusion of experimental thermodynamic data concerning slow or unphysical processes into the parameterization procedure; and (iii) an extensive validation against available experimental data in solution and, to a limited extent, theoretical (quantum-mechanical) data in the gas phase. At present, the 56A(CARBO) force field is restricted to compounds of the elements C, O, and H presenting single bonds only, no oxygen functions other than alcohol, ether, hemiacetal, or acetal, and no cyclic segments other than six-membered rings (separated by at least one intermediate atom). After calibration, this force field is shown to reproduce well the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers. As a result, the 56A(CARBO) force field should be suitable for: (i) the characterization of the dynamics of pyranose ring conformational transitions (in simulations on the microsecond timescale); (ii) the investigation of systems where alternative ring conformations become significantly populated; (iii) the investigation of anomerization or epimerization in terms of free-energy differences; and (iv) the design of simulation approaches accelerating the anomerization process along an unphysical pathway.

Selenoglycosides in silico: ab initio-derived reparameterization of MM4, conformational analysis using histo-blood group ABH antigens and lectin docking as indication for potential of bioactivity

The identification of glycan epitopes such as the histo-blood group ABH determinants as docking sites for bacterial/viral infections and signals in growth regulation fuels the interest to develop non-hydrolysable mimetics for therapeutic applications. Inevitably, the required substitution of the linkage oxygen atom will alter the derivative's topology. Our study addresses the question of the impact of substitution of oxygen by selenium. In order to characterize spatial parameters and flexibility of selenoglycosides, we first performed ab initio calculations on model compounds to refine the MM4 force field. The following application of the resulting MM4R version appears to reduce the difference to ab initio data when compared to using the MM4 estimator. Systematic conformational searches on the derivatives of histo-blood group ABH antigens revealed increased flexibility with acquisition of additional low-energy conformer(s), akin to the behavior of S-glycosides. Docking analysis using the Glide program for eight test cases indicated potential for bioactivity, giving further experimental investigation a clear direction to testing Se-glycosides as lectin ligands.

Molecular simulations of carbohydrates and protein-carbohydrate interactions: motivation, issues and prospects

The characterization of the 3D structure of oligosaccharides, their conjugates and analogs is particularly challenging for traditional experimental methods. Molecular simulation methods provide a basis for interpreting sparse experimental data and for independently predicting conformational and dynamic properties of glycans. Here, we summarize and analyze the issues associated with modeling carbohydrates, with a detailed discussion of four of the most recently developed carbohydrate force fields, reviewed in terms of applicability to natural glycans, carbohydrate-protein complexes and the emerging area of glycomimetic drugs. In addition, we discuss prospectives and new applications of carbohydrate modeling in drug discovery.

Conformational analysis of thioglycoside derivatives of histo-blood group ABH antigens using an ab initio-derived reparameterization of MM4: implications for design of non-hydrolysable mimetics

Histo-blood group ABH antigens serve as recognition sites for infectious microorganisms and tissue lectins in intercellular communication, e.g. in tumor progression. Thus, they are of interest as a starting point for drug design. In this respect, potent non-hydrolysable derivatives such as thioglycosides are of special interest. As prerequisite to enable estimations of ligand properties relative to their natural counterparts, conformational properties of the thioglycosidic derivatives of ABH trisaccharides and their disaccharide units were calculated using systematic and filtered systematic searches with the MM4 force field. Parameters for the glycosidic torsions of thioglycosides were independently derived from ab initio calculations. The resulting energy deviations required a reparameterization of MM4 to a new parameter set called MM4R. The data sets obtained using MM4R reveal that the thioglycosides have somewhat increased levels of flexibility about the major low-energy conformations shared with the corresponding O-glycosides. In the trisaccharides, the thiosubstitution of the Gal[NAc]alpha1-3Gal linkage leads to a preference for a conformation which is the secondary minimum of the natural counterparts. This conformation also generates contacts between the N-acetyl group and the fucose moiety in the blood group A derivative. Calculations further indicate that thiosubstitution of only the Fucalpha1-2Gal linkage does not affect the conformational preferences compared to the natural trisaccharide. Thiosubstitution of both linkages in the trisaccharide results in increased flexibility but the favored conformation of the natural trisaccharides is preferred. The study suggests that thioglycoside derivatives of ABH antigens could have pharmaceutical interest as ligands of lectins and other carbohydrate-binding proteins.

The effect of electrostatic interactions on conformational equilibria of multiply substituted tetrahydropyran oxocarbenium ions

The three-dimensional structures of dioxocarbenium ions related to glycosyl cations were determined by an analysis of spectroscopic, computational, and reactivity data. Hypothetical low-energy structures of the dioxocarbenium ions were correlated with both experimentally determined (1)H NMR coupling constants and diastereoselectivity results from nucleophilic substitution reactions. This method confirmed the pseudoaxial preference of C-3 alkoxy-substituted systems and revealed the conformational preference of the C-5 alkoxymethyl group. Although the monosubstituted C-5 alkoxymethyl substituent preferred a pseudoequatorial orientation, the C-5-C-6 bond rotation was controlled by an electrostatic effect. The preferred diaxial conformer of the trans-4,5-disubstituted tetrahydropyranyl system underscored the importance of electrostatic effects in dictating conformational equilibria. In the 2-deoxymannose system, although steric effects influenced the orientation of the C-5 alkoxymethyl substituent, the all-axial conformer was favored because of electrostatic stabilization.

The important role of lone-pairs in force field (MM4) calculations on hydrogen bonding in alcohols

An expanded treatment of hydrogen bonding has been developed for MM4 force field calculations, which is an extension from the traditional van der Waals-electrostatic model. It adds explicit hydrogen-bond angularity by the inclusion of lone-pair directionality. The vectors that account for this directionality are placed along the hydrogen acceptor and its chemically intuitive electron pairs. No physical lone-pairs are used in the calculations. Instead, an H-bond angularity function, and a lone-pair directionality function, are incorporated into the hydrogen-bond term. The inclusion of the lone-pair directionality results in improved accuracy in hydrogen-bonded geometries and interaction energies. In this work is described hydrogen bonding in alcohols, and also in water and hydrogen fluoride dimer. The extension to other compounds such as aldehydes, ketones, amides, and so on is straightforward and will be discussed in future work. The conformational energies of ethylene glycol are discussed.

Molecular mechanics (MM4) study of amines

The MM4 force field has been extended to include aliphatic amines. About 20 amines have been examined to obtain a set of useful molecular mechanics parameters for this class. The vibrational spectra of seven amines (172 frequencies) calculated by MM4 have an overall rms error of 27 cm(-1), compared with corresponding MM4 value of 24 cm(-1) for alkanes. The rms and signed average errors of the moments of inertia of nine simple amines compared with the experimental data were 0.18% and -0.004%, respectively. The heats of formation of 30 amines were also studied. The MM4 weighted standard deviation is 0.41 kcal/mol, compared with experiment. Electronegativity effects occur in the hydrocarbon portion of an amine from the nitrogen, and are accounted for by including electronegativity induced changes in bond lengths and angles, and induced dipole-dipole interactions in the molecule. Negative hyperconjugation results from the presence of the lone pair of electrons on nitrogen, and leads to the Bohlmann bands in the infrared, and also to strong and unusual geometric changes in the molecules (Bohlmann effect), all of which are fairly well accounted for. The conformational energies in amines appear to be less straightforward than those for most other classes of molecules, apparently because of the Bohlmann effect, and these are probably not yet completely understood. In general, the agreement between the MM4 calculated results and the available data is reasonably good.

Conformation, dynamics, solvation and relative stabilities of selected β-hexopyranoses in water: a molecular dynamics study with the gromos 45A4 force field

The present article reports long timescale (200 ns) simulations of four beta-D-hexopyranoses (beta-D-glucose, beta-D-mannose, beta-D-galactose and beta-D-talose) using explicit-solvent (water) molecular dynamics and vacuum stochastic dynamics simulations together with the GROMOS 45A4 force field. Free-energy and solvation free-energy differences between the four compounds are also calculated using thermodynamic integration. Along with previous experimental findings, the present results suggest that the formation of intramolecular hydrogen-bonds in water is an 'opportunistic' consequence of the close proximity of hydrogen-bonding groups, rather than a major conformational driving force promoting this proximity. In particular, the conformational preferences of the hydroxymethyl group in aqueous environment appear to be dominated by 1,3-syn-diaxial repulsion, with gauche and solvation effects being secondary, and intramolecular hydrogen-bonding essentially negligible. The rotational dynamics of the exocyclic hydroxyl groups, which cannot be probed experimentally, is found to be rapid (10-100 ps timescale) and correlated (flip-flop hydrogen-bonds interconverting preferentially through an asynchronous disrotatory pathway). Structured solvent environments are observed between the ring and lactol oxygen atoms, as well as between the 4-OH and hydroxymethyl groups. The calculated stability differences between the four compounds are dominated by intramolecular effects, while the corresponding differences in solvation free energies are small. An inversion of the stereochemistry at either C(2) or C(4) from equatorial to axial is associated with a raise in free energy. Finally, the particularly low hydrophilicity of beta-D-talose appears to be caused by the formation of a high-occurrence hydrogen-bonded bridge between the 1,3-syn-diaxial 2-OH and 4-OH groups. Overall, good agreement is found with available experimental and theoretical data on the structural, dynamical, solvation and energetic properties of these compounds. However, this detailed comparison also reveals some discrepancies, suggesting the need (and providing a solid basis) for further refinement.

Octa-O-propanoyl-β-maltose: crystal structure, acyl stacking, related structures, and conformational analysis

The crystal structure of beta-maltose octapropanoate (1) was solved to improve understanding of di-, oligo-, and polysaccharide conformations. The O6 and O6' atoms are in gg and gt orientations, respectively. Extrapolation of the coordinates of the non-reducing residue and observed linkage bond and torsion angles of 1 [Formula: see text] yields a left-handed helix similar to amylose triacetate I. The phi and psi values of 1 are also similar to those of other crystalline, acylated maltose compounds as well as some hydroxyl-bearing molecules. Acylated maltose moieties are often stabilized by stacking of the carbonyl groups and alpha-carbons on O3 and O2' as well as by the exo-anomeric effect. The conformation of 1 is within the 1-kcal/mol contour on a hybrid energy map built with a dielectric constant of 7.5, but corresponds to higher energies on maps made with lower dielectric constants. In one region of phi,psi space, both hydroxyl-bearing and derivatized maltose moieties are found but no inter-residue, intramolecular hydrogen-bonding occurs. In another region, only hydroxyl-bearing molecules crystallize and O2'...O3 hydrogen bonds are always found. In agreement with the energy surfaces, amylose helices extrapolated from available linkage geometries were almost all left-handed.

Linkage and pyranosyl ring twisting in cyclodextrins

Acylated beta-cyclodextrins (beta-CDs) were studied to gain perspective on maltose octapropanoate, the crystal structure of which was reported in the preceding paper in this issue. Acylated beta-CDs are distorted so we looked at other CDs and gained increased understanding of distortion in CDs and possibly, shapes in starch. Classic CDs have six to eight glucose residues in a doughnut shape that is stabilized by a ring of inter-residue O3,,,O2' hydrogen bonds. On a phi,psi energy map for a maltose analog that does not form hydrogen bonds, classic CD linkages have higher energies than structures that are stabilized by the exo-anomeric effect. In distorted beta-CDs, which lack hydrogen bonding, some linkages attain low-energies from the exo-anomeric effect and acyl stacking. Those linkages result in left-handed helical geometry so other linkages are forced by the CD macrocycle to have counter-balancing right-handed character. Permethylated gamma-CDs have two 'flipping' linkages as do some larger native CDs. Flipping linkages allow two left-handed segments to join into a macrocycle, thus avoiding the higher-energy, right-handed forms. Some glucose rings in derivatized beta-CDs have substantial positive twists of the pseudo torsion angle O1-C1...C4-O4, adding right-handed character to balance the left-handed linkages. In substituted gamma-CD, all residues have negative twists, giving extra left-handed character to the short, pseudo-helical segments. In non-macrocyclic molecules the twists ranged from -14 degrees to +2 degrees , averaging -6.1 degrees. In these beta- and gamma-CDs, the twists ranged from -22 degrees to +16 degrees for (4)C(1) rings, and the (O)S(2) ring in acetylated beta-CD has a twist of +34 degrees . Glucose residues in other CDs were less twisted.

Molecular Mechanics (MM4) Study of Fluorinated Hydrocarbons

A molecular mechanics study of small saturated hydrocarbons (up to C-6) substituted by up to six fluorines has been carried out with the MM4 force field. A parameter set has been developed for use in the calculation of bond lengths, bond angles, torsion angles, conformational energies, barriers to rotation, dipole moments, moments of inertia, and vibrational frequencies for these compounds. The results are mostly in fair to good agreement with experiment and ab initio calculations. The high electronegativity of fluorine leads to serious geometric consequences in these compounds, but these consequences can be dealt with adequately by suitable cross-terms in the force constant matrix, and by recognizing that some of the reference bond lengths and angles (l(0), theta(0)) and the corresponding stretching and bending constant parameters (k(s), k(theta)) that are usually thought of as constants must in fact be treated as functions of the electronegativity of the substituents. Additionally, the heavy mass of the fluorine (relative to the mass of hydrogen in alkanes) leads to large values for other cross-terms that were found to be unimportant in hydrocarbons. Conformational equilibria for polyfluorinated compounds are affected by the delta-two effect well-known in carbohydrates. A few larger fluorinated and polyfluorinated alkanes, including perfluoropropane, perfluorobutane, and Teflon, have also been studied.

The external-anomeric torsional effect

The rotational barrier for a methyl group at the end of an anomeric system is sometimes lower than we might have anticipated. Thus, in the trans-trans conformation of dimethoxymethane, the barrier to methyl rotation is calculated (B3LYP/6-311++G(2d,2p)) to be 2.22 kcal/mol, just slightly smaller than the corresponding barrier to rotation of the methyl group in methyl propyl ether of 2.32 kcal/mol. However, if the methyl being rotated in dimethoxymethane is placed into a gauche conformation, that rotational barrier is reduced to 1.52 kcal/mol. This substantial (0.80 kcal/mol relative to methyl propyl ether) reduction in barrier height in the latter case is attributed mainly to the change in the bond order of the C-O bond to which the methyl is attached, as a function of conformation, which in turn is a result of the anomeric effect. We have called this barrier lowering the external-anomeric torsional effect. This effect is apparently widespread in carbohydrates, and it results in the changing of conformational energies by up to about 2 kcal/mol. If polysaccharide potential surfaces are to be accurately mapped by molecular mechanics, this effect clearly needs to be accounted for.

Empirical force fields for biological macromolecules: overview and issues

Empirical force field-based studies of biological macromolecules are becoming a common tool for investigating their structure-activity relationships at an atomic level of detail. Such studies facilitate interpretation of experimental data and allow for information not readily accessible to experimental methods to be obtained. A large part of the success of empirical force field-based methods is the quality of the force fields combined with the algorithmic advances that allow for more accurate reproduction of experimental observables. Presented is an overview of the issues associated with the development and application of empirical force fields to biomolecular systems. This is followed by a summary of the force fields commonly applied to the different classes of biomolecules; proteins, nucleic acids, lipids, and carbohydrates. In addition, issues associated with computational studies on "heterogeneous" biomolecular systems and the transferability of force fields to a wide range of organic molecules of pharmacological interest are discussed.

Alcohols, ethers, carbohydrates, and related compounds. IV. Carbohydrates

Ab initio calculations [B3LYP/6-311++G(2d,2p)] have been carried out on 84 conformations of 12 different sugars (hexoses), in both pyranose and furanose forms, with the idea of generating a data base for carbohydrate structural energies that may be used for developing the predictive value of molecular mechanics calculations for carbohydrates. The average value for the apparent gas phase anomeric effect for a series of 31 pairs of pyranose conformations was found to be 1.83 kcal/mol (vs. 2.67 kcal/mol with a smaller basis set used in earlier calculations). In developing MM4 to reproduce these data, it was necessary first to have good energies for simple alcohols and ethers, together with an adequate treatment of hydrogen bonding, and then to include the anomeric effect, and the ethylene glycol type system, as was previously recognized. It was also found that the so-called delta-2 effect, long recognized in carbohydrates, must be explicitly included, in order to obtain acceptable results. When a force field that included all of these items as developed from the small molecules based on the MM4 hydrocarbon force field was applied without any parameter adjustment to the set of hexopyranose and furanose conformations mentioned earlier, the E(beta) - E(alpha) was found to have an average value of 1.88 kcal/mol, versus 1.74 for the quantum calculations. The signed average and RMS deviations of the MM4 from the QM results were +0.15 and 0.87 kcal/mol.