Exploring the inter-molecular interactions in amyloid-β protofibril with molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area free energy calculations

Aggregation of amyloid-β (Aβ) peptides correlates with the pathology of Alzheimer's disease. However, the inter-molecular interactions between Aβ protofibril remain elusive. Herein, molecular mechanics Poisson-Boltzmann surface area analysis based on all-atom molecular dynamics simulations was performed to study the inter-molecular interactions in Aβ(17-42) protofibril. It is found that the nonpolar interactions are the important forces to stabilize the Aβ(17-42) protofibril, while electrostatic interactions play a minor role. Through free energy decomposition, 18 residues of the Aβ(17-42) are identified to provide interaction energy lower than -2.5 kcal/mol. The nonpolar interactions are mainly provided by the main chain of the peptide and the side chains of nine hydrophobic residues (Leu17, Phe19, Phe20, Leu32, Leu34, Met35, Val36, Val40, and Ile41). However, the electrostatic interactions are mainly supplied by the main chains of six hydrophobic residues (Phe19, Phe20, Val24, Met35, Val36, and Val40) and the side chains of the charged residues (Glu22, Asp23, and Lys28). In the electrostatic interactions, the overwhelming majority of hydrogen bonds involve the main chains of Aβ as well as the guanidinium group of the charged side chain of Lys28. The work has thus elucidated the molecular mechanism of the inter-molecular interactions between Aβ monomers in Aβ(17-42) protofibril, and the findings are considered critical for exploring effective agents for the inhibition of Aβ aggregation.     

Can collision‐induced negative‐ion fragmentations of [M–H]– anions be used to identify phosphorylation sites in peptides?

A joint experimental and theoretical investigation of the fragmentation behaviour of energised [M-H](-) anions from selected phosphorylated peptides has confirmed some of the most complex rearrangement processes yet to be reported for peptide negative ions. In particular: pSer and pThr (like pTyr) may transfer phosphate groups to C-terminal carboxyl anions and to the carboxyl anion side chains of Asp and Glu, and characteristic nucleophilic/cleavage reactions accompany or follow these rearrangements. pTyr may transfer phosphate to the side chains of Ser and Thr. The reverse reaction, namely transfer of a phosphate group from pSer or pThr to Tyr, is energetically unfavourable in comparison. pSer can transfer phosphate to a non-phosphorylated Ser. The non-rearranged [M-H](-) species yields more abundant product anions than its rearranged counterpart. If a peptide containing any or all of Ser, Thr and Tyr is not completely phosphorylated, negative-ion cleavages can determine the number of phosphated residues, and normally the positions of Ser, Thr and Tyr, but not which specific residues are phosphorylated. This is in accord with comments made earlier by Lehmann and coworkers.     

Signature of mobile hydrogen bonding of lysine side chains from long-range 15N-13C scalar J-couplings and computation

