Computer simulations of bio-molecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in bio-molecular systems and are therein described by atomic point charges. In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute’s electron density computed with an implicit solvent model and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the AM1-BCC and the MBIS atomic charge methods. The latter includes the solvent polarization and present a root mean square error of 2.0 kcal mol−1 for the 613 organic molecules studied. The largest deviation was observed for phosphorus-containing molecules and the molecules with amide, ester and amine functional groups.

No abstract

Phycocyanin is a phycobiliprotein involved in light harvesting and conduction of light to the reaction centers in cyanobacteria and red algae. The structure of C-phycocyanin from Gracilaria chilensis was solved by X-ray crystallography at 2.0 Å resolution in space group P2 1 . An interaction model between two PC heterohexamers was built, followed by molecular dynamic refinement. The best model showed an inter-hexamer rotation of 23°. The coordinates of a PC heterohexamer (αβ) 6 and of the PC-PC complex were used to perform energy transfer calculations between chromophores pairs using the fluorescence resonance energy transfer approach (FRET). Two main intra PC (

ABSTRACT:Kratzer oscillator is a realistic zero-order model for describing the anharmonic ro-vibrational motion in diatomic molecules. Kratzer oscillator has an energy spectrum containing both discrete and continuum parts. Wavefunctions belonging to the continuum would be useful in the study of transitions to the continuum in molecular dissociation processes. In this article, bound and scattering wavefunctions of the Kratzer oscillator are reviewed and the bound-bound and the bound-free matrix elements are obtained.

A semiempirical methodology to model the intra-phycocyanin and inter-phycocyanin fluorescence resonance energy-transfer (FRET) pathways in the rods of the phycobilisomes (PBSs) from Fremyella diplosiphon is presented. Using the Förster formulation of FRET and combining experimental data and PM3 calculation of the dipole moments of the aromatic portions of the chromophores, transfer constants between pairs of chromophores in the phycocyanin (PC) structure were obtained. Protein docking of two PC hexamers was used to predict the optimal distance and axial rotation angle for the staked PCs in the PBSs' rods. Using the distance obtained by the docking process, transfer constants between pairs of chromophores belonging to different PC hexamers were calculated as a function of the angle of rotation. We show that six preferential FRET pathways within the PC hexameric ring and 15 pathways between hexamers exist, with transfer constants consistent with experimental results. Protein docking predicted the quaternary structure for PCs in rods with inter-phycocyanin distance of 55.6 A and rotation angle of 20.5 degrees . The inter-phycocyanin FRET constant between chromophores at positions beta(155) is maximized at the rotation angle predicted by docking revealing the crucial role of this specific inter-phycocyanin channel in defining the complete set of FRET pathways in the system.