Rapid Prediction of the Hydrogen Bond Cooperativity in N-methylacetamide Chains

Automated Fragmentation Quantum Mechanical Calculation of 15N and 13C Chemical Shifts in a Membrane Protein

In this work, we developed an accurate and cost-effective automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method to calculate the chemical shifts of 15N and 13C of membrane proteins. The convergence of the AF-QM/MM method was tested using Krokinobacter eikastus rhodopsin 2 as a test case. When the distance threshold of the QM region is equal to or larger than 4.0 Å, the results of the AF-QM/MM calculations are close to convergence. In addition, the effects of selected density functionals, basis sets, and local chemical environment of target atoms on the chemical shift calculations were systematically investigated. Our results demonstrate that the predicted chemical shifts are more accurate when important environmental factors including cross-protomer interactions, lipid molecules, and solvent molecules are taken into consideration, especially for the 15N chemical shift prediction. Furthermore, with the presence of sodium ions in the environment, the chemical shift of residues, retinal, and retinal Schiff base are affected, which is consistent with the results of the solid-state nuclear magnetic resonance (NMR) experiment. Upon comparing the performance of various density functionals (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB95, and OPBE), the results show that mPW1PW91 is a suitable functional for the 15N and 13C chemical shift prediction of the membrane proteins. Meanwhile, we find that the improved accuracy of the 13Cβ chemical shift calculations can be achieved by the employment of the triple-ζ basis set. However, the employment of the triple-ζ basis set does not improve the accuracy of the 15N and 13Cα chemical shift calculations nor does the addition of a diffuse function improve the overall prediction accuracy of the chemical shifts. Our study also underscores that the AF-QM/MM method has significant advantages in predicting the chemical shifts of key ligands and nonstandard residues in membrane proteins than most widely used empirical models; therefore, it could be an accurate computational tool for chemical shift calculations on various types of biological systems.

Polarisation effects on the H-bond acceptor properties of secondary amides

H-bonding interactions in networks are stabilised by cooperativity, but the relationship between the chemical structures of the interacting functional groups and the thermodynamic consequences is not well-understood. We have used compounds with an intramolecular H-bond between a pyridine H-bond acceptor and an amide NH group to quantify cooperative effects on the H-bond acceptor properties of the amide carbonyl group. 1H NMR experiments in n-octane confirm the presence of the intramolecular H-bond and show that this interaction is intact in the 1 : 1 complex formed with perfluoro-tert-butanol (PFTB). UV-vis absorption titrations were used to measure the relationship between the association constant for formation of this complex and the H-bond acceptor properties of the pyridine involved in the intramolecular H-bond. Electron-donating substituents on the pyridine increase the strength of the intermolecular H-bond between PFTB and the amide. There is a linear relationship between the H-bond acceptor parameter β measured for the amide carbonyl group and the H-bond acceptor parameter for the pyridine. The cooperativity parameter κ determined from this relationship is 0.2, i.e. β for an amide carbonyl group is increased by one fifth of the value of β of an acceptor that interacts with the NH group. This result is reproduced by DFT calculations of H-bond parameters for the individual molecules in the gas phase, which implies that the observed cooperativity can be understood as polarisation of the electron density in the amide π-system in response to formation of a H-bond. The cooperativity parameter κ measured for the secondary amide H-bond donor and H-bond acceptor is identical, which implies that polarisation of an amide mediates the interaction between an external donor or acceptor in a reciprocal manner.

Force fields optimized against experimental data for large compound families using CombiFF: Validation considering non-target properties and polyfunctional compounds

