Free Enthalpy Difference Research Articles

Structural characterization and keto-enol tautomerization of 4-substituted pyrazolone derivatives with DFT approach

The synthesis of two pyrazolone derivative compounds, PYR-I(4-Acetyl-1-(4-chlorophenyl)-3-isopropyl-1H-pyrazol-5(4H)-one) and PYR-II1-(4-Chlorophenyl))-3-isopropyl-5-oxo-4,5-5-dihydro-1H-pyrazole-4-carbaldehyde, their characterization by FT-IR, NMR, UV–Vis and GC-MS techniques, and the evaluation of the keto-enol tautomerization process of the structures along with the DFT approach and spectral data were reported in this paper. Spectral findings indicated that PYR-I was stable at the keto state. The IR spectrum recorded in solid form showed that the PYR-II structure was stable in the enol state, while the NMR spectrum in the solution medium showed that it was stable in the keto state. DFT-based analyses were realized with the B3LYP hybrid functional and the 6–311++G(d,p) basis set. The modelled keto, transition and enol state molecular geometries of structures were optimized in the gas phase and different solvent media and the total energy and dipole moment values were investigated at the specified theoretical level. The possible keto-enol tautomerism mechanism of the structures was evaluated through some thermodynamic parameters such as the difference in free Gibbs energy (ΔG), enthalpy (ΔH), entropy (ΔS), and predictive tautomeric equilibrium constants (Keq), acidity constants (pKa) and percentages of tautomers at 298.15 K and 1 atm pressure. The results of these analyses based on the DFT approach indicated that the keto-enol tautomer equilibrium heavily favours the keto form for PYR-I and the enol form for PYR-II in all cases. Moreover, natural bond orbital (NBO) analysis was performed for the tautomers, and the chemical reactivity profiles of the most stable tautomers were examined with the values of frontier molecular orbital energy and some reactivity descriptors.

Long-range magnetic order in CePdAl3 enabled by orthorhombic deformation

We investigate the effect of structural deformation on the magnetic properties of orthorhombic CePdAl3 in relation to its tetragonal polymorph. Utilizing x-ray and neutron diffraction, we establish that the crystal structure has the Cmcm space-group symmetry and exhibits pseudotetragonal twinning. According to density functional calculations, the tetragonal-orthorhombic deformation mechanism has its grounds in the relatively small free enthalpy difference between the polymorphs, allowing either phase to be quenched, and fully accounts for the twinned microstructure of the orthorhombic phase. Neutron diffraction measurements show that orthorhombic CePdAl3 establishes long-range magnetic order below TN=5.29(5) K characterized by a collinear, antiferromagnetic arrangement of magnetic moments. Magnetic anisotropies of orthorhombic CePdAl3 arise from strong spin-orbit coupling as evidenced by the crystal-field splitting of the 4f multiplet, fully characterised with neutron spectroscopy. We discuss the potential mechanism of frustration posed by antiferromagnetic interactions between nearest neighbors in the tetragonal phase, which hinders the formation of long-range magnetic order in tetragonal CePdAl3. We propose that orthorhombic deformation releases the frustration and allows for long-range magnetic order. Published by the American Physical Society 2024

Study of influence of anomeric effect on conformational preferences using hybrid density functional theory (DFT) and natural bond orbital (NBO) interpretation

For the relationship between thermodynamic parameters, global hardness, global electronegativity, anomeric effect, structural parameters, and structural properties of dodecahedro-spiro [isochromene-3,2'-pyran] (A), dodecahedro-spiro [isothiocromene-3,2'- thiopyran) (B), dodecahedro-spiro [isochroman-3,2'-thiopyran] (C) and dodecahedro-spiro [isothiochrome-3,2'-pyran] (D) computational methods (B3LYP / 6-311 ++ G **) were used. Results show that the calculated Gibbs free energy and enthalpy differences between the axial- and equatorial- stereoisomers [i.e. ΔG = Geq – Gax, ΔH = Heq- Hax] decrease from compound A to compound B and from compound C to compound D. NBO results show that the anomeric effect associated with the electron delocalization is in benefit of the axial-stereoisomers and with the decrease of the anomeric effect, the values of global hardness (η), and global electronegativity (χ), from compound A to compound B and from compound C to Composition D is also decrease. The correlations between the global hardness, global electronegativity, anomeric effect, structural parameters, and thermodynamic parameters [ΔH, ΔG, and ΔS] of compounds A-D have been investigated. Also, it can be concluded that the stability energy obtained from the numerical values of the anomeric effect and the values calculated theoretically (B3LYP / 6-311 ++ G **) confirmed each other.

