A semi-empirical conformational analysis of the interaction of n-alkanols with dipalmitoylphosphatidylcholine

The interaction of a homologous series of saturated aliphatic n-alkanols (containing 1–13 carbon atoms) with dipalmitoylphosphatidylcholine was studied by a semi-empirical conformational analysis. The minimal conformational energy of the isolated molecule at the hydrocarbon-water interface was calculated as the sum of the contributions resulting from the Van der Waals, torsional, electrostatic and transfer energies. From the conformers of minimal energies were calculated the hydrophilic-hydrophobic balance, the distance between hydrophilic and hydrophobic centres and the energies of interaction between homologous alkanols and between alkanols and lipids. Using these parameters, different modes of conformation, orientation and interaction of n-alkanols and dipalmitoylphosphatidylcholine were described. For methanol, ethanol and n-propanol, the gauche conformers were the most probable interacting only with the lipid polar heads. Only all-trans conformers were obtained for alkanols with longer acyl chains. n-Butanol to n-octanol form clusters in the lipid matrix. Longer n-alkanols are randomly distributed in the lipid layer. However, due to the increase in hydrocarbon chain-length, n-nonanol and higher alkanols have an interaction energy identical or superior to that found in a pure lipid monolayer, leading to a more ordered alkanol-lipid organization.     

Determination of the populations and structures of multiple conformers in an ensemble from NMR data: Multiple-copy refinement of nucleic acid structures using floating weights

A new algorithm is presented for determination of structural conformers and their populations based on NMR data. Restrained Metropolis Monte Carlo simulations or restrained energy minimizations are performed for several copies of a molecule simultaneously. The calculations are restrained with dipolar relaxation rates derived from measured NOE intensities via complete relaxation matrix analysis. The novel feature of the algorithm is that the weights of individual conformers are determined in every refinement step, by the quadratic programming algorithm, in such a way that the restraint energy is minimized. Its design ensures that the calculated populations of the individual conformers are based only on experimental restraints. Presence of internally inconsistent restraints is the driving force for determination of distinct multiple conformers. The method is applied to various simulated test systems. Conformational calculations on nucleic acids are carried out using generalized helical parameters with the program DNAminiCarlo. From different mixtures of A- and B-DNA, minor fractions as low as 10% could be determined with restrained energy minimization. For B-DNA with three local conformers (C2-endo, O4-exo, C3-endo), the minor O4-exo conformer could not be reliably determined using NOE data typically measured for DNA. The other two conformers, C2-endo and C3-endo, could be reproduced by Metropolis Monte Carlo simulated annealing. The behavior of the algorithm in various situations is analyzed, and a number of refinement protocols are discussed. Prior to application of this algorithm to each experimental system, it is suggested that the presence of internal inconsistencies in experimental data be ascertained. In addition, because the performance of the algorithm depends on the type of conformers involved and experimental data available, it is advisable to carry out test calculations with simulated data modeling each experimental system studied.     

The stereochemistry of N-methyl and aryl substituents determine the biological activities of 3-aryl-8-methyl-8-azabicyclo[3.2.1]oct-2,3-enes

Aryl substituted tropanes and their 2,3-ene analogs are highly selective inhibitors of monoamine uptake. The solution structures of a series of aryl tropanes were determined using NMR spectroscopy and molecular modeling to identify conformational preferences that may determine the overall activity. The majority of these analogs undergo nitrogen inversion, and the rate of interconversion between the axial and equatorial N-methyl conformers is fast on the NMR timescale at room temperature but slow between 217 and 243 K allowing us to determine the thermodynamic parameters of interconversion using dynamic and magnetization transfer NMR. The biological activities correlate strongly with the nature and the orientation of the aryl group. The relative orientation of the N-methyl further modulates the activity by directly influencing the ligand interaction in the protein binding pocket and/or by forcing a favorable orientation for the aryl substituent to fit in the binding pocket.     

Conformational study of trinucleoside tetraphosphate d(pCpGpCp): Transition of right-handed form to left-handed form

