Basis set superposition error-counterpoise corrected potential energy surfaces : application to hydrogen peroxide ... X (X=F-, Cl-, Br-, Li+, Na+) complexes

Møller-Plesset (MP2) and Becke-3-Lee-Yang-Parr (B3LYP) calculations have been used to compare the geometrical parameters, hydrogen-bonding properties, vibrational frequencies and relative energies for several X- and X+ hydrogen peroxide complexes. The geometries and interaction energies were corrected for the basis set superposition error (BSSE) in all the complexes (1-5), using the full counterpoise method, yielding small BSSE values for the 6-311 + G(3df,2p) basis set used. The interaction energies calculated ranged from medium to strong hydrogen-bonding systems (1-3) and strong electrostatic interactions (4 and 5). The molecular interactions have been characterized using the atoms in molecules theory (AIM), and by the analysis of the vibrational frequencies. The minima on the BSSE-counterpoise corrected potential-energy surface (PES) have been determined as described by S. Simón, M. Duran, and J. J. Dannenberg, and the results were compared with the uncorrected PES

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

Destabilizing Mutations Alter the Hydrogen Exchange Mechanism in Ribonuclease A

The effect of strongly destabilizing mutations, I106A and V108G of Ribonuclease A (RNase A), on its structure and stability has been determined by NMR. The solution structures of these variants are essentially equivalent to RNase A. The exchange rates of the most protected amide protons in RNase A (35ºC), the I106A variant (35ºC), and the V108G variant (10ºC) yield stability values of 9.9, 6.0, and 6.8 kcal/mol, respectively, when analyzed assuming an EX2 exchange mechanism. Thus, the destabilization induced by these mutations is propagated throughout the protein. Simulation of RNase A hydrogen exchange indicates that the most protected protons in RNase A and the V108G variant exchange via the EX2 regime, whereas those of I106A exchange through a mixed EX1 1 EX2 process. It is striking that a single point mutation can alter the overall exchange mechanism. Thus, destabilizing mutations joins high temperatures, high pH and the presence of denaturating agents as a factor that induces EX1 exchange in proteins. The calculations also indicate a shift from the EX2 to the EX1 mechanism for less protected groups within the same protein. This should be borne in mind when interpreting exchange data as a measure of local stability in less protected regions