A theoretical study on the origin of cooperativity in the formation of 3(10)- and alpha-helices

Stabilization of the Protein Structure by the Many-Body Cooperative Effect in the NH/π Hydrogen-bonding Tryptophan Triad

The indole ring of tryptophan can form NH/π hydrogen bonds, acting both as a hydrogen donor at the NH group in the pyrrole subring and as a hydrogen acceptor at the benzene subring. In the structural core of the trimeric stable protein Pholiota squarrosa lectin (PhoSL), three indoles are symmetrically arranged and form NH/π hydrogen bonds among each other. Here, we conducted quantum chemical calculations on this indole triad by using various methods and basis sets. The analyses revealed cooperativity in triad formation, with the many-body effect contributing approximately -2 kcal mol-1, which significantly stabilizes this protein. Symmetry-adapted perturbation theory ascribed this effect to the induced polarization. The electrostatic potential and atomic charges indeed revealed a charge redistribution through the NH/π hydrogen bond, which was favorable for triad formation.

Cooperativity and Frustration Effects (or Lack Thereof) in Polarizable and Non-polarizable Force Fields

Understanding cooperativity and frustration is crucial for studying biological processes such as molecular recognition and protein aggregation. Force fields have been extensively utilized to explore cooperativity in the formation of protein secondary structures and self-assembled systems. Multiple studies have demonstrated that polarizable force fields provide more accurate descriptions of this phenomenon compared to fixed-charge pairwise nonpolarizable force fields, thanks to the incorporation of polarization effects. In this study, we assess the performance of the AMOEBA polarizable force field and the AMBER and OPLS nonpolarizable pairwise force fields in capturing positive and negative cooperativity recently explored in neutral and charged molecular clusters using density functional theory. Our findings show that polarizable and nonpolarizable force fields qualitatively reproduce the relative cooperativity observed in electron structure calculations. However, AMBER and OPLS fail to describe absolute cooperativity. In contrast, AMOEBA accounts for the absolute cooperativity by considering interactions beyond pairwise interactions. According to the energy decomposition analysis, it is observed that the electrostatic interactions calculated with the AMBER and OPLS force fields seem to play an important and counterintuitive role in reproducing the adiabatic interaction energies calculated with density functional theory. However, it is important to note that these force fields, due to their nature, do not explicitly incorporate many-body effects, which limits their ability to accurately describe cooperativity. On the other hand, frustration in polarizable and nonpolarizable force fields is caused by changes in bond stretching and angle bending terms of the building blocks when they are forming a complex.

Quantum chemical studies on hydrogen bonds in helical secondary structures

We present a brief review of our recent computational studies of hydrogen bonds (H-bonds) in helical secondary structures of proteins, α-helix and 3₁₀-helix, using a Negative Fragmentation Approach with density functional theory. We found that the depolarized electronic structures of the carbonyl oxygen of the ith residue and the amide hydrogen of the (i + 4)th residue cause weaker H-bond in an α-helix than in an isolated H-bond. Our calculations showed that the H-bond energies in the 3₁₀-helix were also weaker than those of the isolated H-bonds. In the 3₁₀-helices, the adjacent N-H group at the (i + 1)th residue was closer to the C=O group of the H-bond pair than the adjacent C=O group in the 3₁₀-helices, whereas the adjacent C=O group at the (i + 1)th residue was close to the H-bond acceptor in α-helices. Therefore, the destabilization of the H-bond is attributed to the depolarization caused by the adjacent residue of the helical backbone connecting the H-bond donor and acceptor. The differences in the change in electron density revealed that such depolarizations were caused by the local electronic interactions in their neighborhood inside the helical structure and redistributed the electron density. We also present the improvements in the force field of classical molecular simulation, based on our findings.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01034-5.

Depolarizing Effects in Hydrogen Bond Energy in 310-Helices Revealed by Quantum Chemical Analysis

