Ion disturbance and clustering in the NaCl water solutions

Ion clustering and the solvation properties in the NaCl solutions are explored by molecular dynamics simulations with several popular force fields. The existence of ions has a negligible disturbance to the hydrogen bond structures and rotational mobility of water beyond the first ion solvation shells, which is suggested by the local hydrogen bond structures and the rotation times of water. The potential of mean force (PMF) of ion pair in the dilute solution presents a consistent view with the populations of ion clusters in the electrolyte solutions. The aggregation level of ions is sensitive to the force field used in the simulations. The ion-ion interaction potential plays an important role in the forming of the contact ion pair. The entropy of water increases as the ion pair approaches each other and the association of ion pair is driven by the increment of water entropy according to the results from the selected force fields. The kinetic transition from the single solvent separated state to the contact ion pair is controlled by the enthalpy loss of solution.

Amino-acid solvation structure in transmembrane helices from molecular dynamics simulations

Understanding the solvation of amino acids in biomembranes is an important step to better explain membrane protein folding. Several experimental studies have shown that polar residues are both common and important in transmembrane segments, which means they have to be solvated in the hydrophobic membrane, at least until helices have aggregated to form integral proteins. In this work, we have used computer simulations to unravel these interactions on the atomic level, and classify intramembrane solvation properties of amino acids. Simulations have been performed for systematic mutations in poly-Leu helices, including not only each amino acid type, but also every z-position in a model helix. Interestingly, many polar or charged residues do not desolvate completely, but rather retain hydration by snorkeling or pulling in water/headgroups--even to the extent where many of them exist in a microscopic polar environment, with hydration levels corresponding well to experimental hydrophobicity scales. This suggests that even for polar/charged residues a large part of solvation cost is due to entropy, not enthalpy loss. Both hydration level and hydrogen bonding exhibit clear position-dependence. Basic side chains cause much less membrane distortion than acidic, since they are able to form hydrogen bonds with carbonyl groups instead of water or headgroups. This preference is supported by sequence statistics, where basic residues have increased relative occurrence at carbonyl z-coordinates. Snorkeling effects and N-/C-terminal orientation bias are directly observed, which significantly reduces the effective thickness of the hydrophobic core. Aromatic side chains intercalate efficiently with lipid chains (improving Trp/Tyr anchoring to the interface) and Ser/Thr residues are stabilized by hydroxyl groups sharing hydrogen bonds to backbone oxygens.     

The role of hydrogen bonding via interfacial water molecules in antigen-antibody complexation. The HyHEL-10-HEL interaction

To study the role of hydrogen bonding via interfacial water molecules in protein-protein interactions, we examined the interaction between hen egg white lysozyme (HEL) and its HyHEL-10 variable domain fragment (Fv) antibody. We constructed three antibody mutants (l-Y50F, l-S91A, and l-S93A) and investigated the interactions between the mutant Fvs and HEL. Isothermal titration calorimetry indicated that the mutations significantly decreased the negative enthalpy change (8-25 kJ mol(-1)), despite some offset by a favorable entropy change. X-ray crystallography demonstrated that the complexes had nearly identical structures, including the positions of the interfacial water molecules. Taken together, the isothermal titration calorimetric and x-ray crystallographic results indicate that hydrogen bonding via interfacial water enthalpically contributes to the Fv-HEL interaction despite the partial offset because of entropy loss, suggesting that hydrogen bonding stiffens the antigen-antibody complex.     

Ab initio molecular dynamics study of proton transfer in a polyglycine analog of the ion channel gramicidin A.