First-Principle Study on the Stability of Lightly Doped (Nb1-x Ti x )C Complex Carbides and Their Verification in 1045 Steel.

The mixing Gibbs free energy and formation enthalpy difference of different Ti-doped (Nb1–xTix)C complex carbides were calculated using the Cambridge Serials Total Energy Package (CASTEP) module of Materials Studio 2019 software. The calculation results predict that (Nb1–xTix)C complex carbides have higher stability than pure NbC and TiC. Therefore, three lightly Ti-doped (Nb1–xTix)C complex carbides with theoretical densities close to that of the 1045 steel were designed for calculations. The calculation results show that the formation energy of (Nb1–xTix)C complex carbides decreases with an increase in the Ti content. These designed (Nb1–xTix)C complex carbides have mechanical stability, and their bulk modulus, shear modulus, Young’s modulus, and hardness are all lower than those of pure NbC. The electronic performance results show that these three structures show good conductivity, and the 3d orbitals of Ti atoms and the 4d orbitals of Nb atoms are strongly hybridized with the 2p orbitals of C atoms. The Nb–C and Ti–C bonds exhibit strong covalent bonds. To verify the stability of the (Nb1–xTix)C complex carbides, the prepared (Nb0.8Ti0.2)C complex carbide was added to the 1045 steel as a refiner. After observing under a transmission electron microscope (TEM), we found that the (Nb0.8Ti0.2)C complex carbide could exist stably as a face-centered cubic structure, which provided a method for the design and synthesis of complex carbides used for refiners.

Development of Efficient Synthetic Method for Tautomeric Dihydropyrimidines and Analysis of Their Functionality

I here present the results of our studies on the synthesis and functional analysis of tautomeric dihydropyrimidines (DPs) and related compounds in two sections. In the first section, we describe our experimental and theoretical studies on the thermodynamics and properties of 2-substituted 1,4- and 1,6-dihydropyrimidine-5-carboxylates by 1H NMR measurements and density functional theory (DFT) calculations, respectively. The concentration ratios of tautomers a/b of DPs 1, 2, and 3 were determined under various conditions to determine the effects of temperature, solvent, and concentration on thermodynamics data. The obtained free energy differences (ΔG), enthalpy differences (ΔH), and entropy differences (ΔS) are discussed in terms of the molecular structures, dipole moments (DM), and electrostatic potential maps obtained by DFT calculations to clarify the nature of DPs 4-8. In the second section, an efficient synthetic method developed for 6-unsubstituted 3,4-dihydropyrimidin-2(1H)-thiones 9 and 2-ones 10 is described. The novelties of the synthesis protocol are as follows: 1) the use of Lewis acid-mediated reaction, 2) good to high yields, and 3) its broad scope of applicability to aldehydes and ureas. Hitherto unavailable 6-unsubstituted 2-amino DP 11 and 2-aryl DP 12 were obtained from 9 by a substitution reaction with the amine and the Liebeskind-Srogl reaction, respectively. The compounds 9, 10, and related 6-methyl derivatives 19-21 were assessed for their antiproliferative effect on the human promyelocytic leukemia cell line HL-60. 4,4-Dipropyl derivative 20 showed relatively strong activity with an IC50 value of 341 nM.

Interface equilibrium modeling of all-solid-state lithium-ion thin film batteries

Dynamic modeling of all-solid-state batteries (SSBs) has gained strong research interest in recent years. However, the intrinsic equilibrium state, important for the exchange current calculation and the interface resistance analysis is hardly addressed. For the first time, we introduce an advanced MPNP-FBV model which takes the electrical double layer (EDL) structure and the unoccupied regular lattice sites (vacancies) into account to illustrate the fundamental phenomena at the equilibrium state. This general model is verified from the chemical perspective and the interface experimental results in our previous investigations, which indicates that the model is able to calculate the exchange current and the charge transfer resistance from fundamental material properties. The influence on concentration and electrostatic potential by the different EDL structures, the state of charge (SOC) and the diffusivity are also investigated. Results show that the total electrostatic potential drop at the equilibrium state is related to the free enthalpy difference between the two materials. Furthermore, the charge transfer resistance with the diffuse double layer structure is higher than that with the compact double layer.

Experimental and theoretical studies on thermodynamics and properties of tautomers of 2-substituted 6(4)-methyl-1,4(1,6)-dihydropyrimidine-5-carboxylates

