Search PubMed⌕ Search

Biomedical subjects

R E Bruccoleri

Publications and source records attributed to R E Bruccoleri.

At least 19 recordsLinked to original sources
PubMed1998

Concordance analysis of microbial genomes.

The set of proteins which are conserved across families of microbes contain important targets of new anti-microbial agents. We have developed a simple and efficient computational tool which determines concordances of putative gene products that show sets of proteins conserved across one set of user specified genomes and not present in another set of user specified genomes. The thresholds and the homology scoring criterion are selectable to allow the user to decide the stringency of the homologies. The system uses a relational database to store protein coding regions from different genomes, and to store the results of a complete comparison of all sequences against all sequences using the FASTA program. Using Web technology, the display of all the related proteins for a given sequence and calculation of multiple sequence alignments (using CLUSTALW) can be performed with the click of a button. The current database holds 97 365 sequences from 19 complete or partial genomes and 8798905 FASTA comparison results. A example concordance is presented which demonstrates that the target of the quinolone antibiotics could have been identified using this tool.

Bacteria↗
PubMed1998

Empirical free energy calculations of phage 434 repressor- and cro-DNA complexes support the 'indirect readout' hypothesis of specificity.

Empirical free energy calculations were done for two X-ray DNA-protein complexes. Complex formation energies were estimated for the phage 434 repressor and the 434 cro protein, both in complex with wild-type versions of their respective nucleic acid binding sites. The resulting quantities were subdivided according to the relative contributions of hydrophobicity, electrostatics and entropy (including both conformational changes and those pertinent to physical association of complex components). The non-linearized finite difference Poisson-Boltzmann equation was solved over a range of salt concentrations to determine magnitudes of coulombic and desolvation effects. Though approximate at best, results confirmed (1) the importance of centrally located charged protein residues that do not make direct contact with the DNA, (2) explicit protein backbone interactions with DNA, and (3) the varied parts of the two proteins used to harness the energy of binding. Although a large difference exists between calculated and experimentally measured free energy values, these can be reconciled by assuming that approximately 1.5 e.u. (about 0.4-0.5 kcal at room temperature) of DNA vibrational and/or conformational entropy is lost for every basepair torsion in contact with the protein. Such a conclusion supports an "indirect readout" ("compound recognition") hypothesis of DNA binding specificity, wherein sequence-dependent differences in protein-DNA interaction are due to the topological peculiarities of the interactions rather than intrinsic chemical differences alone.

Bacteriophage lambda↗
PubMed1997

Empirical free energy calculations: a blind test and further improvements to the method.

Empirical Gibbs functions estimate free energies of non-covalent reactions (deltaG) from atomic coordinates of reaction products (e.g. antibody-antigen complexes). The function previously developed by us has four terms that quantify the effects of hydrophobic, electrostatic and entropy changes (conformational, association) upon complexation. The function was used to calculate delta deltaG of ten lysozyme mutants affecting the stability of the HyHEL-10 antibody-lysozyme complex. The mutants were computer-modeled from the X-ray structure of the wild-type, and free energy calculations produced a correlation coefficient of 0.5 with the experimental delta deltaG data (average absolute error +/-3 kcal). The following changes were then introduced into the Gibbs function: (1) the hydrophobic force was made proportional to the molecular surface, as calculated by the GEPOL93 algorithm, with the scaling constant of 70 cal/mol/A2; (2) calculation of the electrostatics of binding was carried out by the finite difference Poisson-Boltzmann algorithm, which employed uniform grid charging, dielectric boundary smoothing and charge anti-aliasing; and (3) side-chain conformational entropy was estimated from the CONGEN sampling of torsional degrees of freedom. In the new calculations, correlation with experimental data improved to 0.6 or 0.8 if a single outlying mutant, K96M, was neglected. Analysis of the errors remaining in our calculations indicated that molecular mechanics-based modeling of the mutants, rather than the form of our amended Gibbs function, was the main factor limiting the accuracy of the free energy estimates.

Animals↗
PubMed1997

Homology modeling using simulated annealing of restrained molecular dynamics and conformational search calculations with CONGEN: application in predicting the three-dimensional structure of murine homeodomain Msx-1.

