Dispersion interactions govern the strong thermal stability of a protein

Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

London dispersion (LD) interactions are the main contribution of the attractive part of the van der Waals potential. Even though LD effects are the driving force for molecular aggregation and recognition, the role of these omnipresent interactions in structure and reactivity had been largely underappreciated over decades. However, in the recent years considerable efforts have been made to thoroughly study LD interactions and their potential as a chemical design element for structures and catalysis. This was made possible through a fruitful interplay of theory and experiment. This review highlights recent results and advances in utilizing LD interactions as a structural motif to understand and utilize intra- and intermolecularly LD-stabilized systems. Additionally, we focus on the quantification of LD interactions and their fundamental role in chemical reactions.

Exploring London Dispersion and Solvent Interactions at Alkyl-Alkyl Interfaces Using Azobenzene Switches

Interactions on the molecular level control structure as well as function. Especially interfaces between innocent alkyl groups are hardly studied although they are of great importance in larger systems. Herein, London dispersion in conjunction with solvent interactions between linear alkyl chains was examined with an azobenzene-based experimental setup. Alkyl chains in all meta positions of the azobenzene core were systematically elongated, and the change in rate for the thermally induced Z→E isomerization in n-decane was determined. The stability of the Z-isomer increased with longer chains and reached a maximum for n-butyl groups. Further elongation led to faster isomerization. The origin of the intramolecular interactions was elaborated by various techniques, including 1 H NOESY NMR spectroscopy. The results indicate that there are additional long-range interactions between n-alkyl chains with the opposite phenyl core in the Z-state. These interactions are most likely dominated by attractive London dispersion. This work provides rare insight into the stabilizing contributions of highly flexible groups in an intra- as well as an intermolecular setting.

Benchmarking the Effective Fragment Potential Dispersion Correction on the S22 Test Set

The usual modeling of dispersion interactions in density functional theory (DFT) is often limited by the use of empirically fitted parameters. In this study, the accuracies of the popular empirical dispersion corrections and the first-principles derived effective fragment potential (EFP) dispersion correction are compared by computing the DFT-D and HF-D equilibria interaction energies and intermolecular distances of the S22 test set dimers. Functionals based on the local density approximation (LDA) and generalized gradient approximation (GGA), as well as hybrid functionals, are compared for the DFT-D calculations using coupled cluster CCSD(T) at the complete basis set (CBS) limit as the reference method. In general, the HF-D(EFP) method provides accurate equilibrium dimerization energies and intermolecular distances for hydrogen-bonded systems compared to the CCSD(T)/CBS reference data without using any empirical parameters. For dispersion-dominant and mixed systems, the structures and interaction energies obtained with the B3LYP-D(EFP) method are similar to or better than those obtained with the other DFT-D and HF-D methods. Overall, the first-principles derived -D(EFP) correction presents a robust alternative to the empirical -D corrections when used with the B3LYP functional for dispersion-dominant and mixed systems or with Hartree-Fock for hydrogen-bonded systems.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

Implicit treatment of solvent dispersion forces in protein simulations

A model is proposed for the evaluation of dispersive forces in a continuum solvent representation for use in large-scale computer simulations. The model captures the short- and long-range effects of water-exclusion in conditions of partial and anisotropic hydration. The model introduces three parameters, one of which represents the degree of hydration (water occupancy) at any point in the system, which depends on the solute conformation, and two that represent the strength of water-water and water-solute dispersive interactions. The model is optimized for proteins, using hydration data of a suboptimally hydrated binding site and results from dynamics simulations in explicit water. The model is applied to a series of aliphatic-alcohol/protein complexes and a set of binary and ternary complexes of various sizes. Implications for weak and ultra-weak protein-protein association and for simulation in crowded media are discussed.

Hyperthermophile protein behavior: partially-structured conformations of Pyrococcus furiosus rubredoxin monomers generated through forced cold-denaturation and refolding

