Structural insights into steroid hormone binding: the crystal structure of a recombinant anti-testosterone Fab fragment in free and testosterone-bound forms.

The monoclonal anti-testosterone antibody (3-C(4)F(5)) has a relatively high affinity (3 x 10(8) m(-1)) with an overall good specificity profile. However, the earlier characterized binding properties have shown that both the affinity and specificity of this antibody must be improved if it is intended for use in clinical immunoassays. In this paper, the crystal structures of the recombinant anti-testosterone (3-C(4)F(5)) Fab fragment have been determined in the testosterone-bound and free form at resolutions of 2.60 and 2.72 A, respectively. The high affinity binding of the (3-C(4)F(5)) Fab is mainly determined by shape complementarity between the protein and testosterone. Only one direct hydrogen bond is formed between the hydroxyl group of the testosterone D-ring and the main-chain oxygen of Gly100(J)H. The testosterone is deeply bound in a hydrophobic pocket, and the close shape complementarity is mainly formed by the third complementarity-determining regions (CDR) of the heavy and light chain. Comparison of the bound structure with the free structure indicates conformational changes in the protein upon testosterone binding. The conformational changes of the side chains of two residues Glu95H and Tyr99H in the CDR-H3 are particularly essential for the binding. Interesting similarities in the binding of different steroids were also observed upon comparison of the available structures of anti-steroid antibodies.     

Mutagenesis of a stacking contact in the MS2 coat protein-RNA complex.

The thermodynamic contribution of a stacking interaction between Tyr85 in MS2 coat protein and a single-stranded pyrimidine in its RNA binding site has been examined. Mutation of Tyr85 to Phe, His, Cys, Ser and Ala decreased the RNA affinity by 1-3 kcal/mol under standard binding conditions. Since the Phe, His and Cys 85 proteins formed UV photocrosslinks with iodouracil-containing RNA at the same rate as the wild-type protein, the mutant proteins interact with RNA in a similar manner. The pH dependence of KD for the Phe and His proteins differs substantially from the wild-type protein, suggesting that the titration of position 85 contributes substantially to the binding properties. Experiments with specifically substituted phosphorothioate RNAs confirm a hydrogen bond between the hydroxyl group of tyrosine and a phosphate predicted by the crystal structure.     

Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone,androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the binding of testosterone (TES), 5alpha-dihydrotestosterone (5ADHT), androstenedione (AND), and dehydroepiandrosterone sulfate (DHEAS) to the monoclonal antitestosterone antibody 3-C(4)F(5). The relative binding free energy of TES and AND was also calculated with free energy perturbation (FEP) simulations. The antibody 3-C(4)F(5) has a relatively high affinity (3 x 10(8) M(-1)) and on overall good binding profile for testosterone but its cross-reactivity with DHEAS has been the main reason for the failure to use this antibody in clinical immunoassays. The relative binding free energies obtained with the MM-PBSA method were 1.5 kcal/mol for 5ADHT, 3.8 kcal/mol for AND, and 4.3 kcal/mol for DHEAS, as compared to TES. When a water molecule of the ligand binding site, observed in the antibody-TES crystal structure, was explicitly included in MM-PBSA calculations, the relative binding energies were 3.4, 4.9, and 5.4 kcal/mol for 5ADHT, AND, and DHEAS, respectively. The calculated numbers are in correct order but larger than the corresponding experimental energies of 1.3, 1.5, and 2.6 kcal/mol, respectively. The fact that the MM-PBSA method reproduced the relative binding free energies of DHEAS, a steroid having a negatively charged sulfate group, and the neutrally charged TES, 5ADHT, and AND in satisfactory agreement with experiment shows the robustness of the method in predicting relative binding affinities. The 800-ps FEP simulations predicted that the antibody 3-C(4)F(5) binds TES 1.3 kcal/mol tighter than AND. Computational mutagenesis of selected amino acid residues of the ligand binding site revealed that the lower affinities of AND and DHEAS as compared to TES are due to a combined effect of several residues, each contributing a small fraction to the tighter binding of TES. An exception to this is Tyr99H, whose mutation to Ala lowered the binding of DHEAS 0.7 kcal/mol more than the binding of TES. This is probably due to the hydrogen bonding interaction formed between the OH group of Tyr99H and the sulfate group of DHEAS. Computational mutagensis data also showed that the affinity of the steroids to the antitestosterone antibody 3-C(4)F(5) would be enhanced if Trp47H were repositioned so that it would make more extensive contacts with the bound ligands. In addition, the binding of steroids to antitestosterone, antiprogesterone, and antiestradiol antibodies is discussed.     