Hydrogen-bond (H-bond) energies in 3₁₀-helices of short alanine peptides were systematically examined by precise DFT calculations with the negative fragmentation approach (NFA), a modified method based on the molecular tailoring approach. The contribution of each H-bond was evaluated in detail from the 3₁₀-helical conformation of total energies (whole helical model, WH_3-10 model), and the results were compared with the property of H-bond in α-helix from our previous study. The H-bond energies of the WH_3-10 model exhibited tendencies different from those exhibited by the α-helix in that they depended on the helical position of the relevant H-bond pair. H-bond pairs adjacent to the terminal H-bond pairs were observed to be strongly destabilized. The analysis of electronic structures indicated that structural characteristics cause the destabilization of the H-bond in 3₁₀-helices. We also found that the longer the helix length, the more stable the H-bond in the terminal pairs of the WH_3-10 model, suggesting the action of H-bond cooperativity.

Modulation of proton-coupled electron transfer reactions in lysine-containing alpha-helixes: alpha-helixes promoting long-range electron transfer

The proton-coupled electron transfer (PCET) reaction plays an important role in promoting many biological and chemical reactions. Usually, the rate of the PCET reaction increases with an increase in the electron transfer distance because long-range electron transfer requires more free energy barriers. Our density functional theory calculations here reveal that the mechanism of PCET occurring in lysine-containing alpha(α)-helixes changes with an increasing number of residues in the α-helical structure and the different conformations because of the modulation of the excess electron distribution by the α-helical structures. The rate constants of the corresponding PCET reactions are independent of or substantially shallower dependent on the electron transfer distances along α-helixes. This counter-intuitive behavior can be attributed to the fact that the formation of larger macro-cylindrical dipole moments in longer helixes can promote electron transfer along the α-helix with a low energy barrier. These findings may be useful to gain insights into long-range electron transfer in proteins and design α-helix-based electronics via the regulation of short-range proton transfer.

Quantifying Frustrations for Molecular Complexes with Noncovalent Interactions

Molecular systems bound together through noncovalent interactions are involved in a lot of life-essential processes such as molecular recognition, signal transduction, and allosteric regulation. While cooperation as an important effect discovered in these systems focuses on the behavior of system's entirety, we need also examine the behavior of individual parts. In this work, using the distortion energy as the descriptor, we quantify frustration as the energetic loss of individual parts due to the formation of nonadditive molecular complexes. The applicability of our approach has been illustrated by a few simple clusters. Our results show that the frustration effect is smaller than the cooperation effect, but same as cooperativity, it can be both positive and negative. The ultimate benefit of a system made of multiple parts is dictated by the balance between the cooperative behavior among parts and the sacrifice from its individuals. This conflicting yet complementary conceptual pair of cooperation and frustration provides us with a different perspective from the systems' viewpoint for molecular complexes. This new angle of appreciating molecular complexes can be applied in conformational changes, enzymatic catalysis, and many more.

Design of Organo-Peptides As Bipartite PCSK9 Antagonists

Proprotein convertase subtilisin/kexin 9 (PCSK9) has become an important therapeutic target for lipid lowering, since it regulates low-density lipoprotein cholesterol (LDL-c) levels by binding to liver LDL receptors (LDLR) and effecting their intracellular degradation. However, the development of small molecule inhibitors is hampered by the lack of attractive PCSK9 target sites. We recently discovered helical peptides that are able to bind to a cryptic groove site on PCSK9, which is situated in proximity to the main LDLR binding site. Here, we designed potent bipartite PCSK9 inhibitors by appending organic moieties to a helical groove-binding peptide to reach a hydrophobic pocket in the proximal LDLR binding region. The ultimately designed 1-amino-4-phenylcyclohexane-1-carbonyl extension improved the peptide affinity by >100-fold, yielding organo-peptide antagonists that potently inhibited PCSK9 binding to LDLR and preserved cellular LDLR. These new bipartite antagonists have reduced mass and improved potency compared to the first-generation peptide antagonists, further validating the PCSK9 groove as a viable therapeutic target site.

Hybrid polymers bearing oligo-l-lysine(carboxybenzyl)s: synthesis and investigations of secondary structure