Some years ago, we showed that thermo-chemically denatured, partially-unfolded forms of Pyrococcus furiosus triosephosphateisomerase (PfuTIM) display cold-denaturation upon cooling, and heat-renaturation upon reheating, in proportion with the extent of initial partial unfolding achieved. This was the first time that cold-denaturation was demonstrated for a hyperthermophile protein, following unlocking of surface salt bridges. Here, we describe the behavior of another hyperthermophile protein, the small, monomeric, 53 residues-long rubredoxin from Pyrococcus furiosus (PfRd), which is one of the most thermostable proteins known to man. Like PfuTIM, PfRd too displays cold-denaturation after initial thermo-chemical perturbation, however, with two differences: (i) PfRd requires considerably higher temperatures as well as higher concentrations of guanidium hydrochloride (Gdm.HCl) than PfuTIM; (ii) PfRd's cold-denaturation behavior during cooling after thermo-chemical perturbation is incompletely reversible, unlike PfuTIM's, which was clearly reversible (from each different conformation generated). Differential cold-denaturation treatments allow PfRd to access multiple partially-unfolded states, each of which is clearly highly kinetically-stable. We refer to these as ‘Trishanku’ unfolding intermediates (or TUIs). Fascinatingly, refolding of TUIs through removal of Gdm.HCl generates multiple partially-refolded, monomeric, kinetically-trapped, non-native ‘Trishanku’ refolding intermediates (or TRIs), which differ from each other and from native PfRd and TUIs, in structural content and susceptibility to proteolysis. We find that the occurrence of cold denaturation and observations of TUI and TRI states is contingent on the oxidation status of iron, with redox agents managing to modulate the molecule's behavior upon gaining access to PfRd's iron atom. Mass spectrometric examination provides no evidence of the formation of disulfide bonds, but other experiments suggest that the oxidation status of iron (and its extent of burial) together determine whether or not PfRd shows cold denaturation, and also whether redox agents are able to modulate its behavior.

The accuracy of quantum chemical methods for large noncovalent complexes

We evaluate the performance of the most widely used wavefunction, density functional theory, and semiempirical methods for the description of noncovalent interactions in a set of larger, mostly dispersion-stabilized noncovalent complexes (the L7 data set). The methods tested include MP2, MP3, SCS-MP2, SCS(MI)-MP2, MP2.5, MP2.X, MP2C, DFT-D, DFT-D3 (B3-LYP-D3, B-LYP-D3, TPSS-D3, PW6B95-D3, M06-2X-D3) and M06-2X, and semiempirical methods augmented with dispersion and hydrogen bonding corrections: SCC-DFTB-D, PM6-D, PM6-DH2 and PM6-D3H4. The test complexes are the octadecane dimer, the guanine trimer, the circumcoronene…adenine dimer, the coronene dimer, the guanine-cytosine dimer, the circumcoronene…guanine-cytosine dimer, and an amyloid fragment trimer containing phenylalanine residues. The best performing method is MP2.5 with relative root mean square deviation (rRMSD) of 4 %. It can thus be recommended as an alternative to the CCSD(T)/CBS (alternatively QCISD(T)/CBS) benchmark for molecular systems which exceed current computational capacity. The second best non-DFT method is MP2C with rRMSD of 8 %. A method with the most favorable "accuracy/cost" ratio belongs to the DFT family: BLYP-D3, with an rRMSD of 8 %. Semiempirical methods deliver less accurate results (the rRMSD exceeds 25 %). Nevertheless, their absolute errors are close to some much more expensive methods such as M06-2X, MP2 or SCS(MI)-MP2, and thus their price/performance ratio is excellent.

On the importance and origin of aromatic interactions in chemistry and biodisciplines

