Scaled quantum mechanical force field for alanyl-alanine peptide in solution

We have obtained ab initio scale factors and assigned frequencies for the alanine-alanine peptide in water. Calculations were performed on the isolated acidic and basic Ala-Ala structures, two one-water basic Ala-Ala supermolecules, and one two-water acidic and one two-water basic Ala-Ala supermolecules. Force constants were scaled using the experimentally determined Raman and Fourier transform infrared vibrational frequencies of four isotopic species of Ala-Ala in water at pH 13 and pH 1. Most of the 4-31G scale factors were transferable from smaller molecules. All but one scale factor were directly transferable between the pH 1 and pH 13 species for coordinates unchanged by protonation in both the isolated and two-water supermolecule structures. Scale factors for nonpolar coordinates were transferable between all Ala-Ala species with only a few small changes. Good agreement was obtained between the calculated and experimental frequencies for all isotopic species and structures.     

Infrared and Raman spectra and vibrational analyses calculated with Moeller–Plesset perturbation theory of second order of nitrosoethylene and its chloro-derivatives

The conformational and structural stabilities of nitrosoethylene CH₂=CH–N=O, chloronitrosoethylene CH₂=CCl–N=O, and Dichloronitrosoethylene CCl₂=CH–N=O were investigated by ab initio Moeller–Plesset perturbation theory of second order (MP2) calculations using the 6−311+G^** basis set to include electron correlation. From the calculations all three were predicted to exist predominantly in the planar trans structure (C=C and N=O bonds are trans to each other) with high trans-cis rotational barriers of about 9 kcal mol⁻¹ as a result of pronounced conjugation between C=C and N=O bonds. The vibrational frequencies were computed for the three molecules, and also the d ₁ and d ₂ deuterated variants for the parent molecule at the MP2 level. Normal coordinate analyses were carried out and the potential energy distributions (PED), among the symmetry coordinates of the normal modes of the molecule were computed. Complete vibrational assignments were made on the basis of normal coordinate analyses for the molecules. The two chlorinated derivatives of nitrosoethylene were also investigated in the same way. As expected, we then find high Raman and infrared intensities in all modes that contain a high content of chlorine movements because vibrations of C–Cl bonds lead to large changes in polarizability, as well as to a large change in dipole moment. However, modes involving double bonds also have quite large intensities. An appreciable number of modes in these molecules are more or less pure symmetry coordinates.     

Hydration of cis and trans N-methylformamide as revealed by the use of 17O-NMR,molecular mechanics,and ab initio calculations

The solvation of cis and trans N-methylformamide (NMF) by water was investigated using a combination of ¹⁷O-nmr spectroscopy, classical molecular mechanics [MM2(77) and MM2(87)] force field, and ab initio 4-31G* gradient optimization calculations. In dilute aqueous solution, the ¹⁷O-nmr spectra of NMF indicate strong shielding by 66.9 and 66.1 ppm for the cis and trans amide oxygens, respectively, compared to those values obtained in dilute toluene solution. This demonstrates that both isomers are equally solvated by molecules of water, which are further hydrogen bonded to molecules of water of the bulk solvent. Molecular mechanics simulations were carried out for cis and trans NMF in a cluster of water molecules. Radial distribution functions show structural contacts by several water molecules at the amide CO and NH group, which are significantly more pronounced with MM2 (87) calculations. Ab initio 4-31G* gradient optimization calculations on the supermolecule trans NMF-(H₂O)₃ indicates the presence of more than two hydrogen-bond contacts at the carbonyl oxygen. This is in agreement with MM2 calculations and provides further evidence for multiple acceptor properties of the amide oxygen and an out of the amide plane arrangement of the bound molecules of water. Comparison of the integration data to the first radial distribution function (rdf) minima shows that the local solvation of the CO and NH groups is very similar for both cis and trans isomers. The intermolecular geometric parameters of the supermolecule trans NMA–(H₂O)₃ and the first rdf maxima resulting from MM2 (87) and MM2 (77) calculations are compared with distribution of water molecules around the CO and NH groups of peptides and proteins resulting from x-ray and neutron diffraction experiments. The rdfs involving the methyl group of NMF demonstrate the nonrandom distribution of solvent sites with first maxima in reasonable agreement with distribution of water molecules around the apolar side chain of amino acid residues in proteins. © 1995 John Wiley & Sons, Inc.     