Experimental and theoretical studies on the thermodynamics and properties of 2-substituted 6(4)-methyl-1,4(1,6)-dihydropyrimidine-5-carboxylates were undertaken by 1H NMR measurements and DFT (density functional theory) calculations. The ratios of tautomers a/b of dihydropyrimidines (DPs) 1, 2, and 3 were determined under various conditions to reveal the effects of temperature, solvent, and concentration on the thermodynamic data. The observed results, the free energy differences (ΔG), enthalpy differences (ΔH), and entropy differences (ΔS), are discussed in terms of the molecular structures, dipole moments (DM), and the electrostatic potential maps calculated by the DFT to clarify the nature of the DPs.

Estimation of thermodynamic parameters for Au- and Mg-based metallic glasses

The study of temperature dependent thermodynamic parameters; Gibb’s free energy difference (ΔG), entropy difference (ΔS) and enthalpy difference (ΔH) between the undercooled liquid and the corresponding equilibrium solid phases has been proved to be extremely advantageous in the study of the thermodynamic behaviour of Metallic glass (MG) forming melts. In last two decades, Au- and Mg-based alloys were found to form glass phases. In present study, the three thermodynamic parameters viz., ΔG, ΔS and ΔH are calculated theoretically in the entire temperature range Tm (melting temperature) to Tg (glass transition temperature) for both Au- and Mg-based five samples of MGs; Au77Ge13.6Si9.4, Au53.2Pb27.5Sb19.3, Au81.4Si18.6, Mg85.5Cu14.5 and Mg81.6Ga18.4 on the basis of Taylor’s series expansion. A relative study is also made between the present result and the result obtained experimentally as well as on the basis of expressions projected by the earlier researchers. An attempt is also been made to narrate the reduced glass transition temperature with glass forming ability for all five MGs.

Analysis of thermodynamic properties for Pd-based bulk metallic glasses

The analysis of temperature dependent thermodynamic parameters; Gibb's free energy difference (ΔG), entropy difference (ΔS) and enthalpy difference (ΔH) between the undercooled liquid and the corresponding equilibrium solid phases has been proved to be highly beneficial in the study of the thermodynamic behaviour of bulk metallic glass (BMG) forming melts. In last two decades, Pd based alloys were found to form bulk glass phases. In present study, the three thermodynamic parameters viz., ΔG, ΔS and ΔH are calculated theoretically in the entire temperature range Tm (melting temperature) to Tg (glass transition temperature) for Pd-based three samples of BMGs; Pd40Ni40P20, Pd40Cu30Ni10P20 and Pd43Cu27Ni10P20 on the basis of Taylor's series expansion. A comparative study is also made between the present result and the result obtained experimentally as well as on the basis of expressions proposed by the earlier researchers. An attempt is also been made to relate the reduced glass transition temperature with glass forming ability for all three BMGs.

Thermodynamical properties of Zr-based bulk metallic glasses

The temperature dependence of Gibb's free energy difference (ΔG), entropy difference (ΔS) and enthalpy difference (ΔH) between the undercooled melt and the corresponding equilibrium solid phases of bulk metallic glass (BMG) forming melts has been proved to be very useful in the study of their thermodynamical behavior. The present study is made by calculating ΔG, ΔS and ΔH in the entire temperature range Tm (melting temperature) to Tg (glass transition temperature) for three Zr-based samples of BMGs: Zr57Cu15.4Ni12.6Al10Nb5, Zr41.2Ti13.8Ni10Cu12.5Be22.5 and Zr58.5Cu15.6Ni12.8Al10.3Nb2.8. The study is made on the basis of Taylor's series expansion and a comparative study is also performed between the present result and the result obtained in the framework of expansions proposed by earlier workers, and also with the experimental results. An attempt has also been made to study the glass forming ability for BMGs.

Enhanced conformational sampling using enveloping distribution sampling