Using the semiempirical potential functions, conformational energies of the model compounds DMP^?, d(pCp), d(pGp), and d(pCpGpCp) are calculated, and the B → Z transition is discussed along the pseudorotational path of the sugar ring. For dimethylmonophosphate anion, DMP^?, the energy contour map is presented and the stabilities of the phosphodiester conformations discussed. For the sugar ring without the base attached, the minimum energies for each sugar-puckering form are calculated along the pseudorotational path. The energy barrier of the interconversion between the C(3′)-endo form and the C(2′)-endo form is calculated to be about 2.0 kcal/mol. From the conformational energy calculations of the interconversions of mononucleoside diphosphates, d(pCp) and d(pGp), between the C(2′)-endo conformer and the C(3′)-endo conformer, the purine sugar segment is known to be more convertible than the pyrimidine sugar segment. The results also support the finding that the pseudorotational transition occurred with the O(1′)-endo form more easily than with the O(1′)-exo form. Based on the results of conformational studies of DMP^?, d(pCp), and d(pGp), a topological transition of the handedness of the model compound, d(pCpGpCp), is studied. The left-handed Z-form is found to be less stable by about 8.5 kcal/mol than is the right-handed B-form. The energy barrier of the Z → B transition is calculated to be about 17.4 kcal/mol. The contributions of the electrostatic and nonbonded energies to the energy barrier are discussed in connection with the conformation changes of the model compound, d(pCpGpCp).     

Valence orbital response to conformers of n-butane

Individual outer valence orbital responses to rotations of the C–C central bond of butane (C₄H₁₀) are explored on the torsional potential energy surface. Orbital ionization energies, topologies and momentum distributions for the four most significant butane conformation are presented, as snapshots of the conformational variations. The analysis is based on quantum mechanically generated information from coordinate space and momentum space, a technique called dual space analysis (DSA). By comparison with experimental measurements of photo-electron spectra (PES) for energies and of electron momentum spectra (EMS) for energies and Dyson orbitals, we demonstrate that the individual outer valence orbitals of these conformers response differently to the rotations of the central C–C bond of n-butane. Orbital signatures of other higher energy conformations, such as orbitals la ₂ and 5a ₁ of conformation D (C_2v ), are identified. This finding indicates a co-existence of butane conformations, although the global minimum structure of anti-butane, A (C_2h ), is dominant. Orbital topology and electron charges redistribution during the transformation provide useful information on the chemical bonding and related chemical reactions.     

The origin of the generalized anomeric effect: possibility of CH/n and CH/π hydrogen bonds

Ab initio MO calculations were carried out at the MP4/6-311++G(3df,3pd)//MP2/6-311++G(3df,3pd) level to investigate the conformational Gibbs energy of a series of methyl ethers CH₃O-CH₂-X (X = OH, OCH₃, F, Cl, Br, CN, CCH, C₆H₅, CHO). It was found that the Gibbs energy of the gauche conformers is lower in every case than that of the corresponding anti conformers. In the more stable gauche conformers, the interatomic distance between X and the hydrogen atom was shorter than the sum of the van der Waals radii. The natural bonding orbital (NBO) charges of group X were more negative in the gauche conformers than in the anti conformers. We suggest that the CH/n and CH/π hydrogen bonds play an important role in stabilizing the gauche conformation of these compounds.     

Quantum-chemical and empirical calculations on phospholipids IV. Glycero-phosphorylethanolamine in a quasihexagonal planar lattice

The energetically most stable head group conformations of a racemic mixture of diacyl-glycero-phosphorylethanolamine in a planar quasihexagonal lattice were calculated using empirical 1-6-12 atom-atom potential functions for intra- and intermolecular interactions. The results demonstrate that the conformation of phospholipid head groups in bilayer systems is determined by intramolecular interactions as well as by intermolecular interactions with neighbouring phospholipid molecules and with solvent molecules. The most stable conformers are that with a φ2 = guache^? conformation of the phosphodiester group. All conformers with a φ2 = gauche⁽⁺⁾ or trans conformation have energies more than 15 kcal ☆ mol^?1 above that of the global minimum. The calculated torsional angles ?1 and φ1 are in very good agreement with the results of the X-ray diffraction analysis of 1,2-dilauroyl-DL-phosphatidylethanolamine (DLPE) acetic acid single crystals.     

Potentiometric estimation of charges in barnacle muscle fibers under internal perfusion

Summary A statistical mechanical treatment of the fluidity of lipid hydrocarbon chains in phospholipid bilayers is presented, which explicitly takes some account of interchain steric restrictions. With an effective energy separation of 750 cal/mole betweengauche andtrans conformations, it is found possible to account both for the chain dependence of the entropy and enthalpy change at the liquid crystalline transition of saturated lecithins, and also for intensity data in the laser raman spectra of dipalmitoyl lecithin. The method is used to calculate conformational probabilities in the lipid chains, in particular those for 2g 1 kinks. The calculated kink concentrations are found to be in agreement with the molecular permeability theory of H. Träuble (J Membrane Biol. 4:193, 1971).     

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.     

Conformational-energy calculations for oligosaccharides: a comparison of methods and a strategy of calculation