Aromatic systems contain both σ- and π-electrons, which in turn constitute σ- and π-molecular orbitals (MOs). In discussing the properties of these systems, researchers typically refer to the highest occupied and lowest unoccupied MOs, which are π MOs. The characteristic properties of aromatic systems, such as their low ionization potentials and electron affinities, high polarizabilities and stabilities, and small band gaps (in spectroscopy called the N → V1 space), can easily be explained based on their electronic structure. These one-electron properties point to characteristic features of how aromatic systems interact with each other. Unlike hydrogen bonding systems, which primarily interact through electrostatic forces, complexes containing aromatic systems, especially aromatic stacked pairs, are predominantly stabilized by dispersion attraction. The stabilization energy in the benzene dimer is rather small (~2.5 kcal/mol) but strengthens with heteroatom substitution. The stacked interaction of aromatic nucleic acid bases is greater than 10 kcal/mol, and for the most stable stacked pair, guanine and cytosine, it reaches approximately 17 kcal/mol. Although these values do not equal the planar H-bonded interactions of these bases (~29 kcal/mol), stacking in DNA is more frequent than H-bonding and, unlike H-bonding, is not significantly weakened when passing from the gas phase to a water environment. Consequently, the stacking of aromatic systems represents the leading stabilization energy contribution in biomacromolecules and in related nanosystems. Therefore stacking (dispersion) interactions predominantly determine the double helical structure of DNA, which underlies its storage and transfer of genetic information. Similarly, dispersion is the dominant contributor to attractive interactions involving aromatic amino acids within the hydrophobic core of a protein, which is critical for folding. Therefore, understanding the nature of aromatic interactions, which depend greatly on quantum mechanical (QM) calculations, is of key importance in biomolecular science. This Account shows that accurate binding energies for aromatic complexes should be based on computations made at the (estimated) CCSD(T)/complete basis set limit (CBS) level of theory. This method is the least computationally intensive one that can give accurate stabilization energies for all common classes of noncovalent interactions (aromatic-aromatic, H-bonding, ionic, halogen bonding, charge-transfer, etc.). These results allow for direct comparison of binding energies between different interaction types. Conclusions based on lower-level QM calculations should be considered with care.

Towards a better understanding of the specificity of protein-protein interaction

In order to predict interaction interface for proteins, it is crucial to identify their characteristic features controlling the interaction process. We present analysis of 69 crystal structures of dimer protein complexes that provides a basis for reasonable description of the phenomenon. Interaction interfaces of two proteins at amino acids level were localized and described in terms of their chemical composition, binding preferences, and residue interaction energies utilizing Amber empirical force field. The characteristic properties of the interaction interface were compared against set of corresponding intramolecular binding parameters for amino acids in proteins. It has been found that geometrically distinct clusters of large hydrophobic amino acids (leucine, valine, isoleucine, and phenylalanine) as well as polar tyrosines and charged arginines are signatures of the protein-protein interaction interface. At some extent, we can generalize that protein-protein interaction (seen through interaction between amino acids) is very similar to the intramolecular arrangement of amino acids, although intermolecular pairs have generally lower interaction energies with their neighbors. Interfaces, therefore, possess high degree of complementarity suggesting also high selectivity of the process. The utilization of our results can improve interface prediction algorithms and improve our understanding of protein-protein recognition.

On the Structure and Geometry of Biomolecular Binding Motifs (Hydrogen-Bonding, Stacking, X-H···π): WFT and DFT Calculations

The strengths of noncovalent interactions are generally very sensitive to a number of geometric parameters. Among the most important of these parameters is the separation between the interacting moieties (in the case of an intermolecular interaction, this would be the intermolecular separation). Most works seeking to characterize the properties of intermolecular interactions are mainly concerned with binding energies obtained at the potential energy minimum (as determined at some particular level of theory). In this work, in order to extend our understanding of these types of noncovalent interactions, we investigate the distance dependence of several types of intermolecular interactions, these are hydrogen bonds, stacking interactions, dispersion interactions, and X-H···π interactions. There are several methods that have traditionally been used to treat noncovalent interactions as well as many new methods that have emerged within the past three or four years. Here we obtain reference data using estimated CCSD(T) values at the complete basis set limit (using the CBS(T) method); potential energy curves are also produced using several other methods thought to be accurate for intermolecular interactions, these are MP2/cc-pVTZ, MP2/aug-cc-pVDZ, MP2/6-31G*(0.25), SCS(MI)-MP2/cc-pVTZ, estimated MP2.5/CBS, DFT-SAPT/aug-cc-pVTZ, DFT/M06-2X/6-311+G(2df,2p), and DFT-D/TPSS/6-311++G(3df,3pd). The basis set superposition error is systematically considered throughout the study. It is found that the MP2.5 and DFT-SAPT methods, which are both quite computationally intensive, produce potential energy curves that are in very good agreement to those of the reference method. Among the MP2 techniques, which can be said to be of medium computational expense, the best results are obtained with MP2/cc-pVTZ and SCS(MI)-MP2/cc-pVTZ. DFT-D/TPSS/6-311++G(3df,3pd) is the DFT-based method that can be said to give the most well-balanced description of intermolecular interactions.