Proton transfer in biological systems is thought to often proceed through hydrogen-bonded chains of water molecules. The ion channel, gramicidin A (gA), houses within its helical structure just such a chain. Using the density functional theory based ab initio molecular dynamics Car-Parrinello method, the structure and dynamics of proton diffusion through a polyglycine analog of the gA ion channel has been investigated. In the channel, a proton, which is initially present as hydronium (H3O+), rapidly forms a strong hydrogen bond with a nearest neighbor water, yielding a transient H5O2+ complex. As in bulk water, strong hydrogen bonding of this complex to a second neighbor solvation shell is required for proton transfer to occur. Within gA, this second neighbor shell included not only a channel water molecule but also a carbonyl of the channel backbone. The present calculations suggest a transport mechanism in which a priori carbonyl solvation is a requirement for proton transfer.     

Molecular orbital study on the effect of the binding of water to base pair

The effects of binding water to base pairs was studied by means of the CNDO/2 molecular orbital method. The solvation energy is largest when water is bound as a proton donor and is smallest when it is stacked parallel to the plane of the base pair. The effects of two water molecules are nearly additive. The binding of one water molecule to the adenineuracil pair makes one of the two hydrogen bonds stronger and the other weaker. The change in the hydrogen bonding force is explained in terms of electrostatic and charge transfer energies. By the comparison with the adenine-cytosine pair, it is revealed that the binding of water to adenine serves to yield larger solvation energy for the complementary A-U pair than for the non-complementary A-C pair. It was also observed that the solvation energy due to the binding of water to pyrimidine was larger for A-C than for A-U.     

Role of bound water in protein-ligand association processes

The role of water molecules on the protein-ligand interface during macromolecular association has been determined. The free energy of association of insulin has been calculated by the methods of molecular mechanics and continual electrostatics (Poisson-Boltzmann model). The previously developed scheme of the decomposition of association free energy onto contributions from individual interactions has been used to calculate intermolecular interactions, the solvation free energy, and the entropies of the process of macromolecular association. An analysis of the calculated oscillation spectra indicated that the presence of water molecules on the protein-protein interface promotes an increase in the contribution of vibration entropy to the free energy of association due to the enhancement of the flexibility of the complex. It was shown that water molecules involved in the formation of protein-water-ligand hydrogen bond network change the balance of forces in the system.     

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Effective energy function for proteins in solution

A Gaussian solvent-exclusion model for the solvation free energy is developed. It is based on theoretical considerations and parametrized with experimental data. When combined with the CHARMM 19 polar hydrogen energy function, it provides an effective energy function (EEF1) for proteins in solution. The solvation model assumes that the solvation free energy of a protein molecule is a sum of group contributions, which are determined from values for small model compounds. For charged groups, the self-energy contribution is accounted for primarily by the exclusion model. Ionic side-chains are neutralized, and a distance-dependent dielectric constant is used to approximate the charge-charge interactions in solution. The resulting EEF1 is subjected to a number of tests. Molecular dynamics simulations at room temperature of several proteins in their native conformation are performed, and stable trajectories are obtained. The deviations from the experimental structures are similar to those observed in explicit water simulations. The calculated enthalpy of unfolding of a polyalanine helix is found to be in good agreement with experimental data. Results reported elsewhere show that EEF1 clearly distinguishes correctly from incorrectly folded proteins, both in static energy evaluations and in molecular dynamics simulations and that unfolding pathways obtained by high-temperature molecular dynamics simulations agree with those obtained by explicit water simulations. Thus, this energy function appears to provide a realistic first approximation to the effective energy hypersurface of proteins.     

Calculation of Proteins’ Total Side-Chain Torsional Entropy and Its Influence on Protein-Ligand Interactions

Despite the high density within a typical protein fold, the ensemble of sterically permissible side-chain repackings is vast. Here, we examine the extent of this variability that survives energetic biases due to van der Waals interactions, hydrogen bonding, salt bridges, and solvation. Monte Carlo simulations of an atomistic model exhibit thermal fluctuations among a diverse set of side-chain arrangements, even with the peptide backbone fixed in its crystallographic conformation. We have quantified the torsional entropy of this native-state ensemble, relative to that of a noninteracting reference system, for 12 small proteins. The reduction in entropy per rotatable bond due to each kind of interaction is remarkably consistent across this set of molecules. To assess the biophysical importance of these fluctuations, we have estimated side-chain entropy contributions to the binding affinity of several peptide ligands with calmodulin. Calculations for our fixed-backbone model correlate very well with experimentally determined binding entropies over a range spanning more than 80 kJ/(mol·308 K).     

