Transition state docking: a probe for noncovalent catalysis in biological systems. Application to antibody-catalyzed ester hydrolysis

A strategy for pinpointing favorable noncovalent interactions between transition states and active sites of biological catalysts is described. This strategy combines high-level quantum mechanical calculations of transition state geometries with an automated docking procedure using AutoDock. By applying this methodology to antibody-catalyzed hydrolyses of aryl esters (by the 48G7, CNJ206, and 17E8 families of antibodies), varying levels of catalysis are explained in terms of specific hydrogen bonding interactions between combining site residues and transition states. Although these families of antibodies were produced in separate experiments by different researchers using related but different haptens, the mechanism of transition state stabilization appears to be highly conserved. Despite being elicited in response to anionic phosphonate haptens, the best catalysts often utilize hydrogen bond acceptors to stabilize transition states. A mutant of antibody CNJ206, designed based on this observation and predicted to be a better catalyst, is proposed. In the case of antibody 48G7, affinity maturation is shown to produce a catalyst that is highly selective for one of two enantiomeric transition states from a nonselective germline precursor.     

Theoretical investigation of the origins of catalysis of a retro-Diels-Alder reaction by antibody 10F11

The antibody 10F11 catalyzes a retro-Diels-Alder reaction that forms HNO. Deductions about the mechanism of catalysis were made by Reymond, Baumann et al. from X-ray crystal structures and from kinetic measurements for mutated antibodies. We report a study of these reactions with quantum mechanical methods and a study of the substrate and transition state binding to the active site of the antibody 10F11 using density functional theory and empirical docking methods. We have quantitated the likely contributions to catalysis of three residues identified as possible causes of catalysis: Trp H104, Phe H101, and Ser H100. Trp H104 can make a significant contribution to catalysis through dispersive interactions (pi-stacking aromatic-aromatic stabilization). On its own, Phe H101 makes only a small contribution to catalysis. When both aromatic residues are present, they act cooperatively and can make greater contributions to catalysis than expected for each residue alone. Ser H100 and the backbone NH of Phe H101 are expected to act through hydrogen bonding to speed up the reaction, but our calculations suggest that they make only a small contribution to catalysis. Reymond's studies suggest that the hydrogen-bonding network may be mediated through a water molecule in the binding site.     

Modeling the catalysis of anti-cocaine catalytic antibody: competing reaction pathways and free energy barriers

The competing reaction pathways and the corresponding free energy barriers for cocaine hydrolysis catalyzed by an anti-cocaine catalytic antibody, mAb15A10, were studied by using a novel computational strategy based on the binding free energy calculations on the antibody binding with cocaine and transition states. The calculated binding free energies were used to evaluate the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis for each reaction pathway. The free energy barriers for the antibody-catalyzed cocaine hydrolysis were predicted to be the corresponding free energy barriers for the cocaine hydrolysis in water plus the calculated free energy barrier shifts. The calculated free energy barrier shift of -6.87 kcal/mol from the dominant reaction pathway of the cocaine benzoyl ester hydrolysis in water to the dominant reaction pathway of the antibody-catalyzed cocaine hydrolysis is in good agreement with the experimentally derived free energy barrier shift of -5.93 kcal/mol. The calculated mutation-caused shifts of the free energy barrier are also reasonably close to the available experimental activity data. The good agreement suggests that the protocol for calculating the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis may be used in future rational design of possible high-activity mutants of the antibody as anti-cocaine therapeutics. The general strategy of the free energy barrier shift calculation may also be valuable in studying a variety of chemical reactions catalyzed by other antibodies or proteins through noncovalent bonding interactions with the substrates.     

Diels-Alder ribozyme catalysis: a computational approach