The CombiFF scheme is a workflow for the automated calibration of force-field parameters against condensed-phase experimental data considering simultaneously entire classes of organic molecules. The main steps of this scheme are: (i) selection of a molecule family; (ii) enumeration of all isomers; (iii) query for experimental data; (iv) automatic construction of the molecular topologies; (v) iterative refinement of the force-field parameters considering the entire family. In two recent articles, CombiFF was applied to the design of GROMOS-compatible united-atom force fields for the saturated acyclic haloalkanes and for saturated acyclic compounds involving eight common chemical functional groups of oxygen and nitrogen. This calibration and the subsequent initial validation involved two limitations: (i) the experimental data considered was restricted to values for the pure-liquid density ρ_liq and the vaporization enthalpy ΔH_vap of the compounds; (ii) beyond monofunctional compounds, the training set only involved homo-polyhaloalkanes (possibly mixing halogen types) in the first study, and homo-polyfunctional compounds of the considered oxygen or nitrogen functional groups (no mixing of different group types) in the second one. The goal of this article is to further test the accuracy of CombiFF-generated force fields by extending the validation to: (i) nine additional properties that were not used as optimization targets (pure-liquid thermodynamic, dielectric and transport properties, as well as solvation properties); (ii) hetero-polyfunctional molecules that were not included in the calibration and initial validation sets. The results for the nine additional properties show good agreement with experiment, except for the shear viscosity and the dielectric permittivity. There, larger discrepancies are observed, likely due to the united-atom representation adopted for the aliphatic groups and to the implicit treatment of electronic polarization effects. The results for the hetero-polyfunctional molecules also show reasonable agreement with experiment in terms of the monitored properties.

H-bond cooperativity: polarisation effects on secondary amides

Formation of a H-bond with an amide carbonyl oxygen atom increases the strength of subsequent H-bonds formed by the amide NH, due to polarisation of the bond. The magnitude of this effect has been quantified by measuring association constants for the formation of 1 : 1 complexes of 2-hydroxylbenzamides with tri-n-butyl phosphine oxide. In 2-hydroxybenzamides, there is an intramolecular H-bond between the phenol OH group and the carbonyl oxygen atom. Comparison of the association constants measured for compounds with and without the 2-hydroxy group allows direct quantification of the effect of the intramolecular H-bond on the H-bond donor properties of the amide NH group. Substituents were used to modulate the strength of the intramolecular and intermolecular H-bonds. The presence of an intramolecular H-bond increases the strength of the intermolecular H-bond by more than one order of magnitude in n-octane solution. The increase in the H-bond donor parameter used to describe the amide NH group is directly proportional to the H-bond donor parameter of the phenol OH group that makes the intramolecular H-bond. These polarisation effects will lead to substantial cooperativity in complex systems that feature networks of non-covalent interactions, and the measurements described here provide a quantitative basis for understanding such phenomena.

Formation of an intramolecular phenol-amide H-bond leads to a dramatic increase in the H-bond donor strength of the amide NH group. Polarisation of the amide group is directly proportional to the polarity of the phenol H-bond donor.

Modeling cooperative effects in halogen-bonded infinite linear chains

Non-additivity in noncovalent interactions is an important aspect of complex systems that can lead to stronger (cooperative) interactions when three or more molecular units influence each other. The halogen bond (XB) is a highly-directional noncovalent interaction that has been found to be cooperative. Here the strength and nature of cooperativity arising in X-bonded infinite linear chains of cyanogen halides and 4-halopyridines are investigated by means of density functional theory calculations. It is found that cyanogen halide chains are highly cooperative (up to 77%), whereas pyridines are only slightly cooperative (below 21%). It is demonstrated that XB and its non-additivity can be modeled as the sum of a local term, which depends on first nearest-neighbors only, and long-range effective dipole-dipole attractions. It is shown that the local term in cyanogen halides primarily accounts for repulsive short-range screened Coulomb interactions, whereas in 4-halopyridines such a term also includes attractive contributions, which are particularly sizeable in some elongated XB conformations. This outcome reveals differences in the nature of the XBs formed in these molecular systems. Nevertheless, it is shown that both systems behave as effective point dipoles regarding cooperative effects, at any point of the XB dissociation path. As such, these results are useful contributions for the understanding and modeling of non-additive effects of noncovalent interactions.

