Search PubMed⌕ Search

Biomedical subjects

B Tidor

Publications and source records attributed to B Tidor.

At least 19 recordsLinked to original sources
PubMed2001

Side-chain repacking calculations for predicting structures and stabilities of heterodimeric coiled coils.

An important goal in biology is to predict from sequence data the high-resolution structures of proteins and the interactions that occur between them. In this paper, we describe a computational approach that can make these types of predictions for a series of coiled-coil dimers. Our method comprises a dual strategy that augments extensive conformational sampling with molecular mechanics minimization. To test the performance of the method, we designed six heterodimeric coiled coils with a range of stabilities and solved x-ray crystal structures for three of them. The stabilities and structures predicted by the calculations agree very well with experimental data: the average error in unfolding free energies is <1 kcal/mol, and nonhydrogen atoms in the predicted structures superimpose onto the experimental structures with rms deviations <0.7 A. We have also tested the method on a series of homodimers derived from vitellogenin-binding protein. The predicted relative stabilities of the homodimers show excellent agreement with previously published experimental measurements. A critical step in our procedure is to use energy minimization to relax side-chain geometries initially selected from a rotamer library. Our results show that computational methods can predict interaction specificities that are in good agreement with experimental data.

Amino Acid Sequence↗
PubMed2001

Altering dimerization specificity by changes in surface electrostatics.

Arc repressor forms a homodimer in which the subunits intertwine to create a single globular domain. To obtain Arc sequences that fold preferentially as heterodimers, variants with surface patches of excess positive or negative charge were designed. Several but not all oppositely charged sequence pairs showed preferential heterodimer formation. In the most successful design pair, alpha helix B of one subunit contained glutamic acids at positions 43, 46, 47, 48, and 50, whereas the other subunit contained lysines or arginines at these positions. A continuum electrostatic model captures many features of the experimental results and suggests that the most successful designs include elements of both positive and negative design.

Crystallography, X-Ray↗
PubMed2001

Barstar is electrostatically optimized for tight binding to barnase.

We used a novel charge optimization technique to study the small ribonuclease barnase and to analyze its interaction with a natural tight binding inhibitor, the protein barstar. The approach uses a continuum model to explicitly determine the charge distributions that lead to the most favorable electrostatic contribution to binding when competing desolvation and interaction effects are included. Given its backbone fold, barstar is electrostatically optimized for tight binding to barnase when compared with mutants where residues have been substituted with one of the 20 common amino acids. Natural proteins thus appear to use optimization of electrostatic interactions as one strategy for achieving tight binding.

Bacterial Proteins↗
PubMed2001

Optimization of binding electrostatics: charge complementarity in the barnase-barstar protein complex.

Theoretical and experimental studies have shown that the large desolvation penalty required for polar and charged groups frequently precludes their involvement in electrostatic interactions that contribute strongly to net stability in the folding or binding of proteins in aqueous solution near room temperature. We have previously developed a theoretical framework for computing optimized electrostatic interactions and illustrated use of the algorithm with simplified geometries. Given a receptor and model assumptions, the method computes the ligand-charge distribution that provides the most favorable balance of desolvation and interaction effects on binding. In this paper the method has been extended to treat complexes using actual molecular shapes. The barnase-barstar protein complex was investigated with barnase treated as a target receptor. The atomic point charges of barstar were varied to optimize the electrostatic binding free energy. Barnase and natural barstar form a tight complex (K(d) approximately 10(-14) M) with many charged and polar groups near the interface that make this a particularly relevant system for investigating the role of electrostatic effects on binding. The results show that sets of barstar charges (resulting from optimization with different constraints) can be found that give rise to relatively large predicted improvements in electrostatic binding free energy. Principles for enhancing the effect of electrostatic interactions in molecular binding in aqueous environments are discussed in light of the optima. Our findings suggest that, in general, the enhancements in electrostatic binding free energy resulting from modification of polar and charged groups can be substantial. Moreover, a recently proposed definition of electrostatic complementarity is shown to be a useful tool for examining binding interfaces. Finally, calculational results suggest that wild-type barstar is closer to being affinity optimized than is barnase for their mutual binding, consistent with the known roles of these proteins.

Algorithms↗
PubMed2000

Rational modification of protein stability by the mutation of charged surface residues.