We have developed an automatic approach for homology modeling using restrained molecular dynamics and simulated annealing procedures, together with conformational search algorithms available in the molecular mechanics program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168). The accuracy of the method is validated by "predicting" structures of two homeodomain proteins with known three-dimensional structures, and then applied to predict the three-dimensional structure of the homeodomain of the murine Msx-1 transcription factor. Regions of the unknown protein structure that are highly homologous to the known template structure are constrained by "homology distance constraints," whereas the conformations of nonhomologous regions of the unknown protein are defined only by the potential energy function. A full energy function (excluding explicit solvent) is employed to ensure that the calculated structures have good conformational energies and are physically reasonable. As in NMR structure determinations, information on the consistency of the structure prediction is obtained by superposition of the resulting family of protein structures. In this paper, our homology modeling algorithm is described and compared with related homology modeling methods using spatial constraints derived from the structures of homologous proteins. The software is then used to predict the DNA-bound structures of three homeodomain proteins from the X-ray crystal structure of the engrailed homeodomain protein (Kissinger CR et al., 1990, Cell 63:579-590). The resulting backbone and side-chain conformations of the modeled yeast Mat alpha 2 and D. melanogaster Antennapedia homeodomains are excellent matches to the corresponding published X-ray crystal (Wolberger C et al., 1991, Cell 67:517-528) and NMR (Billeter M et al., 1993, J Mol Biol 234:1084-1097) structures, respectively. Examination of these structures of Msx-1 reveals a network of highly conserved surface salt bridges that are proposed to play a role in regulating protein-protein interactions of homeodomains in transcription complexes.

Algorithms↗
PubMed1996

Simulated annealing with restrained molecular dynamics using CONGEN: energy refinement of the NMR solution structures of epidermal and type-alpha transforming growth factors.

The new functionality of the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168; Bassolino-Klimas D et al., 1996, Protein Sci 5:593-603) has been applied for energy refinement of two previously determined solution NMR structures, murine epidermal growth factor (mEGF) and human type-alpha transforming growth factor (hTGF alpha). A summary of considerations used in converting experimental NMR data into distance constraints for CONGEN is presented. A general protocol for simulated annealing with restrained molecular dynamics is applied to generate NMR solution structures using CONGEN together with real experimental NMR data. A total of 730 NMR-derived constraints for mEGF and 424 NMR-derived constraints for hTGF alpha were used in these energy-refinement calculations. Different weighting schemes and starting conformations were studied to check and/or improve the sampling of the low-energy conformational space that is consistent with all constraints. The results demonstrate that loosened (i.e., "relaxed") sets of the EGF and hTGF alpha internuclear distance constraints allow molecules to overcome local minima in the search for a global minimum with respect to both distance restraints and conformational energy. The resulting energy-refined structures of mEGF and hTGF alpha are compared with structures determined previously and with structures of homologous proteins determined by NMR and X-ray crystallography.

Animals↗
PubMed1996

Simulated annealing with restrained molecular dynamics using a flexible restraint potential: theory and evaluation with simulated NMR constraints.

A new functional representation of NMR-derived distance constraints, the flexible restraint potential, has been implemented in the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168) for molecular structure generation. In addition, flat-bottomed restraint potentials for representing dihedral angle and vicinal scalar coupling constraints have been introduced into CONGEN. An effective simulated annealing (SA) protocol that combines both weight annealing and temperature annealing is described. Calculations have been performed using ideal simulated NMR constraints, in order to evaluate the use of restrained molecular dynamics (MD) with these target functions as implemented in CONGEN. In this benchmark study, internuclear distance, dihedral angle, and vicinal coupling constant constraints were calculated from the energy-minimized X-ray crystal structure of the 46-amino acid polypeptide crambin (ICRN). Three-dimensional structures of crambin that satisfy these simulated NMR constraints were generated using restrained MD and SA. Polypeptide structures with extended backbone and side-chain conformations were used as starting conformations. Dynamical annealing calculations using extended starting conformations and assignments of initial velocities taken randomly from a Maxwellian distribution were found to adequately sample the conformational space consistent with the constraints. These calculations also show that loosened internuclear constraints can allow molecules to overcome local minima in the search for a global minimum with respect to both the NMR-derived constraints and conformational energy. This protocol and the modified version of the CONGEN program described here are shown to be reliable and robust, and are applicable generally for protein structure determination by dynamical simulated annealing using NMR data.

Computer Simulation↗
PubMed1996

Factors influencing accuracy of computer-built models: a study based on leucine zipper GCN4 structure.

