Force Field Modeling of Amino Acid Conformational Energies

Effect of cation complexation on the structure of a conformationally flexible multiply charged anion: stabilization of excess charge in the Na+ x adenosine 5'-triphosphate dianion ion-pair complex

We report a computational study of the conformationally and tautomerically flexible cation-dianion complex of Na(+) with doubly deprotonated adenosine 5'-triphosphate (ATP) using a hierarchical selection method. The method uses molecular dynamics to generate initial conformeric structures, followed by a classification process that groups conformers into five "families" to ensure that a representative sample of structures is retained for further analysis, while very similar conformational structures are eliminated. Hierarchical ab initio calculations (DFT and MP2) of typical conformers of the families are then performed to identify the lowest-energy conformeric structures. The procedure described should provide a useful methodology for conducting higher-level ab initio calculations of medium-sized gas-phase biological molecules for interpreting contemporary laser spectroscopy measurements. For Na(+) x [ATP-2H](2) (considering tautomers where the phosphate chain of ATP is doubly deprotonated), the calculations reveal that the sodium cation interacts directly with the negatively charged phosphates (maximum distance = 2.54 A) in all of the low-energy conformers, while a number of the structures also display close cation-adenine interactions producing compact ball-like structures. These compact structures generally correspond to the lowest-energy conformers. The structural variation between the bare [ATP-2H](2-) molecular ion (Burke et al. J. Phys. Chem. A 2005 , 109 , 9775-9785) and the Na(+) x [ATP-2H](2-) cluster is discussed in detail, including the effect of sodiation on the intramolecular hydrogen-bonding network within ATP in a gas-phase environment.

Effect of polarization on the opsin shift in rhodopsins. 2. Empirical polarization models for proteins

The explicit treatment of polarization as a many-body interaction in condensed-phase systems represents a current problem in empirical force-field development. Although a variety of efficient models for molecular polarization have been suggested, polarizable force fields are still far from common use nowadays. In this work, we consider interactive polarization models employing Thole's short-range damping scheme and assess them for application on polypeptides. Despite the simplicity of the model, we find mean polarizabilities and anisotropies of amino acid side chains in excellent agreement with MP2/cc-pVQZ benchmark calculations. Combined with restrained electrostatic potential (RESP) derived atomic charges, the models are applied in a quantum-mechanical/molecular-mechanical (QM/MM) approach. An iterative scheme is used to establish a self-consistent mutual polarization between the QM and MM moieties. This ansatz is employed to study the influence of the protein polarizability on calculated optical properties of the protonated Schiff base of retinal in rhodopsin (Rh), bacterio-rhodopsin (bR), and pharaonis sensory rhodopsin II (psRII). The shifts of the excitation energy due to the instantaneous polarization response of the protein to the charge transfer on the retinal chromophore are quantified using the high level ab initio multireference spectroscopy-oriented configuration interaction (SORCI) method. The results are compared with those of previously published QM1/QM2/MM models for bR and psRII.