A computational comparison of the Diels-Alder reaction of a maleimide and an anthracene in water and the active site of the ribozyme Diels-Alderase is reported. During the course of the catalyzed reaction, the maleimide is held in the hydrophobic pocket while the anthracene approaches to the maleimide through the back passage of the active site. The active site is so narrow that the anthracene has to adopt a tilted approach angle toward maleimide. The conformation of the active site changes marginally at different states of the reaction. Active site dynamics contribution to catalysis has been ruled out. The active site stabilizes the product more than the transition state (TS). The reaction coordinates of the ribozyme reaction in TS, RC1-CD1 and RC4-CD2, are 2.35 and 2.33 A, respectively, compared to 2.37 and 2.36 A in water. The approach angle of anthracene toward maleimide is twisted by 18 degrees in the TS structure of ribozyme reaction while no twisted angle is found in TS of the reaction in water. The free energy barriers for reactions in both ribozyme and water were obtained by umbrella sampling combined with SCCDFTB/MM. The calculated free energy barriers for the ribozyme and water reactions are in good agreement with the experimental values. As expected, Mulliken charges of the atoms involved in the ribozyme reaction change in a similar manner as that of the reaction in water. The proficiency of the Diels-Alder ribozyme reaction originates from the active site holding the two reactants in reactive conformations, in which the reacting atoms are brought together in van der Waals distances and reactants approach to each other at an appropriate angle.     

Catalysis on the coastline: theozyme,molecular dynamics,and free energy perturbation analysis of antibody 21D8 catalysis of the decarboxylation of 5-nitro-3-carboxybenzisoxazole

Antibody 21D8 catalyzes the decarboxylation of 5-nitro-3-carboxybenzisoxazole. The hapten used was designed to induce an antibody binding site with anion binders for the carboxylate, plus a nonpolar environment to accelerate decarboxylation. A recent X-ray crystal structure of 21D8 has shown that the binding pocket contains an array of both polar and charged residues. Nevertheless, 21D8 is able to catalyze a reaction that involves a decrease in polarity from reactant to transition state. The origins of this phenomenon were explored using various computational strategies-quantum mechanics, theozyme models, docking, molecular dynamics, free energy perturbation, and linear interaction energy-the combination of which has produced a consistent picture of catalysis. By partially desolvating the charged carboxylate, 21D8 manages to effect "catalysis on the coastline," without burying the carboxylate in a nonpolar region of the binding pocket. The results have implications for that broad class of enzyme and antibody catalyzed reactions that involve the conversion of a substrate with a relatively localized charge into a transition state with a highly dispersed charge.     

Origin of selectivity in the antibody 20F10-catalyzed Yang cyclization

The first antibody-catalyzed Yang (Norrish type II) cyclization has been achieved with antibodies that were elicited against cis- and trans-2,3-diaryloxetanes. The photocyclization of 1,4-diarylbutan-1-one produced a single stereoisomer of cis-1,2-diarylcyclobutanol with very high enantioselectivity. The background photochemical reaction in the absence of the antibody yielded only fragmentation products. The antibody 20F10-catalyzed reaction was studied in detail, exploring its selectivity, substituent effects, substrate and hapten binding, kinetic parameters and irradiation wavelength dependence. Quantum mechanical calculations suggest that the activation enthalpy of fragmentation pathway is favored by 7.9 kcal/mol over cyclization pathway. Hapten, substrate, and transition state docking studies on a homology based modeled antibody binding site indicate that the trans hapten, substrate and the cyclization transition state have similar binding modes. By contrast, the fragmentation transition state is bound in a different way, not easily accessible within the lifetime of the bound substrate excited state. Several side chain residues were identified that can act as local sensitizers to enhance the cyclization process.     

Just a near attack conformer for catalysis (chorismate to prephenate rearrangements in water,antibody, enzymes,and their mutants)

Standard free energies for formation of ground-state reactive conformers (DeltaGN degrees ) and transition states (DeltaG) in the conversion of chorismate to prephenate in water, B. subtilis mutase, E. coli mutase, and their mutants, as well as a catalytic antibody, are related by DeltaG = DeltaGN degrees + 16 kcal/mol. Thus, the differences in the rate constants for the water reaction and catalysts reactions reside in the mole fraction of substrate present as reactive conformers (NACs). These results, and knowledge of the importance of transition state stabilization in other cases, suggest a proposal that enzymes utilize both NAC and transition state stabilization in the mix required for the most efficient catalysis.     