Predicting melting temperature directly from protein sequences

Proteins of both hyperthermophilic and mesophilic microorganisms generally constitute from the same 20 amino acids; however, the extent of thermal tolerance of any given protein is an inherent property of its amino acid sequence. The present study is the first to report a rapid method for predicting Tm (melting temperature), the temperature at which 50% of the protein is unfolded, directly from protein sequences (the Tm Index program is available at http://tm.life.nthu.edu.tw/). We examined 75 complete microbial genomes using the Tm Index, and the analysis clearly differentiated hyperthermophilic from mesophilic microorganisms on this global genomic basis. These results are consistent with the previous hypothesis that hyperthermophiles express a greater number of high Tm proteins compared with mesophiles. The Tm Index will be valuable for modifying existing proteins (enzymes, protein drugs and vaccines) or designing novel proteins having a desired melting temperature.

Molecular design of protein-based nanocapsules for stimulus-responsive characteristics

Hsp16.5, a small heat-shock protein (sHSP) from hyperthermophilic archaeon, forms a homogeneous complex comprised of 24 subunits with a molecular mass of 400 kDa. This complex self-organizes under physiological conditions, and the structure of the complex is a nanoscale spherical capsule with small pores. Furthermore, this natural nanocapsule exhibits very high thermal stability. In this paper, we functionalized the nanocapsule to control the structure in response to external stimuli such as a protease signal and temperature. For this purpose, several mutations (Mut1-10) to create a cleavage site for a specific protease, Factor Xa, were introduced on the outer surface of the nanocapsule using a genetic engineering strategy. The resulting mutants were expressed to high levels in Escherichia coli. One of these mutants, Mut6, which has the most accessible cleavage site located at the triangular pore on the surface of the capsule, formed a spherical assembly similar to that observed for the wild-type protein. Mut6 showed the highest sensitivity to Factor Xa, and the structure of the protease digested Mut6 disassembled irreversibly after heating. In contrast, the nanocapsule comprising the wild-type Hsp16.5 was not influenced by the dual stimuli. These results suggest that Mut6 acts as a stimulus-responsive nanocapsule. Such a characteristic of the protein-based nanocapsule has attractive potential as a versatile intelligent system.

Identifying stabilizing key residues in proteins using interresidue interaction energy matrix

We are proposing an interresidue interaction energy map (IEM)--a new tool for protein structure analysis and protein bioinformatics. This approach employs the sum of pair-wise interaction energies of a particular residue as a measure of its structural importance. We will show that the IEM can serve as a means for identifying key residues responsible for the stability of a protein. Our method can be compared with the interresidue contact map but has the advantage of weighting the contacts by the stabilization energy content which they bring to the protein structure. For the theoretical adjustment of the proposed method, we chose the Trp-cage mini protein as a model system to compare a spectrum of computational methods ranging from the ab initio MP2 level through the DFT method to empirical force-field methods. The IEM method correctly identifies Tryptophane 6 as the key residue in the Trp-cage. The other residues with the highest stabilizing contributions correspond to the structurally important positions in the protein. We have further tested our method on the Trp2Cage miniprotein--a P12W mutant of the Trp-cage and on two proteins from the rubredoxin family that differ in their thermostability. Our method correctly identified the thermodynamically more stable variants in both cases and therefore can also be used as a tool for the relative measurement of protein stability. Finally, we will point out the important role played by dispersion energy, which contributes significantly to the total stabilization energy and whose role in aromatic pairs is clearly dominant. Surprisingly, the dispersion energy plays an even more important role in the interaction of prolines with aromatic systems.