Molecular dynamics study of the solvation of calcium carbonate in water

We performed molecular dynamics simulations of diluted solutions of calcium carbonate in water. To this end, we combined and tested previous polarizable models. The carbonate anion forms long-living hydrogen bonds with water and shows an amphiphilic character, in which the water molecules are expelled in a region close to its C(3) symmetry axis. The calcium cation forms a strongly bound ion pair with the carbonate. The first hydration shell around the CaCO(3) pair is found to be very similar to the location of the water molecules surrounding CaCO(3) in ikaite, the hydrated mineral.     

Side-chain dynamics of an alpha-helical polypeptide monitored by fluorescence

The "blob" model, developed to analyze the fluorescence decays of polymers randomly labeled with pyrene, has been applied to a series of pyrene-labeled poly(glutamic acid)s (PyPGA) in DMF and carbonated buffer solutions at pH 9. Poly(glutamic acid) (PGA) exists in the ionized form in the buffer solutions as poly(sodium glutamate) (PGNa). PGA adopts an alpha-helical conformation in DMF, whereas in aqueous solution PGNa is a random coil. Fluorescence, UV-vis absorption, and circular dichroism measurements indicate that in our studies pyrene pendants attached themselves along PGA in a clustered manner. Simulations were carried out to establish that the geometry of the PGA alpha-helix induces the high level of pyrene clustering. Since the level of pyrene clustering decreased with lower pyrene content, information about naked PGA was retrieved by extrapolating the trends obtained by fluorescence to zero pyrene content. Analysis of the fluorescence decays demonstrated that during its lifetime an excited pyrene probes a 32 amino acid section of the PGA alpha-helix. This result was supported by molecular mechanics optimizations. This study establishes that the blob model, originally used to monitor the encounters between pyrenes attached randomly onto a polymer adopting a random coil conformation, can also be applied to study the dynamics of the side chains of structured proteins. Since the blob model helps in monitoring the encounters between amino acids in the initial state (i.e., random coil) and in the final state (i.e., structured protein) of the folding pathway of a protein, it could be applicable to the study of protein folding.     

Molecular dynamics simulation of a hydrated diphytanol phosphatidylcholine lipid bilayer containing an alpha-helical bundle of four transmembrane domains of the influenza A virus M2 protein

An alpha-helical bundle composed of four transmembrane portions of the M2 protein from the Influenza A virus has been studied in a hydrated diphytanol phosphatidylcholine bilayer using molecular dynamics (MD) calculations. Experimentally, the sequence utilized is known to aggregate as a four-helix bundle and act as a pH-gated proton-selective ion channel, which is blocked by the drug amantadine hydrochloride. In the presented simulation, the ion channel was initially set up as a parallel four-helix bundle. The all-atom simulation consisted of almost 16,000 atoms, described classically, using a forcefield from the CHARMM22 database. Bilayers with and without the bundle were shown to be stable throughout the nanosecond timescale of the MD simulation. Structural and dynamical properties of the bilayer both with and without the transmembrane protein are reported.     

The dynamics of carbene solvation: an ultrafast study of p-biphenylyltrifluoromethylcarbene

Ultrafast photolysis (lambda(ex) = 308 nm) of p-biphenylyltrifluoromethyl diazomethane (BpCN2CF3) releases singlet p-biphenylyltrifluoromethylcarbene (BpCCF3) which absorbs strongly at 385 nm in cyclohexane, immediately after the 300 fs laser pulse. The initial absorption maximum shifts to longer wavelengths in coordinating solvents (nitrile, ether, and alcohol). In low viscosity coordinating solvents, the initial absorption maximum further red shifts between 2 and 10 ps after the laser pulse. Similar effects are observed upon ultrafast photolysis of 2-fluorenyltrifluoromethyl diazomethane (FlCN2CF3) and therefore cannot be associated with torsional motion around the two phenyl rings of the biphenyl compound. Instead, the effect is attributed to the dynamics of solvation of the singlet carbene. The time constant of solvation in normal alcohols lengthens with solvent viscosity in a linear manner. Furthermore, the time constants of the red shift in methanol-O-d (16 ps), ethanol-O-d (26 ps), 2-propanol-OD (40 ps), and 2,2,2-trifluoroethanol-O-d (14 ps) are longer than those recorded in methanol (9.6 ps, KIE = 1.7), ethanol (14.3 ps, KIE = 1.8), 2-propanol (28 ps, KIE = 1.4), and 2,2,2-trifluoroethanol (4.4 ps, KIE = 3.2), which indicates that the solvent reorganization involves formation of hydrogen bonds. The kinetic data are consistent with motion of the solvent to achieve a specific interaction with the carbene, with the creation of a new hydrogen bond. The solvated carbene reacts with the solvent over tens, hundreds, and thousands of ps, depending upon the solvent.     

