Molecular dynamics study of hydration in ethanol-water mixtures using a polarizable force field

The abnormal physicochemical characteristics of ethanol solvation in water are commonly attributed to the phenomenon of hydrophobic hydration. To investigate the structural organization of hydrophobic hydration in water-ethanol mixtures, we use molecular dynamics simulations based on detailed atomic models. Induced polarization is incorporated into the potential function on the basis of the classical Drude oscillator model. Water-ethanol mixtures are simulated at 11 ethanol molar fractions, from 0.05 to 0.9. Although the water and ethanol models are parametrized separately to reproduce the vaporization enthalpy, static dielectric constant, and self-diffusion constant of neat liquids at ambient conditions, they also reproduce the energetic and dynamical properties of the mixtures accurately. Furthermore, the calculated dielectric constant for the various water-alcohol mixtures is in excellent agreement with experimental data. The simulations provide a detailed structural characterization of the mixtures. A depletion of water-water hydrogen bonding in the first hydration shell of ethanol is compensated by an enhancement in the second hydration shell. The structuring effect from the second solvation shell gives rise to a net positive hydrogen-bonding excess for ethanol molar fractions up to approximately 0.5. For larger molar fractions, the second hydration shell is not sufficiently populated to overcome the net H-bond depletion from the first shell.     

Water concentration profiles in membranes measured by ESEEM of spin-labeled lipids

Electron spin-echo envelope modulation (ESEEM) spectroscopy of phospholipids spin-labeled systematically down the sn-2 chain was used to detect the penetration of water (D2O) into bilayer membranes of dipalmitoyl phosphatidylcholine with and without 50 mol % cholesterol. Three-pulse stimulated echoes allow the resolution of two superimposed 2H-ESEEM spectral components of different widths, for spin labels located in the upper part of the lipid chains. Quantum chemical calculations (DFT) and ESEEM simulations assign the broad spectral component to one or two D2O molecules that are directly hydrogen bonded to the N-O group of the spin label. Classical ESEEM simulations establish that the narrow spectral component arises from nonbonded water (D2O) molecules that are free in the hydrocarbon chain region of the bilayer membrane. The amplitudes of the broad 2H-ESEEM spectral component correlate directly with those of the narrow component for spin labels at different positions down the lipid chain, reflecting the local H-bonding equilibria. The D2O-ESEEM amplitudes decrease with position down the chain toward the bilayer center, displaying a sigmoidal dependence on position that is characteristic of transmembrane polarity profiles established by other less direct spin-labeling methods. The midpoint of the sigmoidal profile is shifted toward the membrane center for membranes without cholesterol, relative to those with cholesterol, and the D2O-ESEEM amplitude in the outer regions of the chain is greater in the presence of cholesterol than in its absence. For both membrane types, the D2O amplitude is almost vanishingly small at the bilayer center. The water-penetration profiles reverse correlate with the lipid-chain packing density, as reflected by 1H-ESEEM intensities from protons of the membrane matrix. An analysis of the H-bonding equilibria provides essential information on the binding of water molecules to H-bond acceptors within the hydrophobic interior of membranes. For membranes containing cholesterol, approximately 40% of the nitroxides in the region adjacent to the lipid headgroups are H bonded to water, of which ca. 15% are doubly H bonded. Corresponding H-bonded populations in membranes without cholesterol are ca. 20%, of which ca. 6% are doubly bonded.     

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.     

Intercalation of Cyclic Amines into α-Titanium Phosphate

The intercalation of some amines (aniline, benzylamine, cyclohexylamine,piperidine, pyridine, pyrazine and piperazine) into -titaniumphosphate, Ti(HPO₄)_2×H₂O,has been investigated by the batch method and/or by exposing the host to thevapour of the amines. The changes in the interlayer distance of the solidduring the intercalation process was followed by X-ray powder diffraction.The new intercalates were characterised by chemical and thermal analysis.Materials with a monolaminar and/or bilaminar arrangement of amine moleculesin the phosphate interlayer region are obtained depending on the nature ofthe amine. Due to steric hindrance, saturated phases are not obtained forall amines studied. The thermal decomposition of the intercalates (nitrogenatmosphere), takes place in three stages: dehydration, removal of amines andcondensation of the hydrogenphosphate to pyrophosphate groups.     

Final elucidation of the absolute configuration of the signal metabolite hormaomycin