Dynamic Water-Mediated Hydrogen Bonding in a Collagen Model Peptide

In the canonical (G-X-Y)(n) sequence of the fibrillar collagen triple helix, stabilizing direct interchain hydrogen bonding connects neighboring chains. Mutations of G can disrupt these interactions and are linked to connective tissue diseases. Here we integrate computational approaches with nuclear magnetic resonance (NMR) to obtain a dynamic view of hydrogen bonding distributions in the (POG)(4)(-)(POA)-(POG)(5) peptide, showing that the solution conformation, dynamics, and hydrogen bonding deviate from the reported X-ray crystal structure in many aspects. The simulations and NMR data provide clear evidence of inequivalent environments in the three chains. Molecular dynamics (MD) simulations indicate direct interchain hydrogen bonds in the leading chain, water bridges in the middle chain, and nonbridging waters in the trailing chain at the G → A substitution site. Theoretical calculations of NMR chemical shifts using a quantum fragmentation procedure can account for the unusual downfield NMR chemical shifts at the substitution sites and are used to assign the resonances to the individual chains. The NMR and MD data highlight the sensitivity of amide shifts to changes in the acceptor group from peptide carbonyls to water. The results are used to interpret solution NMR data for a variety of glycine substitutions and other sequence triplet interruptions to provide new connections between collagen sequences, their associated structures, dynamical behavior, and their ability to recognize collagen receptors.

Rapid evaluation of the interaction energies for carbohydrate-containing hydrogen-bonded complexes via the polarizable dipole–dipole interaction model combined with NBO or AM1 charge

The polarizable dipole–dipole interaction model has been developed to rapidly and accurately estimate the hydrogen bond distances and interaction energies for carbohydrate-containing hydrogen-bonded complexes.

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model

We have performed a density functional theory (DFT) calculation of the amide proton NMR chemical shift in proteins using a recently developed automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. Systematic investigation was carried out to examine the influence of explicit solvent molecules, cooperative hydrogen bonding effects, density functionals, size of the basis sets, and the local geometry of proteins on calculated chemical shifts. Our result demonstrates that the predicted amide proton ((1)HN) NMR chemical shift in explicit solvent shows remarkable improvement over that calculated with the implicit solvation model. The cooperative hydrogen bonding effect is also shown to improve the accuracy of (1)HN chemical shifts. Furthermore, we found that the OPBE exchange-correlation functional is the best density functional for the prediction of protein (1)HN chemical shifts among a selective set of DFT methods (namely, B3LYP, B3PW91, M062X, M06L, mPW1PW91, OB98, OPBE), and the locally dense basis set of 6-311++G**/4-31G* is shown to be sufficient for (1)HN chemical shift calculation. By taking ensemble averaging into account, (1)HN chemical shifts calculated by the AF-QM/MM approach can be used to validate the performance of various force fields. Our study underscores that the electronic polarization of protein is of critical importance to stabilizing hydrogen bonding, and the AF-QM/MM method is able to describe the local chemical environment in proteins more accurately than most widely used empirical models.

Rapid evaluation of the binding energies in hydrogen-bonded amide-thymine and amide-uracil dimers in gas phase

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 48 hydrogen-bonded amide-thymine and amide-uracil dimers are evaluated from the analytic potential energy function established in our lab recently. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6-311+G** calculations by including the BSSE correction. For all the 48 dimers, the analytic potential energy function yields the binding energies of the MP2/6-311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 46 dimers, only two differences are larger than 0.50 kcal/mol and the largest one is only 0.60 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6-311+G** with BSSE correction within the error limits of 0.050 Å for all the 48 dimers. The analytic potential energy function is further applied to four more complicated hydrogen-bonded amide-base systems involving amino acid side chain and β-sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen-bonded amide-base dimers quickly and accurately.