Molecular dynamics study of 2rotaxanes: influence of solvation and cation on co-conformation

The conformational preference of a [2]rotaxane system has been examined by molecular dynamics simulations. The rotaxane wheel consists of two bridged binding components: a cis-dibenzo-18-crown-6 ether and a 1,3-phenyldicarboxamide, and the penetrating axle consists of a central isophthaloyl unit with phenyltrityl capping groups. The influence of solvation on the co-conformation of the [2]rotaxane was evaluated by comparing the conformational flexibility in two solvents: chloroform and dimethyl sulfoxide. Attention was also paid to the effect of cation binding on the dynamical properties of the [2]rotaxane. The conformational stability of the [2]rotaxane was calculated using a MM/PB-SA strategy, and the occurrence of specific motions was examined by essential dynamics analysis. The changes in the co-conformational properties in the two solvents and upon cation binding are discussed in light of the available NMR data. The results indicate that in chloroform solution the [2]rotaxane system exists as a mixture of co-conformational states including some that have hydrogen bonds between axle C=O and wheel NH groups. Analysis of the simulations allow us to hypothesize that the [2]rotaxane's circumrotation motion can occur as the result of a dynamic process that combines a preliminary axle sliding step that breaks these hydrogen bonds and a conformational change in the ester group more distant from the wheel. In contrast, no hydrogen-bonded co-conformation was found in dimethyl sulfoxide, which appears to be due to the preferential formation of hydrogen bonds between the wheel NH groups with solvent molecules. Moreover, the axle experiences notable changes in anisotropic shielding, which would explain why the NMR signals are broadened in this solvent. Insertion of a sodium cation into the crown ether reduces co-conformational flexibility due to an interaction of the axle with the cation. Overall, the results reveal how both solvent and ionic atmosphere can influence the co-conformational preferences of rotaxanes.     

Molecular dynamics simulation study on the transient response of solvation structure during the translational diffusion of solute

The transient response function of the density profile of the solvent around a solute during the translational diffusion of the solute is formulated based on the generalized Langevin formalism. The resultant theory is applied to both neat Lennard-Jones fluids and cations in liquid water, and the response functions are obtained from the analysis of the molecular dynamics simulations. In the case of the self-diffusion of Lennard-Jones fluids, the responses of the solvation structures are in harmony with conventional pictures based on the mode-coupling theory, that is, the binary collision in the low-density fluids, the backflow effect from medium to high density fluids, and the backscatter effect in the liquids near the triple point. In the case of cations in water, the qualitative behavior is strongly dependent on the size of cations. The pictures similar to simple dense liquids are obtained for the large ion and the neutral molecule, while the solvent waters within the first solvation shell of small ions show an oscillatory response in the short-time region. In particular, the oscillation is remarkably underdumped for lithium ion. The origin of the oscillation is discussed in relation to the theoretical treatment of the translational diffusion of ions in water.     

Heterocyclic peptide backbone modifications in an alpha-helical coiled coil

In this paper, we present 1,2,3-triazole epsilon2-amino acids incorporated as a dipeptide surrogate at three positions in the sequence of a known alpha-helical coiled coil. Biophysical characterization indicates that the modified peptides retain much of the helical structure of the parent sequence, and that the thermodynamic stability of the coiled coil depends on the position of the incorporation of the epsilon-residue. Crystal structures obtained for each peptide give insight into the chemical behavior and conformational preferences of the non-natural amino acid and show that the triazole ring can participate in the backbone hydrogen bonding of the alpha-helix as well as template an interhelical crossing between chains in the bundle.     