Experimental determination of the absolute enantioselectivity of an antibody-catalyzed Diels-Alder reaction and theoretical explorations of the origins of stereoselectivity

The exo and endo Diels-Alder adducts of p-methoxycarbonylbenzyl trans-1,3-butadiene-1-carbamate and N,N-dimethylacrylamide have been synthesized, and the absolute configurations of resolved enantiomers have been determined. On the basis of this information, the absolute enantioselectivities of the Diels-Alder reaction catalyzed by antibodies 13G5 and 4D5 as well as other catalytic antibodies elicited in the same immunizations have been established. The effects of different arrangements of catalytic residues on the structure and energetics of the possible Diels-Alder transition states were modeled quantum mechanically at the B3LYP/6-311++G**//B3LYP/6-31+G** level of theory. Flexible docking of these enantiomeric transition states in the antibody active site followed by molecular dynamics on the resulting complexes provided a prediction of the transition-state binding modes and an explanation of the origin of the observed enantioselectivity of antibody 13G5.     

Nonspecific medium effects versus specific group positioning in the antibody and albumin catalysis of the base-promoted ring-opening reactions of benzisoxazoles

The mechanisms by which solvents, antibodies, and albumins influence the rates of base-catalyzed reactions of benzisoxazoles have been explored theoretically. New experimental data on substituent effects and rates of reactions in several solvents, in an antibody, and in an albumin are reported. Quantum mechanical calculations were carried out for the reactions in water and acetonitrile, and docking of the transition state into a homology model of antibody 34E4 and an X-ray structure of human serum albumin was accomplished. A microenvironment made up of catalytic polar groups (glutamate in antibody 34E4 and lysine in human serum albumin) surrounded by relatively nonpolar groups is present in both catalytic proteins.     

Metal-macrocycle framework (MMF): supramolecular nano-channel surfaces with shape sorting capability

Hollow nanostructures for the functional assembly of chemical groups with inner surface geometry and regulable stoichiometry enable steric design of interior reaction centers. Herein we report a metal-macrocycle framework (MMF) that forms single-crystalline nanochannels with five distinct enantiomeric pairs of guest binding pockets. During crystal-soaking experiments, the MMF crystals can encapsulate aromatic molecules with high site selectivity. First, constitutional isomers of dibromobenzene are captured and sorted into different binding pockets. Second, each of the optical isomers of (1R/1S)-1-(3-chlorophenyl)ethanol is included diastereoselectively into one of an enantiomeric pair of binding pockets. An advantage of this strategy is that the interior walls can be "repainted" via replacement of the trapped molecules with alternatives. Such guest uptake behaviors would allow highly regioselective or stereoselective reactions within the nanochannel.     

Molecular dynamics simulation study of the negative correlation in antibody AZ28-catalyzed oxy-cope rearrangement

The oxy-Cope rearrangement reaction in the antibody AZ28 is investigated using ab initio molecular orbital calculations and molecular mechanical molecular dynamics simulations. This antibody, AZ28, is known as one of the few systems where the mature catalytic antibody shows a negative correlation between the transition state analogue (TSA) binding affinity and the catalytic rate of the oxy-Cope rearrangement compared to the germ line catalytic antibody. The ab initio optimized structure shows that the transition state structure has a more planar configuration than the TSA. The favorable electrostatic interactions between AZ28 and the transition state analogue overcome the unfavorable van der Waals interactions; thus, AZ28 shows higher binding affinity for the TSA than the germ line. However, the AZ28 is not flexible enough to accept the relatively planar transition state structure. Because the lower flexibility causes poorer antibody-hapten interaction energies, the activation free energy of the oxy-Cope rearrangement becomes larger in the mature antibody than the germ line. We show that the differences in flexibility between the germ line and the mature form and the differences in structure between TSA and the transition state are the origin of the negative correlation in AZ28-catalyzed oxy-Cope rearrangement. The mutation of residue 34 of the light chain, 34(L), affects the binding free energies through the interresidue interaction because it is the closest to the hapten among the six mutatable residues. However, it does not affect the negative correlation.     

