Department of Chemistry, University of Kuopio, PO Box 1627, FIN-70211 Kuopio, Finland
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Keywords: antibody/free energy/molecular dynamics/mutation/steroid
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Characterization of antibodies by protein X-ray crystallography and other biochemical methods has provided a lot of data on structure and function of antibodies (Arevalo et al., 1994; Wilson and Stanfield, 1994
; Charbonnier et al., 1997
; Gigant et al., 1997
; Wedemayer et al., 1997
). The structures of several steroid-binding antibodies have been determined by X-ray crystallography, providing an atomic level understanding of steroid binding (Arevalo et al., 1994
; Jeffrey et al., 1993
; Trinh et al., 1997
). The X-ray structure of the monoclonal anti-progesterone (anti-PRG) antibody DB3, which is investigated in this study, in an unliganded form and complexed with PRG and several PRG-like steroids have been determined (Arevalo et al, 1993a
,b
,1994
). DB3 binds PRG with high affinity (Ka = 1x109 M) and cross-reacts with several PRG-like steroids with nanomolar affinities (Arevalo et al., 1993b
).
The two steroids studied, PRG and 5ß-androstane-3,17-dione (5AD, Figure 1), are of quite different shape. PRG is a flat molecule with an angle of 25° between the planes defined by the ring carbon atoms of the steroid A ring and B, C and D rings. In the case of 5AD the corresponding angle is 110°. Furthermore, the X-ray structures revealed that the two molecules bind to the ligand-binding site of DB3 in different orientations: the D ring of both molecules binds to the same binding pocket but the faces of the steroids point in the opposite directions (Figures 2 and 3
). In spite of these differences, the binding affinity of PRG (1 nM) is close to that of 5AD (8 nM, Arevalo et al., 1993b
).
|
|
|
![]() |
Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The Cornell et al. (Cornell et al., 1995) force field, which was used in the MD simulations of this work, does not have some bond, angle and torsion angle parameters needed for the steroids. The missing parameters compatible with the Cornell et al. force field were developed by following the guidelines of Fox and Kollman (Fox and Kollman, 1998
). The equilibrium bond distances were taken from the structures optimized at the HF/6-31G* and the corresponding bond stretching force constants were interpolated from similar type of bonds of the Cornell et al. force field. The angle parameters were taken from the corresponding angles of the Cornell et al. force field. The missing torsion parameters were developed using small model molecules. Torsion profiles were calculated for the models at the MP2/6-31G**//HF/6-31G* level and the torsion angle parameters were fitted to reproduce the quantum mechanical (QM) energy profile. Atomic point charges were calculated with the two-stage RESP method (Cornell et al., 1993
) at the HF/6-31G* level for the steroids and the model compounds using the geometries optimized at the same level. The parameters added to the parm96 parameter file of the AMBER5 program (Case et al., 1997), the QM and MM energy profiles and steroid charges are available as supplementary material.
Molecular dynamics simulations
The initial coordinates for the antibody complexes were obtained from the X-ray crystal structures of the Fab' fragment of the anti-PRG monoclonal antibody DB3 complexed with PRG (Protein Data Bank code: 1DBB; Arevalo et al., 1993a) and 5AD (Protein Data Bank code:1DBK; Arevalo et al., 1993b
) determined to 2.7 and 3.0 Å resolution, respectively. For the MD simulations the ligand-binding sites of the antibodyligand complexes were solvated by adding a sphere of TIP3P (Jorgensen et al., 1983
) water molecules with a 25 Å radius from the mass-centre of the ligands with the use of the cap option of the LEAP program (Schafmeister et al., 1995
). This resulted in 848 water molecules for the DB3PRG and 874 for the DB35AD complex. In the simulations, residues having atoms within 12 Å from the atoms of the ligand molecules were allowed to move. This resulted in 82 moving residues for the PRG and 92 residues for the 5AD simulation. The water molecules of the solvated systems were first energy-minimized for 1000 steps with the use of the conjugate gradient algorithm and then equilibrated by a 30 ps MD simulation at a constant temperature of 300 K with the use of the Berendsen temperature coupling algorithm (Berendsen et al., 1984
). After that, the water molecules and the moving part of the protein were energy-minimized for 1000 steps and equilibrated for 100 ps at 300 K. The production simulations of 500 ps were then started. In the simulations the non-bonded cut-off of 12 Å and a time step of 2.0 fs were used. The SHAKE algorithm (Ryckaert et al., 1977
) was used to constrain bond distances to their equilibrium values. During the 500 ps simulations structures were saved every 5 ps for further analyses. All MD simulations were done using the AMBER5 program with the Cornell et al. force field and parameters developed in this work.
Free energy analyses
The binding free energies were calculated using the MM-PBSA method (Kollman et al., 2000). In the method free energies are calculated for snapshot structures taken from the MD trajectory of the system studied. The average binding free energy (
GBIND) is calculated from the average molecular mechanical gas-phase energies (EMM), solvation free energies (
Gsolv) and entropy contributions (TS):
![]() | (1) |
![]() | (2) |
![]() | (3) |
In this work the energy contribution from entropy changes upon ligand binding was not included. This was justified by the fact that it is likely that entropy does not contribute much to the relative binding free energies of the two steroids to the same receptor. More importantly, there is no straightforward way to quantitatively calculate the entropy contribution to binding. The normal mode analysis, which is often used to estimate entropy changes, gives only qualitative estimates (Kollman et al., 2000).
The MM-PBSA method was used for the computational alanine scanning (Reyes and Kollman, 2000a) of selected binding-site residues. In the scanning, the structures of the MD trajectory of the proteinligand complex are mutated by truncating the side chains of the residues one at a time and binding free energies are calculated for each of the mutated trajectories. Binding free energies and computational mutagenesis of selected binding-site residues were calculated using single proteinligand trajectories. This means that the snapshot structures for the energy calculations of the proteinsteroid complex and separated protein and steroid were taken from the MD trajectory of the proteinsteroid complex. It has been observed that the single trajectory method provides fairly good estimates for the relative binding energies (Kollman et al., 2000
). From the 500 ps antibodyligand MD trajectories 50 snapshots were taken at even intervals for the binding energy analyses. The binding free energies reported are averages from the 50 snapshots.
![]() |
Results and discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The starting structures for the MD simulations of the DB3PRG and DB35AD complexes were obtained from the corresponding X-ray structures. In the simulations only the residues within 12 Å from the steroids were allowed to move and the moving part was solvated with a 25 Å sphere of water molecules. The production simulations of 500 ps carried out for these systems were stable on the basis of the total and potential energies of the systems (data not show) and the root-mean square (RMS) deviations from the X-ray structures (Figure 4). The magnitudes of the atomic fluctuations of the two systems are similar but the structure of the antibody5AD complex has an average RMS deviation of 1.1 Å as compared to the 0.6 Å of the antibodyPRG complex. This difference is most probably due to the larger number of moving amino acid residues in the 5AD (92) than the PRG simulation (82). The medium resolution antibodysteroid X-ray structures and use of a water sphere to solvate the simulation systems may also be responsible for the different behaviour of the systems. However, the deviations are comparable to those observed earlier for similar systems (Chong et al., 1999
). Visual inspection of the simulation trajectories showed that the residues of the ligand-binding pocket and the ligands stayed close to the positions they have in the X-ray structures. The X-ray structures of the complexes between the DB3 antibody and the steroids studied in this work and other similar steroids have been analysed in detail by Arevalo et al. (Arevalo et al., 1993a
,b
,1994
).
|
Free energies were calculated for the binding of PRG and 5AD to the DB3 antibody with the MM-PBSA method (Tables IIII). The calculated binding free energy, without the contribution of entropy, was 119.3 kJ/mol for PRG and 110.5 kJ/mol for 5AD. The experimental binding free energies are 51.4 and 46.2 kJ/mol for PRG and 5AD, respectively (Arevalo et al., 1993b
). The calculated relative binding free energy, 8.8 kJ/mol, is in fair agreement with the relative experimental energy, 5.4 kJ/mol. Note that the standard deviations of the calculated binding energies of PRG and 5AD are of the same magnitude as the relative energy. The relatively large statistical error in calculated binding energy is an inherent feature of the MM-PBSA method because the energy is calculated from averages of several large numbers. In spite of this, good agreement between the calculated and experimental relative binding energies have been obtained in several cases with this method (Kollman et al., 2000
). The inclusion of change in entropy upon association, T
S, would bring the calculated absolute energies closer to the experimental ones. The entropy contribution has been calculated to be 77.0 and 80.8 kJ/mol for a phosphonate hapten binding to the germ line and mature forms of the antibody 48G7 using normal-mode analysis in the gas-phase (Chong et al., 1999
). In addition, T
S has been calculated to be 4284 kJ/mol for the binding of a set of nine ligands to avidin (Kuhn and Kollman, 2000
). Since the calculation of entropy contribution requires several approximations and provides only rough estimates, especially in the case of a simulation in which only a part of the protein is moving, it was not calculated here.
|
|
|
The sum of the electrostatic proteinligand interaction energy and electrostatic solvation energy (GPB,elec) is 17.3 kJ/mol less favourable for the binding of PRG than 5AD. However, because of the 25.2 kJ/mol more favourable
GvdW, PRG is bound tighter to DB3 than 5AD. The more favourable van der Waals and proteinligand electrostatic interaction energies of PRG are in agreement with the larger number of van der Waals and hydrogen-bonding interactions observed in the X-ray structures of PRG than 5AD. There are 59 such interactions in the antibodyPRG complex but only 20 in the 5AD complex (Arevalo et al., 1993b
). The number of DB3 antibodysteroid interactions has been found to correlate with the slightly higher affinities of the steroids with a 5ß substituent or a double bond to carbon 5 (like PRG) compared to steroids with a 5
substituent (like 5AD; Arevalo et al., 1993b
). On the basis of the data for 5AD, the binding of other 5ß derivatives to the DB3 antibody is also favoured by solvation energies.
Computational alanine scanning of the ligand-binding residues
To study the role of individual amino acid residues on the cross-reactivity of DB3 with PRG and 5AD, selected residues of the ligand-binding pocket were mutated to alanines (Figure 5.). The residues mutated and the binding-site structures of the antibodyPRG and antibody5AD complexes are shown in Figures 2 and 3
. Arevalo et al. have made a detailed crystal structure analysis of the steroid binding to the DB3 antibody (Arevalo et al., 1994
). They divided the binding site into four different compartments (pockets) named P1, P2, P3 and P3'. The P1 compartment is the most hydrophobic cavity of the binding site and interacts with the steroid D ring. The side chain of NH35 located at the bottom of the P1 compartment provides a hydrogen-bond donor for the bound steroids. This interaction is important for the antibody specificity by enforcing a hydrogen-bond acceptor as the functional group of the steroid D ring. There is a greater degree of flexibility and less antigen complementarity in the P2 compartment. It makes interactions with the B, C and D rings of steroids. The P3 and P3' compartments are available for the binding of the two different conformations of the A ring.
|
There are two compartments in the ligand-binding cavity where the steroid skeleton can be placed. The walls of these compartments are mainly formed by WH50 and WH100 located on the opposite sides of the binding cavity. The two steroids are bound to the binding site in different orientations. The -face (the face with no CH3 groups) of PRG is stacked face-to-face with WH100, whereas the ß-face (the face with CH3 groups) of 5AD makes this interaction. These interactions have large contributions to the binding free energies. Mutation WH100 decreases the binding of PRG by 23.0 kJ/mol and that of 5AD by 20.1 kJ/mol. Mutation of WH50A has a 7.9 kJ/mol larger effect on the binding of 5AD than that of PRG. In addition to the different binding orientations of PRG and 5AD, a slight twist about the long axis of the steroid skeleton differentiates the two binding modes. This rotation allows extensive interactions between the
-faces of the steroids and WH50 (in the case of 5AD) or WH100 (in the case of PRG). Mutation WH47A decreases the binding energy of 5AD 4.6 kJ/mol more than that of PRG. WH47 is <4.5 Å away from the C ring of both steroids. However, due to the bent geometry of 5AD, its A ring is <5 Å from WH47, which does not interact with the A ring of PRG. In addition to this, the A ring of 5AD makes several additional interactions with the antibody. The A ring of 5AD is bound nicely to the P3' binding pocket defined by residues WH47, WH50, VL94, PL95 and PL96. There is a tight contact between 5AD and VL94 with distances of 3.7 Å between the methyl carbons of VL94 and the A ring of 5AD. In contrast, only one of the methyl groups of VL94 interacts with PRG. These differences are seen in the calculated binding free energies. Mutation VL94A decreases the binding of 5AD 5.9 kJ/mol more than that of PRG. The carbonyl group of PRG's A ring makes a hydrogen bond with N
H of HL27d (part of the P3 compartment). This interaction is absent for 5AD, whose carbonyl group of the A ring is solvated by water molecules. Therefore, mutation HL27dA decreases the binding of PRG by 8.4 kJ/mol but that of 5AD only by 3.8 kJ/mol. In addition to the hydrogen bond with HL27d, the A and B rings of PRG make several hydrogen-bonding type interactions with the main-chain carbonyl carbons of SL91 and HL93. The heavy atom distances between the carbonyl carbons and C4, C6 and C7 of PRG are 34 Å. In contrast, only the C19 of 5AD has a short distance (3.9 Å) to SL91.
The P3 and P3' compartments are responsible for the observed steroid cross-reactivity to the DB3 antibody. The comparison of the antibodysteroid structures shows that the compartments have different types of interactions with the bound steroids: the P3 compartment operates mainly by hydrogen-bonding and electrostatic interactions, whereas the binding to the P3' compartment is mainly hydrophobic in nature. These differences are seen in the free energy components. Although the calculated electrostatic solvation free energies of PRG and 5AD are similar (51.9 and 54.0 kJ/mol, respectively), the solvation free energy for binding is 47.7 kJ/mol less unfavourable for 5AD. This reflects the more hydrophobic nature of the P3' compartment. That the P3 compartment has hydrogen-bonding like interactions with PRG, is seen in the larger electrostatic proteinligand interaction energies (Eelec) for PRG (73.4 kJ/mol) than 5AD (41.9 kJ/mol).
To conclude, the MM-PBSA calculations of the present study reproduced the relative binding free energy of PRG and 5AD, which bind to slightly different binding pockets of anti-PRG antibody DB3, in fair agreement with experimental values (8.8 versus 5.5 kJ/mol). Analyses of the binding energetics showed that in agreement with earlier computational ligand-binding studies the van der Waals energy is the most important complex formation favouring energy term. In addition, electrostatic antibodysteroid interactions were found to be important for the binding of PRG, whereas van der Waals and hydrophobic interactions were important for 5AD. These differences are mostly due to the binding of the A rings of the steroids to different binding pockets of DB3.
![]() |
Notes |
---|
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Arevalo,J.H., Hassig,C.A., Stura,E.A., Sims,M.J., Taussig,M.J. and Wilson,I.A. (1994) J. Mol. Biol., 241, 663690.[ISI][Medline]
Arevalo,J.H., Stura,E.A., Taussig,M.J. and Wilson,I.A. (1993a) J. Mol. Biol., 231, 103118.[ISI][Medline]
Arevalo,J.H., Taussig,M.J. and Wilson,I.A. (1993b) Nature, 365, 859863.[ISI][Medline]
Arnold,F.H. and Volkov,A.A. (1999) Curr. Opin. Chem. Biol., 3, 5459.[ISI][Medline]
Berendsen,H.J.C., Postma,J.P.M., van Gunsteren,W.F., DiNola,A.D. and Haak,J.R. (1984) J. Chem. Phys., 81, 36843690.[ISI]
Case,D.A. et al. (1997) AMBER 5.0. University of California, San Francisco.
Charbonnier,J.-B., Golinelli-Pimpaneau,B., Gigant,B., Tawfik,D.S., Chap,R., Schindler,D.G., Kim,S.-H., Green,B.S., Eshhar,Z. and Knossow,M. (1997) Science, 275, 11401142.
Chong,L.T., Duan,Y., Wang,L., Massova,I. and Kollman,P.A. (1999) Proc. Natl Acad. Sci. USA, 96, 1433014335.
Cornell,W.D., Cieplak,P., Bayly,C.I. and Kollman,P.A. (1993) J. Am. Chem. Soc., 115, 96209631.[ISI]
Cornell,W.D., Cieplak,P., Bayly,C.I., Gould,I.R., Merz,K.M.,Jr, Ferguson,D.M., Spellmeyer,D.C., Fox,T., Caldwell,J.W. and Kollman,P.A. (1995) J. Am. Chem. Soc., 117, 51795197.[ISI]
Fox,T. and Kollman,P.A. (1998) J. Phys. Chem. B, 102, 80708079.[ISI]
Gigant,B., Charbonnier,J.-B., Eshhar,Z., Green,B.S. and Knossow,M. (1997) Proc. Natl Acad. Sci. USA, 94, 78577861.
Hellinge,H.W. (1998) Nat. Struct. Biol., 5, 525527.[ISI][Medline]
Hill,R.B., Raleigh,D.P. and Lombardi,W.F. (2000) Acc. Chem. Res., 33, 745754.[ISI][Medline]
Jayaram,B., Sprous,D., Young,M.A. and Beveridge,D.L. (1998) J. Am. Chem. Soc., 120, 1062910633.[ISI]
Jeffrey,P.D., Strong,R.K., Sieker,L.C., Chang,C.Y.Y., Campbell,R.L., Petsko,G.A., Haber,E., Margolies,M.N. and Sheriff,S. (1993) Proc. Natl Acad. Sci. USA, 90, 1031010314.[Abstract]
Jorgensen,W.L., Chandrasekhar,J., Madura,J., Impley,R.W. and Klein,M.L. (1983) J. Chem. Phys., 79, 926935.[ISI]
Kollman,P.A. et al. (2000) Acc. Chem. Res., 33, 889897.[ISI][Medline]
Kuhn,B. and Kollman,P.A. (2000) J. Med. Chem., 43, 37863791.[ISI][Medline]
Lee,M.R., Duan,Y. and Kollman,P.A. (2000) Proteins, 39, 309316.[ISI][Medline]
Liu,D.R. and Schultz,P.G. (1999) Angew. Chem., Int. Ed. Engl., 38, 3654.[ISI]
Mayo,S.L. and Gordon,D.B. (2000) J. Mol. Biol., 299, 789803.[ISI][Medline]
Miyamoto,S. and Kollman,P.A. (1993) Proc. Natl Acad. Sci. USA, 90, 84028406.
Nicholls,A., Sharp,K.A. and Honig,B. (1990) Delphi. Department of Biochemistry and Molecular Biophysics, Columbia University, New York.
Reyes,C.M. and Kollman,P.A. (2000a) J. Mol. Biol., 295, 16.[ISI][Medline]
Reyes,C.M. and Kollman,P.A. (2000b) J. Mol. Biol., 297, 11451158.[ISI][Medline]
Ryckaert,J.P., Ciccotti,G. and Berendsen,H.J.C. (1977) J. Comput. Phys., 23, 327341.[ISI]
Sanner,M.F., Olson,A.J. and Spehner,J.C. (1996) Biopolymers, 38, 305320.[ISI][Medline]
Schafmeister,C.E.A.F., Ross,W.S. and Romanovski,V. (1995) Leap. University of California, San Francisco.
Schultz,P.G. and Lerner,R.A. (1995) Science, 269, 18351842.[ISI][Medline]
Sinha,S.C., Barbas,C.F.I. and Lerner,R.A. (1998) Proc. Natl Acad. Sci. USA, 95, 1460314608.
Sitkoff,D., Sharp,K.A. and Honig,B. (1994) J. Phys. Chem., 98, 19781988.[ISI]
Srinivasan,J., Cheatham,T.E., Cieplak,P., Kollman,P.A. and Case,D. (1998) J. Am. Chem. Soc., 120, 94019409.[ISI]
Trinh,C.H., Hemmington,S.D., Verhoeyen,M.E. and Phillips,S.E.V. (1997) Science, 7, 937948.
Vorobjev,Y.N., Almagro,J.C. and Hermans,J. (1998) Proteins, 32, 399413[ISI][Medline]
Wedemayer,G.J., Patten,P.A., Wang,L.H., Schultz,P.G. and Stevens,R.C. (1997) Science, 276, 16651669.
Wilson,I.A. and Stanfield,R.L. (1994) Curr. Opin. Struct. Biol., 4, 857867.[ISI][Medline]
Received March 13, 2001; revised July 24, 2001; accepted July 31, 2001.