A three-dimensional model of the leucine zipper GCN4 built from its amino acid sequence had been reported previously by us. When the two alternative x-ray structures of the GCN4 dimer became available, the root mean square (r.m.s.) shifts between our model and the structures were determined as approximately 2.7 A on all atoms. These values are similar to the r.m.s. shift of 2.8 A between the two GCN4 structures in the different crystal forms (C2 and P2(1)2(1)2(1)). CONGEN conformational searches were run to better understand the conditions that may determine the preference of different conformers in different environments and to test the sensitivity of our current modeling techniques. With a judicious choice of CONGEN search parameters, the backbone r.m.s. deviation improved to 0.8 A and 2.5 A on all atoms. The side-chain conformations of Val and Leu at the helical interface were well reproduced (1.2 A r.m.s.), and the large side-chain misplacements occurred with only a small number of charged amino acids and a tyrosine. Inclusion of the crystal environment (C2 symmetry), as a passive background, into the side-chain conformational search further improved the accuracy of the model to an r.m.s. deviation of 2.1 A. Conformational searches carried out in the two different crystal environments and employing the AMBER protein/DNA forcefield, as implemented in CONGEN, gave the r.m.s. values of 2.2 A (for the C2 symmetry) and 2.5 A (for the P2(1)2(1)2(1) symmetry). In the C2 symmetry crystal, as much as 40% of the surface of each dimer was involved in crystal contacts with other dimers, and the charged residues on the surface often interacted with immobilized water molecules. Thus, occasional large r.m.s. deviations between the model and the x-ray side chains were due to specific conditions that did not occur in solution.

Amino Acid Sequence↗
PubMed1995

High-resolution solution structure of siamycin II: novel amphipathic character of a 21-residue peptide that inhibits HIV fusion.

The 21-amino acid peptides siamycin II (BMY-29303) and siamycin I (BMY-29304), derived from Streptomyces strains AA3891 and AA6532, respectively, have been found to inhibit HIV-1 fusion and viral replication in cell culture. The primary sequence of siamycin II is CLGIGSCNDFAGCGYAIVCFW. Siamycin I differs by only one amino acid; it has a valine residue at position 4. In both peptides, disulfide bonds link Cys1 with Cys13 and Cys7 with Cys19, and the side chain of Asp9 forms an amide bond with the N-terminus. Siamycin II, when dissolved in a 50:50 mixture of DMSO and H2O, yields NOESY spectra with exceptional numbers of cross peaks for a peptide of this size. We have used 335 NOE distance constraints and 13 dihedral angle constraints to generate an ensemble of 30 siamycin II structures; these have average backbone atom and all heavy atom rmsd values to the mean coordinates of 0.24 and 0.52 A, respectively. The peptide displays an unusual wedge-shaped structure, with one face being predominantly hydrophobic and the other being predominantly hydrophilic. Chemical shift and NOE data show that the siamycin I structure is essentially identical to siamycin II. These peptides may act by preventing oligomerization of the HIV transmembrane glycoprotein gp41, or by interfering with interactions between gp41 and the envelope glycoprotein gp120, the cell membrane or membrane-bound proteins [Frèchet, D. et al. (1994) Biochemistry, 33, 42-50]. The amphipathic nature of siamycin II and siamycin I suggests that a polar (or apolar) site on the target protein may be masked by the apolar (or polar) face of the peptide upon peptide/protein complexation.

Amino Acid Sequence↗
PubMed1994

Contribution of a single heavy chain residue to specificity of an anti-digoxin monoclonal antibody.

Two distinct spontaneous variants of the murine anti-digoxin hybridoma 26-10 were isolated by fluorescence-activated cell sorting for reduced affinity of surface antibody for antigen. Nucleotide and partial amino acid sequencing of the variant antibody variable regions revealed that 1 variant had a single amino acid substitution: Lys for Asn at heavy chain position 35. The second variant antibody had 2 heavy chain substitutions: Tyr for Asn at position 35, and Met for Arg at position 38. Mutagenesis experiments confirmed that the position 35 substitutions were solely responsible for the markedly reduced affinity of both variant antibodies. Several mutants with more conservative position 35 substitutions were engineered to ascertain the contribution of Asn 35 to the binding of digoxin to antibody 26-10. Replacement of Asn with Gln reduced affinity for digoxin 10-fold relative to the wild-type antibody, but maintained wild-type fine specificity for cardiac glycoside analogues. All other substitutions (Val, Thr, Leu, Ala, and Asp) reduced affinity by at least 90-fold and caused distinct shifts in fine specificity. The Ala mutant demonstrated greatly increased relative affinities for 16-acetylated haptens and haptens with a saturated lactone. The X-ray crystal structure of the 26-10 Fab in complex with digoxin (Jeffrey PD et al., 1993, Proc Natl Acad Sci USA 90:10310-10314) reveals that the position 35 Asn contacts hapten and forms hydrogen bonds with 2 other contact residues. The reductions in affinity of the position 35 mutants for digoxin are greater than expected based upon the small hapten contact area provided by the wild-type Asn. We therefore performed molecular modeling experiments which suggested that substitution of Gln or Asp can maintain these hydrogen bonds whereas the other substituted side chains cannot. The altered binding of the Asp mutant may be due to the introduction of a negative charge. The similarities in binding of the wild-type and Gln-mutant antibodies, however, suggest that these hydrogen bonds are important for maintaining the architecture of the binding site and therefore the affinity and specificity of this antibody. The Ala mutant eliminates the wild-type hydrogen bonding, and molecular modeling suggests that the reduced side-chain volume also provides space that can accommodate a congener with a 16-acetyl group or saturated lactone, accounting for the altered fine specificity of this antibody.