Peptide folding driven by Van der Waals interactions

Contrary to the widespread view that hydrogen bonding and its entropy effect play a dominant role in protein folding, folding into helical and hairpin-like structures is observed in molecular dynamics (MD) simulations without hydrogen bonding in the peptide-solvent system. In the widely used point charge model, hydrogen bonding is calculated as part of the interaction between atomic partial charges. It is removed from these simulations by setting atomic charges of the peptide and water to zero. Because of the structural difference between the peptide and water, van der Waals (VDW) interactions favor peptide intramolecular interactions and are a major contributing factor to the structural compactness. These compact structures are amino acid sequence dependent and closely resemble standard secondary structures, as a consequence of VDW interactions and covalent bonding constraints. Hydrogen bonding is a short range interaction and it locks the approximate structure into the specific secondary structure when it is included in the simulation. In contrast to standard molecular simulations where the total energy is dominated by charge-charge interactions, these simulation results will give us a new view of the folding mechanism.     

Kinetics of rearrangement of solvation sheath of an excited molecule of the fluorescent probe 4"-dimethylaminochalcone

Processes accompanying the quenching of the fluorescent probe 4"-dimethylaminochalcone by hydroxyl groups of the proton-donor solvent 1-butanol have been studied. The kinetics of the deactivation of the excited state of 4"-dimethylaminochalcone has been monitored from the transition absorption spectra at a time resolution of 50 fs and fluorescence decay at a time resolution of 30 ps. The data obtained allow thinking that the next picture occurs in 1-butanol. At first stage, the 4"-dimethylaminochalcone molecule in its ground state forms a hydrogen bond with an alcohol molecule. At the second stage, the absorption of light quantum and corresponding rise of the dipole moment of 4"-dimethylaminochalcone take place, the initially existing hydrogen bond is retained. The third stage consists in the rearrangement of the 4"-dimethylaminochalcone solvation shell formed by alcohol dipole molecules due to an increase of the dipole of moment 4"-dimethylaminochalcone; this rearrangement takes an energy of about 24 kJ/mol, the arrangement time constant is close to 40 ps; the initial hydrogen bond is retained. The fourth stage involves processes that lead to fluorescence quenching; their time constant is about 200 ps. Taking into account that the quenching is a much slower process than the relaxation of the solvation shell, it was supposed that the quenching is not a direct consequence of the solvation shell relaxation or the existence of the hydrogen bond formed prior to excitation. Then the fluorescence quenching of 4"-dimethylaminochalcone can be accomplished through some other processes that are observed in other fluorescent molecules: (a) rearrangement of the initial hydrogen bond from a conformation that cannot quench the fluorescence of 4"-dimethylaminochalcone to a more "effective" conformation, (b) charge transfer between the excited of molecule 4"-dimethylaminochalcone and alcohol, or (c) solvent-induced twist of the 4"-dimethylaminochalcone amino group (its withdrawal from the molecule plane) by the action of the solvent.     

Thermodynamics of oligosaccharide binding to a monoclonal antibody specific for a Salmonella O-antigen point to hydrophobic interactions in the binding site.