Hybrid polymers of peptides resembling (partially) folded protein structures are promising materials in biomedicine, especially in view of folding-interactions between different segments. In this study polymers bearing repetitive peptidic folding elements, composed of N-terminus functionalized bis-ω-ene-functional oligo-l-lysine(carboxybenzyl(Z))s (Lys_n) with repeating units (n) of 3, 6, 12, 24 and 30 were successfully synthesized to study their secondary structure introduced by conformational interactions between their chains. The pre-polymers of ADMET, narrowly dispersed Lys_ns, were obtained by ring opening polymerization (ROP) of N-carboxyanhydride (NCA) initiated with 11-amino-undecene, following N-terminus functionalization with 10-undecenoyl chloride. The resulting Lys_ns were subsequently polymerized via ADMET polymerization by using Grubbs’ first generation (G1) catalyst in 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) generating the ADMET polymers (A-[Lys_n]_m) (m = 2–12) with molecular weights ranging from 3 to 28 kDa, displaying polydispersity (Đ) values in the range of 1.5–3.2. After chemical analyses of Lys_ns and A-[Lys_n]_ms by ¹H-NMR, GPC and MALDI-ToF MS, secondary structural investigations were probed by CD spectroscopy and IR spectroscopy in 2,2,2-trifluoroethanol (TFE). In order to study A-[Lys_n]_ms with defined molecular weights and low polydispersity values (Đ = 1.03–1.48), the ADMET polymers A-[Lys_n=3]_m=3 and A-[Lys_n=24]_m=4 were fractionated by preparative GPC, and subsequently analysed by ¹H-NMR, analytical GPC, MALDI-ToF MS and CD spectroscopy. We can demonstrate the influence of chain length of the generated polymers on the formation of secondary structures by comparing Lys_ns with varying n values to the ADMET-polymers with the help of spectroscopic techniques such as CD and FTIR-spectroscopy in a helicogenic solvent.

We demonstrate the influence of chain length of segmented polymers bearing dynamic folding elements onto the formation of secondary structures with the help of spectroscopic techniques such as CD and FTIR-spectroscopy in a helicogenic solvent.

The relative stability of trpzip1 and its mutants determined by computation and experiment

The single-point mutations of tprzip1 are indicated at left, and their relative energetics are compared at right.

Homogeneous Molecular Systems are Positively Cooperative, but Charged Molecular Systems are Negatively Cooperative

Molecular systems bound together through noncovalent interactions are ubiquitous in nature, many of which are involved in essential life processes, yet little is known about the principles governing their structure, stability, and function. Cooperativity as one of the intrinsic properties in these systems plays a key role. In this work, on the basis of our recent quantification scheme of the cooperativity effect, we present a general pattern to identify which systems are positively cooperative and which are negatively cooperative. We show that cooperativity in homogeneous molecular systems is positive, but cooperativity in charged molecular systems is negative. We also employ analytical tools from energetics and information perspectives to appreciate the origin of the cooperativity effect. We find that positive cooperativity is dominated by the exchange-correlation interaction and steric effect, whereas negative cooperativity is governed by the electrostatic interaction. Our results should have strong implications for better understanding molecular recognition, protein folding, signal transduction, allosteric regulation, and other processes.

Understanding the R882H mutation effects of DNA methyltransferase DNMT3A: a combination of molecular dynamics simulations and QM/MM calculations

The AML-related high-frequent R882H mutation of DNA (cytosine-5)-methyltransferase 3A (DNMT3A), a key enzyme forde novoepigenetic methylation in human beings, was characterized by a disturbing conformation ofS-adenosylmethionine (SAM).

Quantification and origin of cooperativity: insights from density functional reactivity theory