To lessen the problem of insufficient conformational sampling in biomolecular simulations is still a major challenge in computational biochemistry. In this article, an application of the method of enveloping distribution sampling (EDS) is proposed that addresses this challenge and its sampling efficiency is demonstrated in simulations of a hexa-β-peptide whose conformational equilibrium encompasses two different helical folds, i.e., a right-handed 2.7(10∕12)-helix and a left-handed 3(14)-helix, separated by a high energy barrier. Standard MD simulations of this peptide using the GROMOS 53A6 force field did not reach convergence of the free enthalpy difference between the two helices even after 500 ns of simulation time. The use of soft-core non-bonded interactions in the centre of the peptide did enhance the number of transitions between the helices, but at the same time led to neglect of relevant helical configurations. In the simulations of a two-state EDS reference Hamiltonian that envelops both the physical peptide and the soft-core peptide, sampling of the conformational space of the physical peptide ensures that physically relevant conformations can be visited, and sampling of the conformational space of the soft-core peptide helps to enhance the transitions between the two helices. The EDS simulations sampled many more transitions between the two helices and showed much faster convergence of the relative free enthalpy of the two helices compared with the standard MD simulations with only a slightly larger computational effort to determine optimized EDS parameters. Combined with various methods to smoothen the potential energy surface, the proposed EDS application will be a powerful technique to enhance the sampling efficiency in biomolecular simulations.

Refinement of the application of the GROMOS 54A7 force field to β-peptides

In this study, a hexa-β-peptide whose conformational equilibrium encompasses two different helical folds, a right-handed 2.7(10/12)-helix and a left-handed 3(14)-helix, is simulated using different GROMOS force-field parameter sets. When applying the recently developed GROMOS 54A7 force field, a significant destabilization effect on the 2.7(10/12)-helix of the peptide is observed, and the agreement with the experimental NOE distance bounds is much worse compared with the ones using previous versions of the GROMOS force field. This led us to investigate the free enthalpy difference between the two helices as a function of a variation of different subsets of force-field parameters. Both long time molecular dynamics simulations and one-step perturbation predictions suggest that the disagreement with the experimental NMR data when using the 54A7 force field is caused by the use for β-peptides of the new backbone φ-/ψ-torsional-angle energy terms introduced in this force field which were based on conformational fitting of backbone φ/ψ angles for a large set of proteins. This means that these parameters of backbone φ- and ψ-torsional-angle terms should not be applied to non-α-peptides such as β-peptides. This modified assignment of torsional-angle energy terms and parameters is denoted as 54A7_β. It corrects the wrong description of the conformational ensemble of the hexa-β-peptide obtained using the previous assignment and yields as good agreement with NMR data for other β-peptides that adopt a single helical or a hairpin fold.

Intramolecular Charge Transfer with Crystal Violet Lactone in Acetonitrile As a Function of Temperature: Reaction Is Not Solvent-Controlled

Intramolecular charge transfer (ICT) with crystal violet lactone (CVL) in the excited singlet state takes place in solvents more polar than n-hexane, such as ethyl acetate, tetrahydrofuran, and acetonitrile (MeCN). In these solvents, the fluorescence spectrum of CVL consists of two emission bands, from a locally excited (LE) and an ICT state. The dominant deactivation channel of the lowest excited singlet state is internal conversion, as the quantum yields of fluorescence (0.007) and intersystem crossing (0.015) in MeCN at 25 °C are very small. CVL is a weakly coupled electron donor/acceptor (D/A) molecule, similar to an exciplex (1)(A(-)D(+)). A solvatochromic treatment of the LE and ICT emission maxima results in the dipole moments μe(LE) = 17 D and μe(ICT) = 33 D, much larger than those previously reported. This discrepancy is attributed to different Onsager radii and spectral fluorimeter calibration. The LE and ICT fluorescence decays of CVL in MeCN are double exponential. As determined by global analysis, the LE and ICT decays at 25 °C have the times τ2 = 9.2 ps and τ1 = 1180 ps, with an amplitude ratio of 35.3 for LE. From these parameters, the rate constants ka = 106 × 10(9) s(-1) and kd = 3.0 × 10(9) s(-1) of the forward and backward reaction in the LE ⇄ ICT equilibrium are calculated, resulting in a free enthalpy difference ΔG of -8.9 kJ/mol. The amplitude ratio of the ICT fluorescence decay equals -1.0, which signifies that the ICT state is not prepared by light absorption in the S0 ground state, but originates exclusively from the directly excited LE precursor. From the temperature dependence of the fluorescence decays of CVL in MeCN (-45 to 75 °C), activation energies E(a) = 3.9 kJ/mol (LE → ICT) and E(d) = 23.6 kJ/mol (ICT → LE) are obtained, giving an enthalpy difference ΔH (= E(a) - E(d)) of -19.7 kJ/mol, and an entropy difference ΔS = -35.5 J mol(-1) K(-1). These data show that the ICT reaction of CVL in MeCN is not barrierless. The ICT reaction time of 9.2 ps is much longer than the mean solvent relaxation time of MeCN (0.26 ps), indicating, in contrast with earlier reports in the literature, that the reaction is not solvent controlled. This conclusion is supported by the observation of double exponential LE and ICT fluorescence with the same decay times.