Amino Acid Sequence↗
PubMed1993

Heavy chain position 50 is a determinant of affinity and specificity for the anti-digoxin antibody 26-10.

Antibody produced by a variant of the murine antidigoxin hybridoma 26-10 has reduced affinity for digoxin but enhanced recognition of the digoxin 12-hydroxyl due to a Tyr to His substitution at heavy chain position 50 (Schildbach, J. F., Panka, D. J., Parks, D. R., Jager, G. C., Novotny, J., Herzenberg, L. A., Mudgett-Hunter, M., Bruccoleri, R. E., Haber, E., and Margolies, M. N. (1991) J. Biol. Chem. 266, 4640-4647). Consistent with these data, the 26-10 Fab-digoxin x-ray crystal structure (Jeffrey, P. D., Strong, R. K., Sieker, L. C., Chang, C. Y., Campbell, R. L., Petsko, G. A., Haber, E., Margolies, M. N., and Sheriff, S. (1993) Proc. Natl. Acad. Sci. U. S. A., in press) reveals that Tyr-50 contacts a region of digoxin that includes the hapten-12 carbon. To determine the effects of other heavy chain position 50 substitutions, mutant antibodies were engineered, and their affinities for digoxin and digoxin analogues were measured. The affinity of the mutant antibodies for digoxin roughly correlates with the size of the position 50 side chain. Substitutions of Trp or Phe have no effect on affinity, whereas substitutions of Asn, His, Leu, Ala, Gly, and Asp confer progressively lower affinities. Although Trp and Phe mutants exhibit wild-type specificity, Asn and Asp mutants have improved affinity for digoxin relative to digitoxin (12-deshydroxydigoxin). Leu, Ala, and Gly mutants have improved affinity for 12-acetyldigoxin relative to digoxin as compared with 26-10. These results indicate that position 50 is a determinant of both antibody affinity and fine specificity for antibody 26-10 and that single-amino acid substitutions can alter antibody fine specificity. Models of the mutants were computationally constructed, and haptens were docked into the modeled binding sites. The results suggest that 12-acetyldigoxigenin occupies different orientations in the 26-10 and in the Ala mutant binding sites, resulting in altered binding.

Amino Acids↗
PubMed1993

Modulation of antibody affinity by a non-contact residue.

Antibody LB4, produced by a spontaneous variant of the murine anti-digoxin monoclonal antibody 26-10, has an affinity for digoxin two orders of magnitude lower than that of the parent antibody due to replacement of serine with phenylalanine at position 52 of the heavy chain variable region (Schildbach, J.F., Panka, D.J., Parks, D.R., et al., 1991, J. Biol. Chem. 266, 4640-4647). To examine the basis for the decreased affinity, a panel of engineered antibodies with substitutions at position 52 was created, and their affinities for digoxin were measured. The antibody affinities decreased concomitantly with increasing size of the substituted side chains, although the shape of the side chains also influenced affinity. The crystal structure of the 26-10 Fab complexed with digoxin (P.D.J., R.K. Strong, L.C. Sieker, C. Chang, R.L. Campbell, G.A. Petsko, E.H., M.N.M., & S.S., submitted for publication) shows that the serine at heavy chain position 52 is not in contact with hapten, but is adjacent to a tyrosine at heavy chain position 33 that is a contact residue. The mutant antibodies were modeled by applying a conformational search procedure to position side chains, using the 26-10 Fab crystal structure as a starting point. The results suggest that each of the substituted side chains may be accommodated within the antibody without substantial structural rearrangement, and that none of these substituted side chains are able to contact hapten. These modeling results are consistent with the substituents at position 52 having only an indirect influence upon antibody affinity.(ABSTRACT TRUNCATED AT 250 WORDS)

