Molecular recognition in proteins. Simulation analysis of substrate binding by a tyrosyl-tRNA synthetase mutant.

Alchemical molecular dynamics simulations are performed to determine the difference in the free energy of binding of the tyrosine substrate between the wild type of tyrosyl-tRNA synthetase (TyrRS) from Bacillus stearothermophilus and the mutant Tyr169-->Phe. The results are of general interest because the Tyr169 hydroxyl group interacts with the ammonium group of the substrate in a manner corresponding to that found in other amino acid binding proteins (e.g. the Asp receptor of the chemotactic bacterium Salmonella typhimurium and class I major histocompatibility complex molecules). The calculated free-energy change due to the Tyr169-->Phe mutation is 3.4 kcal/mol (the statistical error is +/- 0.5 kcal/mol) in satisfactory agreement with the experimental value of 3(+/- 0.5) kcal/mol. By use of thermodynamic integration, the contribution of the different terms to the free energy change are estimated. The path dependence of such a decomposition is discussed and it is suggested that the alchemical choice is of primary interest for understanding the interactions involved. There are large protein contributions to the alchemical free energy difference of the bound and free enzyme that cancel in the overall result. Due to this cancellation, the essential interactions contributing to the free-energy change are those between the OH group of Tyr169 and water in the free enzyme and those between the OH group of Tyr169 and the ammonium group of the substrate in the bound system. The results thus support simple models based on a balance of hydrogen bonding interactions.

All-atom empirical potential for molecular modeling and dynamics studies of proteins.

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent-solvent, solvent-solute, and solute-solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volumes, and crystal pressures and structures were used in the optimization. The resulting protein parameters were tested by applying them to noncyclic tripeptide crystals, cyclic peptide crystals, and the proteins crambin, bovine pancreatic trypsin inhibitor, and carbonmonoxy myoglobin in vacuo and in crystals. A detailed analysis of the relationship between the alanine dipeptide potential energy surface and calculated protein φ, χ angles was made and used in optimizing the peptide group torsional parameters. The results demonstrate that use of ab initio structural and energetic data by themselves are not sufficient to obtain an adequate backbone representation for peptides and proteins in solution and in crystals. Extensive comparisons between molecular dynamics simulations and experimental data for polypeptides and proteins were performed for both structural and dynamic properties. Energy minimization and dynamics simulations for crystals demonstrate that the latter are needed to obtain meaningful comparisons with experimental crystal structures. The presented parameters, in combination with the previously published CHARMM all-atom parameters for nucleic acids and lipids, provide a consistent set for condensed-phase simulations of a wide variety of molecules of biological interest.

Dissection of the effector-binding site and complementation studies of Escherichia coli phosphofructokinase using site-directed mutagenesis.

A systematic study by site-directed mutagenesis has been conducted on the effector site of phosphofructokinase from Escherichia coli to delineate the role of side chains in binding the allosteric activator, GDP, and inhibitor, PEP, and to search for key residues in the allosteric transtion. Target residues were identified from the crystal structure of the enzyme-nucleoside diphosphate complex. It is found that both activator and inhibitor bind to the same set of amino acid side chains. Deletion of positively charged groups (Arg21, Arg25, Arg54, Arg154, and Lys213 mutated to alanine) weakens binding of both effectors by 2-3 kcal/mol, consistent with the disruption of charged hydrogen bonds. Residue Glu187, which is known from the crystal structure to bind the coordinated Mg2+ ion of GDP, is found to have a unique behavior on mutation and appears to be crucial in triggering the allosteric transition. All other residues mutated simply weaken binding of both PEP and GDP in a parallel manner. However, mutation of Glu----Ala187 reverses the roles of GDP and PEP, causing GDP to become an allosteric inhibitor and PEP an activator. Mutation of Glu----Gln187 has only a small effect on the binding of PEP, and both PEP and GDP are inhibitors. Studies are described in which mutations in different subunits of a tetrameric complex complement each other. The effector site is composed of residues from two subunits. In particular, Arg21 and Lys213 in each site are from different subunits. Mutations of either one of these residues abolishes activation by GDP of the homotetramer.(ABSTRACT TRUNCATED AT 250 WORDS)

Conversion of allosteric inhibition to activation in phosphofructokinase by protein engineering.

Many enzymes are subject to allosteric control, often with inhibitors and activators binding to the same effector site. Phosphofructokinase in Escherichia coli is such an enzyme, being inhibited by phosphoenolpyruvate (PEP) and activated by ADP and GDP. How do individual interactions with effectors affect the balance between activation and inhibition, especially when both ligands share aspects of the same binding site? We find that mutation of a single residue in the effector site, Glu----Ala 187, leads to PEP being an activator rather than an inhibitor. With low concentrations of the substrate fructose-6-phosphate, the mutant enzyme is more than one hundred times more active than wild-type enzyme at millimolar concentrations of PEP. The classical Monod-Wyman-Changeux two-state model is too simple to account for the properties of the mutant enzyme.

Insertion of an elastase-binding loop into interleukin-1 beta.

The protease-binding sequence EAIPMSIPPE from alpha 1-antitrypsin has been inserted into the cytokine interleukin-1 beta, replacing residues 50-53. The resulting mutant protein was cleaved specifically at a single site by elastase and chymotrypsin, but not by trypsin. The cleavage by elastase was shown to be between Met and Ser of the inserted loop. In contrast, wild-type interleukin is not susceptible to cleavage by any of these enzymes. The mutant protein acts as an inhibitor of elastase, with a KI of approximately 30 microM. The wild type displays no such inhibitory activity. The overall structure of the mutant, as demonstrated by CD, appears to be indistinguishable from that of the wild type. These results indicate that the protease-binding region of alpha 1-antitrypsin can be recognized and is active even within the context of an entirely different protein structure. Given that interleukin-1 beta binds to, and is internalized by, many types of cells, this hybrid protein also demonstrates the feasibility of using interleukin-1 beta as a delivery system for useful therapeutic agents.

Sequential assignment of the 1H nuclear magnetic resonance spectrum of barnase.

Two-dimensional nuclear magnetic resonance spectroscopy has been used to study the bacterial ribonuclease barnase (MW 12,382). Resonance assignments have been made for protons in all of the 110 residues. Analysis of medium- and long-range contacts in NOESY spectra has demonstrated that the major elements of secondary structure in barnase in solution are essentially identical with those found in the crystal structure.

Site-directed mutagenesis in the effector site of Escherichia coli phosphofructokinase.

A new vector for the expression of phosphofructokinase (pfk-1) was constructed with pEMBL, which allows reliable, inducible, high-expression, and facile mutagenesis of the gene. Two mutants in the effector site of the enzyme were produced by site-specific mutagenesis of residue Tyr-55 to assess the role of its side chain in binding an allosteric inhibitor, phosphoenolpyruvate (PEP), and an activator, guanosine 5'-diphosphate (GDP): Tyr-55----Phe-55 and Try-55----Gly-55. The dissociation constant of PEP from the T state is unaffected by the mutations. Mutation of Tyr-55----Phe-55 only slightly increases the dissociation constant of GDP from the R state, indicating a minimal involvement of the hydroxyl group in binding. A 5.5-fold increase in the dissociation constant of GDP on the mutation of Tyr-55----Gly-55 suggests a small hydrophobic interaction of the aromatic ring of the tyrosine residue with guanine of GDP.