Continuum methods were used to calculate the electrostatic contributions of charged and polar side chains to the overall stability of a small 41-residue helical protein, the peripheral subunit-binding domain. The results of these calculations suggest several residues that are destabilizing, relative to hydrophobic isosteres. One position was chosen to test the results of these calculations. Arg8 is located on the surface of the protein in a region of positive electrostatic potential. The calculations suggest that Arg8 makes a significant, unfavorable electrostatic contribution to the overall stability. The experiments described in this paper represent the first direct experimental test of the theoretical methods, taking advantage of solid-phase peptide synthesis to incorporate approximately isosteric amino acid substitutions. Arg8 was replaced with norleucine (Nle), an amino acid that is hydrophobic and approximately isosteric, or with alpha-amino adipic acid (Aad), which is also approximately isosteric but oppositely charged. In this manner, it is possible to isolate electrostatic interactions from the effects of hydrophobic and van der Waals interactions. Both Arg8Nle and Arg8Aad are more thermostable than the wild-type sequence, testifying to the validity of the calculations. These replacements led to stability increases at 52.6 degrees C, the T(m) of the wild-type, of 0.86 and 1.08 kcal mol(-)(1), respectively. The stability of Arg8Nle is particularly interesting as a rare case in which replacement of a surface charge with a hydrophobic residue leads to an increase in the stability of the protein.

Acetyltransferases↗
PubMed1999

Electrostatic interactions in the GCN4 leucine zipper: substantial contributions arise from intramolecular interactions enhanced on binding.

The GCN4 leucine zipper is a peptide homodimer that has been the subject of a number of experimental and theoretical investigations into the determinants of affinity and specificity. Here, we utilize this model system to investigate electrostatic effects in protein binding using continuum calculations. A particularly novel feature of the computations made here is that they provide an interaction-by-interaction breakdown of the electrostatic contributions to the free energy of docking that includes changes in the interaction of each functional group with solvent and changes in interactions between all pairs of functional groups on binding. The results show that (1) electrostatic effects disfavor binding by roughly 15 kcal/mol due to desolvation effects that are incompletely compensated in the bound state, (2) while no groups strongly stabilize binding, the groups that are most destabilizing are charged and polar side chains at the interface that have been implicated in determining binding specificity, and (3) attractive intramolecular interactions (e.g., backbone hydrogen bonds) that are enhanced on binding due to reduced solvent screening in the bound state contribute significantly to affinity and are likely to be a general effect in other complexes. A comparison is made between the results obtained in an electrostatic analysis carried out calculationally and simulated results corresponding to idealized data from a scanning mutagenesis experiment. It is shown that scanning experiments provide incomplete information on interactions and, if overinterpreted, tend to overestimate the energetic effect of individual side chains that make attractive interactions. Finally, a comparison is made between the results available from a continuum electrostatic model and from a simpler surface-area dependent solvation model. In this case, although the simpler model neglects certain interactions, on average it performs rather well.

DNA-Binding Proteins↗
PubMed1998

High-resolution protein design with backbone freedom.

Recent advances in computational techniques have allowed the design of precise side-chain packing in proteins with predetermined, naturally occurring backbone structures. Because these methods do not model protein main-chain flexibility, they lack the breadth to explore novel backbone conformations. Here the de novo design of a family of alpha-helical bundle proteins with a right-handed superhelical twist is described. In the design, the overall protein fold was specified by hydrophobic-polar residue patterning, whereas the bundle oligomerization state, detailed main-chain conformation, and interior side-chain rotamers were engineered by computational enumerations of packing in alternate backbone structures. Main-chain flexibility was incorporated through an algebraic parameterization of the backbone. The designed peptides form alpha-helical dimers, trimers, and tetramers in accord with the design goals. The crystal structure of the tetramer matches the designed structure in atomic detail.

Amino Acid Sequence↗
PubMed1998

AmbiPack: a systematic algorithm for packing of macromolecular structures with ambiguous distance constraints.