The thermodynamic characteristics of oligosaccharide binding to an antibody binding site that is dominated by aromatic amino acids suggest that the hydrophobic effect contributes substantially to complex formation as well as hydrogen bonding and van der Waals interactions. A detailed titration microcalorimetric study on the temperature dependence of the binding of a trisaccharide, representing the epitope of a Salmonella O-antigen, showed that its maximum binding to the monoclonal antibody Se155-4 occurs just below room temperature and both enthalpy and entropy changes are strongly dependent on temperature in a mutually compensating manner. The heat capacity change also shows an unusually strong temperature dependence being large and negative above room temperature and positive below. van't Hoff analysis of the temperature dependence of the binding constant yielded a biphasic curve with two apparent intrinsic enthalpy estimations (approximately -100 kJ mol-1 above 18 degrees C and approximately +100 kJ mol-1 below), each very different from the calorimetrically determined enthalpies (ranging from about -60 kJ mol-1 to -20 kJ mol-1). This was interpreted as being due to large enthalpy contributions from concomitant reactions, most notably changes in solvation. Linear plots, -delta H0 versus -T delta S0, observed for temperature-dependent measurements mirror the behavior seen for a series of functional group replacements, suggesting that the molecular and physical origin of these phenomena are closely related and linked to the role of water in complex formation. The thermodynamic results are compared to the mode of binding determined from a 2.05-A resolution structure of the Fab-oligosaccharide complex, and with literature data for the heat capacities of sugars in aqueous solution and for the thermodynamics of carbohydrate binding to transport proteins and lectins.     

Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.     

Energetics of protein folding

The energetics of protein folding determine the 3D structure of a folded protein. Knowledge of the energetics is needed to predict the 3D structure from the amino acid sequence or to modify the structure by protein engineering. Recent developments are discussed: major factors are reviewed and auxiliary factors are discussed briefly. Major factors include the hydrophobic factor (burial of non-polar surface area) and van der Waals interactions together with peptide hydrogen bonds and peptide solvation. The long-standing model for the hydrophobic factor (free energy change proportional to buried non-polar surface area) is contrasted with the packing-desolvation model and the approximate nature of the proportionality between free energy and apolar surface area is discussed. Recent energetic studies of forming peptide hydrogen bonds (gas phase) are reviewed together with studies of peptide solvation in solution. Closer agreement is achieved between the 1995 values for protein unfolding enthalpies in vacuum given by Lazaridis-Archontis-Karplus and Makhatadze-Privalov when the solvation enthalpy of the peptide group is taken from electrostatic calculations. Auxiliary factors in folding energetics include salt bridges and side-chain hydrogen bonds, disulfide bridges, and propensities to form alpha-helices and beta-structure. Backbone conformational entropy is a major energetic factor which is discussed only briefly for lack of knowledge.     

Acceptable protein and solvent behavior in primary hydration shell simulations of hen lysozyme

The "primary hydration shell" method in molecular dynamics simulations uses a two- to three-layer thick shell of explicitly represented water molecules as the solvent around the protein of interest. We show that despite its simplicity, this computationally cheap model is capable of predicting acceptable water and protein behavior using the CHARMM22/CMAP potential function. For protein dynamics, comparisons are made with Lipari-Szabo order parameters. These have been derived from NMR relaxation parameters for pico-nano second motions of the NH groups in the main-chain and NH(2) groups in Asn/Gln side chains in hen lysozyme. It is also shown that an even simpler, and therefore faster, water-shell model leads to results in similarly good agreement with experiments, and also compared with simulations using a full box of water with periodic boundary conditions or with an implicit solvation model. Thus, the primary hydration shell method should be useful in making larger systems accessible to extensive simulations.     

Enthalpies of ligand binding to bovine neurophysins