Monte Carlo Simulation of Hydration of the Nucleic Acid Fragments

Abstract

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of “water economy”. Hydration shell patterns of NH₂-groups of AU and GC helical stacks differ significantly.     

Conformational analysis of perezone and dihydroperezone using vibrational circular dichroism

The vibrational circular dichroism (VCD) spectra of perezone and dihydroperezone measured from CDCl₃ solutions were quite similar, suggesting analogous conformations for both molecules. Their absolute configurations were confirmed by comparison of the experimental VCD spectrum of each compound with curves generated from theoretical calculations using density functional theory (DFT) at the B3LYP/DGDZVP level of theory taking into account their conformational mobility. Conformational analysis of the 8-(R) enantiomer showed 19 low energy conformers in a 2.4 kcal/mol energy range, while for 8-(R), with the saturated side alkyl chain, 34 conformers were considered in the first 2 kcal/mol. Initial analyses were carried out using a Monte Carlo searching with the MMFF94 molecular mechanics force field, all MMFF94 conformers were geometrically optimized using DFT at the B3LYP/6-31G(d) level of theory, followed by reoptimization and calculations of their vibrational frequencies at the B3LYP/DGDZVP level. Good agreement between the theoretical 8-(R) enantiomers and experimental VCD curves were observed for both.     

Structure of guanylyl-3',5'-cytidine monophosphate. II. Description of the molecular and crystal structure of the calcium derivative in space group P2(1).

The structural features of calcium guanosine-3′,5′-cytidine monophosphate (GpC) have been elucidated by X-ray diffraction analysis. The molecule was crystallized in space group P2₁ with cell constants of a = 21.224 Å, b = 34.207 Å, c = 9.327 Å, and β = 90.527°, Z = 8. The hydration of the crystal is 21% by weight with 72 water molecules in the unit cell. The four GpC molecules in the asymmetric unit occur as two Watson-Crick hydrogen-bonded dimers related by a pseudo-C face centering. Each dimer consists of two independent GpC molecules whose bases are hydrogen bonded to each other in the traditional Watson-Crick fashion. Each dimer possesses a pseudo twofold axis broken by a calcium ion and associated solvent. The four molecules are conformationally similar to helical RNA, but are not identical to it or to each other. Instead, values of conformational angles reflect the intrinsic flexibility of the molecule within the range of basic helical conformations. All eight bases are anti, sugars are all C3′-endo, and the C4′-C5′ bond rotations are gauche-gauche. The R factor is 12.6% for 2918 observed reflections at 1.2-Å resolution.     

An X-ray diffraction study of poly-L-homoarginine hydrochloride.

An X-ray study of the synthetic polypeptide poly(L -homoarginine hydrochloride) has been made to investigate whether, like the chemically related polypeptides poly(L -lysine hydrochloride), poly(L -arginine hydrochloride), and poly(L -ornithine hydrobromide), it can undergo conformational transitions merely from variations in its degree of hydration. X-ray photographs of powder and oriented specimens containing one to 15 molecules of water per L -homoarginine hydrochloride residue showed that this polymer forms only a β-pleated-sheet structure. The pleated sheets, formed by antiparallel polypeptide chains hydrogen-bonded to each other, are piled up along the b axis in an alternating sequence (“sandwich structure”). This structure did not appreciably change with variations of the degree of hydration, and the observed reflections at 56% relative humidity (1.8 molecules of water per residue) could be indexed satisfactorily in terms of a monoclinic unit cell, of space group P2₁, with a = 9.34 Å, b = 40.07 Å, c = 6.94 Å, and γ = 106°. These dimensions are shown by models to be compatible with the proposed structure, and the calculated density of 1.27 g/cm³ agrees well with the experimental value of 1.29 g/cm³. Removal of the last molecule of water results in a very diffuse pattern, while specimens containing 20 molecules of water per residue show only reflections due to water.     