A theoretical conformational analysis of dimethoxymethane, 2-methoxytetrahydropyran, cellobiose, and maltose has been performed. The validity of several commonly used classical approaches to conformational energy, assuming non-bonded interactions, torsional terms, and the exo-anomeric contribution, and the MM2CARB method (a modified version of the MM2 force-field program using the Jeffrey-Taylor parameters) was tested against available experimental data or previous quantum-chemical calculations. The MM2CARB method correctly reproduces the energies and the variations in bond lengths and bond angles for conformers of dimethoxymethane and 2-methoxytetrahydropyran. Prediction of the observed conformers with simple potential functions appears to be less satisfactory. In particular, calculations that take into account non-bonded interactions and the exo-anomeric contribution based on dimethoxymethane give predicted energy differences that are 2–3 times higher than the experimental values. The general shapes of the (Φ, Ψ) potential-energy surfaces for cellobiose and maltose provided by potential-function calculations suggest the presence of several minima whose energies depend, to a great extent, on the choice of molecular geometry. The MM2CARB-calculated structures of seven cellobiose and five maltose conformers demonstrate clearly the variation of disaccharide geometry with change of conformation around the glycosidic linkage. The relative energies calculated by simple methods differ from MM2CARB energies and indicate that the simple potential-functions methods give only a qualitative estimate of oligosaccharide conformers. Based on these results, we propose a general strategy and two different approaches for the investigation of conformational properties of oligosaccharides.     

Conformational characteristics of the trinucleoside diphosphate d(ApApA) from energy-minimization studies

Energy-minimization studies were carried out on the trinucleoside diphosphate d(ApApA). The potential energy contributions from nonbonded, electrostatic, hydrogen-bonding, and torsional interactions were minimized by treating the 13 relevant dihedral angles as simultaneous variables. For the C(3′)-endo trimer, 14 low-energy conformations are within 10 kcal/mol above the lowest energy found, compared to only 3 in the case of the C(2′)-endo trimer. This result shows the flexible character of the C(3′)-endo unit. The hairpin-type, loop-promoting conformer with (ω′,ω) = (101°, 59°) was found to be the most favored structure at the 3′-terminus of d(ApApA). The predicted U- and L-type bend conformers were found to lie within 5 kcal/mol, compared to the lowest energy B-DNA structure. The A-DNA and Watson-Crick DNA types of helical conformers also lie within very small energy barriers. The phosphate group at the 5′-end of the nucleotide residue has a definite influence on the base of the corresponding nucleotide, keeping it in the normal anti-region, and hence on the base-stacking property. The results are compared with the presently available experimental data, mainly with the tRNA^Phe crystal.     

Mass Spectrometry Study and Infrared Spectroscopy of the Complex Between Camphor and the Two Enantiomers of Protonated Alanine: The Role of Higher‐Energy Conformers in the Enantioselectivity of the Dissociation Rate Constants

The properties of the protonated complexes built from S camphor and R or S alanine were studied in a Paul ion trap at room temperature by collision‐induced dissociation (CID) and infrared multiple‐photon dissociation spectroscopy (IRMPD), as well as molecular dynamics and ab initio calculations. While the two diastereomer complexes display very similar vibrational spectra in the fingerprint region, in line with similar structures, and almost identical calculated binding energies, their collision‐induced dissociation rates are different. Comparison of the IRMPD results to computed spectra shows that the SS and SR complexes both contain protonated alanine strongly hydrogen‐bonded to the keto group of camphor. The floppiness of this structure around the NH⁺…O = C hydrogen bond results in a complex potential energy surface showing multiple minima. Calculating the dissociation rate constant within the frame of the transition state theory shows that the fragmentation rate larger for the heterochiral SR complex than the homochiral SS complex can be explained in terms of two almost isoenergetic low‐energy conformers in the latter that are not present for the former. Chirality 25:436‐443, 2013. © 2013 Wiley Periodicals, Inc.     

Conformational characteristics of the RNA subunit ApA from energy minimization studies