催化2-(对甲氧苯氧基)-四氢吡喃缩醛键水解的抗体研究

以2-[对-(3-羧基)-丙基氧苯]-亚氨基哌啶(2)为半抗原, 联接载体蛋白后免疫动物, 并经杂交瘤技术获得十株单克隆抗体, 其中一株单抗1D1在pH=5.6, 37℃条件下能催化2-(对-甲氧苯氧基)-四氢吡喃(1)缩醛键水解, 加速37倍。     

Bicyclic proline analogues as organocatalysts for stereoselective aldol reactions: an in silico DFT study

Density functional theory has been employed in investigating the efficiency of a series of bicyclic analogues of proline as stereoselective organocatalysts for the aldol reaction. Three classes of conformationally restricted proline analogues, as part of either a [2.2.1] or [2.1.1] bicyclic framework, have been studied. Transition states for the stereoselective C-C bond formation between enamines derived from [2.2.1] and [2.1.1] bicyclic amino acids and p-nitrobenzaldehyde, leading to enantiomeric products, have been identified. Analysis of the transition state geometries revealed that the structural rigidity of catalysts, improved transition state organization as well as other weak interactions influence the relative stabilities of diastereomeric transition states and help contribute to the overall stereoselectivity in the aldol reaction. These bicyclic catalysts are predicted to be substantially more effective in improving the enantiomeric excess than the widely used organocatalyst proline. Enantiomeric excesses in the range 82-95% are predicted for these bicyclic catalysts when a sterically unbiased substrate such as p-nitrobenzaldehyde is employed for the asymmetric aldol reaction. More interestingly, introduction of substituents, as simple as a methyl group, at the ortho position of the aryl aldehyde bring about an increase in the enantiomeric excess to values greater than 98%. The reasons behind the vital energy separation between diastereomeric transition states has been rationalized with the help of a number of weak interactions such as intramolecular hydrogen bonding and Coulombic interactions operating on the transition states. These predictions could have wider implications for the rational design of improved organocatalysts for stereoselective carbon-carbon bond-forming reactions.     

Concerted and stepwise mechanisms in heterogeneous acid catalysis

Conclusion The above developed ideas should not be regarded as criticism of the carbenium ionic conception in heterogeneous acid catalysis. They are aimed rather at its modernization.We consider the present form of this theory as oversimplified. Owing to the absence of solvation effects, surface intermediates of high temperature catalytic conversions of hydrocarbons are of much less ionic character than it is usually believed. Therefore, the carbenium ionic properties are characteristic of transition states but not stable intermediates of these reactions.Another feature of heterogeneous acid catalysis that is not usually taken into account is the bifunctional nature of its active sites. The acidic part protonates adsorbed molecules and the basic one either stabilizes intermediate structures or enables a reversible proton abstraction from transition states. In the case of substrates with a low proton affinity, this results in concerted mechanisms, where simultaneously with a proton transfer a new bond with the surface is formed. Deprotonation of intermediates in the case of catalytic transformations of hydrocarbons with a higher proton affinity supplies them with some features of adsorbed carbenes or ylides. This makes it possible to suggest new mechanisms of such acidic catalytic reactions as synthesis of hydrocarbons from methanol on high silica containing zeolites, olefin oligomerization, etc.     

The origins of noncovalent catalysis of intermolecular Diels-Alder reactions by cyclodextrins,self-assembling capsules,antibodies, and RNAses

The catalysis of Diels-Alder reactions by noncovalent binding by synthetic, protein, and nucleic acid hosts has been surveyed and compared. These catalysts consist of binding cavities that form complexes containing both the diene and the dienophile; the cycloaddition reaction occurs in the cavity. The binding requires no formation of covalent bonds and is driven principally by the hydrophobic (or solvophobic) effect. A molecular mechanics and dynamics study of the cyclodextrin catalysis of a Diels-Alder reaction is used to exemplify and probe this form of catalysis. Detailed kinetic data is available for catalysis by antibodies, RNA, cyclodextrins, and Rebek's tennis ball capsules. Some of these catalysts stabilize the reactants more than the transition state and consequently will only have catalytic effect under conditions of low substrate-to-catalyst ratios. None of the hosts achieve significant specific binding of transition states that is the hallmark of enzyme catalysis.     