Picosecond relaxations in model substances for proteins: A millimeter-wave investigation on crystalline alkyl amides

The dielectric absorption at millimeter-wave (mm-wave) frequencies (50–150 GHz) of N-methylacetamide (NMA), N,N-dimethylacetamide (DiNMA), and N,N-dimethylacrylamide (DiNMAcry) is measured. Measurements are performed using the oversized-cavity technique in the temperature range from liquid helium to room temperature. Additionally, a mm-wave interferometeric measurement at room temperature is made. NMA and DiNMAcry exhibit monotonic increases of the absorption coefficient with temperature as well as with frequency. For DiNMA a monotonic increase of the absorption coefficient with frequency is also found, while the absorption coefficient as a function of temperature shows a pronounced maximum at approximately 30 K. At this maximum the absorption coefficient of DiNMA exceeds those of NMA and DiNMAcry by about two orders of magnitude. The dielectric behavior of the three substances can be described by relaxation processes in asymmetric double-well potentials. For the low-temperature relaxation in DiNMA the double well could be established by two possible positions of the molecule in the crystal that are separated by a rotational movement. Hydrogen bonds and long side chains may hinder these relaxational movements in NMA and DiNMAcry, respectively, and thereby account for their comparatively lower absorption. The results are compared with similar results recently obtained on proteins and synthetic biopolymers.     

Ab Initio Calculations of Vibrational Frequencies,Potential Functions of Internal Rotations and Vibrational Infrared and Raman Spectra for 3,3,3-Trifluoropropanal

The conformational behavior and structure of 3,3,3,-trifluoropropal have been investigated by utilizing ab initio calculations with the 6-31G^** basis set (valence double zeta basis with polarization functions on all atoms) at the restricted Hartree Fock (RHF), second-order Møller-Plesset perturbation (MP2), and Density Functional (B3LYP) levels. The molecule is predicted to have a cis Û gauche conformational equilibrium. Full optimization of the transition states was performed and the rotational barriers of both the CHO and CF₃ rotors were calculated. Vibrational frequencies were computed at the three levels and the zero-point corrections were included into the calculated asymmetric CHO rotational barrier. Complete vibrational assignments were made on the basis of normal coordinate calculations for both stable conformers of the molecule.     

Chiroptical properties of 2,2’‐bioxirane

The two enantiomers of 2,2′‐bioxirane were synthesized, and their chiroptical properties were thoroughly investigated in various solvents by polarimetry, vibrational circular dichroism (VCD), and Raman optical activity (ROA). Density functional theory (DFT) calculations at the B3LYP/aug‐cc‐pVTZ level revealed the presence of three conformers (G₊, G_?, and cis) with Gibbs populations of 51, 44, and 5% for the isolated molecule, respectively. The population ratios of the two main conformers were modified for solvents exhibiting higher dielectric constants (G_? form decreases whereas G₊ form increases). The behavior of the specific optical rotation values with the different solvents was correctly reproduced by time‐dependent DFT calculations using the polarizable continuum model (PCM), except for the benzene for which explicit solvent model should be necessary. Finally, VCD and ROA spectra were perfectly reproduced by the DFT/PCM calculations for the Boltzmann‐averaged G₊ and G_? conformers.     

Zwitterionic conformers of pyrrolysine and their interactions with metal ions—a theoretical study