Amino acid side chains involved in hydrogen bonds and electrostatic interactions are crucial for protein function. However, detailed investigations of such side chains in solution are rare. Here, through the combination of long-range (15)N-(13)C scalar J-coupling measurements and an atomic-detail molecular dynamics (MD) simulation, direct insight into the structural dynamic behavior of lysine side chains in human ubiquitin has been gained. On the basis of (1)H/(13)C/(15)N heteronuclear correlation experiments selective for lysine NH(3)(+) groups, we analyzed two different types of long-range (15)N-(13)C J-coupling constants: one between intraresidue (15)Nζ and (13)Cγ nuclei ((3)J(NζCγ)) and the other between (15)Nζ and carbonyl (13)C' nuclei across a hydrogen bond ((h3)J(NζC')). The experimental (3)J(NζCγ) data confirm the highly mobile nature of the χ(4) torsion angles of lysine side chains seen in the MD simulation. The NH(3)(+) groups of Lys29 and Lys33 exhibit measurable (h3)J(NζC') couplings arising from hydrogen bonds with backbone carbonyl groups of Glu16 and Thr14, respectively. When interpreted together with the (3)J(NζCγ)-coupling constants and NMR-relaxation-derived S(2) order parameters of the NH(3)(+) groups, they strongly suggest that hydrogen bonds involving NH(3)(+) groups are of a transient and highly dynamic nature, in remarkably good agreement with the MD simulation results.     

Solid–solvent molecular interactions observed in crystal structures of β-chitin complexes

Three β-chitin structures [anhydrous, di-hydrate, mono-ethylenediamine (EDA)] recently determined by synchrotron X-ray and neutron fiber diffraction were reviewed from the viewpoint of molecular interactions. Both water and EDA molecules interact with the chitin chains through multiple hydrogen bonds. When water complexes with chitin, the hydrogen bonding pattern rearranges with the replacement of an intrachain chitin hydrogen bond by a stronger hydrogen bond between chitin and water, with an associated reduction in the degrees of freedom; the water oxygen is a much stronger acceptor than the O5 ring atom. The behavior of hydrogen exchange by deuterium supports this interpretation. EDA-molecules change the conformation of hydroxymethyl group from gg to gt, accompanied by changes in hydrogen bonds due to the strong accepting ability of the EDA nitrogen atoms. Some important interactions are in common with experimental crystallographic results of cellulosic crystals and of molecular dynamics studies. These new insights into solid–solvent interactions are valuable in understanding molecular interactions in other polysaccharides-solvents system in solution or on surface.     

The Nature of Activated Non‐classical Hydrogen Bonds: A Case Study on Acetylcholinesterase–Ligand Complexes

Molecular recognition events in biological systems are driven by non‐covalent interactions between interacting species. Here, we have studied hydrogen bonds of the CH???Y type involving electron‐deficient CH donors using dispersion‐corrected density functional theory (DFT) calculations applied to acetylcholinesterase–ligand complexes. The strengths of CH???Y interactions activated by a proximal cation were considerably strong; comparable to or greater than those of classical hydrogen bonds. Significant differences in the energetic components compared to classical hydrogen bonds and non‐activated CH???Y interactions were observed. Comparison between DFT and molecular mechanics calculations showed that common force fields could not reproduce the interaction energy values of the studied hydrogen bonds. The presented results highlight the importance of considering CH???Y interactions when analysing protein–ligand complexes, call for a review of current force fields, and opens up possibilities for the development of improved design tools for drug discovery.     

Side-chain and backbone ordering in homopolymers

In order to study the relation between backbone and side-chain ordering in proteins, we have performed multicanonical simulations of deka-peptide chains with various side groups. Glu(10), Gln(10), Asp(10), Asn(10), and Lys(10) were selected to cover a wide variety of possible interactions between the side chains of the monomers. All homopolymers undergo helix-coil transitions. We found that peptides with long side chains that are capable of hydrogen bonding, i.e., Glu(10), and Gln(10), exhibit a second transition at lower temperatures connected with side-chain ordering. This occurs in the gas phase as well as in solvent, although the character of the side-chain structure is different in each case. However, in polymers with short side chains capable of hydrogen bonding, i.e., Asp(10) and Asn(10), side-chain ordering takes place over a wide temperature range and exhibits no phase transition-like character. Moreover, non-backbone hydrogen bonds show enhanced formation and fluctuations already at the helix-coil transition temperature, indicating competition between side-chain and backbone hydrogen bond formation. Again, these results are qualitatively independent of the environment. Side-chain ordering in Lys(10), whose side groups are long and polar, also takes place over a wide temperature range and exhibits no phase transition-like character in both environments. Reasons for the observed chain length threshold and consequences from these results for protein folding are discussed.     

First-second shell interactions in metal binding sites in proteins: a PDB survey and DFT/CDM calculations

The role of the second shell in the process of metal binding and selectivity in metalloproteins has been elucidated by combining Protein Data Bank (PDB) surveys of Mg, Mn, Ca, and Zn binding sites with density functional theory/continuum dielectric methods (DFT/CDM). Peptide backbone groups were found to be the most common second-shell ligand in Mg, Mn, Ca, and Zn binding sites, followed (in decreasing order) by Asp/Glu, Lys/Arg, Asn/Gln, and Ser/Thr side chains. Aromatic oxygen- or nitrogen-containing side chains (Tyr, His, and Trp) and sulfur-containing side chains (Cys and Met) are seldom found in the second coordination layer. The backbone and Asn/Gln side chain are ubiquitous in the metal second coordination layer as their carbonyl oxygen and amide hydrogen can act as a hydrogen-bond acceptor and donor, respectively, and can therefore partner practically every first-shell ligand. The second most common outer-shell ligand, Asp/Glu, predominantly hydrogen bonds to a metal-bound water or Zn-bound histidine and polarizes the H-O or H-N bond. In certain cases, a second-shell Asp/Glu could affect the protonation state of the metal ligand. It could also energetically stabilize a positively charged metal complex more than a neutral ligand such as the backbone and Asn/Gln side chain. As for the first shell, the second shell is predicted to contribute to the metal selectivity of the binding site by discriminating between metal cations of different ionic radii and coordination geometries. The first-shell-second-shell interaction energies decay rapidly with increasing solvent exposure of the metal binding site. They are less favorable but are of the same order of magnitude as compared to the respective metal-first-shell interaction energies. Altogether, the results indicate that the structure and properties of the second shell are dictated by those of the first layer. The outer shell is apparently designed to stabilize/protect the inner-shell and complement/enhance its properties.     

Influence of hydrogen bonds on charge distribution and conformation of L-arginine

The molecular recognition of adenosine-5'-triphosphate (ATP) with L-arginine (Arg) through hydrogen bonding interactions has been found using 1H NMR, H-H NOESY, acidity titration and fluorescence spectra techniques. The interactions could influence charge distribution in Arg and induce Arg conformational variation. It is realized that Arg conformation change from a partly folded state to an extended state through the rotation of CC single bonds of Arg side chain during the molecular recognition process.     

Stabilizing interactions between aromatic and basic side chains in alpha-helical peptides and proteins. Tyrosine effects on helix circular dichroism

Here we investigate the structures and energetics of interactions between aromatic (Phe or Tyr) and basic (Lys or Arg) amino acids in alpha-helices. Side chain interaction energies are measured using helical peptides, by quantifying their helicities with circular dichroism at 222 nm and interpreting the results with Lifson-Roig-based helix/coil theory. A difficulty in working with Tyr is that the aromatic ring perturbs the CD spectrum, giving an incorrect helicity. We calculated the effect of Tyr on the CD at 222 nm by deriving the intensities of the bands directly from the electronic and magnetic transition dipole moments through the rotational strengths corresponding to each excited state of the polypeptide. This gives an improved value of the helix preference of Tyr (from 0.48 to 0.35) and a correction to the helicity for the peptides containing Tyr. We find that Phe-Lys, Lys-Phe, Phe-Arg, Arg-Phe, and Tyr-Lys are all stabilizing by -0.10 to -0.18 kcal.mol-1 when placed i, i + 4 on the surface of a helix in aqueous solution, despite the great difference in polarity between these residues. Interactions between these side chains have previously been attributed to cation-pi bonds. A survey of protein structures shows that they are in fact predominantly hydrophobic interactions between the CH2 groups of Lys or Arg and the aromatic rings.     

Molecular dynamics study of the solvation of an alpha-helical transmembrane peptide by DMSO

A 10-ns molecular dynamics study of the solvation of a hydrophobic transmembrane helical peptide in dimethyl sulfoxide (DMSO) is presented. The objective is to analyze how this aprotic polar solvent is able to solvate three groups of amino acid residues (i.e., polar, apolar, and charged) that are located in a stable helical region of a transmembrane peptide. The 25-residue peptide (sMTM7) used mimics the cytoplasmic proton hemichannel domain of the seventh transmembrane segment (TM7) from subunit a of H(+)-V-ATPase from Saccharomyces cerevisiae. The three-dimensional structure of peptide sMTM7 in DMSO has been previously solved by NMR spectroscopy. The radial and spatial distributions of the DMSO molecules surrounding the peptide as well as the number of hydrogen bonds between DMSO and the side chains of the amino acid residues involved are extracted from the molecular dynamics simulations. Analysis of the molecular dynamics trajectories shows that the amino acid side chains are fully embedded in DMSO. Polar and positively charged amino acid side chains have dipole-dipole interactions with the oxygen atom of DMSO and form hydrogen bonds. Apolar residues become solvated by DMSO through the formation of a hydrophobic pocket in which the methyl groups of DMSO are pointing toward the hydrophobic side chains of the residues involved. The dual solvation properties of DMSO cause it to be a good membrane-mimicking solvent for transmembrane peptides that do not unfold due to the presence of DMSO.     

Competition between intramolecular hydrogen bonds and solvation in phosphorylated peptides: simulations with explicit and implicit solvent

The atomic-level mechanisms of protein regulation by post-translational phosphorylation remain poorly understood, except in a few well-studied systems. Molecular mechanics simulations can in principle be used to help understand and predict the effects of protein phosphorylation, but the accuracy of the results will of course depend on the quality of the force field parameters for the phosphorylated residues as well as the quality of the solvent model. The phosphorylated residues typically carry a -2 charge at physiological pH; however, the effects of phosphorylation can sometimes be mimicked by substituting Asp or Glu for the phosphorylated residue. Here we examine the suitability of explicit and implicit solvent models for simulating phospho-serine in both the -1 and -2 charge states. Specifically, we simulate a capped phosphorylated peptide, Ace-Gly-Ser-pSer-Ser-Nme, and compare the results to each other and to experimental observables from an NMR experiment. The first major conclusion is that explicit water models (TIP3P, TIP4P and SPC/E) and a Generalized Born implicit solvent model provide reasonable agreement with the experimental observables, given appropriate partial charges for the phosphate group. The Generalized Born results, however, show greater hydrogen bonding propensity than the explicit solvent results. Distance dependent dielectric treatments perform poorly. The second major conclusion is that many ensemble-averaged properties obtained for the phosphopeptide in the -1 and -2 charge states are strikingly similar; the -1 species has a slightly higher propensity to form internal hydrogen bonds. All of the results can be rationalized by quantifying the strength of the P-O/H-N hydrogen bond, which depends on a sensitive balance between strongly favorable charge/dipole and dipole/dipole interactions and strongly unfavorable desolvation.     

Weak Intermolecular Hydrogen Bonds with Fluorine: Detection and Implications for Enzymatic/Chemical Reactions,Chemical Properties,and Ligand/Protein Fluorine NMR Screening

It is known that strong hydrogen‐bonding interactions play an important role in many chemical and biological systems. However, weak or very weak hydrogen bonds, which are often difficult to detect and characterize, may also be relevant in many recognition and reaction processes. Fluorine serving as a hydrogen‐bond acceptor has been the subject of many controversial discussions and there are different opinions about it. It now appears that there is compelling experimental evidence for the involvement of fluorine in weak intramolecular or intermolecular hydrogen bonds. Using established NMR methods, we have previously characterized and measured the strengths of intermolecular hydrogen‐bond complexes involving the fluorine moieties CH₂F, CHF₂, and CF₃, and have compared them with the well‐known hydrogen‐bond complex formed between acetophenone and the strong hydrogen‐bond donor p‐fluorophenol. We now report evidence for the formation of hydrogen bonds involving fluorine with significantly weaker donors, namely 5‐fluoroindole and water. A simple NMR method is proposed for the simultaneous measurement of the strengths of hydrogen bonds between an acceptor and a donor or water. Important implications of these results for enzymatic/chemical reactions involving fluorine, for chemical and physical properties, and for ligand/protein ¹⁹F NMR screening are analyzed through experiments and theoretical simulations.     

Hydrogen bonding mediated by key orbital interactions determines hydration enthalpy differences of phosphate water clusters

Electronic structure calculations have been carried out to provide a molecular interpretation for dihydrogen phosphate stability in water relative to that of metaphosphate. Specifically, hydration enthalpies of biologically important metaphosphate and dihydrogen phosphate with one to three waters have been computed with second-order M?ller-Plesset perturbation and density functional theory (B3LYP) with up to the aug-cc-pvtz basis set and compared to experiment. The inclusion of basis set superposition error corrections and supplemental diffuse functions are necessary to predict hydration enthalpies within experimental uncertainty. Natural bond orbital analysis is used to rationalize underlying hydrogen bond configurations and key orbital interactions responsible for the experimentally reported difference in hydration enthalpies between metaphosphate and dihydrogen phosphate. In general, dihydrogen phosphate forms stronger hydrogen bonds compared to metaphosphate due to a greater charge transfer or enhanced orbital overlap between the phosphoryl oxygen lone pairs, n(O), and the antibonding O-H bond of water. Intramolecular distal lone pair repulsion with the donor n(O) orbital of dihydrogen phosphate distorts symmetric conformations, which improves n(O) and sigma*(O-H) overlap and ultimately the hydrogen bond strength. Unlike metaphosphate, water complexed to dihydrogen phosphate can serve as both a hydrogen bond donor and a hydrogen bond acceptor, which results in cooperative charge transfer and a reduction of the energy gap between n(O) and sigma*(O-H), leading to stronger hydrogen bonds. This study offers insight into how orbital interactions mediate hydrogen bond strengths with potential implications on the understanding of the kinetics and mechanism in enzymatic phosphoryl transfer reactions.     

Determination of aliphatic side-chain conformation using cross-correlated relaxation: application to an extraordinarily stable 2'-aminoethoxy-modified oligonucleotide triplex

The structural basis for the extraordinary stability of a triple-stranded oligonucleotide in which the third strand contains 2'-aminoethoxy-substituted riboses is investigated by NMR spectroscopy. The enhanced stability of the modified triplex in comparison to the unmodified DNA triplex of the same sequence can be attributed to strong interactions of the aminoethoxy groups of the third strand with the phosphate groups of the purine strand. In molecular dynamics calculations the aminoethoxy side chain was found to be rather flexible, allowing for the presence of hydrogen bonds between the aminoethoxy group of the third strand and two different phosphates of the backbone of the second strand. To investigate the conformational preference of the aminoethoxy side chain a new NMR method has been developed which relies on CH-CH dipolar-dipolar cross-correlated relaxation rates. The results indicate that the aminoethoxy side chains adopt mainly a gauche(+) conformation, for which only one of the two hydrogen bonds inferred by NMR and molecular dynamics simulations is possible. This demonstrates a highly specific interaction between the amino group of the third strand and one of the phosphate groups of the purine strand.     

Theoretical studies of the tautomeric equilibria for five-member N-heterocycles in the gas phase and in solution

Tautomeric equilibria have been studied for five-member N-heterocycles and their methyl derivatives in the gas phase and in different solvents with dielectric constants of epsilon = 4.7-78.4. The free energy changes differently for tautomers upon solvation as compared to the gas phase, resulting in a shift of the equilibrium constant in solution. Solvents with increasing dielectric constant produce more negative solute-solvent interaction energies and increasing internal energies. The methyl-substituted imidazole and pyrrazole form delicate equilibria between two tautomeric forms. Depending on the solvent, the methyl-substituted triazoles and tetrazole have one or two major tautomers in solution. When estimating the relative solvation free energies by means of an explicit solvent model and using the FEP/MC method, one observes that the preferred tautomers differ in several cases from those predicted by the continuum solvent model. The 1,2-prototropic shift, as an intramolecular tautomerization path, requires about 50 kcal/mol activation energy for imidazole in the gas phase, and this route is also disfavored in a solution. The calculated activation free energy along the intramolecular path is 48-50 kcal/mol in chloroform and water as compared to a literature value of 13.6 kcal/mol for pyrrazole in DMSO. A molecular dynamics computer experiment favors the formation of an imidazole chain in chloroform, making the 1,3-tautomerization feasible along an intermolecular path in nonprotic solvents. In aqueous solution, one strong N-H...Ow hydrogen bond is formed for each species, whereas all other nitrogens in the ring form weaker, N...HwOw type hydrogen bonds. The tetrahydrofuran solvent acts as a hydrogen bond acceptor and forms N-H...Oether bonds. Molecules of the dichloromethane solvent are in favorable dipole-dipole interactions with the solute. The results obtained are useful in the design of N-heterocyclic ligands forming specified hydrogen bonds with protein side chains.     

Examination of the folding of a short alanine-based helical peptide with salt bridges using molecular dynamics simulation

A molecular dynamics simulation of the folding of a short alanine-based helical peptide of 17 residues with three Glu...Lys (i, i + 4) salt bridge pairs, referred to as the AEK17 peptide, was carried out. The simulation gave an estimated simulation folding time of 2.5 ns, shorter than 12 ns for an alanine-based peptide of 16 residues with three Lys residues only, referred to as the AK16 peptide, simulated previously. After folded, the AEK17 peptide had a helical content of 77%, in excellent agreement with the experimentally determined value of 80%. An examination of the folding pathways of AEK17 indicated that the peptide proceeded via three-turn helix conformations more than the helix-turn-helix conformation in the folding pathways. An analysis of interactions indicated that the formation of hydrogen bonds between Lys residue side chains and backbone carbonyls is a major factor in the abundant conformation of the three-turn helix intermediate. The substitution of three Ala with Glu residues reduces the extent of hydrophobic interaction in alanine-based AK peptides with the result that the breaking of the interactions of Lys epsilon-NH3+(side chain)...C=O(backbone) is a major activation action for the AEK17 to achieve a complete fold, in contrast to the AK16 peptide, in which breaking non-native hydrophobic interaction is the rate-determining step.     

Chemistry of aqueous silica nanoparticle surfaces and the mechanism of selective peptide adsorption

Control over selective recognition of biomolecules on inorganic nanoparticles is a major challenge for the synthesis of new catalysts, functional carriers for therapeutics, and assembly of renewable biobased materials. We found low sequence similarity among sequences of peptides strongly attracted to amorphous silica nanoparticles of various size (15-450 nm) using combinatorial phage display methods. Characterization of the surface by acid base titrations and zeta potential measurements revealed that the acidity of the silica particles increased with larger particle size, corresponding to between 5% and 20% ionization of silanol groups at pH 7. The wide range of surface ionization results in the attraction of increasingly basic peptides to increasingly acidic nanoparticles, along with major changes in the aqueous interfacial layer as seen in molecular dynamics simulation. We identified the mechanism of peptide adsorption using binding assays, zeta potential measurements, IR spectra, and molecular simulations of the purified peptides (without phage) in contact with uniformly sized silica particles. Positively charged peptides are strongly attracted to anionic silica surfaces by ion pairing of protonated N-termini, Lys side chains, and Arg side chains with negatively charged siloxide groups. Further, attraction of the peptides to the surface involves hydrogen bonds between polar groups in the peptide with silanol and siloxide groups on the silica surface, as well as ion-dipole, dipole-dipole, and van-der-Waals interactions. Electrostatic attraction between peptides and particle surfaces is supported by neutralization of zeta potentials, an inverse correlation between the required peptide concentration for measurable adsorption and the peptide pI, and proximity of cationic groups to the surface in the computation. The importance of hydrogen bonds and polar interactions is supported by adsorption of noncationic peptides containing Ser, His, and Asp residues, including the formation of multilayers. We also demonstrate tuning of interfacial interactions using mutant peptides with an excellent correlation between adsorption measurements, zeta potentials, computed adsorption energies, and the proposed binding mechanism. Follow-on questions about the relation between peptide adsorption on silica nanoparticles and mineralization of silica from peptide-stabilized precursors are raised.     

Toward assessing the position-dependent contributions of backbone hydrogen bonding to beta-sheet folding thermodynamics employing amide-to-ester perturbations

An amide-to-ester backbone substitution in a protein is accomplished by replacing an alpha-amino acid residue with the corresponding alpha-hydroxy acid, preserving stereochemistry, and conformation of the backbone and the structure of the side chain. This substitution replaces the amide NH (a hydrogen bond donor) with an ester O (which is not a hydrogen bond donor) and the amide carbonyl (a strong hydrogen bond acceptor) with an ester carbonyl (a weaker hydrogen bond acceptor), thus perturbing folding energetics. Amide-to-ester perturbations were used to evaluate the thermodynamic contribution of each hydrogen bond in the PIN WW domain, a three-stranded beta-sheet protein. Our results reveal that removing a hydrogen bond donor destabilizes the native state more than weakening a hydrogen bond acceptor and that the degree of destabilization is strongly dependent on the location of the amide bond replaced. Hydrogen bonds near turns or at the ends of beta-strands are less influential than hydrogen bonds that are protected within a hydrophobic core. Beta-sheet destabilization caused by an amide-to-ester substitution cannot be directly related to hydrogen bond strength because of differences in the solvation and electrostatic interactions of amides and esters. We propose corrections for these differences to obtain approximate hydrogen bond strengths from destabilization energies. These corrections, however, do not alter the trends noted above, indicating that the destabilization energy of an amide-to-ester mutation is a good first-order approximation of the free energy of formation of a backbone amide hydrogen bond.     

The Effects of Various Alcohols on the Stability of Melittin: A Molecular Dynamics Study

In this study, 200 ps molecular dynamics simulations were conducted to investigate the effects of various alcohols on the structural stability of melittin. The averaged helicity of melittin remained 80% in pure butanol, whereas it was below 60% both in pure water and in pure methanol. The α‐helix propensity of melittin increased with the aliphatic chain length of the alcohol. Charge‐charge interaction between Lys21 and Arg24 and polar‐nonpolar interaction between Trp19 and Arg22 are probably responsible for the higher structural integrity of the C‐terminal α‐helix over the N‐terminal one. The weaker dielectric constant of longer aliphatic chain length of alcohol possibly reduces the hydrogen bonding between amide protons and surrounding solvent molecules and simultaneously promotes the intramolecular hydrogen bonding in melittin and therefore stabilizes the secondary structure of melittin. The effect of various alcohols on stabilizing melittin is most likely due to their ability to form clusters on the surface of melittin effectively, favoring the formation of intramolecular hydrogen bonds instead of intermolecular hydrogen bonds and promoting the formation of stable α‐helices.     

Molecular dynamics study on the solvent dependent heme cooling following ligand photolysis in carbonmonoxy myoglobin

The time scale and mechanism of vibrational energy relaxation of the heme moiety in myoglobin was studied using molecular dynamics simulation. Five different solvent models, including normal water, heavy water, normal glycerol, deuterated glycerol and a nonpolar solvent, and two forms of the heme, one native and one lacking acidic side chains, were studied. Structural alteration of the protein was observed in native myoglobin glycerol solution and native myoglobin water solution. The single-exponential decay of the excess kinetic energy of the heme following ligand photolysis was observed in all systems studied. The relaxation rate depends on the solvent used. However, this dependence cannot be explained using bulk transport properties of the solvent including macroscopic thermal diffusion. The rate and mechanism of heme cooling depends upon the detailed microscopic interaction between the heme and solvent. Three intermolecular energy transfer mechanisms were considered: (i) energy transfer mediated by hydrogen bonds, (ii) direct vibration-vibration energy transfer via resonant interaction, and (iii) energy transfer via vibration-translation or vibration-rotation interaction, or in other words, thermal collision. The hydrogen bond interaction and vibration-vibration interaction between the heme and solvent molecules dominates the energy transfer in native myoglobin aqueous solution and native myoglobin glycerol solutions. For modified myoglobin, the vibration-vibration interaction is also effective in glycerol solution, different from aqueous solution. Thermal collisions form the dominant energy transfer pathway for modified myoglobin in water solution, and for both native myoglobin and modified myoglobin in a nonpolar environment. For native myoglobin in a nonpolar solvent solution, hydrogen bonds between heme isopropionate side chains and nearby protein residues, absent in the modified myoglobin nonpolar solvent solution, are key interactions influencing the relaxation pathways.