Atoms-in-molecules study of the genetically encoded amino acids. II. Computational study of molecular geometries

The geometries of the 20 genetically encoded amino acids were optimized at the restricted Hartree-Fock level of theory using the 6-31+G* basis set. A detailed comparison showed the calculated geometries to be in excellent agreement with those determined by X-ray crystallography. The study demonstrated that the geometric parameters for the main-chain group and for the bonds and common functional groups of the side-chains exhibit a high degree of transferability among the members of this set of molecules. This geometric transferability is a necessary prerequisite for the corresponding transferability of their electron density distributions and hence of their bond and atomic properties. The transferability of the electron distributions will be demonstrated and exploited in the following paper of this series, which uses the topology of the electron density to define an atom within the quantum theory of atoms in molecules. Particular features of the geometries of the amino acids are discussed. It has been shown, for example, how the apparent anomaly of the Calpha-N bond length in a peptide being shorter than in the charged species Calpha-NH3+ is resolved when the charge separation is gauged by the differences in the charges of the Calpha and N atoms as opposed to the use of formal charges. A compilation of literature sources on experimental geometries covering each member of the 20 amino acids is presented. A set of rules for labeling the atoms and bonds, complementing the generally accepted IUPAC-IUB rules, is proposed to uniquely identify every atom and bond in the amino acids.     

Molecular dynamics study of the primary charge separation reactions in Photosystem I: Effect of the replacement of the axial ligands to the electron acceptor A0

Molecular dynamics (MD) calculations, a semi-continuum (SC) approach, and quantum chemistry (QC) calculations were employed together to investigate the molecular mechanics of ultrafast charge separation reactions in Photosystem I (PS I) of Thermosynechococcus elongatus. A molecular model of PS I was developed with the aim to relate the atomic structure with electron transfer events in the two branches of cofactors. A structural flexibility map of PS I was constructed based on MD simulations, which demonstrated its rigid hydrophobic core and more flexible peripheral regions. The MD model permitted the study of atomic movements (dielectric polarization) in response to primary and secondary charge separations, while QC calculations were used to estimate the direct chemical effect of the A_0A/A_0B ligands (Met or Asn in the 688/668 position) on the redox potential of chlorophylls A_0A/A_0B and phylloquinones A_1A/A_1B. A combination of MD and SC approaches was used to estimate reorganization energies λ of the primary (λ₁) and secondary (λ₂) charge separation reactions, which were found to be independent of the active branch of electron transfer; in PS I from the wild type, λ₁ was estimated to be 390 ± 20 mV, while λ₂ was estimated to be higher at 445 ± 15 mV. MD and QC approaches were used to describe the effect of substituting Met688_PsaA/Met668_PsaB by Asn688_PsaA/Asn668_PsaB on the energetics of electron transfer. Unlike Met, which has limited degrees of freedom in the site, Asn was found to switch between two relatively stable conformations depending on cofactor charge. The introduction of Asn and its conformation flexibility significantly affected the reorganization energy of charge separation and the redox potentials of chlorophylls A_0A/A_0B and phylloquinones A_1A/A_1B, which may explain the experimentally observed slowdown of secondary electron transfer in the M688N_PsaA variant. This article is part of a Special Issue entitled: Photosynthesis research for sustainability: Keys to produce clean energy.