Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Stabilization of microtubule dynamics by estramustine by binding to a novel site in tubulin: A possible mechanistic basis for its antitumor action

The cellular targets for estramustine, an antitumor drug used in the treatment of hormone-refractory prostate cancer, are believed to be the spindle microtubules responsible for chromosome separation at mitosis. Estramustine only weakly inhibits polymerization of purified tubulin into microtubules by binding to tubulin (Kd, ≈30 μM) at a site distinct from the colchicine or the vinblastine binding sites. However, by video microscopy, we find that estramustine strongly stabilizes growing and shortening dynamics at plus ends of bovine brain microtubules devoid of microtubule-associated proteins at concentrations substantially below those required to inhibit polymerization of the microtubules. Estramustine strongly reduced the rate and extent both of shortening and growing, increased the percentage of time the microtubules spent in an attenuated state, neither growing nor shortening detectably, and reduced the overall dynamicity of the microtubules. Significantly, the combined suppressive effects of vinblastine and estramustine on the rate and extent of shortening and dynamicity were additive. Thus, like the antimitotic mechanisms of action of the antitumor drugs vinblastine and taxol, the antimitotic mechanism of action of estramustine may be due to kinetic stabilization of spindle microtubule dynamics. The results may explain the mechanistic basis for the benefit derived from combined use of estramustine with vinblastine or taxol, two other drugs that target microtubules, in the treatment of hormone-refractory prostate cancer.

Theoretical investigation of the [1,2]-sigmatropic hydrogen migration in the mechanism of oxidation of 2-aminobenzoyl-CoA by 2-aminobenzoyl-CoA monooxygenase/reductase

The flavin hydroperoxide at the active site of the mixed-function oxidase 2-aminobenzoyl-CoA monooxygenase/reductase (Azoarcus evansii) transfers an oxygen to the 5-position of the 2-aminobenzoyl-CoA substrate to provide the alkoxide intermediate II−. Hydrogen migration from C5 to C6 follows this monooxygenation. The nature of the monooxygenation intermediate and plausible competing reactions leading to hydrogen migration have been considered. Ab initio molecular orbital theory has been used to calculate structures and electron distributions in intermediate and transition state structures. Electrostatic potential surface calculations establish that the transition state and product, associated with the C5 to C6 hydrogen transfer, are stabilized by electron distribution to the benzoyl-CoA thioester carbonyl oxygen. This is not so for the transition state and product associated with hydrogen transfer from C5 to C4. The activation energy for the 5,6-shift is 2.5 kcal/mol lower than that for the 5,4-shift. In addition, the product of the hydrogen 5,6-shift is more stable than is the product of the hydrogen 5,4-shift, by ≈6 kcal/mol. These results explain why only the shift of hydrogen from C5 to C6 is observed experimentally. Oxygen transfer and hydrogen migration almost coincide in the gas phase (activation energy of ≈0.6 kcal/mol, equivalent to a single bond vibration). Enzymatic formation of alkoxide II− requires its stabilization; thus, the rate constant for its breakdown would be slower than in the gas phase.

Inhibiting transthyretin amyloid fibril formation via protein stabilization

Transthyretin (TTR) amyloid fibril formation is observed systemically in familial amyloid polyneuropathy and senile systemic amyloidosis and appears to be the causative agent in these diseases. Herein, we demonstrate conclusively that thyroxine (10.8 μM) inhibits TTR fibril formation efficiently in vitro and does so by stabilizing the tetramer against dissociation and the subsequent conformational changes required for amyloid fibril formation. In addition, the nonnative ligand 2,4,6-triiodophenol, which binds to TTR with slightly increased affinity also inhibits TTR fibril formation by this mechanism. Sedimentation velocity experiments were employed to show that TTR undergoes dissociation (linked to a conformational change) to form the monomeric amyloidogenic intermediate, which self-assembles into amyloid in the absence, but not in the presence of thyroxine. These results demonstrate the feasibility of using small molecules to stabilize the native fold of a potentially amyloidogenic human protein, thus preventing the conformational changes, which appear to be the common link in several human amyloid diseases. This strategy and the compounds resulting from further development should prove useful for critically evaluating the amyloid hypothesis—i.e., the putative cause-and-effect relationship between TTR amyloid deposition and the onset of familial amyloid polyneuropathy and senile systemic amyloidosis.

Physical reasons for the unusual α-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH2 in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

Stabilization of tetanus and diphtheria toxoids against moisture-induced aggregation.

The progress toward single-dose vaccines has been limited by the poor solid-state stability of vaccine antigens within controlled-release polymers, such as poly(lactide-co-glycolide). For example, herein we report that lyophilized tetanus toxoid aggregates during incubation at 37 degrees C and elevated humidity--i.e., conditions relevant to its release from such systems. The mechanism and extent of this aggregation are dependent on the moisture level in the solid protein, with maximum aggregation observed at intermediate moisture contents. The main aggregation pathway is consistent with formaldehyde-mediated cross-linking, where reactive electrophiles created and stored in the vaccine upon formalinization (exposure to formaldehyde during vaccine preparation) react with nucleophiles of a second vaccine molecule to form intermolecular cross-links. This process is inhibited by the following: (i) succinylating the vaccine to block reactive amino groups; (ii) treating the vaccine with sodium cyanoborohydride, which presumably reduces Schiff bases and some other electrophiles created upon formalinization; and (iii) addition of low-molecular-weight excipients, particularly sorbitol. The moisture-induced aggregation of another formalinized vaccine, diphtheria toxoid, is also retarded by succinylation, suggesting the generality of this mechanism for formalinized vaccines. Hence, mechanistic stability studies of the type described herein may be important for the development of effective single-dose vaccines.