Characterization of the solvation dynamics of an ionic liquid via molecular dynamics simulation

The solvation dynamics of ionic liquids have been the subject of intense experimental study but remain poorly understood. We present the results of molecular dynamics simulations of the solvation dynamics of the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate in response to photoexcitation of the fluorescent dye coumarin-153. We reproduce the time-resolved fluorescence Stokes shift using linear response theory, then use novel statistical techniques to analyze cation and anion contributions to the signal. We find that the solvation dynamics are dominated by collective ionic motion and characterize the time scale for various features of the collective response. Further, we use the Steele analysis [Mol. Phys. 61, 1031 (1987)] to characterize the contributions to the observed Stokes shift made by translational and rovibrational degrees of freedom. Our results indicate that in contrast to molecular liquids, the rovibrational response is trivial and the observed fluorescence response arises almost entirely from ionic translation. Our results resolve previously open questions in the literature about the nature of the rapid dynamics in room-temperature ionic liquids and offer insight into the physical principles governing ionic liquid behavior on longer time scales.     

Molecular fibers based on the honeycomb-like self-assembly of an alpha-helical polypeptide

We have succeeded in fabricating well-grown molecular fibers of a polypeptide on substrates by using a conventional solution spin-coating and drying process. These molecular fibers were found to consist of a honeycomb-like molecular assembly formed via the hexagonal close packing of the polypeptide chains in the alpha-helix conformation.     

Femtosecond study of partially folded states of cytochrome C by solvation dynamics