A total of 16 pyrrolysine conformers in their zwitterionic forms are studied in gas and simulated aqueous phase using a polarizable continuum model (PCM). These conformers are selected on the basis of our study on the intrinsic conformational properties of non-ionic pyrrolysine molecule in gas phase [Das and Mandal (2013) J Mol Model 19:1695?1704]. In aqueous phase, the stable zwitterionic pyrrolysine conformers are characterized by full geometry optimization and vibrational frequency calculations using B3LYP/6-311++G(d,p) level of theory. Single point calculations are also carried out at MP2/6-311++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. The calculated relative energy range of the conformers at B3LYP/6-311++G(d,p) level is 5.19 kcal mol^?1 whereas the same obtained by single point calculations at MP2/6-311++G(d,p) level is 4.58 kcal mol^?1. A thorough analysis reveals that four types of intramolecular H-bonds are present in the conformers; all of which play key roles in determining the energetics and in imparting the observed conformations to the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of the H-bonds. This study also points out that conformers with diverse structural motifs may differ in their thermodynamical stability by a narrow range of relative energy. The effects of metal coordination on the relative stability order and structural features of the conformers are examined by complexing five zwitterionic conformers of pyrrolysine with Cu⁺² through their carboxylate groups. The interaction enthalpies and Gibbs energies, rotational constants, vibrational frequencies and dipole moments of the metal complexes calculated at B3LYP level are also reported. The zwitterionic conformers of pyrrolysine are not stable in gas phase; after geometry optimization they are converted to the non-ionic forms.     

Effects of the position and manner of hydration on the stability of solvated N-methylacetamides and the strength of binding between N-methylacetamide and water clusters: a computational study

This work mainly studies the effects of the position (there are two possible hydrated sites) and the manner (i.e., whether water acts as a proton donor or acceptor) of hydration by various numbers of water molecules on the stability of 14 solvated N-methylacetamide structures, NMA-(H₂O) _n (n = 1–3), as well as the binding strength between the NMA and the water cluster, using molecular dynamics (MD) and B3LYP methods. Natural bond orbital (NBO) analysis is used to explore the origin of these effects. Some novel observations are obtained from the work. Our results show that monohydration at the carbonyl site favors stability and binding strength compared to monohydration at the amino site. Similarly, the preferred hydration at the carbonyl site is observed for dihydrated NMAs when the second water is added as a proton donor to the C=O group or the first water is H-bonded to the C=O group. However, unfavorable hydration at the C=O site occurs if the second water acts as a proton acceptor. Trihydration by a ring cluster of three water molecules at either the carbonyl site or the amino one yields relatively stable complexes, but significantly disfavors binding strength. The other trihydrated NMAs show similar behavior to dihydrated NMAs. In addition, our results show that the C=O and N–H frequencies can still be utilized to examine the H-bond effects of the water cluster.     

On the possible influence of different pathways of binding of amides with water on “Pyramidalization” of valence bonds of peptide group nitrogen

With the aim to study solvation effects in peptide structure organization, the behavior of the energy of different types of hydration in simple amines and amides has been analyzed. On the example of quantum-chemical DFT and PM3 calculations of amino derivatives of composition CH₃-(CH₂)₃)-NH₂, (CH3)₂-NH, CH₃-NH₂, NH₃, CH₂=CH-NH₂, H-CC-NH₂, O=C(CH)₃-N(CH₃)₂, O=C(CH₃)-NH(CH₃), O=C(CH₃)-NH₂, O=CH-N(CH₃)H, and O=CH-NH₂ it has been shown that: (1) in the given set of molecules, the proton acceptor N…H-O variant of hydrogen bonding of NH₂ group with a water molecule is dominating only for the simplest amines. Being primordially weaker, the proton donor N-H…OH variant of water H-bonding gradually increases in energy in the given set as the basicity of the compound decreases, and for the case of amides of carboxylic acids it becomes already a significant channel of the hydration; (2) the intermolecular N-H…O=C bonding of trans-N-methylacetamides, which models the peptide hydrogen bonds in proteins, induces “planarization” of its initially nonplanar O=C-NH fragments. However, the addition of water molecules to the complex through the proton acceptor N…H-O variant of binding of N atom not only restores but even strengthens the “pyramidalization” of valence bonds of peptide groups.     

2-(Phenylethyl)ammonium hydrogensquarate hemihydrate: crystal structure,solid-state IR-spectroscopic and theoretical characterization

The title compound, 2-(phenylethyl)ammonium hydrogensquarate hemihydrate, was synthesized and structurally and spectroscopically characterized by a single crystal X-ray diffraction and solid-state polarized IR spectroscopy of oriented colloids in a nematic host. The crystal structure consists of two crystallographically independent 2-(phenylethyl)ammonium cations, joined in a 2D hydrogen-bonded network with hydrogensquarate anions and solvent water molecules. Surprisingly, the crystallographically non-equivalent cations exhibit differing pseudo T and G trans configurations.     