Contribution of the hydrogen-bond network involving a tyrosine triad in the active site to the structure and function of a highly proficient ketosteroid isomerase from Pseudomonas putida biotype B

Delta(5)-3-Ketosteroid isomerase from Pseudomonas putida biotype B is one of the most proficient enzymes catalyzing an allylic isomerization reaction at rates comparable to the diffusion limit. The hydrogen-bond network (Asp99... Wat504...Tyr14...Tyr55...Tyr30) which links the two catalytic residues, Tyr14 and Asp99, to Tyr30, Tyr55, and a water molecule in the highly apolar active site has been characterized in an effort to identify its roles in function and stability. The DeltaG(U)(H2O) determined from equilibrium unfolding experiments reveals that the elimination of the hydroxyl group of Tyr14 or Tyr55 or the replacement of Asp99 with leucine results in a loss of conformational stability of 3.5-4.4 kcal/mol, suggesting that the hydrogen bonds of Tyr14, Tyr55, and Asp99 contribute significantly to stability. While decreasing the stability by about 6.5-7.9 kcal/mol, the Y55F/D99L or Y30F/D99L double mutation also reduced activity significantly, exhibiting a synergistic effect on k(cat) relative to the respective single mutations. These results indicate that the hydrogen-bond network is important for both stability and function. Additionally, they suggest that Tyr14 cannot function efficiently alone without additional support from the hydrogen bonds of Tyr55 and Asp99. The crystal structure of Y55F as determined at 1.9 A resolution shows that Tyr14 OH undergoes an alteration in orientation to form a new hydrogen bond with Tyr30. This observation supports the role of Tyr55 OH in positioning Tyr14 properly to optimize the hydrogen bond between Tyr14 and C3-O of the steroid substrate. No significant structural changes were observed in the crystal structures of Y30F and Y30F/Y55F, which allowed us to estimate approximately the interaction energies mediated by the hydrogen bonds Tyr30...Tyr55 and Tyr14...Tyr55. Taken together, our results demonstrate that the hydrogen-bond network provides the structural support that is needed for the enzyme to maintain the active-site geometry optimized for both function and stability.     

Bacillus circulans MH-K1 chitosanase: amino acid residues responsible for substrate binding

To identify the amino acids responsible for the substrate binding of chitosanase from Bacillus circulans MH-K1 (MH-K1 chitosanase), Tyr148 and Lys218 of the chitosanase were mutated to serine and proline, respectively, and the mutated chitosanases were characterized. The enzymatic activities of Y148S and K218P were found to be 12.5% and 0.16% of the wild type, respectively. When the (GlcN)3 binding ability to the chitosanase was evaluated by fluorescence spectroscopy and thermal unfolding experiments, the binding abilities of both mutant enzymes were markedly reduced as compared with the wild type enzyme. The affinity of the enzyme for the trisaccharide decreased by 1.0 kcal/mol of binding free energy for Y148S, and 3.7 kcal/mol for K218P. The crystal structure of K218P revealed that Pro218 forms a cis-peptide bond and that the state of the flexible loop containing the 218th residue is considerably affected by the mutation. Thus, we conclude that the flexible loop containing Lys218 plays an important role in substrate binding, and that the role of Tyr148 is less critical, but still important, due to a stacking interaction or hydrogen bond.     