Influence of the physical state of phospholipid monolayers on protein binding

Langmuir monolayers were used to characterize the influence of the physical state of phospholipid monolayers on the binding of protein Retinis Pigmentosa 2 (RP2). The binding parameters of RP2 (maximum insertion pressure (MIP), synergy and ΔΠ(0)) in monolayers were thus analyzed in the presence of phospholipids bearing increasing fatty acyl chain lengths at temperatures where their liquid-expanded (LE), liquid-condensed (LC), or solid-condensed (SC) states can be individually observed. The data show that a larger value of synergy is observed in the LC/SC states than in the LE state, independent of the fatty acyl chain length of phospholipids. Moreover, both the MIP and the ΔΠ(0) increase with the fatty acyl chain length when phospholipids are in the LC/SC state, whereas those binding parameters remain almost unchanged when phospholipids are in the LE state. This effect of the phospholipid physical state on the binding of RP2 was further demonstrated by measurements performed in the presence of a phospholipid monolayer showing a phase transition from the LE to the LC state at room temperature. The data collected are showing that very similar values of MIP but very different values of synergy and ΔΠ(0) are obtained in the LE (below the phase transition) and LC (above the phase transition) states. In addition, the binding parameters of RP2 in the LE (below the phase transition) as well as in the LC (above the phase transition) states were found to be indistinguishable from those where single LC and LE states are respectively observed. The preference of RP2 for binding phospholipids in the LC state was then confirmed by the observation of a large modification of the shape of the LC domains in the phase transition. Therefore, protein binding parameters can be strongly influenced by the physical state of phospholipid monolayers. Moreover, measurements performed with the α/β domain of RP2 strongly suggest that the β helix of RP2 plays a major role in the preferential binding of this protein to phospholipids in the LC state.     

Positional ordering of reacting groups contributes significantly to the efficiency of proton transfer at an antibody active site

Catalytic antibody 34E4 accelerates the conversion of benzisoxazoles to salicylonitriles with surprising efficiency, exploiting a carboxylate base with an elevated pKa for proton abstraction. Mutagenesis of this antibody, produced as a chimeric Fab, confirms the prediction of a homology model that GluH50 is the essential catalytic residue. Replacement of this residue by glutamine, alanine, or glycine reduces catalytic activity by more than 2.6 x 10(4)-fold. By comparing the chemical proficiencies of the parent antibody with the chemical proficiencies of acetate and the mutants, the effective concentration of the catalytic side chain was estimated to be >51 000 M. The 2.1 kcal/mol destabilization of the transition state observed when GluH50 is replaced by aspartate suggests that positional ordering imposed by the antibody active site contributes significantly to the efficiency of proton transfer. The observation that the GluH50Ala and GluH50Gly variants could not be chemically rescued by exogenous addition of high concentrations of formate or acetate further underscores the advantage the antibody derives from covalently fixing its base at the active site. Although medium effects also play an important role in 34E4, for example in enhancing the reactivity of the carboxylate side chain through desolvation, comparison of 34E4 with less proficient antibodies shows that positioning a carboxylate in a hydrophobic binding pocket alone is insufficient for efficient general base catalysis. Our results demonstrate that structural complementarity between the antibody and its substrate in the transition state is an important and necessary component of 34E4's high activity. By harnessing an additional catalytic group that could serve as a general acid to stabilize developing negative charge in the leaving group, overall efficiencies rivaling those of highly evolved enzymes should be accessible.     

A theoretical study of Favorskii reaction stereochemistry. Lessons in torquoselectivity