Free Energy Profile and Mechanism of Self-Assembly of Peptide Amphiphiles Based on a Collective Assembly Coordinate

By combining targeted molecular dynamics (TMD) simulations, umbrella sampling, and the weighted histogram analysis method (WHAM), we have calculated the potential of mean force (PMF) for the transition between the bound and free states of 90 peptide amphiphiles (PAs) in aqueous solution, with the bound state corresponding to a cylindrical micelle fiber. We specifically consider a collective reaction coordinate, the radius of gyration of the PAs, to describe assembly in this work. It is found that the free energy, enthalpy, and entropy differences between the free and bound states are -126 kcal/mol, -185 kcal/mol, and -190 cal/(mol K), respectively, for the self-assembly process. This indicates that the driving force to form the micelle structure is enthalpic. The enthalpic driving forces originate from several factors, including the conformational energy of PAs and the electrostatic and van der Waals interaction energy between solvent molecules and between solvent and PAs. Among these interactions, the solvent electrostatic interaction is the dominating one, contributing 54% of the total driving force. The PMF profile can be recognized as involving two stages of assembly: (1) PAs initially approach each other in mostly random configurations and loosely aggregate, resulting in significant desolvation and initiation of head-tail conformational reorganization; (2) PAs undergo a conformational disorder-to-order transition, including forming secondary structures that include more β-sheets and fewer random coils, along with tail-head core-shell alignment and condensation that leads to total exclusion of water from the core. The PMF decreases slowly in the first stage, but rapidly in the second. This study demonstrates a hierarchy of assembly steps in which PA structural changes, solvation, and redistribution of solvent molecules play significant roles in the PA self-assembly process.

Ab Initio Calculations of the Conformational Preferences of 1,3-Oxathiane S-Oxide and its Analogs Containing S and SE Atoms—Evidence for Stereoelectronic Interactions Associated with the Anomeric Effects

Stereoelectronic interactions associated with the AE and also the conformational and structural properties of 1,3-oxathiane S-oxide (1), 1,3-dithiane S-oxide (2), 1,3-thiaselenane S-oxide (3), 1,3-oxaselenane Se-oxide (4), 1,3-thiaselenane Se-oxide (5), and 1,3-diselenane Se-oxide (6) were investigated using quantum mechanical methods. These compounds were fully optimized at the B3LYP/6 – 311 + G** and HF/6 – 311 + G** levels of theory. The Gibbs Free Energy, Enthalpy, and Entropy differences (i.e., ΔG, ΔH, and ΔS) between the axial and equatorial conformations were calculated at the B3LYP/6 – 311 + G** and HF/6 – 311 + G** levels of theory. The decrease of the AE is corresponding to the decreases of calculated Δ(G ax – G eq) value of the above mentioned compounds. The calculated AE values are more significant for the justification of the conformational dominances of the compounds than steric effects. In this work, the relations between the Anomeric Effects, donor and acceptor orbital energies, occupancies, structural parameters, dipole–dipole interactions, and conformational behavior of the compounds have been studied.

Relative Free Enthalpies for Point Mutations in Two Proteins with Highly Similar Sequences but Different Folds

Enveloping distribution sampling was used to calculate free-enthalpy changes associated with single amino acid mutations for a pair of proteins, GA95 and GB95, that show 95% sequence identity yet fold into topologically different structures. Of the L → A, I → F, and L → Y mutations at positions 20, 30, and 45, respectively, of the 56-residue sequence, the first and the last contribute the most to the free-enthalpy difference between the native and non-native sequence-structure combinations, in agreement with the experimental findings for this protein pair. The individual free-enthalpy changes are almost sequence-independent in the four-strand/one-helix structure, the stable form of GB95, while in the three-helix bundle structure, the stable form of GA95, an interplay between residues 20 and 45 is observed.

The effect of branched side chains on the relative stability of α- and π-helices: a combination of the enveloping distribution sampling and one-step perturbation methods†