Flow microcalorimetry and batch microcalorimetry have been used to survey the energetics of ligand binding by bovine neurophysins I and II. Calorimetry studies were supplemented by van't Hoff analyses of binding constants determined by circular dichroism. Free energies of binding of a series of di- and tripeptides that bind to the strong hormone binding site of neurophysin were partitioned into their enthalpic and entropic components. The results indicate that, at 25 degrees C, the binding of most peptides is an enthalpy-driven reaction associated with negative entropy and heat capacity changes. Studies elsewhere, supported by evidence here, indicate that the principal component of the negative enthalpy change does not arise from the increase in neurophysin dimerization associated with peptide binding. Accordingly, the negative enthalpy change is attributed to direct bonding interactions with peptide and possibly also to peptide-induced changes in tertiary or quaternary organization. Comparison of the binding enthalpies of different peptides indicated two types of bonding interactions that contribute to the negative enthalpy change of peptide ligation. Substitution of an aromatic- or sulfur-containing side chain for an aliphatic side chain in position 1 of bound peptides led to increases in negative enthalpy of from 1 to 6 kcal/mol, demonstrating that interactions typically classified as hydrophobic can have a significant exothermic component at 25 degrees C. Similarly, loss of hydrogen bonding potential in the peptide decreased the enthalpy change upon binding, in keeping with the expected enthalpic contribution of hydrogen bonds. In particular, the data suggested that the peptide backbone between residues 2 and 3 and the phenolic hydroxyl group in position 2 participate in hydrogen bonding.(ABSTRACT TRUNCATED AT 250 WORDS)     

Heat Effects of Dehydration of Human Serum Albumin in Hydrophilic Organic Solvents

A thermochemical model for describing the transfer of water from the protein phase to the organic solvent liquid phase and for determining how the solvation ability of organic solvents affects this process was developed. Enthalpy changes on the interaction of dried and hydrated human serum albumin (HSA) with hydrophilic organic solvents (dimethyl sulfoxide, formamide, ethanol, methanol and acetic acid) and water were measured by isothermal calorimetry at 25 °C. The initial hydration level of human serum albumin was varied in the entire water content range from 0–30 % [g water/g HSA]. The dependence of the interaction enthalpies on the initial water content is complex. The interaction enthalpies of the dried HSA with organic solvents are exothermic. At low water contents (less than 0.1 g/g), there is a sharp increase in the interaction enthalpy values. At the highest water contents (more than 0.2 g/g), the interaction enthalpies are endothermic for acetic acid and formamide and exothermic for DMSO, methanol, and ethanol. These thermochemical data were analyzed in conjunction with the results for the water adsorption in organic solvents to calculate the molar enthalpies of dehydration of HSA in organic liquids. It was found that the dehydration enthalpy changes may be endothermic or exothermic depending on the initial water content and the water solvation enthalpy value. From the results obtained, it can be concluded that: (i) only the solvation of water by hydrophilic organic solvent determines the changes in the dehydration enthalpy values, and (ii) the data for the enthalpies of solvation of water by the solvent at infinite dilution reflect this effect.     

Self-association behaviour of atactic polymethacrylic acid in aqueous solution investigated by atomistic molecular dynamics simulations

The self-association behaviour of atactic poly(methacrylic acid) (a-PMA) in water was investigated by atomistic molecular dynamics (MD) simulations. Simulations show that interchain association of a-PMA occurs only in its un-neutralised form, by hydrogen bonding between –COOH groups, which is in agreement with the experimental observation. Chain conformations, dihedral angle distributions, hydration behaviour, scattering structure factor and enthalpy-of-hydration (i.e. aqueous solvation) were analysed as a function of concentration for un-neutralised PMA, across dilute to concentrated regimes. The average 〈R_g〉 of the chain remains unaffected in solution and also for amorphous undissolved a-PMA phase, confirming the occurrence of the approximate theta-solution condition for the first time, as revealed by simulations, in a polar hydrogen-bonding polymer aqueous solution. Chain hydration behaviour and scattering structure factor show significant changes in concentrated regime. Scattering intensity collapse occurs in concentrated PMA solution, due to the existence of the swollen regime captured for the first time by explicit-MD-simulations. The hydration of PMA is driven by H-bonding, specifically between H atoms of the COOH groups and O atoms of water molecules in the closest coordination shell. The enthalpy of hydration of PMA is dominated by PMA–water interactions (charges and H-bonding). The thermodynamic contributions of PMA–PMA and PMA–water interactions towards the electrostatics as well as the dispersion components of the total solvation-enthalpy become more favourable than water–water interactions.