In an attempt to better our understanding of the conformational stabilities in RNAs, an intensive theoraticl study has been carried out on one of its dimeric subunits, ApA, using an improved set of atom-atom interaction energy parameters and an improved version of energy-minimization technique. The C(3′)0endo and the C(2′)-endo sugar ApA units were sperately considered and 38 probable conformations have been analyzed in each case. The total potential energy, comprising nonbonded, electrostatic, and torsional contributions, was minimized by varying all seven relevant dihedral angles simumtaneously. The result reveal that 17 conformations in the case of C(3′)-endo sugar ApA and 7 confomations in the case of C(2′)-endo sugar ApA unit, the lowest energy conformation corresponds to a nonhelical structure and the A-RNA and the Watson-Crick-yype conformations lie at energy levels of about 0.5 and 1.0 Kcal/mo., respectively, above the lowest energy found. For ApA with the lops of different types in the backbone and they all differ in energies by about 3.5 Kcal/mol with refrence to the lowest energy founs. It is noted that the order ofmprefrence of the base stacking is observed in the A-RNA and the Watson-Crick type conformers. The ApA unit with C(2′)-endo sugar is forced to assume phosphodiester conformations with large deviations fom the expected staggered conformations compared to the ApA unit with C(3′)-endo sugar. The result obtained for ApA are discussed with refrence to those previously obtained for the dApdA unit. Te theoretical predictions are compared with the experimental data on the tRNA^Phe crystal, as well as those on fibrous RNAs and RNa subunitlike crystal structures. This study brings out many important aspects of the conformational stability of ApA which have been missed by studies made by others on this system.     

Conformational Changes of 3,5,3′-Triiodo <Emphasis Type="SmallCaps">L</Emphasis>-Thyronine Induced by Interactions with Phospholipid: Physiological Speculations

The conformational changes of 3,5,3′-triiodo L-thyronine induced by interaction with phospholipids were analyzed by Raman spectroscopy. The spectra were interpreted in terms of two conformers of this hormone in equilibrium in the lipid medium, depending on the orientation of the 3′-iodine with respect to the ring α. Theoretical geometry optimizations on both conformers in vacuo and in different solvents, together with the respective calculated energies support the experimental results. The presence of only one iodine atom in the phenolic ring allows assumption of a higher flexibility of 3,5,3′-triiodo L-thyronine and a better accommodation into the lipid medium compared to 3,5,3′,5′-tetraiodo L-thyronine. The possible physiological implications of structural differences that appear in membrane models between 3,5,3′-triiodo L-thyronine and 3,5,3′,5′-tetraiodo L-thyronine are discussed.     

VCD study of α‐methylbenzyl amine derivatives: Detection of the unchanged chiral motif

Chiral α‐methylbenzyl amine is a well known and often used chiral auxiliary, e.g., in the resolution of racemates or asymmetric catalysis. In this work, α‐methylbenzyl amine and its derivatives N,α‐dimethylbenzyl amine, N,N,α‐trimethylbenzyl amine, and bis[α‐methylbenzyl] amine were investigated by vibrational circular dichroism (VCD) spectroscopy and density functional theory (DFT). For all compounds, stable low energy conformers were obtained by the DFT calculations and based on those, the theoretical vibrational absorption (VA) and VCD spectra were calculated and compared with experimental spectra. Hence, the absolute configurations and conformational preferences were determined. A qualitative comparison of all the experimental VCD spectra of the investigated chiral molecules supported by the calculated ones is given which clearly shows similarities between the spectra of the different chiral amines. These can be assigned to vibrations of the unchanged chiral center. Chirality 2010. © 2010 Wiley‐Liss, Inc.     

Configuration-dependent properties of randomly coiling polynucleotide chains. II. The role of the phosphodiester linkage

The dependence of the unperturbed dimensions of randomly coiling polynucleotides on the rotations about the phosphodiester linkages of the chain has been examined in order to understand the conformational discrepancies, set forth in paper I, regarding these angles (ω′ and ω). Large values of the characteristic ratio 〈r²〉₀/nl² , which agree with the experimental behavior of the chain, are obtained only if a sizeable proportion of the polymer residues have trans ω′ values. The asymmetric torsional potential that is believed to arise from gauche effects associated with the P-O bonds has been approximated using a hard core model. The calculated characteristic ratio exhibits a strong dependence upon the magnitude of this torsional barrier (separating trans and gauche conformations) and shows agreement with experimental values for polyribonucleotides only if this energy difference is 1 kcal/mol or less.     

Solvent‐dependent inversion of circular dichroism signal in naproxen: An unusual effect!

The electronic circular dichroism (ECD) spectra of naproxen enantiomers were studied as a function of solvents using experimental (circular dichroism) and theoretical (time‐dependent density functional theory) approaches. The (R)‐ and (S)‐naproxen enantiomers presented an unusual inversion in their ECD signals in the presence of ethanol and water when compared with polar aprotic solvents such as acetonitrile. From a practical point of view, these findings deserve great attention because these solvents are widely used for high‐performance liquid chromatography analysis in quality control of chiral pharmaceutical drugs. This is particularly relevant to naproxen because the (S)‐naproxen has anti‐inflammatory properties, whereas (R)‐naproxen is hepatotoxic. A time‐dependent density functional theory computer simulation was conducted to investigate the signal inversion using the solvation model based on density, a reparameterization of polarized continuum model. Electronic circular dichroism signals of conformers were calculated by computer simulation and their contribution to the combined spectra obtained according to Boltzmann weighting. It was found that the experimentally observed ECD signal inversion can be associated with the minor or major contribution of different conformers of naproxen.     