The crystal and molecular structure of 1,6-anhydro-β-d-mannopyranose

The crystal structure of 1,6-anhydro-β-d-mannopyranose, C₆H₁₀O₅, is orthorhombic, P2₁2₁2₁, with a = 10.971(2), b = 13.935(3), c = 9.012(1) Å, V = 1377.76 »³ (MoKα, λ = 0.7107 Å), Z = 8, D_x = 1.563 M.gm⁻³, D_m = 1.565 M.gm⁻³. the structure was solved by MULTAN and refined to R(F) = 0.043 for 2355 reflections. The two symmetry-independent molecules in the unit cell have similar conformations, except for the orientation of one of the three hydroxyl groups. The conformation of the pyranose rings is ¹C₄ distorted towards E_o, and that of the anhydro rings is E. There are significant differences between the two molecules in two of the four C---O bond-lengths. These occur where there are important differences in the hydrogen-bonding environment of the oxygen atoms. The molecules are hydrogen-bonded by three linear and three bifurcated O---H···O interactions which form four-membered loops linked into infinite chains. Empirical force-field calculations with MMI-CARB reproduced the geometry of the molecules within the variations observed experimentally between the two molecules, except for a C---O bond in one of the molecules. The effect of excluding the anomeric effect from the theoretical calculations was not significant. Calculations for an intramolecularly hydrogen-bonded molecule were also carried out as a model for the molecules in a non-polar solvent.     

The effect of hydration on the mobility of phospholipids in the gel state. A proton nuclear magnetic resonance spin echo study

Proton nuclear magnetic resonance (NMR) dipolar echo studies are presented for the gel state of dipalmitoylglycerophosphocholine (dipalmitoyl-GPC) — heavy water dispersions. The mobility and the mean order of the chains and the head group of dipalmitoyl-GPC were determined for different water concentrations and temperatures. For smaller than 5 mol D₂O per mol dipalmitoyl-GPC the molecule undergoes temperature- and hydration-dependent restricted rotational oscillations about the long axis of the molecule. For hydration numbers equal or larger than 5 mol D₂O per mol dipalmitoyl-GPC the molecules rotate effectively about their long axes and intermolecular dipolar interactions between proton groups of neighbour molecules are averaged. The onset of the lateral diffusion of dipalmitoyl-GPC is observed which averages out all intermolecular dipolar interactions. Deviations of the individual segments of the chains from the all-trans state have to be considered. The widely accepted model that the dipalmitoyl-GPC molecules rotate about their long axes with stiff all-trans chains should be modified. The polar head groups of dipalmitoyl-GPC effectively rotate about the bilayer-normal and restricted rotations about single bonds in the head group are allowed. An order parameter of about 0.6 for the head group was obtained for fully hydrated dipalmitoyl-GPC molecules at ambient temperature.     

On the conformation of luteinizing hormone-releasing hormone, nuclear Overhauser observations.

Observations of proton nuclear Overhauser effects in the molecule Luteinizing Hormone-Releasing Hormone indicate that a high population of a particular set of conformers exists in water solution. The results can be interpreted as two distinct conformers in which the pGlu¹ ring is in close proximity to aromatic residues further along the sequence of the linear structure. The observed nuclear Overhauser effects were in agreement with the enhancements calculated from models obtained by conformational energy calculations.     

Peptides and Peptoids - A Systematic Structure Comparison

A systematic analysis of the conformational space of the basic structure unit of peptoids in comparison to the corresponding peptide unit was performed based on ab initio MO theory and complemented by molecular mechanics (MM) and molecular dynamics (MD) calculations both in the gas phase and in aqueous solution.The calculations show three minimum conformations denoted as C_7ß, a_D and a that do not correspond to conformers on the gas phase peptide potential energy hypersurface. The influence of aqueous solvation was estimated by means of continuum models. The MD simulations indicate the a_D form as the preferred conformation in solution both in cis and trans peptide bond orientations.