Animals↗
PubMed1992

Solution conformation of a cyclic pentapeptide endothelin antagonist. Comparison of structures obtained from constrained dynamics and conformational search.

The structure of a cyclic pentapeptide, cyclo-(D-Trp-D-Asp-L-Pro-D-Val-L-Leu), that has high selectively for the endothelin ETA receptor has been determined by NMR spectroscopy using constrained molecular dynamics and conformational search procedures. Structures obtained using two methods of refinement, namely (i) constrained molecular dynamics; and (ii) systematic searches of conformational space for optimal satisfaction of distance constraints, were compared to those obtained from systematic searches of conformational space without NMR data. The two different procedures of refinement produce similar conformations that are consistent with the NMR distance constraints. Conformational searches for optimal energy without any NMR distance constraints produced several low-energy structures, two of which have essentially the same backbone as those structures derived from distance-constrained procedures and one of these even reproduces several side-chain positions well. The pentapeptide backbone consists of a linked gamma- and beta-turn conformation, with the leucine and tryptophan as corner residues of the type II beta-turn. The side chains are highly ordered both in aqueous solvent and in dimethyl sulfoxide. In aqueous media the leucine side chain is directed towards the indole ring, presumably to reduce the non-polar surface exposure, producing unusual upfield shifts for the methyls (and particularly H gamma). This structural feature was reproduced in one of the structures obtained from conformational searches performed without NMR data. Exhaustive conformational searches appear to provide an alternative method for structure generation for cyclic peptides.

Endothelins↗
PubMed1992

Modeling the antigen combining site of an anti-dinitrophenyl antibody, ANO2.

A model structure has been constructed for a monoclonal anti-dinitrophenyl antibody. The antibody, ANO2, has been sequenced and cloned (Anglister, J., Frey, T., & McConnell, H.M., 1984, Biochemistry 23, 1138-1142). Its amino acid sequence shows striking homology with the anti-lysozyme Fab fragments HyHel5 (83%) and HyHel10 (73%). Based on this homology, a model for the ANO2 variable heavy and variable light chain framework was constructed using a hybrid of the HyHel5 light chain and the HyHel10 heavy chain backbone, omitting the hypervariable loops. These coordinates were used as scaffolds for the model building of ANO2. The CONGEN conformational sampling algorithm (Bruccoleri, R.E. & Karplus, M., 1987, Biopolymers 26, 127-196) was used to model the six hypervariable loops that contain the antigen-combining site. All the possible conformations of the loop backbones were constructed and the best loop structures were selected using a combination of the CHARMM potential energy function and evaluation of the solvent-accessible surface area of the conformers. The order in which the loops were searched was carried out based on the relative locations of the loops with reference to the framework of the beta-barrel, namely, L2-H1-L3-H2-H3-L1. The model structures thus obtained were compared to the high resolution X-ray structure (Brünger, A.T., Leahy, D.J., Hynes, T.R., & Fox, R.O., 1991, J. Mol. Biol. 221, 239-256).

Algorithms↗
PubMed1992

Application of a directed conformational search for generating 3-D coordinates for protein structures from alpha-carbon coordinates.

A directed conformational search algorithm using the program CONGEN (ref. 3), which samples backbone conformers, is described. The search technique uses information from the partially built structures to direct the search process and is tested on the problem of generating a full set of backbone Cartesian coordinates given only alpha-carbon coordinates. The method has been tested on six proteins of known structure, varying in size and classification, and was able to generate the original backbone coordinates with RMSs ranging from 0.30-0.87A for the alpha-carbons and 0.5-0.99A RMSs for the backbone atoms. Cis peptide linkages were also correctly identified. The procedure was also applied to two proteins available with only alpha-carbon coordinates in the Brookhaven Protein Data Bank; thioredoxin (SRX) and triacylglycerol acylhydrolase (TGL). All-atom models are proposed for the backbone of both these proteins. In addition, the technique was applied to randomized coordinates of flavodoxin to assess the effects of irregularities in the data on the final RMS. This study represents the first time a deterministic conformational search was used on such a large scale.

Algorithms↗
PubMed1991

Altered hapten recognition by two anti-digoxin hybridoma variants due to variable region point mutations.

