PubMed Health⌕ Search

Biomedical subjects

B Tidor

Publications and source records attributed to B Tidor.

At least 19 recordsLinked to original sources
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↗
PubMed1993

Probing the roles of residues at the e and g positions of the GCN4 leucine zipper by combinatorial mutagenesis.

Combinatorial mutagenesis with an alphabet limited to alanine, glutamic acid, lysine, and threonine was used to probe the role of interactions involving surface residues in stabilizing a short alpha-helical coiled coil. The residues at eight e and g positions in the leucine zipper of the Saccharomyces cerevisiae transcription factor GCN4 were randomized to these four residues in a lambda repressor-leucine zipper fusion protein, resulting in 65,536 possible residue combinations. Roughly three-fourths of these combinations allowed stable coiled-coil formation as assayed by DNA binding by the fusion protein. To understand the basis for the activity differences, functional and non-functional mutants were sequenced and statistical tests were applied to identify structure/function correlations. Helix-forming propensity and favorable intrasubunit and intersubunit charge-charge interactions were positively correlated with activity. These studies suggest that the identities of surface side chains at the e and g positions of coiled coils contribute modestly to stability; by comparison with previous work, however, the e and g positions are far less critical than residues at the a and d positions, which form the hydrophobic core of the dimer interface.

Amino Acid Sequence↗
PubMed1993

The contribution of cross-links to protein stability: a normal mode analysis of the configurational entropy of the native state.

The vibrational entropy of native BPTI, with three disulfide bonds, was determined by use of normal mode calculations and compared with that of folded variants having either one less disulfide bond or lacking a peptide bond at the trypsin-reactive site. Favorable contributions to the free energy of 2.5-5.1 kcal/mol at 300 K were calculated for the reduction of disulfide bonds in the folded state, whereas no favorable contribution was found for the hydrolysis of the peptide bond cleaved by trypsin. This is on the order of the effect of disulfides in the unfolded state. The implications of these results for the stabilization of a folded protein by the introduction of cross-links are discussed.

Aprotinin↗
PubMed1991

Contribution of the hydrophobic effect to protein stability: analysis based on simulations of the Ile-96----Ala mutation in barnase.

Molecular dynamics simulations have been used to compute the difference in the unfolding free energy between wild-type barnase and the mutant in which Ile-96 is replaced by alanine. The simulations yield results (-3.42 and -5.21 kcal/mol) that compare favorably with experimental values (-3.3 and -4.0 kcal/mol). The major contributions to the free energy difference arise from bonding terms involving degrees of freedom of the mutated side chain and from nonbonded interactions of that side chain with its environment in the folded protein. By comparison with simulations of an extended peptide in the absence of solvent, used as a reference state, hydration effects are shown to play a minor role in the overall free energy balance for the Ile----Ala transformation. The implications of these results for our understanding of the hydrophobic effect and its contribution to protein stability are discussed.

Alanine↗
PubMed1991

Arginine-mediated RNA recognition: the arginine fork.

Short peptides that contain the basic region of the HIV-1 Tat protein bind specifically to a bulged region in TAR RNA. A peptide that contained nine arginines (R9) also bound specifically to TAR, and a mutant Tat protein that contained R9 was fully active for transactivation. In contrast, a peptide that contained nine lysines (K9) bound TAR poorly and the corresponding protein gave only marginal activity. By starting with the K9 mutant and replacing lysine residues with arginines, a single arginine was identified that is required for specific binding and transactivation. Ethylation interference experiments suggest that this arginine contacts two adjacent phosphates at the RNA bulge. Model building suggests that the arginine eta nitrogens and the epsilon nitrogen can form specific networks of hydrogen bonds with adjacent pairs of phosphates and that these arrangements are likely to occur near RNA loops and bulges and not within double-stranded A-form RNA. Thus, arginine side chains may be commonly used to recognize specific RNA structures.

Amino Acid Sequence↗
PubMed1991

Simulation analysis of the stability mutant R96H of T4 lysozyme.

Free energy simulation methods are used to analyze the effects of the mutation Arg 96----His on the stability of T4 lysozyme. The calculated stability change and the lack of significant structural rearrangement in the folded state due to the mutation are in agreement with experimental studies [Kitamura, S., & Sturtevant, J. M. (1989) Biochemistry 28, 3788-3792; Weaver, L. H., et al. (1989) Biochemistry 28, 3793-3797]. By use of thermodynamic integration, the contributions of specific interactions to the free energy change are evaluated. It is shown that a number of contributions that stabilize the wild type or the mutant partially cancel in the overall free energy difference; some of these involve the unfolded state. Comparison of the results with conclusions based on structural and thermodynamic data leads to new insights into the origin of the stability difference between wild-type and mutant proteins. Of particular interest is the importance of the contributions of more distant residues, solvent water, and the covalent linkage of the mutated amino acid. Also, the analysis of the interactions of Arg/His 96 with the C-terminal end of a helix (residues 82-90) makes it clear that the nearby carbonyl groups (Tyr 88 and Asp 89) make the dominant contribution, that the amide groups do not contribute significantly, and that the helix-dipole model is inappropriate for this case.

Arginine↗
PubMed1991

Simulation analysis of the stability mutants R96H of bacteriophage T4 lysozyme and I96A of barnase.

Free energy simulation methods are used to analyse the effects of the mutation Arg-96----His on the stability of bacteriophage T4 lysozyme and of Ile-96----Ala on the stability of barnase. By use of thermodynamic integration, the contributions of specific interactions to the free energy change are evaluated. It is shown that a number of contributions that stabilize the wild-type or the mutant partially cancel in the overall free energy difference; some of these involve the unfolded state. Comparison of the results with conclusions based on structural and thermodynamic data leads to new insights into the origin of the stability difference between wild-type and mutant proteins. For the charged-to-charged amino acid mutation in T4 lysozyme, the importance of the contributions of more distant residues, solvent water and the covalent linkage involving the mutated amino acid are of particular interest. Also, the analysis of the Arg-96 to His mutation with respect to the interactions with the C-terminal end of a helix (residues 82-90) indicates that the nearby carbonyl groups (Tyr-88 and Asp-89) make the dominant contribution, that the amide groups do not contribute significantly and that the helix dipole model is inappropriate for this case. For the non-polar-to-non-polar amino acid mutation in barnase, the solvent contribution is unimportant, and covalent terms are shown to be significant because they do not cancel between the folded and unfolded state.

Bacterial Proteins↗
PubMed1990

Free energy of sickling: A simulation analysis.

Molecular dynamics simulations were performed to calculate the difference between the dimerization free energies of normal human deoxyhemoglobin (HbA) and the mutant sickle-cell deoxyhemoglobin HbS (Glu-beta 6----Val) for one of the lateral contacts in the HbS x-ray structure. The simulations yield a value of--15 kcal/mol. Although there is no quantitative experimental value for comparison, this is in qualitative agreement with the experimental result that HbS self-assembles into multistranded fibers that are responsible for erythrocyte sickling, while HbA does not. The free-energy difference was decomposed into enthalpic and entropic terms, both of which are significant, and the contributions of individual protein residues and of the solvent were examined. Electrostatic effects play the dominant role in favoring dimerization of HbS compared with HbA; van der Waals interactions make a negligible contribution to the difference. Both differential solvation and protein-protein interactions are important. Interactions within the donor tetramer (i.e., that containing the Glu-beta 6 mutation site), as well as those with the acceptor tetramer, contribute to the preferential free energy of dimerization of HbS.

Calorimetry↗