The method of enveloping distribution sampling (EDS) to efficiently obtain free enthalpy differences between different molecular systems from a single simulation can also be used to compute free enthalpy differences between two different conformations of a system. A combination (EDS–OSP) of EDS and the one-step perturbation (OSP) method that yields many free enthalpy differences from a single simulation allows for rapid prediction of conformational free enthalpy differences for many not too different molecular systems from a single simulation. Here, we applied EDS–OSP to predict the free enthalpy differences between a π-helix and an α-helix for a set of four deca-peptides with as fifth residue Ala, Val, Leu, or Ile in aqueous solution. First, an EDS simulation of a designed soft-core reference-state Hamiltonian was carried out to sample both helices in a single simulation. Then the soft-core atoms are perturbed into real atoms of each of the four peptides. Thus, the free enthalpy differences between the π-helix and the α-helix for all four perturbed-state, real peptides can be predicted with only one simulation using EDS–OSP. EDS–OSP and the original EDS method gave very similar free enthalpy values, i.e., the average absolute deviation between the two methods for the four peptides is 0.8 kJ mol−1, while EDS–OSP required almost four times less computational effort. Side chains branched at the C β -position were found to slightly decrease the stability of the α-helical conformation with respect to the π-helical one.

Multistate reweighting and configuration mapping together accelerate the efficiency of thermodynamic calculations as a function of molecular geometry by orders of magnitude

We present an approach to calculate free energy and other thermodynamic property differences between molecules which have very little or no overlap in configuration space, but where a one-to-one mapping between the molecule geometries exists. The approach combines multistate reweighting with remapping of phase space between simulated states. We apply this method to calculate the free energy differences between non-overlapping, truncated harmonic oscillators, the free energy, enthalpy, and entropy differences between different parameterizations of rigid water, and differences in free energy of solvation between dipoles of different lengths. Previously difficult or impossible problems become either trivially easy or are improved in efficiency by two to five orders of magnitude.

Steered Molecular Dynamics Studies of the Potential of Mean Force for Peptide Amphiphile Self-Assembly into Cylindrical Nanofibers

Steered molecular dynamics (SMD) simulations were applied to determine the potential of mean force for the self-assembly of peptide amphiphile (PA) nanofibers, specifically considering a single PA adding to a growing cylindrical nanofiber at 310 K. It is found that the free energy, enthalpy, and entropy differences for this assembly process are -67 kcal/mol, -71.5 kcal/ml, and -14.5 cal/(mol K), respectively, and therefore that enthalpy provides the driving force for self-assembly to form a fiber. A pairwise interaction analysis shows that both electrostatic and van der Waals interactions play important roles in the self-assembly process, with the van der Waals interaction being the larger effect. The mechanistic picture that emerges from this work is that as the PA is pulled from the fiber, the interaction evolves through three stages: (1) initially electrostatic interactions between the charged head of the pulled PA and other PAs, and between the pulled PA and solvent are dominant, (2) after the charged head emerges, the rest of the peptide comes out, with both PA-solvent electrostatic interactions and van der Waals interactions being significant, and (3) in the last step, the alkane tail emerges, dominated by van der Waals interactions with either peptide or solvent.

Free Enthalpy Differences between α-, π-, and 310-Helices of an Atomic Level Fine-Grained Alanine Deca-Peptide Solvated in Supramolecular Coarse-Grained Water.

Atomistic molecular dynamics simulations of peptides or proteins in aqueous solution are still limited to the multi-nanosecond time scale and multi-nanometer range by computational cost. Combining atomic solutes with a supramolecular solvent model in hybrid fine-grained/coarse-grained (FG/CG) simulations allows atomic detail in the region of interest while being computationally more efficient. We used enveloping distribution sampling (EDS) to calculate the free enthalpy differences between different helical conformations, i.e., α-, π-, and 310-helices, of an atomic level FG alanine deca-peptide solvated in a supramolecular CG water solvent. The free enthalpy differences obtained show that by replacing the FG solvent by the CG solvent, the π-helix is destabilized with respect to the α-helix by about 2.5 kJ mol(-1), and the 310-helix is stabilized with respect to the α-helix by about 9 kJ mol(-1). In addition, the dynamics of the peptide becomes faster. By introducing a FG water layer of 0.8 nm around the peptide, both thermodynamic and dynamic properties are recovered, while the hybrid FG/CG simulations are still four times more efficient than the atomistic simulations, even when the cutoff radius for the nonbonded interactions is increased from 1.4 to 2.0 nm. Hence, the hybrid FG/CG model, which yields an appropriate balance between reduced accuracy and enhanced computational speed, is very suitable for molecular dynamics simulation investigations of biomolecules.