Salt- and sequence-dependence of the secondary structure of DNA in Solution by 31P-nmr spectroscopy

We examined three sonicated, specific-seqiemce polydeoxynucleotides in solution over a wide range of concentrations of several salts by ¹³P-nmr spectroscopy, and we found that the alternating copolymer poly(dAdT)·poly(dAdT) exhibits a dinucleotide repeat unit in all five salts and at all concentrations studied, as indicated by the presence of a doubled in its ³¹P-nmr spectra. The two components of the doublet show selective shift effects. The upfield component is assigned to dApdT in the gauche^?-gauche^? conformation and shifts upfield in all four monovalent salts used, relative to a single-stranded oligonucleotide control. The downfield component is assigned to dTpdA in the trans-gauche^? conformation and shifts downfield with increasing CsF concentration but remains essentially constant in LiCl, NaCl, and CsCl. These changes indicate a fast noncooperative transition for poly(dAdT)·poly-(dAdT) from a presumed right-handed dinucleotide-repeat B-form to another conformation with a dinucleotide-repeat structure, via a continuum of structures that may differ in the extent of the winding of the double helix. Ethanol causes the upfield component to collapse into the other component, indicating conversion to a structure with a mononucleotide repeat unit and a trans-gauche^? conformation. Up to 1M Mg²⁺ appears to have no significant effect on the phosphodiester conformations of poly(dAdT)·poly(dAdT). By contrast, poly-(dGdC)·poly(dGdC) gives a slow cooperative transition from what is considered to be a right-handed regular B-form to a left-handed Z-form on increasing MgCl₂ and NaCl concentrations, although we observed no changes in chemical shifts below the transition points. The homopolymer poly(dA)·poly(dT) exhibits no unusual shift effects or transitions upon the addition of salts when compared to the oligonucleotide control and is considered to be a regular B-form with a gauche^?-gauche^? phosphodiester backbone conformation. These differences emphasize the distinct secondary structures of DNAs of different sequences and their selective responses to changes in solution conditions.     

Linear free energy relationship for the anomeric effect: MP2, DFT and ab initio study of 2-substituted-1,4-dioxanes

The anomeric effect of 2-substituted 1,4-dioxane derivatives was calculated and compared with the values for substituted cyclohexane. The bond lengths, bond angles, torsion angles, and relative energies of axial and equatorial conformers of 2-substituted 1,4-dioxanes were calculated by the second-order Møller–Plesset (MP2), density functional theory (DFT/B3LYP), and Hartree–Fock (HF) methods using 6-31G^∗ basis set. The energy differences between the axial and equatorial conformers, endo and exo-anomeric effects, repulsive non-bond and H-bonding interactions were investigated. A linear free energy relationship (LFER) between calculated (MP2/6-31G^∗) anomeric effect and inductive substituent constants (σ_I) was obtained for 2-substituted-1,4-dioxanes (slope = 6.19 and r² = 0.967). The calculated energy differences indicate lower equatorial orientation for 2-substituted-1,4-dioxanes compared to the 2-substituted-tetrahydropyrans. The contribution of resonance, hyperconjugation, inductive, steric, hydrogen bonding, electrostatic interaction, and level of theory influences the anomeric effect.     

Parallel nucleic acid helices with hoogsteen base pairing: Symmetry and structure

Molecular structures for parallel DNA and RNA double helices with Hoogsteen pairing are proposed for the first time. The DNA helices have sugars in the C2′-endo region and the phosphodiester conformations are (trans, gauche^?), and the RNA helices have sugars in the C3′-endo region and the phosphodiester conformations are (gauche^?, gauche^?). A pseudorotational symmetry relates the two parallel strands of DNA helices and a screw symmetry relates the two strands of RNA helices, which have an associated tilt of the The conformational space of parallel helices with Hoogsteen base pairing, unlike the Watson-Crick duplex, is highly restricted due to the unique positioning of the symmetry axis in the former case. The features of the parallel double helix with Hoogsteen pairing are compared with the Watson-Crick duplex and the corresponding triple helix. © 1994 John Wiley & Sons, Inc.