The determination of structures of multimers presents interesting new challenges. The structure(s) of the individual monomers must be found and the transformations to produce the packing interfaces must be described. A substantial difficulty results from ambiguities in assigning intermolecular distance measurements (from nuclear magnetic resonance, for example) to particular intermolecular interfaces in the structure. Here we present a rapid and efficient method to solve the packing and the assignment problems simultaneously given rigid monomer structures and (potentially ambiguous) intermolecular distance measurements. A promising application of this algorithm is to couple it with a monomer searching protocol such that each monomer structure consistent with intramolecular constraints can be subsequently input to the current algorithm to check whether it is consistent with (potentially ambiguous) intermolecular constraints. The algorithm AmbiPack uses a hierarchical division of the search space and the branch-and-bound algorithm to eliminate infeasible regions of the space. Local search methods are then focused on the remaining space. The algorithm generally runs faster as more constraints are included because more regions of the search space can be eliminated. This is not the case for other methods, for which additional constraints increase the complexity of the search space. The algorithm presented is guaranteed to find all solutions to a predetermined resolution. This resolution can be chosen arbitrarily to produce outputs at various level of detail. Illustrative applications are presented for the P22 tailspike protein (a trimer) and portions of beta-amyloid (an ordered aggregate).

Algorithms↗
PubMed1998

Molecular-modeling calculations of enzymatic enantioselectivity taking hydration into account.

A new molecular-modeling methodology has been applied to explain enzymatic enantioselectivity in water. This methodology, which combines vacuum molecular mechanics and the continuum solvation method, should provide a more realistic view of the solvent-enzyme and solvent-substrate interactions than the heretofore used approaches involving the vacuum molecular mechanics only. The methodology described herein has been validated using the experimental data on alpha-chymotrypsin's enantioselectivity in the hydrolysis of four chiral substrates. The reasons why the vacuum molecular mechanics, although not taking hydration into account, still in most cases provide a satisfactory approximation of reality are discussed.

Chymotrypsin↗
PubMed1998

Computation of electrostatic complements to proteins: a case of charge stabilized binding.

Recent evidence suggests that the net effect of electrostatics is generally to destabilize protein binding due to large desolvation penalties. A novel method for computing ligand-charge distributions that optimize the tradeoff between ligand desolvation penalty and favorable interactions with a binding site has been applied to a model for barnase. The result is a ligand-charge distribution with a favorable electrostatic contribution to binding due, in part, to ligand point charges whose direct interaction with the binding site is unfavorable, but which make strong intra-molecular interactions that are uncloaked on binding and thus act to lessen the ligand desolvation penalty.

Bacterial Proteins↗
PubMed1998

Effects of salt bridges on protein structure and design.

Theoretical calculations (Hendsch ZS & Tidor B, 1994, Protein Sci 3:211-226) and experiments (Waldburger CD et al., 1995, Nat Struct Biol 2:122-128; Wimley WC et al., 1996, Proc Natl Acad Sci USA 93:2985-2990) suggest that hydrophobic interactions are more stabilizing than salt bridges in protein folding. The lack of apparent stability benefit for many salt bridges requires an alternative explanation for their occurrence within proteins. To examine the effect of salt bridges on protein structure and stability in more detail, we have developed an energy function for simple cubic lattice polymers based on continuum electrostatic calculations of a representative selection of salt bridges found in known protein crystal structures. There are only three types of residues in the model, with charges of -1, 0, or + 1. We have exhaustively enumerated conformational space and significant regions of sequence space for three-dimensional cubic lattice polymers of length 16. The results demonstrate that, while the more highly charged sequences are less stable, the loss of stability is accompanied by a substantial reduction in the degeneracy of the lowest-energy state. Moreover, the reduction in degeneracy is greater due to charges that pair than for lone charges that remain relatively exposed to solvent. We have also explored and illustrated the use of ion-pairing strategies for rational structural design using model lattice studies.

Models, Molecular↗
PubMed1996

Protein stabilization by removal of unsatisfied polar groups: computational approaches and experimental tests.

The role of polar and charged side chains in partially buried protein environments has been probed in a variant of Arc repressor (MYL) in which hydrophobic interactions between Met31, Tyr36, and Leu40 replace the wild-type salt-bridge interactions between Arg31, Glu36, and Arg40. In the absence of this salt-bridge triad, three additional side chains were identified by continuum electrostatic calculations as incurring larger desolvation penalties during folding than were recovered in favorable electrostatic interactions in the folded state. These side chains (Asn29, Ser44, and Glu48) were mutated singly and collectively to alanine in the MYL background, and the thermodynamic stabilities of the resulting mutant proteins were found to be increased by 0.1 to 1.3 kcal/mol of dimer. All of the mutants displayed cooperative thermal melts and appeared to have well-packed hydrophobic cores by near-UV circular dichroism spectroscopy, indicating that conformational specificity is maintained. The Arc variant (MYL-NA29/SA44/EA48) in which the entire six-residue polar network is replaced by nonpolar groups is 5.1 kcal/mol of dimer more stable than wild type, indicating that the strategy of replacing buried or partially buried charged and polar side chains with hydrophobic residues can lead to substantial stabilization.

