From Structure to Function with Binding Free Energy Calculations for Codon Reading, Riboswitches and Lectins

Abstract: Molecular association is part of many important processes in living cells. Computational methods for calculating binding free energies allows for a quantitative examination of biomolecular structures and hypotheses drawn from biochemical experiments. Here, binding free energy calculations for tRNAs and release factors binding to mRNA codons on the ribosome, sugars binding to lectins and purine analogs binding to the purine riboswitch are presented.The relative affinities between cognate and non-cognate tRNAs for different states involved in codon reading on the ribosome were determined. The calculations show that tRNA discrimination varies between different conformations of the 30S subunit, where the existence of both low and high selectivity states provides an efficient common mechanism for initial selection and proofreading. The simulations reveal a desolvation mechanism for the 30S conformational switch with which the accuracy of peptide bond formation can be amplified.When an mRNA stop codon (UAA, UAG or UGA) is located in the ribosomal A-site release factors bind to the ribosome and the synthesized protein is released. RF1 is specific for UAA and UAG whereas RF2 is specific for UAA and UGA. The free energy calculations and an analysis of the performed simulations show the mechanisms for how RF1 and RF2 are able to read the stop codons with different specificities. Also mitochondrial release factors were investigated. Vertebrate mitochondria have four stop codons, UAA, UAG, AGA and AGG and two release factors mtRF1 and mtRF1a. The calculations show how the specificities of both mtRF1 and mtRF1a agree with RF1 and that none of them are likely to read the non-standard stop codons AGA and AGG.The linear interaction energy method has also been examined for the RSL and PA-IIL lectins and for the purine riboswitch. The standard parameterization of the method works well for RSL, but fails for PA-IIL and the purine riboswitch due to compositions of the active sites in these systems. The development of new parameterizations to overcome these problems leads to a better understanding of both the method and the binding mechanisms in these systems.