The mechanisms of the chloroenolate-->cyclopropanone step of the "normal" Favorskii rearrangement have been investigated in detail using high-level ab initio calculations. A series of simple alpha-chloroenolates, based on chloroacetone (6), all monomethyl derivatives (7-9), a dimethyl analogue (10), and 1-acetyl-1-chlorocyclohexane (11) was first used to explore and define the basic features of the mechanism, which include the finding of both an "inversion" and a "retention" transition state and that in most cases these arise from separate ground-state conformations of the chloroenolate. These theoretical studies were then extended to an isomeric pair of chloroenolates 1 and 2, the cis- and trans-2-methyl derivatives of 11, which are the reactive intermediates involved in a well-known experimental study carried out by Stork and Borowitz (S-B). Finally, three alpha-chlorocyclohexanone enolate systems 12-14 were studied, since these intermediates have a more restricted enolate geometry. The "inversion" mechanism has been described as an SN2 process but the present results, while supporting a concerted process, is better described as an oxyallyl structure undergoing concerted ring closure. The "retention" mechanism has been described as SN1-like, but the calculations show that this process is also concerted, although much less so, and again involves oxyallyl-like transition-states. The model systems 6-8, 10, and 11 with a potential plane of symmetry have two enantiomeric transition states for inversion and another two for retention of configuration (at the C-Cl center). With 9 and the S-B models 1 and 2, with no symmetry plane, there are a calculated total of four diastereomeric transition states for cyclopropanone ring closure in each case, two for inversion and two for retention. While the transition-state energies calculated for simple chloroenolates favor the inversion process, the S-B models 1 and 2 have almost equal inversion-retention transition-state energies. Solvation simulation calculations of ground states and transition states suggest that the retention mechanism becomes relatively more favored in polar solvents, in agreement with some experimental results. In the chloroenolates 12-14, both inversion and retention mechanisms were also located, these arising from two different ground-state ring conformations of the enolate. In these models, one also finds similar inversion and retention transition-state energies, but again with a small preference for the inversion process.     

Mechanism and functional role of antibody catalysis

The light (L) chain of a model antibody (Ab) was deduced to contain a serine protease-like catalytic site capable of cleaving peptide bonds. The catalytic site is encoded by a germline V_L gene. The catalytic activity can potentially be improved by somatic sequence diversification and pairing of the L chain with the appropriate heavy chain. Autoimmune disease is associated with increased synthesis of antigen (Ag)-specific Abs, but the reasons for this phenomenon are not known. Only recently has attention turned to the functional role of the catalytic function. Preliminary studies confirm that the catalytic cleavage of peptide bonds is a more potent means to achieve Ag neutralization, compared to reversible Ag binding. Administration of a monoclonal Ab to VIP in experimental animals induces an inflammatory response in the airways, suggesting that catalytic autoantibodies to this peptide found in airway disease and lupus are capable of causing airway dysfunction. The phenomenon of autoantibody catalysis can potentially be applied to isolate efficient catalysts directed against tumor or microbial Ags by exposing the autoimmune repertoire to such Ags or their analogs capable of recruiting the germline V_L gene encoding the catalytic site.     

Complex formation of a calcium-dependent antibody: a thermodynamical consideration

The elution of FLAG-fusions (an octapeptide with the sequence DYKDDDDK) from immobilized anti-FLAG antibody M1 cannot be explained by a switch of the equilibrium binding constant to a lower value. To get a further insight into thermodynamics, the binding of anti-FLAG antibody M1 to the FLAG peptide was studied by real-time biosensor technology at seven different temperatures in the range from 5 to 35 degrees C. Binding studies were performed in the presence and absence of calcium. Thermodynamic parameters such as change in Gibbs free energy (deltaG), enthalpy (deltaH) and entropy (deltaS) were evaluated from the corresponding equilibrium data applying the integrated Van't Hoff equation. In contrast to similar kinetic data obtained, the contribution of deltaH and deltaS to deltaG in the presence or absence of calcium results in a different conformation of the antibody-antigen complex under binding and non-binding conditions. Therefore, complex dissociation with EDTA must be effected during a transition state of complex formation and dissociation.