Using femtosecond time-resolved fluorescence spectroscopy, it is shown that the solvation dynamics in the two partially folded states (IS' and IS' ') of a protein, cytochrome C, are very different. In the case of IS' (formed by the addition of 2 mM sodium dodecyl sulfate, SDS) almost the entire dynamic solvent shift of coumarin 153 (C153) is captured in a picosecond setup and the contribution of the ultrafast component (0.5 ps) is very small (5%). Solvation dynamics of IS' ' (formed by 2 mM SDS and 5 M urea) displays a major component (47%) of 1.3 ps. This indicates that the structure of IS' ' is much more open and exposed compared to that of IS'. The difference in the dynamics of IS' and IS' ' is attributed to differences in their structure, particularly near the heme region, and the presence of urea in IS' '.     

Molecular dynamics study of the vaporization of an ionic drop

The melting of a microcrystal in vacuum and subsequent vaporization of a drop of NaCl were studied through molecular dynamics simulations with the Born-Mayer-Huggins-Tosi-Fumi rigid-ion effective potential. The vaporization was studied for a single isochor at increasing temperatures until the drop completely vaporized, and gaseous NaCl formed. Examination of the vapor composition shows that the vapor of the ionic drop and gaseous NaCl are composed of neutral species, the most abundant of which, ranging from simple NaCl monomers (ion pairs) to nonlinear polymers, (Na(n)Cl(n))(n=2-4). The enthalpies of sublimation, vaporization, and dissociation of the different vapor species are found to be in reasonable agreement with available experimental data. The decrease of the enthalpy of vaporization of the vapor species, with the radius of the drop decrease, accounts for a larger fraction of trimers and tetramers than that inferred from experiments. Further, the rhombic dimer is significantly more abundant than its linear isomer although the latter increases with the temperature. The present results suggest that both trimers and linear dimers may be important to explain the vapor pressure of molten NaCl at temperatures above 1500 K.     

Mechanisms of the interaction of alpha-helical transmembrane peptides with phospholipid bilayers

The synthetic peptide acetyl-K(2)-G-L(24)-K(2)-A-amide (P(24)) and its analogs have been successfully utilized as models of the hydrophobic transmembrane alpha-helical segments of integral membrane proteins. The central polyleucine region of these peptides was designed to form a maximally stable, very hydrophobic alpha-helix which will partition strongly into the hydrophobic environment of the lipid bilayer core, while the dilysine caps were designed to anchor the ends of these peptides to the polar surface of the lipid bilayer and to inhibit the lateral aggregation of these peptides. Moreover, the normally positively charged N-terminus and the negatively charged C-terminus have both been blocked in order to provide a symmetrical tetracationic peptide, which will more faithfully mimic the transbilayer region of natural membrane proteins and preclude favorable electrostatic interactions. In fact, P(24) adopts a very stable alpha-helical conformation and transbilayer orientation in lipid model membranes. The results of our recent studies of the interaction of this family of alpha-helical transmembrane peptides with phospholipid bilayers are summarized here.     

Molecular dynamics simulations of peptide carboxylate hydration

Aqueous solvation of carboxylate groups, as present in the glycine zwitterion and the dipeptide aspartylalanine, is studied employing a force-field that includes distributed multipole electrostatics and induction contributions (Amoebapro: P. Ren and J. W. Ponder, J. Comput. Chem., 2002, 23, 1497; P. Ren and J. W. Ponder, J. Phys. Chem. B, 2003, 107, 5933; J. W. Ponder and D. A. Case, Adv. Protein Chem., 2003, 66, 27). Radial and orientation distribution functions, as well as hydration numbers, are calculated and compared with existing simulation data derived from Car-Parrinello molecular dynamics (CPMD), and also distributed-charge force-fields. Connections are also made with experimental data for solvation of carboxylates in water. Our findings show that Amoebapro yields carboxylate solvation properties in very good agreement with CPMD results, significantly closer agreement than can be obtained from traditional force-fields. We also demonstrate that the influence of solvation on the conformation of the dipeptide is markedly different using Amoebapro compared with the other force-fields.     

Glyco-helix: parallel lactose-lactose interactions stabilize an alpha-helical structure of multi-glycosylated peptide

Glyco-helix is designed as a novel model system to investigate cis carbohydrate-carbohydrate interactions. Adhesive Lac-Lac interactions stabilize alpha-helix of Lac-peptide in the presence of fluorinated alcohols, but no such an interaction was observed in GlcNAc-peptide.     

Femtosecond infrared study of the dynamics of solvation and solvent caging.

The ultrafast reaction dynamics following 295-nm photodissociation of Re2CO10 were studied experimentally with 300-fs time resolution in the reactive, strongly coordinating CCl4 solution and in the inert, weakly coordinating hexane solution. Density-functional theoretical (DFT) and ab initio calculations were used to further characterize the transient intermediates seen in the experiments. It was found that the quantum yield of the Re-Re bond dissociation is governed by geminate recombination on two time scales in CCl4, approximately 50 and approximately 500 ps. The recombination dynamics are discussed in terms of solvent caging in which the geminate Re(CO)5 pair has a low probability to escape the first solvent shell in the first few picoseconds after femtosecond photolysis. The other photofragmentation channel resulted in the equatorially solvated dirhenium nonacarbonyl eq-Re2(CO)9(solvent). Theoretical calculations indicated that a structural reorganization energy cost on the order of 6-7 kcal/mol might be required for the unsolvated nonacarbonyl to coordinate to a solvent molecule. These results suggest that for Re(CO)5 the solvent can be treated as a viscous continuum, whereas for the Re2(CO)9 the solvent is best described in molecular terms.     

Capillary electrophoresis of amphipathic alpha-helical peptide diastereomers

We have made a rigorous assessment of the ability of capillary electrophoresis to resolve peptide diastereomers through its application to the separation of a series of synthetic 18-residue, amphipathic alpha-helical monomeric peptide analogues, where a single site in the centre of the hydrophobic face of the alpha-helix is substituted by 19 L- or D-amino acids. Such L- and D-peptide pairs have the same mass-to-charge ratio, amino acid sequence and intrinsic hydrophobicity, varying only in the stereochemistry of one residue. CE approaches assessed in their ability to separate diastereomeric peptide pairs included capillary zone electrophoresis (uncoated capillary), micellar electrokinetic chromatography (uncoated capillary in the presence of 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate, CHAPS), open-tubular capillary electrochromatography (C(8)-coated capillary in the presence of 25% 2,2,2-trifluoroethanol (TFE) or 25% ethanol). Overall, the OT-CEC methods were the most effective at separating the most peptide pairs, particularly for those containing hydrophilic side chains. However, the MEKC approach proved most effective for separation of peptide pairs containing hydrophobic or aromatic side chains.     

A comparative study of solvation dynamics in room-temperature ionic liquids

The solvation dynamics of ionic liquids have been the subject of many experimental and theoretical studies but remain poorly understood. We analyze these dynamics by modeling the time-resolved fluorescence response of coumarin 153 in two room-temperature ionic liquids: 1-butyl-1-methylpyrrolidinium bromide and 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide. Our results demonstrate that phenomena such as electrostatic screening operate significantly differently in the two liquids, and the relative importance of translational and rovibrational components of the ionic response depends significantly on the character of the ions involved. However, collective motion dominates the response of both ionic liquids, and the qualitative features of this collective behavior are strikingly similar in both cases.     

Molecular dynamics study of the beta amyloid peptide of Alzheimer's disease and its divalent copper complexes

The Abeta1-42 monomer structure was assessed with a 790 ns molecular dynamics (MD) simulation, and the results were compared with the NMR experiment on Abeta10-35 and Abeta1-40. Previous theoretical work in a model of the His13-His14 region of Abeta defined the possible Cu(II) binding geometries at this site (Raffa et al. J. Biol. Inorg. Chem. 2005, 10, 887-902). MD simulations totalling almost 2 micros were also carried out on Cu(II)/Abeta1-42 systems, using the ab initio structures as templates for the copper binding site. This work finds that the copper-free Abeta1-42 system may stabilize after approximately 350 ns into a collapsed coil conformation, and we find good agreement with some, but not all, of the structural features determined experimentally for the Abeta10-35 and Abeta1-40 peptides. The results of the Cu(II)/Abeta1-42 systems are compared to the Cu(II)-free Abeta1-42 simulation.     

The study of the interaction of a model alpha-helical peptide with lipid bilayers and monolayers

We studied the interaction of the alpha-helical peptide acetyl-Lys(2)-Leu(24)-Lys(2)-amide (L(24)) with tethered bilayer lipid membranes (tBLM) and lipid monolayers formed at an air-water interface. The interaction of L(24) with tBLM resulted in adsorption of the peptide to the surface of the bilayer, characterized by a binding constant K(c)=2.4+/-0.6 microM(-1). The peptide L(24) an induced decrease of the elasticity modulus of the tBLM in a direction perpendicular to the membrane surface, E(radial). The decrease of E(radial) with increasing peptide concentration can be connected with a disordering effect of the peptide to the tBLM structure. The pure peptide formed a stable monolayer at the air/water interface. The pressure-area isotherms were characterized by a transition of the peptide monolayer, which probably corresponds of the partial intercalation of the alpha-helixes at higher surface pressure. Interaction of the peptide molecules with lipid monolayers resulted in an increase of the mean molecular area of phospholipids both in the gel and liquid crystalline states. With increasing peptide concentration, the temperature of the phase transition of the monolayer shifted toward lower temperatures. The analysis showed that the peptide-lipid monolayer is not an ideally miscible system and that the peptide molecules form aggregates in the monolayer.     

Molecular probes of solvation phenomena

The properties of the molecules present in any chemical or biological system are dependent on interactions with the environment, and a quantitative understanding of solvation phenomena remains a major challenge. Molecular recognition probes provide a new approach to quantitatively measure the properties of solvents. Traditionally, solvent polarity scales have been based on spectroscopic probes that provide insight into the nature of solvent-solute interactions. This review compares the solvent polarity parameters obtained from the wavelengths of UV/Visible absorption maxima with solute H-bond parameters obtained from the free energies of solution equilibria. The similarity of the solvent and solute H-bond scales leads to a general H-bond scale that uses the same parameters to describe both solvent and solute. The general H-bond scale provides a framework for understanding the relationship between local intermolecular interactions and the properties of the bulk medium. Intermolecular interactions are sensitive to solvation equilibria, so molecular recognition probes provide fundamentally different information from spectroscopic probes that are sensitive to the populations of different solvation states of the solute. Studies of mixed solvents demonstrate the potential of molecular recognition probes for providing new insights into solvation phenomena.