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.

The ribosome is a macromolecular machine that produces proteins in all kingdoms of life. The proteins, in turn, control the biochemical processes within the cell. It is thus of extreme importance that the machine that makes the proteins works with high precision. By using three dimensional structures of the ribosome and homology modelling, we have applied molecular dynamics simulations and free-energy calculations to study the codon specificity of protein synthesis in initiation and termination on an atomistic level. In addition, we have examined the binding of small molecules to riboswitches, which can change the expression of an mRNA. The relative affinities on the ribosome between the eukaryotic initiator tRNA to the AUG start codon and six near-cognate codons were determined. The free-energy calculations show that the initiator tRNA has a strong preference for the start codon, but requires assistance from initiation factors 1 and 1A to uphold discrimination against near-cognate codons. When instead a stop codon (UAA, UGA or UAG) is positioned in the ribosomal A-site, a release factor binds and terminates protein synthesis by hydrolyzing the nascent peptide chain. However, vertebrate mitochondria have been thought to have four stop codons, namely AGA and AGG in addition to the standard UAA and UAG codons. Furthermore, two release factors have been identified, mtRF1 and mtRF1a. Free-energy calculations were used to determine if any of these two factors could bind to the two non-standard stop codons, and thereby terminate protein synthesis. Our calculations showed that the mtRF’s have similar stop codon specificity as bacterial RF1 and that it is highly unlikely that the mtRF’s are responsible for terminating at the AGA and AGG stop codons. The eukaryotic release factor 1, eRF1, on the other hand, can read all three stop codons singlehandedly. We show that eRF1 exerts a high discrimination against near-cognate codons, while having little preference for the different cognate stop codons. We also found an energetic mechanism for avoiding misreading of the UGG codon and could identify a conserved cluster of hydrophobic amino acids which prevents excessive solvent molecules to enter the codon binding site. The linear interaction energy method was used to examine binding of small molecules to the purine riboswitch and the FEP method was employed to explicitly calculate the LIE b-parameters. We show that the purine riboswitches have a remarkably high degree of electrostatic preorganization for their cognate ligands which is fundamental for discriminating against different purine analogs.