Calorimetry↗
PubMed1995

Repacking protein cores with backbone freedom: structure prediction for coiled coils.

Progress in homology modeling and protein design has generated considerable interest in methods for predicting side-chain packing in the hydrophobic cores of proteins. Present techniques are not practically useful, however, because they are unable to model protein main-chain flexibility. Parameterization of backbone motions may represent a general and efficient method to incorporate backbone relaxation into such fixed main-chain models. To test this notion, we introduce a method for treating explicitly the backbone motions of alpha-helical bundles based on an algebraic parameterization proposed by Francis Crick in 1953 [Crick, F. H. C. (1953) Acta Crystallogr. 6, 685-689]. Given only the core amino acid sequence, a simple calculation can rapidly reproduce the crystallographic main-chain and core side-chain structures of three coiled coils (one dimer, one trimer, and one tetramer) to within 0.6-A root-mean-square deviations. The speed of the predictive method [approximately 3 min per rotamer choice on a Silicon Graphics (Mountain View, CA) 4D/35 computer] permits it to be used as a design tool.

Protein Structure, Secondary↗
PubMed1994

The contribution of vibrational entropy to molecular association. The dimerization of insulin.

The association reaction of two molecules to form a single complex must overcome a large entropic barrier due to the loss of translational and rotational degrees of freedom; estimates of the T delta S term are of the order of 30 kcal/mol for proteins. The approach of Chandler and Pratt is used to provide a statistical mechanical formulation for the connection between the gas-phase and solution binding free energies. This makes possible a clear separation of the vibrational contribution to the gas-phase binding enthalpy and entropy from the solvation terms. Further, it suggests that the calculated gas-phase result, should be a good approximation in solution for many systems. To illustrate the formulation, a harmonic dynamics model is used to study the dimerization of insulin. The vibrational entropy increase in the dimer complex, relative to the two separate monomers, is 23 entropy units. This contributes -7.2 kcal/mol to the dimerization free energy. It is not possible to identify a small number of specific dimer modes that give rise to this entropy contribution. Instead, small alterations in the frequencies of many modes below 400 to 600 cm-1 are found to contribute. The relative importance of vibrational and other effects in macromolecule-macromolecule and macromolecule-small molecule associations is discussed.

Binding Sites↗
PubMed1994

Do salt bridges stabilize proteins? A continuum electrostatic analysis.

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).

Crystallography, X-Ray↗
PubMed1994

Helix-capping interaction in lambda Cro protein: a free energy simulation analysis.

The stability mutant Tyr-26-->Asp was studied in the Cro protein from bacteriophage lambda using free energy molecular dynamics simulations. The mutant was calculated to be more stable than the wild type by 3.0 +/- 1.7 kcal/mol/monomer, in reasonable agreement with experiment (1.4 kcal/mol/monomer). Moreover, the aspartic acid in the mutant was found to form a capping interaction with the amino terminus of the third alpha-helix of Cro. The simulations were analyzed to understand better the source of the stability of this helix-capping interaction and to examine the results in light of previous explanations of stabilizing helix caps--namely, a model of local unsatisfied hydrogen bonds at the helix termini and the helix macrodipole model. Analysis of the simulations shows that the stabilizing effect of this charged helical cap is due both to favorable hydrogen bonds with backbone NH groups at the helix terminus and to favorable electrostatic interactions (but not hydrogen bonds) with their carbonyls (effectively the next row of local dipoles in the helix). However, electrostatic interactions are weak or negligible with backbone dipolar groups in the helix further away from the terminus. Moreover, the importance of other local electrostatic interactions with polar side chains near the helix terminus, which are neglected in most treatments of this effect, are shown to be important. Thus, the results support a model that is intermediate between the two previous explanations: both unsatisfied hydrogen bonds at the helix terminus and other, local preoriented dipolar groups stabilize the helix cap. These findings suggest that similar interactions with preoriented dipolar groups may be important for cooperativity in other charge-dipole interactions and may be employed to advantage for molecular design.

Aspartic Acid↗