Cooperativity is a widely used chemical concept whose existence is ubiquitous in chemical and biological systems but whose quantification is still controversial and origin much less appreciated. In this work, using the interaction energy of a molecular system, which is composed of multiple copies of a building block, we propose a quantitative measurement to evaluate the cooperativity effect. This quantification approach is then applied to six molecular systems, i.e., water cluster, argon cluster, protonated water cluster, zinc atom cluster, water cluster on top of a graphene sheet, and alpha helix of glycine amino acids, each with up to 20 copies of the building block. Cooperativity is seen in all these systems. Both positive and negative cooperativity effects are observed. Employing the two energy partition schemes in density functional theory and the information-theoretic quantities such as Shannon entropy, Fisher information, information gain, etc., we then examine the origin of the cooperativity effect for these systems. Strong linear correlations between the cooperativity measure and some of these theoretical quantities have been unveiled. With these correlations, we are able to quantitatively account for their origin of cooperativity. Our results show that the interactions governing the existence and validity of the cooperativity effect are complicated. An opposite mechanism in enthalpy-entropy compensation for positive and negative cooperativity has been unveiled. These results should provide new insights and understandings from a different viewpoint about the nature and origin of cooperativity to appreciate this vastly important chemical concept.

Structural and functional aspects of the interaction partners of the small heat-shock protein in Synechocystis

The canonical function of small heat-shock proteins (sHSPs) is to interact with proteins destabilized under conditions of cellular stress. While the breadth of interactions made by many sHSPs is well-known, there is currently little knowledge about what structural features of the interactors form the basis for their recognition. Here, we have identified 83 in vivo interactors of the sole sHSP in the cyanobacterium Synechocystis sp. PCC 6803, HSP16.6, reflective of stable associations with soluble proteins made under heat-shock conditions. By performing bioinformatic analyses on these interactors, we identify primary and secondary structural elements that are enriched relative to expectations from the cyanobacterial genome. In addition, by examining the Synechocystis interactors and comparing them with those identified to bind sHSPs in other prokaryotes, we show that sHSPs associate with specific proteins and biological processes. Our data are therefore consistent with a picture of sHSPs being broadly specific molecular chaperones that act to protect multiple cellular pathways.

Tuning peptide self-assembly by an in-tether chiral center

The self-assembly of peptides into ordered nanostructures is important for understanding both peptide molecular interactions and nanotechnological applications. However, because of the complexity and various self-assembling pathways of peptide molecules, design of self-assembling helical peptides with high controllability and tunability is challenging. We report a new self-assembling mode that uses in-tether chiral center-induced helical peptides as a platform for tunable peptide self-assembly with good controllability. It was found that self-assembling behavior was governed by in-tether substitutional groups, where chirality determined the formation of helical structures and aromaticity provided the driving force for self-assembly. Both factors were essential for peptide self-assembly to occur. Experiments and theoretical calculations indicate long-range crystal-like packing in the self-assembly, which was stabilized by a synergy of interpeptide π-π and π-sulfur interactions and hydrogen bond networks. In addition, the self-assembled peptide nanomaterials were demonstrated to be promising candidate materials for applications in biocompatible electrochemical supercapacitors.

Charge-patching method for the calculation of electronic structure of polypeptides

Based on the CPM method, the charge densities of polypeptides can be generated and their electronic structure can be further calculated.

Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field

An accurate potential energy model is crucial for biomolecular simulations. Despite many recent improvements of classical protein force fields, there are remaining key issues: much weaker temperature dependence of folding/unfolding equilibrium and overly collapsed unfolded or disordered states. For the latter problem, a new water model (TIP4P-D) has been proposed to correct the significantly underestimated water dispersion interactions. Here, using TIP4P-D, we reveal problems in current force fields through failures in folding model systems (a polyalanine peptide, Trp-cage, and the GB1 hairpin). By using residue-specific parameters to achieve better match between amino acid sequences and native structures and adding a small H-bond correction to partially compensate the missing many-body effects in α-helix formation, the new RSFF2+ force field with the TIP4P-D water model can excellently reproduce experimental melting curves of both α-helical and β-hairpin systems. The RSFF2+/TIP4P-D method also gives less collapsed unfolded structures and describes well folded proteins simultaneously.

Utilizing Electron Spin Echo Envelope Modulation To Distinguish between the Local Secondary Structures of an α-Helix and an Amphipathic 310-Helical Peptide

