Selectivity and specificity of substrate binding in methionyl-tRNA synthetase

The accuracy of in vivo incorporation of amino acids during protein biosynthesis is controlled to a significant extent by aminoacyl-tRNA synthetases (aaRS). This paper describes the application of the HierDock computational method to study the molecular basis of amino acid binding to the Escherichia coli methionyl tRNA synthetase (MetRS). Starting with the protein structure from the MetRS cocrystal, the HierDock calculations predict the binding site of methionine in MetRS to a root mean square deviation in coordinates (CRMS) of 0.55 A for all the atoms, compared with the crystal structure. The MetRS conformation in the cocrystal structure shows good discrimination between cognate and the 19 noncognate amino acids. In addition, the calculated binding energies of a set of five methionine analogs show a good correlation (R(2) = 0.86) to the relative free energies of binding derived from the measured in vitro kinetic parameters, K(m) and k(cat). Starting with the crystal structure of MetRS without the methionine (apo-MetRS), the putative binding site of methionine was predicted. We demonstrate that even the apo-MetRS structure shows a preference for binding methionine compared with the 19 other natural amino acids. On comparing the calculated binding energies of the 20 natural amino acids for apo-MetRS with those for the cocrystal structure, we observe that the discrimination against the noncognate substrate increases dramatically in the second step of the physical binding process associated with the conformation change in the protein.     

Test of the Binding Threshold Hypothesis for olfactory receptors: explanation of the differential binding of ketones to the mouse and human orthologs of olfactory receptor 912-93

We tested the Binding Threshold Hypothesis (BTH) for activation of olfactory receptors (ORs): To activate an OR, the odorant must bind to the OR with binding energy above some threshold value. The olfactory receptor (OR) 912-93 is known experimentally to be activated by ketones in mouse, but is inactive to ketones in human, despite an amino acid sequence identity of approximately 66%. To investigate the origins of this difference, we used the MembStruk first-principles method to predict the tertiary structure of the mouse OR 912-93 (mOR912-93), and the HierDock first-principles method to predict the binding site for ketones to this receptor. We found that the strong binding of ketones to mOR912-93 is dominated by a hydrogen bond of the ketone carbonyl group to Ser105. All ketones predicted to have a binding energy stronger than EBindThresh = 26 kcal/mol were observed experimentally to activate this OR, while the two ketones predicted to bind more weakly do not. In addition, we predict that 2-undecanone and 2-dodecanone both bind sufficiently strongly to activate mOR912-93. A similar binding site for ketones was predicted in hOR912-93, but the binding is much weaker because the human ortholog has a Gly at the position of Ser105. We predict that mutating this Gly to Ser in human should lead to activation of hOR912-93 by these ketones. Experimental substantiations of the above predictions would provide further tests of the validity of the BTH, our predicted 3D structures, and our predicted binding sites for these ORs.     

Making sense of olfaction through predictions of the 3-D structure and function of olfactory receptors

We used the MembStruk first principles computational technique to predict the three-dimensional (3-D) structure of six mouse olfactory receptors (S6, S18, S19, S25, S46 and S50) for which experimental odorant recognition profiles are available for a set of 24 odorants (4-9 carbons aliphatic alcohols, acids, bromo-acids and diacids). We used the HierDock method to scan each predicted OR structure for potential odorant binding site(s) and to calculate binding energies of each odorant in these binding sites. The calculated binding affinity profiles are in good agreement with experimental activation profiles, validating the predicted 3-D structures and the predicted binding sites. For each of the six ORs, the binding site is located between trans-membrane domains (TMs) 3-6, with contributions from extracellular loops 2 and 3. In particular, we find six residue positions in TM3 and TM6 to be consistently involved in the binding modes of the odorants. Indeed, the differences in the experimental recognition profiles can be explained on the basis of these critical residues alone. These predictions are also consistent with mutation data on ligand binding for catecholamine receptors and sequence hypervariability studies for ORs. Based on this analysis, we defined amino acid patterns associated with the recognition of short aliphatic alcohols and mono-acids. Using these two sequence fingerprints to probe the alignment of 869 OR sequences from the mouse genome, we identified 34 OR sequences matching the fingerprint for aliphatic mono-acids and 36 corresponding to the recognition pattern for aliphatic alcohols. We suggest that these two sets of ORs might function as basic arrays for uniquely recognizing aliphatic alcohols and acids. We screened a library of 89 additional molecules against the six ORs and found that this set of ORs is likely to respond to aldehydes and esters with longer carbon chains than their currently known agonists. We also find that compounds associated with the flavor in foods are often among the best calculated binding affinities. This suggests that physiologic ligands for these ORs may be found among aldehydes and esters associated with flavor.