Protein:Ligand binding free energies: A stringent test for computational protein design.

A computational protein design method is extended to allow Monte Carlo simulations where two ligands are titrated into a protein binding pocket, yielding binding free energy differences. These provide a stringent test of the physical model, including the energy surface and sidechain rotamer definition. As a test, we consider tyrosyl-tRNA synthetase (TyrRS), which has been extensively redesigned experimentally. We consider its specificity for its substrate l-tyrosine (l-Tyr), compared to the analogs d-Tyr, p-acetyl-, and p-azido-phenylalanine (ac-Phe, az-Phe). We simulate l- and d-Tyr binding to TyrRS and six mutants, and compare the structures and binding free energies to a more rigorous "MD/GBSA" procedure: molecular dynamics with explicit solvent for structures and a Generalized Born + Surface Area model for binding free energies. Next, we consider l-Tyr, ac- and az-Phe binding to six other TyrRS variants. The titration results are sensitive to the precise rotamer definition, which involves a short energy minimization for each sidechain pair to help relax bad contacts induced by the discrete rotamer set. However, when designed mutant structures are rescored with a standard GBSA energy model, results agree well with the more rigorous MD/GBSA. As a third test, we redesign three amino acid positions in the substrate coordination sphere, with either l-Tyr or d-Tyr as the ligand. For two, we obtain good agreement with experiment, recovering the wildtype residue when l-Tyr is the ligand and a d-Tyr specific mutant when d-Tyr is the ligand. For the third, we recover His with either ligand, instead of wildtype Gln.

Sequence requirements for terminators and antiterminators in the T box transcription antitermination system: disparity between conservation and functional requirements.

The T box transcription termination control system is used in Gram-positive bacteria to regulate expression of aminoacyl-tRNA synthetase and other amino acid-related genes. Readthrough of a transcriptional terminator located in the leader region of the target gene is dependent on a specific interaction between the nascent leader transcript and the cognate uncharged tRNA. This interaction is required for formation of an antiterminator structure in the leader, which prevents formation of a competing transcriptional terminator stem-loop. The antiterminators and terminators of genes in this family are highly conserved in both secondary structure and primary sequence; the antiterminator contains the T box sequence, which is the most highly conserved leader element. These conserved features were investigated by phylogenetic and mutational analysis. Changes at highly conserved positions in the bulge and in the helix above the bulge reduced function, while alteration of other positions that were as much as 96% conserved did not have a major effect. The disparity between sequence conservation and function may be due to the requirement for maintaining base pairing in both the antiterminator and terminator structures.