Two spontaneous variants of the murine anti-digoxin antibody-producing hybridoma cell line 26-10 were isolated by two-color fluorescence-activated cell sorting on the basis of altered hapten binding. The variable region sequences of the antibodies produced by the mutant lines revealed that each contains a single amino acid change in the heavy chain second complementarity determining region. A Tyr to His change at position 50 leads to a 40-fold reduction in affinity for digoxin. A Ser to Phe mutation at position 52 results in a 300-fold reduction in affinity for digoxin. A competition assay involving 33 digoxin analogues was used to examine the specificity of hapten binding of 26-10 and the two mutant antibodies. The position 50 mutant has a distinct specificity change; it exhibits a preference for digoxin congeners containing a hydroxyl group at the steroid 12 position, whereas the 26-10 parent does not. The affinities of all three antibodies for hapten are progressively lowered by substitutions of increasing size at the digoxin steroid D ring 16 position. Although 26-10 binds digoxin and its genin form equally, 12 and 16 steroid position substitutions which lower affinity also confer a preference for a sugar at the steroid 3 position. These results suggest that position 50 contributes to specificity of the antibody and that alterations of the hapten can lead to differences in recognition, possibly through a shift in hapten orientation within the binding site.

Amino Acid Sequence↗
PubMed1991

Stabilities of leucine zipper dimers estimated by an empirical free energy method.

The leucine zipper motif is a characteristic amino acid sequence found in dimeric DNA-binding proteins. Computer-generated models for leucine zippers were constructed as alpha-helical coiled dimers with leucine repeated every seventh residue. An empirical Gibbs free energy, delta G, function which incorporates hydrophobic force, electrostatic interactions, and conformational entropy loss as the major intermolecular interactions was used to estimate the delta G of dimer formation in fos, jun, and GCN4 zipper sequences. The calculations showed that complexes known to form stable homo- or heterodimers have favorable (negative) delta G, while other less stable complexes have unfavorable (positive) delta G. Leucines in position d of the coiled coil contribute large hydrophobic stabilization energies while residues in the a position contribute less to dimer stability. Hydrophobic contributions show little sequence specificity, however, and do not contribute significantly to homo/heterodimer preference. Charged residues in the e and g positions, on the other hand, determine homo/heterodimer specificity. In GCN4 homodimers, residues GLU el, Glu b2, Lys g2, and Lys e4 greatly contribute to dimer stability. The preferential stability of fos-jun heterodimer over the jun-jun and fos-fos homodimers is primarily due to the side chains Asp b1, Glu g1, Asp b2, Glu e2, Glu g2, Glu g3, and Lys a5 of the fos helix, and Arg c1, Lys g1, Lys b2, Lys e2, Arg e4, and Glu g4 of the jun helix.

Leucine Zippers↗
PubMed1990

A V kappa-J kappa junctional change in an antidigoxin recombinant antibody destroys digoxin-binding activity.

A set of high affinity antidigoxin antibodies were previously identified with high homologous V kappa 1A L chain sequences but were associated with two entirely different VH regions and two dramatically different specificities for digoxin analogs. Antibodies 40-20, 40-60, 40-90, and 40-100 displayed similar binding specificities but differed from that of antibody 26-10. In a previous study using somatic cell fusion for Ig chain recombination we demonstrated that a recombinant antibody consisting of the H chain of antibody 26-10 and the L chain of antibody 40-20 retained digoxin binding and the 26-10 Id, but displayed a binding specificity pattern dominated by the 26-10 H chain donor. In the present study we produced three additional chain recombinant antibodies that contain the 26-10 H chain recombined with each of the L chains of antibodies 40-60, 40-90, and 40-100. All four recombinants expressed the 26-10 Id indistinguishably from the 26-10 antibody. Two of the recombinants (using the 40-60 and 40-90 L chains) bind digoxin; however, the recombinant using the 40-100 L chain failed to bind digoxin. Complete sequence analyses of the 40-20, 40-60, 40-90, and 40-100 VH and VL regions were performed. Antibodies 40-90 and 40-100 have identical VH region sequences but differed only in their L chains at position 96 (proline/leucine). This single difference at the VK-JK junction abolished digoxin binding in the context of one H chain (26-10), but does not cause a significant change in binding in association with the "normal" parental chains 40-90 and 40-100. Thus, structurally closely related VL regions can recombine with different VH regions to form digoxin binding sites of different specificity; in one binding site the identity of a L chain junctional residue is critical whereas in the second binding site that residue is unimportant. Molecular modeling studies revealed major differences between calculated binding site structures for 26-10 when leucine is substituted for proline at position 96 in the 26-10 VL region.

Amino Acid Sequence↗