Electron spin echo envelope modulation (ESEEM) spectroscopy was used to distinguish between the local secondary structures of an α-helix and a 3₁₀-helix. Previously, we have shown that ESEEM spectroscopy in combination with site-directed spin labeling (SDSL) and ²H-labeled amino acids (i) can probe the local secondary structure of α-helices, resulting in an obvious deuterium modulation pattern, where i+4 positions generally show larger ²H ESEEM peak intensities than i+3 positions. Here, we have hypothesized that due to the unique turn periodicities of an α-helix (3.6 residues per turn with a pitch of 5.4 Å) and a 3₁₀-helix (3.1 residues per turn with a pitch of 5.8-6.0 Å), the opposite deuterium modulation pattern would be observed for a 3₁₀-helix. In this study, ²H-labeled d₁₀-leucine (Leu) was substituted at a specific Leu residue (i) and a nitroxide spin label was positioned 2, 3, and 4 residues away (denoted i+2 to i+4) on an amphipathic model peptide, LRL₈. When LRL₈ is solubilized in trifluoroethanol (TFE), the peptide adopts an α-helical structure, and alternatively, forms a 3₁₀-helical secondary structure when incorporated into liposomes. Larger ²H ESEEM peaks in the FT frequency domain data were observed for the i+4 samples when compared to the i+3 samples for the α-helix whereas the opposite pattern was revealed for the 3₁₀-helix. These unique patterns provide pertinent local secondary structural information to distinguish between the α-helical and 3₁₀-helical structural motifs for the first time using this ESEEM spectroscopic approach with short data acquisition times (∼30 min) and small sample concentrations (∼100 μM) as well as providing more site-specific secondary structural information compared to other common biophysical approaches, such as CD.

Proton Dynamics in Protein Mass Spectrometry

Native electrospray ionization/ion mobility-mass spectrometry (ESI/IM-MS) allows an accurate determination of low-resolution structural features of proteins. Yet, the presence of proton dynamics, observed already by us for DNA in the gas phase, and its impact on protein structural determinants, have not been investigated so far. Here, we address this issue by a multistep simulation strategy on a pharmacologically relevant peptide, the N-terminal residues of amyloid-β peptide (Aβ(1-16)). Our calculations reproduce the experimental maximum charge state from ESI-MS and are also in fair agreement with collision cross section (CCS) data measured here by ESI/IM-MS. Although the main structural features are preserved, subtle conformational changes do take place in the first ∼0.1 ms of dynamics. In addition, intramolecular proton dynamics processes occur on the picosecond-time scale in the gas phase as emerging from quantum mechanics/molecular mechanics (QM/MM) simulations at the B3LYP level of theory. We conclude that proton transfer phenomena do occur frequently during fly time in ESI-MS experiments (typically on the millisecond time scale). However, the structural changes associated with the process do not significantly affect the structural determinants.

Density functional study of molecular interactions in secondary structures of proteins

Proteins play diverse and vital roles in biology, which are dominated by their three-dimensional structures. The three-dimensional structure of a protein determines its functions and chemical properties. Protein secondary structures, including α-helices and β-sheets, are key components of the protein architecture. Molecular interactions, in particular hydrogen bonds, play significant roles in the formation of protein secondary structures. Precise and quantitative estimations of these interactions are required to understand the principles underlying the formation of three-dimensional protein structures. In the present study, we have investigated the molecular interactions in α-helices and β-sheets, using ab initio wave function-based methods, the Hartree-Fock method (HF) and the second-order Møller-Plesset perturbation theory (MP2), density functional theory, and molecular mechanics. The characteristic interactions essential for forming the secondary structures are discussed quantitatively.

Generalized energy-based fragmentation approach and its applications to macromolecules and molecular aggregates