The roles of Tyr391 and Tyr619 in RB69 DNA polymerase replication fidelity

In the family-B DNA polymerase of bacteriophage RB69, the conserved aromatic palm-subdomain residues Tyr391 and Tyr619 interact with the last primer-template base-pair. Tyr619 interacts via a water-mediated hydrogen bond with the phosphate of the terminal primer nucleotide. The main-chain amide of Tyr391 interacts with the corresponding template nucleotide. A hydrogen bond has been postulated between Tyr391 and the hydroxyl group of Tyr567, a residue that plays a key role in base discrimination. This hydrogen bond may be crucial for forcing an infrequent Tyr567 rotamer conformation and, when the bond is removed, may influence fidelity. We investigated the roles of these residues in replication fidelity in vivo employing phage T4 rII reversion assays and an rI forward assay. Tyr391 was replaced by Phe, Met and Ala, and Tyr619 by Phe. The Y391A mutant, reported previously to decrease polymerase affinity for incoming nucleotides, was unable to support DNA replication in vivo, so we used an in vitro fidelity assay. Tyr391F/M replacements affect fidelity only slightly, implying that the bond with Tyr567 is not essential for fidelity. The Y391A enzyme has no mutator phenotype in vitro. The Y619F mutant displays a complex profile of impacts on fidelity but has almost the same mutational spectrum as the parental enzyme. The Y619F mutant displays reduced DNA binding, processivity, and exonuclease activity on single-stranded DNA and double-stranded DNA substrates. The Y619F substitution would disrupt the hydrogen bond network at the primer terminus and may affect the alignment of the 3' primer terminus at the polymerase active site, slowing chemistry and overall DNA synthesis.     

Contribution of a tyrosine side chain to ribonuclease A catalysis and stability.

An intricate architecture of covalent bonds and noncovalent interactions appear to position the side chain of Lys 41 properly within the active site of bovine pancreatic ribonuclease A (RNase A). One of these interactions arises from Tyr 97, which is conserved in all 41 RNase A homologues of known sequence. Tyr 97 has a solvent-inaccessible side chain that donates a hydrogen bond to the main-chain oxygen of Lys 41. Here, the role of Tyr 97 was examined by replacing Tyr 97 with a phenylalanine, alanine, or glycine residue. All three mutant proteins have diminished catalytic activity, with the value of Kcat being perturbed more significantly than that of Km. The free energies with which Y97F, Y97A, and Y97G RNase A bind to the rate-limiting transition state during the cleavage of poly(cytidylic acid) are diminished by 0.74, 3.3, and 3.8 kcal/mol, respectively. These results show that even though Tyr 97 is remote from the active site, its side chain contributes to catalysis. The role of Tyr 97 in the thermal stability of RNase A is large. The conformational free energies of native Y97F, Y97A, and Y97G RNase A are decreased by 3.54, 12.0, and 11.7 kcal/mol, respectively. The unusually large decrease in stability caused by the Tyr-->Phe mutation could result from a decrease in the barrier to isomerization of the Lys 41-Pro 42 peptide bond.     

Oxidation and electronic state dependence of proton transfer in the enzymatic cycle of cytochrome P450eryF