The complete absolute configuration of hormaomycin 1 a has been established by HPLC and HPLC/MS experiments with appropriately derivatized 4-propylprolines, (2S,4S)-6 and (2R,4R)-6, as well as 4-(Z)-propenylprolines, cis-5 and trans-5, and also feeding experiments with enantiomerically pure samples of the deuterium-labeled 3-(2'-nitrocyclopropyl)alanine, (2S)-3,3-[D2]15 and (2S)-2,2'-[D2]15, and 4-(Z)-propenylproline 2',4-[D2]-(2S,4R)-5. The latter five amino acids were prepared for the first time and allowed one to unequivocally assign the hitherto unknown absolute configurations of the last four stereocenters in hormaomycin 1 a. As a bonus, some new information about the biosynthesis of this molecule has also been gathered.     

Orientation-selected 15N-HYSCORE detection of weakly coupled nitrogens around the archaeal rieske [2Fe-2S] center

The weakly coupled 15N atoms around a reduced Rieske [2Fe-2S] cluster of the uniformly 15N-labeled, hyperthermostable archaeal Rieske protein appear to produce readily observable cross-peaks in the HYSCORE spectra, with the well-resolved couplings of 0.3-0.4 MHz for the Nepsilon and 1.1 MHz for the peptide backbone nitrogens, in addition to the contributions from the coordinated Ndelta atoms. These features can be used for structure-mechanism studies of the biological redox protein system involving the weakly coupled nitrogens in coupled electron-proton transfer reactions.     

Photodissociation of heme distal methionine in ferrous cytochrome C revealed by subpicosecond time-resolved resonance Raman spectroscopy

Cytochrome c (cyt c) is an electron-transfer heme protein that also binds nitric oxide (NO). In resting cyt c, two endogenous ligands of the heme iron are histidine-18 (His) and methionine-80 (Met) side chains, and NO binding requires the cleavage of one of the axial bonds. Previous femtosecond transient absorption studies suggested the photolysis of either Fe-His or Fe-Met bonds. We aimed at unequivocally identifying the internal side chain that is photodissociated in ferrous cyt c and at monitoring heme structural dynamics, by means of time-resolved resonance Raman (TR3) spectroscopy with approximately 0.6 ps time resolution. The Fe-His stretching mode at 216 cm-1 has been observed in photoproduct TR3 spectra for the first time for a c-type heme. The same transient mode was observed for a model ferrous cyt c N-fragment (residues 1-56) ligated with two His in the resting state. Our TR3 data reveal that upon ferrous cyt c photoexcitation, (i) distal Met side chain is instantly released, producing a five-coordinated domed heme structure, (ii) proximal His side chain, coupled to the heme, exhibits distortion due to strain exerted by the protein, and (iii) alteration in heme-cysteine coupling takes place along with the relaxation of the protein-induced deformations of the heme macrocycle.     

Vibrational spectrum of BiH3: Six-dimensional variational calculations on high-level ab initio potential energy surfaces

We report a theoretical study of the ground electronic state of BiH₃. The potential energy surface (PES) is obtained from coupled cluster CCSD(T) calculations with a large basis set (289 contracted Gaussian functions). The previously available quartic force field (P4) is extended by adding the dominant quintic and sextic stretching terms to yield improved potential functions in symmetry coordinates (P6) and Morse-type coordinates (M4). Second-order rovibrational perturbation calculations on the P4-PES and full variational calculations on the P6-PES and M4-PES yield almost identical vibrational term values which is rationalized by considering the local mode behavior of BiH₃ and the Morse-type character of the M4-PES. The remaining deviations between the computed and observed vibrational term values must thus be caused by imperfections in the CCSD(T) surface. A refinement of this ab initio surface by a restrained fit to experimental data allows an essentially perfect reproduction of the observed vibrational term values. Variational calculations on this refined surface provide predictions for several overtone and combination bands that have not yet been observed. Dedicated to Hermann Stoll on the occasion of his 60th birthday     

Ring-closing olefin metathesis on ruthenium carbene complexes: model DFT study of stereochemistry

Ring-closing metathesis (RCM) is the key step in a recently reported synthesis of salicylihalamide and related model compounds. Experimentally, the stereochemistry of the resulting cycloolefin (cis/trans) depends strongly on the substituents that are present in the diene substrate. To gain insight into the factors that govern the observed stereochemistry, density functional theory (DFT) calculations have been carried out for a simplified dichloro(2-propylidene)(imidazole-2-ylidene)ruthenium catalyst I, as well as for the real catalyst II with two mesityl substituents on the imidazole ring. Four model substrates are considered, which are closely related to the systems studied experimentally, and in each case, two pathways A and B are possible since the RCM reaction can be initiated by coordination of either of the two diene double bonds to the metal center. The first metathesis yields a carbene intermediate, which can then undergo a second metathesis by ring closure, metallacycle formation, and metallacycle cleavage to give the final cycloolefin complex. According to the DFT calculations, the stereochemistry is always determined in the second metathesis reaction, but the rate-determining step may be different for different catalysts, substrates, and pathways. The ancillary N-heterocyclic carbene ligand lies in the Ru-Cl-Cl plane in the simplified catalyst I, but is perpendicular to it in the real catalyst II, and this affects the relative energies of the relevant intermediates and transition states. Likewise, the introduction of methyl substituents in the diene substrates influences these relative energies appreciably. Good agreement with the experimentally observed stereochemistry is only found when using the real catalyst II and the largest model substrates in the DFT calculations.