Conspectus The generalized energy-based fragmentation (GEBF) approach provides a very simple way of approximately evaluating the ground-state energy or properties of a large system in terms of ground-state energies of various small "electrostatically embedded" subsystems, which can be calculated with any traditional ab initio quantum chemistry (X) method (X = Hartree-Fock, density functional theory, and so on). Due to its excellent parallel efficiency, the GEBF approach at the X theory level (GEBF-X) allows full quantum mechanical (QM) calculations to be accessible for systems with hundreds and even thousands of atoms on ordinary workstations. The implementation of the GEBF approach at various theoretical levels can be easily done with existing quantum chemistry programs. This Account reviews the methodology, implementation, and applications of the GEBF-X approach. This method has been successfully applied to optimize the structures of various large systems including molecular clusters, polypeptides, proteins, and foldamers. Such investigations could allow us to elucidate the origin and nature of the cooperative interaction in secondary structures of long peptides or the driving force of the self-assembly processes of aromatic oligoamides. These GEBF-based QM calculations reveal that the structures and stability of various complex systems result from a subtle balance of many types of noncovalent interactions such as hydrogen bonding and van der Waals interactions. The GEBF-based ab initio molecular dynamics (AIMD) method also allows the investigation of dynamic behaviors of large systems on the order of tens of picoseconds. It was demonstrated that the conformational dynamics of two model peptides predicted by GEBF-based AIMD are noticeably different from those predicted by the classical force field MD method. With the target of extending QM calculations to molecular aggregates in the condensed phase, we have implemented the GEBF-based multilayer hybrid models, which could provide satisfactory descriptions of the binding energies between a solute molecule and its surrounding waters and the chain-length dependence of the conformational changes of oligomers in aqueous solutions. A coarse-grained polarizable molecular mechanics model, furnished with GEBF-X dipole moments of subsystems, exhibits some advantages of treating the electrostatic polarization with reduced computational costs. We anticipate that the GEBF approach will continue to develop with the ultimate goal of studying complicated phenomena at mesoscopic scales and serve as a practical tool to elucidate the structure and dynamics of chemical and biological systems.

A B3LYP and MP2(full) theoretical investigation into cooperativity effects, aromaticity and thermodynamic properties in the Na(+)⋯benzonitrile⋯H2O ternary complex

The cooperativity effects between H-bonding and Na(+)⋯π or Na(+)⋯σ interactions in Na(+)⋯benzonitrile⋯H2O complexes were investigated using the B3LYP and MP2(full) methods with 6-311++G(2d,p) and aug-cc-pVTZ basis sets. The thermodynamic cooperativity and the influence of this cooperativity on aromaticity was evaluated by nucleus-independent chemical shifts (NICS). The results showed that the influence of the Na(+)⋯σ or Na(+)⋯π interaction on the hydrogen bond is more pronounced than that of the latter on the former. The cooperativity effect appeared in the Na(+)⋯σ interaction complex while the anti-cooperativity effect tended to be in the Na(+)⋯π system. The change in enthalpy is the major factor driving cooperativity. Thermodynamic cooperativity is not in accordance with the cooperativity effect evaluated by the change of interaction energy. The ring aromaticity of is weakened while the bond dissociation energy (BDE) of the C-CN bond increases upon ternary complex formation. The cooperativity effect (E coop) correlates with R c (NICS(1)ternary/NICS(1)binary) and ΔΔδ (Δδ ternary - Δδ binary) involving the ring and C ≡ N bond, as well as R BDE(C-CN) [BDE(C-CN)ternary/BDE(C-CN)binary], respectively. AIM (atoms in molecules) analysis confirms the existence of cooperativity.

A peptide loop and an α-helix N-terminal serving as alternative electron hopping relays in proteins

This work presents a density functional theory calculational study for clarifying that peptide loops (-[peptide](n)-) including the N-terminal and the C-terminal oligopeptides and the α-helix N-terminal can serve as an intriguing kind of relay elements, as an addition to the known relay stations served by aromatic amino acids for electron hopping migration. For these protein motifs, an excess electron generally prefers to reside at the -NH(3)(+) group in a Rydberg state for the N-terminal peptides, or at the -COOH group in a dipole-bound state for the C-terminal peptides, and at the N-terminal in a dipole-bound π*-orbital state for the peptide loops and α-helices. The electron binding ability can be effectively enhanced by elongation for the α-helix N-terminal, and by bending, twisting, and even β-turning for the peptide chains. The relay property is determined by the local dipole instead of the total dipole of the peptide chains. Although no direct experiment supports this hypothesis, a series of recent studies regarding charge hopping migration associated with the peptide chains and helices could be viewed as strong evidence. But, further studies are still needed by considering the effects from relative redox potential between the donor and acceptor sites, protein environment, and structure water molecules.