Hydrogen bond networks, consisting of hydrogen bonded waters anchored by polar/acidic amino acid sidechains, are often present in the vicinity of the oxygen binding clefts of P450s. Density functional and quantum dynamics calculations of a O(2) binding cleft network model of cytochrome P450eryF(CYP107A1) indicate that such structural motifs facilitate ultrafast proton transfer from network waters to the dioxygen of the reduced oxyferrous species via a multiple proton translocation mechanism with barriers of 7-10 kcal/mol on its doublet ground state, and that the energies of the proton transfer reactant and constrained proton transfer products have an electronic and oxidation state dependence [J. Am. Chem. Soc. 124 (2002) 1430]. In the present study, the origin of the oxidation state dependence is shown to have its roots in differential proton affinities while the electronic state dependence of the reduced oxyferrous heme has its origins in subtle differences in network topologies near the transition state of the initial proton transfer event. Relaxed potential surface scans and unconstrained proton transfer product optimizations indicate that the proton transfer product in both the singlet oxyferrous heme and the reduced oxyferrous heme species in a quartet state are not viable stable (bound) states relative to the reactant form. While the proton affinity of H(3)O(+) is sufficient for it to protonate both the oxyferrous and the reduced oxyferrous heme species, hydrogen bond network stabilized water is only capable of protonating the reduced oxyferrous form. This interpretation is substantiated by study of the NO bound reduced ferrous heme of P450nor, which is isoelectronic with the oxyferrous heme and has a similar proton affinity. Density functional calculations on a more extensive O(2) binding cleft model support the multiple proton translocation mechanism of transfer but indicates that the significant negative charge density on the bound dioxygen of the reduced oxyferrous heme species, in its doublet ground state, polarizes the associated hydrogen bond network sufficiently so as to result in short, strong, low-barrier hydrogen bonds. The computed O-H-O bond distances are less than 2.55 A and have a near degeneracy of the proton transfer reactant and initial (sudden) proton transfer products. These low-barrier hydrogen bond features, in addition to the finding of a (zero point uncorrected) barrier of 1.3 kcal/mol, indicate that proton transfer from water to the distal oxygen should be rapid, facile and may not require large curvature tunneling as originally suggested by use of a smaller model. An initial assessment of protonation of the reduced oxyferrous heme distal oxygen by a model of 6-deoxyerythronolide B (6-DEB) indicates it to be low barrier (3.8 kcal/mol) and exothermic (-2.9 kcal/mol). The combined results indicate the plausibility of simultaneous diprotonation of the distal oxygen of the reduced oxyferrous heme, leading to O-O bond scission, using the combined water network and 6-DEB substrate protonation agents.     

Hsp90 chaperones have an energetic hot‐spot for binding inhibitors

Although Hsp90‐family chaperones have been extensively targeted with ATP‐competitive inhibitors, it is unknown whether high affinity is achieved from a few highly stabilizing contacts or from many weaker contacts within the ATP‐binding pocket. A large‐scale analysis of Hsp90α:inhibitor structures shows that inhibitor hydrogen‐bonding to a conserved aspartate (D93 in Hsp90α) stands out as most universal among Hsp90 inhibitors. Here we show that the D93 region makes a dominant energetic contribution to inhibitor binding for both cytosolic and organelle‐specific Hsp90 paralogs. For inhibitors in the resorcinol family, the D93:inhibitor hydrogen‐bond is pH‐dependent because the associated inhibitor hydroxyl group is titratable, rationalizing a linked‐protonation event previously observed by the Matulis group. The inhibitor hydroxyl group pK_a associated with the D93 hydrogen‐bond is therefore critical for optimizing the affinity of resorcinol derivatives, and we demonstrate that spectrophotometric measurements can determine this pK_a value. Quantifying the energetic contribution of the D93 hotspot is best achieved with the mitochondrial Hsp90 paralog, yielding 3–6 kcal/mol of stabilization (35–60% of the total binding energy) for a diverse set of inhibitors. The Hsp90 Asp93?Asn substitution has long been known to abolish nucleotide binding, yet puzzlingly, native sequences of structurally similar ATPases, such as Topoisomerasese II, have an asparagine at this same crucial site. While aspartate and asparagine sidechains can both act as hydrogen bond acceptors, we show that a steric clash prevents the Hsp90 Asp93?Asn sidechain from adopting the necessary rotamer, whereas this steric restriction is absent in Topoisomerasese II.     