Direct evaluation of individual hydrogen bond energy in situ in intra- and intermolecular multiple hydrogen bonds system

The results of evaluating the individual hydrogen bond (H-bond) strength are expected to be helpful for the rational design of new strategies for molecular recognition or supramolecular assemblies. Unfortunately, there is few obvious and unambiguous means of evaluating the energy of a single H-bond within a multiple H-bonds system. We present a local analytic model, ABEEMσπ H-bond energy (HBE) model based on ab initio calculations (MP2) as benchmark, to directly and rapidly evaluate the individual HBE in situ in inter- and intramolecular multiple H-bonds system. This model describes the HBE as the sum of electrostatic and van der Waals (vdW) interactions which all depend upon the geometry and environment, and the ambient environment of H-bond in the model is accounted fairly. Thus, it can fairly consider the cooperative effect and secondary effect. The application range of ABEEMσπ HBE model is rather wide. This work has discussed the individual H-bond in DNA base pair and protein peptide dimers. The results indicate that the interactions among donor H atom, acceptor atom as well as those atoms connected to them with 1,2 or 1,3 relationships are all important for evaluating the HBE, although the interaction between the donor H atom and the acceptor atom is large. Furthermore, our model quantitatively indicates the polarization ability of N, O, and S in a new style, and gives the percentage of the polarization effect in HBE, which can not be given by fixed partial charge force field.

Polarization of intraprotein hydrogen bond is critical to thermal stability of short helix

Simulation result for protein folding/unfolding is highly dependent on the accuracy of the force field employed. Even for the simplest structure of protein such as a short helix, simulations using the existing force fields often fail to produce the correct structural/thermodynamic properties of the protein. Recent research indicated that lack of polarization is at least partially responsible for the failure to successfully fold a short helix. In this work, we develop a simple formula-based atomic charge polarization model for intraprotein (backbone) hydrogen bonding based on the existing AMBER force field to study the thermal stability of a short helix (2I9M) by replica exchange molecular dynamics simulation. By comparison of the simulation results with those obtained by employing the standard AMBER03 force field, the formula-based atomic charge polarization model gave the helix melting curve in close agreement with the NMR experiment. However, in simulations using the standard AMBER force field, the helix was thermally unstable at the temperature of the NMR experiment, with a melting temperature almost below the freezing point. The difference in observed thermal stability from these two simulations is the effect of backbone intraprotein polarization, which was included in the formula-based atomic charge polarization model. The polarization of backbone hydrogen bonding thus plays a critical role in the thermal stability of helix or more general protein structures.

N-Protonated isomers as gateways to peptide ion fragmentation

According to the popular "mobile proton model" for peptide ion fragmentation in tandem mass spectrometry, peptide bond cleavage is typically preceded by intramolecular proton transfer from basic sites to an amide nitrogen in the backbone. If the intrinsic barrier to dissociation is the same for all backbone sites, the fragmentation propensity at each amide bond should reflect the stability of the corresponding N-protonated isomer. This hypothesis was tested by using ab initio and force-field computations on several polyalanines and Leu-enkephalin. The results agree acceptably with experimental reports, supporting the hypothesis. It was found that backbone N-protonation is most favorable near the C-terminus. The preference for C-terminal N-protonation, which is stronger for longer polyalanines, may be understood in terms of the well known "helix macrodipole" in the corresponding helical conformations. The opposite stability trend is found for peptides constrained to be linear, which is initially surprising but turns out to be consistent with the reversed direction of the macrodipole in the linear conformation.