Identification of functionally important residues in the pyridoxal-5'-phosphate-dependent catalytic antibody 15A9

Antibody 15A9 is unique in its ability to catalyze the transamination reaction of hydrophobic D-amino acids with pyridoxal-5'-phosphate (PLP). Both previous chemical modification studies and a three dimensional (3-D) homology model indicated the presence of functionally important tyrosine residues in the antigen-binding cavity of antibody 15A9. To gain further insight into the hapten, ligand binding, and catalytic mechanism of 15A9, all tyrosine residues in the complementarity-determining regions (CDRs) and the single arginine residue in CDR3 of the light chain were subject to an alanine scan. Substitution of Tyr(H33), Tyr(L94), or Arg(L91) abolished the catalytic activity and reduced the affinity for PLP and N(a)-(5'-phosphopyridoxyl)-amino acids, which are close analogues of covalent PLP-substrate adducts. The Tyr(H100b)Ala mutant possessed no detectable catalytic activity, while its affinity for each ligand was essentially the same as that of the wild-type antibody. The binding affinity for the hapten was drastically reduced by a Tyr(L32)Ala mutation, suggesting that the hydroxyphenyl group of Tyr(L32) participates in the binding of the extended side chain of the hapten. The other Tyr --> Ala substitutions affected both binding and catalytic activity only to a minor degree. On the basis of the information obtained from the mutagenesis study, we docked N(alpha)-(5'-phosphopyridoxyl)-D-alanine into the antigen-binding site. According to this model, Arg(L91) binds the alpha-carboxylate group of the amino acid substrate and Tyr(H100b) plays an essential role in the catalytic mechanism of antibody 15A9 by facilitating the Calpha/C4' prototropic shift. In addition, the catalytic apparatus of antibody 15A9 revealed several mechanistic features that overlap with those of PLP-dependent enzymes.     

Kinetics and equilibria for the reactions of coenzymes with wild type and the Y70F mutant of Escherichia coli aspartate aminotransferase.

The Y70F mutant of aspartate aminotransferase has reduced affinity for coenzymes compared to the wild type. The equilibrium dissociation constants for pyridoxamine phosphate (PMP) holoenzymes, KPMPdiss, were determined from the association and dissociation rate constants to be 1.3 nM and 30 nM for the wild type and mutant, respectively. This increase in KPMPdiss for Y70F is due to a 27-fold increase in the dissociation rate constant. Pyridoxal phosphate (PLP) association kinetics are complex, with three kinetic processes detectable for wild type and two for Y70F. A directly determined, accurate value of KPLPdiss for wild type enzyme has been difficult to obtain because of the low value of this constant. The values of KPLPdiss for the holoenzymes were determined indirectly through the measured values for KPMPdiss, glutamate-alpha-ketoglutarate half-reaction equilibrium constants, and the equilibrium constant for the transamination of PLP by glutamate catalyzed by Y70F. The values of KPLPdiss obtained by this procedure are 0.4 pM for wild type and 40 pM for Y70F. The increases in KPMPdiss and KPLPdiss for Y70F correspond to delta delta G values of 1.9 and 2.7 kcal/mol, respectively, and are directly attributed to the loss of the hydrogen bond from the phenolic hydroxyl group of Tyr70 to the coenzyme phosphate. The delta G for association of PLP with wild type enzyme is 4.7 kcal/mol more favorable than that for PMP.     

Trans-substitution of the proximal hydrogen bond in myoglobin: II. Energetics, functional consequences, and implications for hemoglobin allostery

The trans-substituted histidine to glycine mutant of sperm whale myoglobin (H93G Mb) is used to study energetics of proximal hydrogen bonding, proximal ligand-heme interactions, and coupling to distal ligand binding. Comparison of mono- and dimethylimidazole structural isomers shows that the hydrogen bond between the proximal ligand and the neighboring Ser92 hydroxyl (position F7) is stabilizing. The range of hydrogen bond stabilities measured here for different distal ligand complexes ranges from -0.7 kcal/mol (monomethylimidazole isomers to MbCO) to -4.1 kcal/mol (dimethylimidazole isomers to MbCN). This range of hydrogen bond stabilities, which is similar to that seen in protein mutagenesis unfolding studies, demonstrates the high sensitivity of the hydrogen bond to modest structural perturbations. The degree to which the 2-methyl group destabilizes proximal ligand binding is found to depend inversely on the total electronic spin. For monomethylimidazole proximal ligands, distal ligand binding weakens the proximal hydrogen bond compared to deoxyMb. Surprisingly, this trend is largely reversed for the dimethylimidazole proximal ligands. These results demonstrate strong coupling between the proximal protein matrix and distal ligand binding. These results provide an explanation for the strong avoidance of hydrogen bonding residues at position F7 in hemoglobin sequences.     

Molecular basis of in vitro affinity maturation and functional evolution of a neutralizing anti-human GM-CSF antibody

X-ray structure analysis of 4 antibody Fab fragments, each in complex with human granulocyte macrophage colony stimulating factor (GM-CSF), was performed to investigate the changes at the protein-protein binding interface during the course of in vitro affinity maturation by phage display selection. The parental antibody MOR03929 was compared to its derivatives MOR04252 (CDR-H2 optimized), MOR04302 (CDR-L3 optimized) and MOR04357 (CDR-H2 and CDR-L3 optimized). All antibodies bind to a conformational epitope that can be divided into 3 sub-epitopes. Specifically, MOR04357 binds to a region close to the GM-CSF N-terminus (residues 11–24), a short second sub-epitope (residues 83–89) and a third at the C-terminus (residues 112–123). Modifications introduced during affinity maturation in CDR-H2 and CDR-L3 led to the establishment of additional hydrogen bonds and van der Waals contacts, respectively, providing a rationale for the observed improvement in binding affinity and neutralization potency. Once GM-CSF is complexed to the antibodies, modeling predicts a sterical clash with GM-CSF binding to GM-CSF receptor α and β chain. This predicted mutually exclusive binding was confirmed by a GM-CSF receptor α chain ligand binding inhibition assay. Finally, high throughput sequencing of clones obtained after affinity maturation phage display pannings revealed highly selected consensus sequences for CDR-H2 as well for CDR-L3, which are in accordance with the sequence of the highest affinity antibody MOR04357. The resolved crystal structures highlight the criticality of these strongly selected residues for high affinity interaction with GM-CSF.     

Ser45 plays an important role in managing both the equilibrium and transition state energetics of the streptavidin-biotin system

The contribution of the Ser45 hydrogen bond to biotin binding activation and equilibrium thermodynamics was investigated by biophysical and X-ray crystallographic studies. The S45A mutant exhibits a 1,700-fold greater dissociation rate and 907-fold lower equilibrium affinity for biotin relative to wild-type streptavidin at 37 degrees C, indicating a crucial role in binding energetics. The crystal structure of the biotin-bound mutant reveals only small changes from the wild-type bound structure, and the remaining hydrogen bonds to biotin retain approximately the same lengths. No additional water molecules are observed to replace the missing hydroxyl, in contrast to the previously studied D128A mutant. The equilibrium deltaG degrees, deltaH degrees, deltaS degrees, deltaC degrees(p), and activation deltaG++ of S45A at 37 degrees C are 13.7+/-0.1 kcal/mol, -21.1+/-0.5 kcal/mol, -23.7+/-1.8 cal/mol K, -223+/-12 cal/mol K, and 20.0+/-2.5 kcal/mol, respectively. Eyring analysis of the large temperature dependence of the S45A off-rate resolves the deltaH++ and deltaS++ of dissociation, 25.8+/-1.2 kcal/mol and 18.7+/-4.3 cal/mol K. The large increases of deltaH++ and deltaS++ in the mutant, relative to wild-type, indicate that Ser45 could form a hydrogen bond with biotin in the wild-type dissociation transition state, enthalpically stabilizing it, and constraining the transition state entropically. The postulated existence of a Ser45-mediated hydrogen bond in the wild-type streptavidin transition state is consistent with potential of mean force simulations of the dissociation pathway and with molecular dynamics simulations of biotin pullout, where Ser45 is seen to form a hydrogen bond with the ureido oxygen as biotin slips past this residue after breaking the native hydrogen bonds.     

Protein-carbohydrate interactions defining substrate specificity in Bacillus 1,3-1,4-beta-D-glucan 4-glucanohydrolases as dissected by mutational analysis

The carbohydrate-binding site of Bacillus macerans 1,3-1, 4-beta-D-glucan 4-glucanohydrolase has been analyzed through a mutational analysis to probe the role of protein-carbohydrate interactions defining substrate specificity. Amino acid residues involved in substrate binding were proposed on the basis of a modeled enzyme-substrate complex [Hahn, M., Keitel, T., and Heinemann, U. (1995) Eur. J. Biochem. 232, 849-859]. The effects of the mutations at 15 selected residues on catalysis and binding were determined by steady-state kinetics using a series of chromogenic substrates of different degree of polymerization to assign the individual H-bond and hydrophobic contributions to individual subsites in the binding site cleft. The glucopyranose rings at subsites -III and -II are tightly bound by a number of H-bond interactions to Glu61, Asn24, Tyr92, and Asn180. From k(cat)/K(M) values, single H-bonds account for 1.8-2.2 kcal mol(-)(1) transition-state (TS) stabilization, and a charged H-bond contributes up to 3.5 kcal mol(-)(1). Glu61 forms a bidentated H-bond in subsites -III and -II, and provides up to 6.5 kcal mol(-)(1) TS stabilization. With a disaccharide substrate that fills subsites -I and -II, activation kinetics were observed for the wild-type and mutant enzymes except for mutations on Glu61, pointing to an important role of the bidentate interaction of Glu61 in two subsites. Whereas removal of the hydroxyl group of Tyr121, initially proposed to hydrogen-bond with the 2OH of Glcp-I, has essentially no effect (Y121F mutant), side-chain removal (Y121A mutant) gave a 100-fold reduction in k(cat)/K(M) and a 10-fold lower K(I) value with a competitive inhibitor. In subsite -IV, only a stacking interaction with Tyr22 (0.7 kcal mol(-)(1) TS stabilization) is observed.     

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.     

Free-energy simulations reveal that both hydrophobic and polar interactions are important for influenza hemagglutinin antibody binding

Antibodies binding to conserved epitopes can provide a broad range of neutralization to existing influenza subtypes and may also prevent the propagation of potential pandemic viruses by fighting against emerging strands. Here we propose a computational framework to study structural binding patterns and detailed molecular mechanisms of viral surface glycoprotein hemagglutinin (HA) binding with a broad spectrum of neutralizing monoclonal antibody fragments (Fab). We used rigorous free-energy perturbation (FEP) methods to calculate the antigen-antibody binding affinities, with an aggregate underlying molecular-dynamics simulation time of several microseconds (～2 μs) using all-atom, explicit-solvent models. We achieved a high accuracy in the validation of our FEP protocol against a series of known binding affinities for this complex system, with <0.5 kcal/mol errors on average. We then introduced what to our knowledge are novel mutations into the interfacial region to further study the binding mechanism. We found that the stacking interaction between Trp-21 in HA2 and Phe-55 in the CDR-H2 of Fab is crucial to the antibody-antigen association. A single mutation of either W21A or F55A can cause a binding affinity decrease of ΔΔG > 4.0 kcal/mol (equivalent to an ～1000-fold increase in the dissociation constant K_d). Moreover, for group 1 HA subtypes (which include both the H1N1 swine flu and the H5N1 bird flu), the relative binding affinities change only slightly (< ±1 kcal/mol) when nonpolar residues at the αA helix of HA mutate to conservative amino acids of similar size, which explains the broad neutralization capability of antibodies such as F10 and CR6261. Finally, we found that the hydrogen-bonding network between His-38 (in HA1) and Ser-30/Gln-64 (in Fab) is important for preserving the strong binding of Fab against group 1 HAs, whereas the lack of such hydrogen bonds with Asn-38 in most group 2 HAs may be responsible for the escape of antibody neutralization. These large-scale simulations may provide new insight into the antigen-antibody binding mechanism at the atomic level, which could be essential for designing more-effective vaccines for influenza.     

Crystal structure of Fab198, an efficient protector of the acetylcholine receptor against myasthenogenic antibodies.

The crystal structure of the Fab fragment of the rat monoclonal antibody 198, with protective activity for the main immunogenic region of the human muscle acetylcholine receptor against the destructive action of myasthenic antibodies, has been determined and refined to 2.8 A resolution by X-ray crystallographic methods. The mouse anti-lysozyme Fab D1.3 was used as a search model in molecular replacement with the AMORE software. The complementarity determining regions (CDR)-L2, CDR-H1 and CDR-H2 belong to canonical groups. Loops CDR-L3, CDR-H2 and CDR-H3, which seem to make a major contribution to binding, were analyzed and residues of potential importance for antigen-binding are examined. The antigen-binding site was found to be a long crescent-shaped crevice. The structure should serve as a model in the rational design of very high affinity humanized mutants of Fab198, appropriate for therapeutic approaches in the model autoimmune disease myasthenia gravis.     

Molecular modeling of cardiac glycoside binding by the human sequence monoclonal antibody 1B3

The amino acid sequences of the heavy- and light-chain variable regions of the high-affinity human sequence antidigoxin monoclonal antibody 1B3 (mAb 1B3) were determined, and a structural model for the mAb's variable region was developed by homology modeling techniques. The structural model provided the basis for computationally docking digoxin and eight related cardiac glycosides into the putative binding site of mAb 1B3. Analysis of the consensus binding mode obtained for digoxin showed that the cardenolide moiety of digoxin is deeply embedded in a predominantly hydrophobic, narrow cavity, whereas the terminal, gamma-carbohydrate group is solvent-exposed. The docking results indicated that the primary driving forces for digoxin binding by mAb 1B3 are hydrophobic interactions with the digoxin steroid ring system and hydrogen bonds with the digitoxose groups. The binding model accounts for the experimentally observed variations in mAb 1B3 binding affinity for various structural analogs of digoxin used previously to develop a 3D structure-activity relationship model of drug binding (Farr CD, Tabet MR, Ball WJ Jr, Fishwild DM, Wang X, Nair AC, Welsh WJ. Three-dimensional quantitative structure-activity relationship analysis of ligand binding to human sequence antidigoxin monoclonal antibodies using comparative molecular field analysis. J Med Chem 2002;45:3257-3270). In particular, the hydrogen bond pattern is consistent with the unique sensitivity of mAb 1B3's binding affinity to the number of sugar residues present in a cardiac glycoside. The hydrophobic environment about the steroid moiety of digoxin is compatible with the mAb's reduced affinity for ligands that possess hydrophilic hydroxyl and acetyl group modifications in this region. The model also indicated that most of the amino acid residues in contact with the ligand reside in or about the three complementarity determining regions (CDRs) of the heavy chain and the third CDR of the light chain. A comparison of the 1B3 binding model with the crystal structures of two murine antidigoxin mAbs revealed similar binding patterns used by the three mAbs, such as a high frequency of occurrence of aromatic, hydrophobic residues in the CDRs and a dominant role of the heavy chain CDR3 in antigen binding.