Skip to main content

Synergy between trastuzumab and pertuzumab for human epidermal growth factor 2 (Her2) from colocalization: an in silicobased mechanism



Human epidermal growth factor 2 (Her2), a receptor tyrosine kinase, is overexpressed in breast cancers. It has been successfully targeted by small molecule kinase inhibitors and by antibodies. Recent clinical data show a synergistic response in patients when two antibodies, trastuzumab and pertuzumab, are given in combination.


This unexpected effect is rationalized through computer models and molecular dynamic simulations by hypothesizing that the two antibodies can co-localize on the same molecule of the Her2 extracellular domain.


Simulations suggest that the clinical synergism observed for the two antibodies arises partly from enhanced affinity that originates in cooperative interactions between these two antibodies when they are co-localized on Her2 and "clamp" it; this may inhibit dimerization and possibly higher oligomerizations with neighboring receptors. In the presence of trastuzumab, the receptor becomes highly plastic, especially domains I to III, and this appears to promote increased association with pertuzumab. Further, the presence of pertuzumab evokes novel interactions between the receptor and trastuzumab. Indeed, splicing out of this region in silico results in a big reduction in the interactions of the antibody with the receptor.


If validated, these findings will bring about a new direction in the design of antibodies whereby different epitopes on the same antibody may be targeted to lead to synergistic/cooperative inhibition and contribute to generate more potent therapeutics and to increase clinical efficacy.


ErbB receptors are the prototypical founders of the growth factor receptor tyrosine kinase (RTK) family. They are activated by the binding of different ligands and are involved in the transmission of signals from the extracellular space to the cytoplasm and nucleus of a cell, thus orchestrating biological processes [1]. Among the members of this family, which also includes the epidermal growth factor receptor (EGFR or ErbB1), ErbB3, ErbB4, Her2 or ErbB2 is homologous to, but distinct from the others, since it is not activated upon ligand binding. This ligand-independent activation makes Her2 the universal heterodimeric partner for each of the other ErbB family members [2]. The architecture of these receptors [3, 4] reveals an extracellular domain (made up of four subdomains), a single transmembrane helix, and an intracellular domain (consisting of a juxtamembrane region, a tyrosine kinase domain, and a C-terminal tail harboring autophosphorylation sites that serve as docking sites for adaptor proteins [5, 6]). The importance of the tight regulation of these receptors and in particular Her2 is signified in human breast cancers, where Her2 is overexpressed by approximately 20 to 30% and this is normally associated with more aggressive tumours and a poorer prognosis [7, 8]. Due to this active role in human cancers, a number of therapeutic approaches are currently under development to block the effects of Her2 overexpression, including kinase inhibitors (such as lapatinib) and monoclonal antibodies (trastuzumab and pertuzumab) [9]. The anti-Her2 monoclonal antibody trastuzumab binds to domain IV of Her2, a region that is not involved in receptor dimerization and is thought to work in a complex manner [10, 11]. In contrast, pertuzumab binds to Her2 near the centre of the domain II dimerization arm (Figure 1). The overlap between the pertuzumab epitope and the probable heterodimer interface suggests steric occlusion by physically blocking the formation of Her2-containing heterodimers. Recent clinical observations have demonstrated that combining trastuzumab and pertuzumab together yields surprisingly synergistic results in tumour inhibition (that is, the combined effect of using the two molecules was far greater in effect than when using either alone) during the treatment of Her2 positive breast cancers that progressed during prior trastuzumab therapy [12]. This suggests a cooperative mechanism of inhibition that can lead to clinical improvement in the treatment of these tumours.

Figure 1

Representation of the time evolution of the Her2 and Abs complex and their interactions. Representation of the in silico pertuzumab (in orange) -induced trastuzumab (in black) epitope that emerges during the course of the MD simulations of Her2 ectodomain (in red-blue electrostatic surface where red represents regions of acidity and blue represents regions of basicity) in presence of trastuzumab and/or pertuzumab. (A) Her2:trastuzumab complex structure at the beginning of the simulation(time = 0 ns) showing no contact between trastuzumab and the new epitope and (B) snapshot taken during the MD simulation (time = 15 ns) showing interactions between trastuzumab and the new epitope; (C) dynamic Cα-Cα contact matrix (representing the distance between each pair of residues comprising the three dimensional structure of protein complex), the data from both pertuzumab and trastuzumab co-localized on Her2 are shown below the diagonal; the data from the Her2-trastuzumab complex are shown above the diagonal. Trastuzumab has been coloured in black; while pertuzumab is shown in orange. This plot shows the occurrence of interactions between the heavy chain of trastuzumab (see vertical axis) and the region of the ectodomain around 310 to 330 (see the horizontal axis). During the course of simulation, when pertuzumab is bound to Her2:trastuzumab these contacts are formed, in contrast to the simulation where only Her2:trastuzumab are present (comparison with the same region above the diagonal) where these contacts are absent.

We have used molecular modelling to develop a mechanism that may partly explain this cooperative effect. A simple model is that cooperativity could arise from the physical co-localization of the two antibodies on the same Her2 molecule, as has been previously demonstrated for the insulin-like growth factor receptor [13]. Models of the binary structures of Her2:trastuzumab and Her2:pertuzumab were combined to create a model of the ternary complex of Her2 and the two Fab fragments of these antibodies. This system was subjected to structural/energetic refinement followed by molecular dynamics (MD) simulations (Additional file 1).


Building the models

The binary associations of Her2 and the two antibodies have been experimentally characterized and are available in the Protein Data Base (pdb entries: 1n8z[3] and 1s78 [14]). The electronic density of the C-terminal region of the receptor is the Her2:pertuzumab complex is missing, due probably to the high inherent flexibility of this part of the structure. For the coordinates of the residues in this region, we have assumed that this stretch will adopt a conformation similar to the one found in the complex with trastuzumab and have been so modeled. A few other missing residues, located in loops, have been modeled based on homologue structures. The ternary model of the complex Her2 and the two Fab domains of these antibodies, trastuzumab and pertuzumab, were manually constructed using the coordinates of the previous models, after superposition on the backbone of the receptor ectodomain.

MD simulations and binding energy calculations

We have used seven different systems for this study (Additional file 1); the complete 3D model structure of the complex (Her2-P-T), the receptor bound to trastuzumab (Her2-T) and pertuzumab (Her2-P), the apo-receptor (Her2), the two independent antibodies (pert and trast) and both antibodies co-localized as they are found in the complex (Abs). All these systems have been refined using several energy minimization and short MD simulation cycles in order to equilibrate the structures after which they were subjected to the production phases of MD simulation (around 18 to 20 ns each) in explicit water, as detailed. The AMBER (Assisted Model Building with Energy Refinement) 9 package [15] together with the ff99 force field [16] was used. The protonation states of the residues were assigned using PDB2PQR [17], after which the systems were surrounded with a box of TIP3 or Transferrable Interaction Potential 3 water molecules that result in an 8 Å layer of water molecules around the proteins and neutralized by the necessary numbers of Cl-or Na+ counter-ions using the t-leap module of AMBER 9. The solvated systems were then subjected to a minimization phase, mainly by the steepest descend method and small number of conjugate gradient steps, in order to remove bad contacts. After that a 50 ps heating phase was carried out, where the proteins were heated from 0 to 300 K, followed by another 50 ps phase density equilibration at a constant pressure of 1 atm. A final 2 ns unconstrained MD simulation was carried out. The production runs were extended to around 20 ns at constant temperature using a Langevin thermostat [18] and constant pressure (1 atm) by coupling to a Berendsen barostat These two queries do not need to be address since these are mathematical methods and not devices [19]. The particle mesh Ewald (PME) method [20] was used with a direct-space non-bonded cutoff of 8 Å to account for long range interactions. Covalent bond lengths to hydrogen atoms were constrained using the SHAKE algorithm [21], which allowed a time step of 2 fs for all the MD simulations. Snapshots in the interval of 8 to 16 ns were extracted for further processing. The analysis was carried out using functionalities in the ptraj module of AMBER 9, the module Bio3D [22] of R and Dyndom [23]. The PyMOL molecular viewer, a software package [24] was used for visualizing and generating the figures.

The binding free energy ΔGbind of the different systems was calculated according to the so-called Molecular Mechanics/Generalized-Born Surface Area (MM/GBSA) procedure [25]. Here ΔGbind denotes the contributions from the molecular mechanics (MM) energies of the molecules, and originates in contributions from internal energies, electrostatics (ΔGELE) and van der Waals (ΔGVDW) contributions. These terms have been calculated using the mm_pbsa module of AMBER 9. The polar component was evaluated with the Generalized-Born (GB) approach. In this approach, the non-polar contribution to the solvation energy was calculated as ΔGNP = γ (SASA) + β, in which γ = 0.00542 kcal/Å, β = 0.92 kcal/mol and SASA is the solvent-accessible surface area estimated with LCPO (linear combination of pairwise overlaps) method [26] as implemented in AMBER. While the real free energy of a system also includes contributions from the entropies, given the large sizes of the systems studied here, computations of entropies were prohibitive and hence could not be carried out. However, the free energies computed here have been averaged over several conformations generated through MD simulations, and hence the effects of entropy to some extent are implicitly included in the current calculations.

Essential dynamics analysis and Dyndom analysis of Her2 and its antibodies

Essential dynamics [27, 28] have been widely used for filtering large scale concerted fluctuations from an ensemble of conformations, like those obtained from MD simulations [29]. The method is comparable to a principal component analysis, and it is based on the construction of the diagonalized covariance matrix containing the atomic fluctuations. We have performed a standard essential dynamics (or Principal Component Analysis, PCA) analysis in order to identify the most relevant motions occurring in the receptor in the different bound states considered on the MD simulations carried out in this study (Additional file 2); all the covariance matrices have been calculated taking the apo receptor as a reference structure. This yields a set of eigenvectors in a 3N-dimensional space, and eigenvalues, reflecting the magnitude of the fluctuations along the corresponding eigenvectors. Each eigenvector describing a collective motion can be represented by two structures, one at the maximum value, the other at the minimum value of the collective variable when the trajectory is projected onto it. These two structures thus generated for the first eigenvector of each trajectory were used as input for the program DynDom [23], using default parameters, in order to determine dynamic domains, hinge axes and hinge-bending regions.

Results and discussion

Conformational space of Her2 in different bound forms

Simulations show that the ectodomain of Her2 adopts various conformations, depending on its bound state. The binary states (when complexed to trastuzumab or pertuzumab) are more stable compared to the apo or the ternary complex (Additional file 3). This could stem from the high intrinsic flexibility of the apo receptor or/and the ternary complex. Although a complete understanding of the behaviour of such a receptor demands the full-length structure, we believe that the models used here represent valid starting structures for the type of analyses carried out in this study. As expected, the presence of an antibody immobilizes its corresponding epitope (Additional file 4). The largest fluctuations in the pertuzumab binding region (ranging from 240 to 300) and in the trastuzumab epitope (around 550 to 600) are seen in the apo state or in the Her2:trastuzumab or Her2:pertuzumab complexes respectively. Interestingly the dimerization arm, a β hairpin in domain II that protrudes from the rest of Her2, together with the spatially adjacent regions, undergoes enhanced fluctuations when trastuzumab binds. This may reinforce the fact that the Her2:trastuzumab complex can actually promote the association of Her2 with other members of the ErbB family.

We have also investigated the flexibility of the antibodies and the stability of their binding interfaces with Her2. In general, the complementary determining regions (CDR) exhibit smaller fluctuations in comparison to the rest of the antibody; whereas regions far from the receptor's influence remain extremely flexible. In the ternary complex, the C-terminal half of the light chain of pertuzumab shows an increased flexibility, induced by the presence of trastuzumab (Additional file 4). Trastuzumab shows a more complicated pattern. When only trastuzumab is present in the complex, the whole antibody shows a high flexibility, with the exception of the CDR for the light chain. In the ternary complex, the light chain undergoes the largest fluctuations in the region interacting with the receptor binding epitope; however the most relevant difference is an attenuation in the mobility of the N-terminal half of the heavy chain due to the proximity of one of the disulphide-bonded modules of domain II, which in turn shows lower flexibility (residues 310 to 325; see section below on new epitopes). We can conclude that the dynamics of pertuzumab seem to be regulated mainly by the binding to the receptor; while for trastuzumab, the co-binding of the second antibody has a significant effect. When the two antibodies are analysed in the absence of the receptor, trastuzumab per se seems more flexible than pertuzumab, which can be the origin of the more complex dynamics observed for the systems including this antibody. These intricate changes and regulation between the different components of the system reflect the importance of flexibility and could point to "allosteric" changes in explaining the in vivo synergistic effect.

Uncovering new epitopes in Her2

The analysis of the interactions between Her2 and the two antibodies reveals that the perturbation of the crystallographically observed interfacial region of Her2 is small during the simulations of the binary complexes (Additional file 5). However, the co-localization of the two Abs on Her2 surprisingly triggers the formation of additional and novel contacts of Her2 with trastuzumab; these lie outside of its epitopes (Figure 1). This β-hairpin loop region, termed the pertuzumab-induced trastuzumab epitope, is located in domain II and is in proximity to the dimerization arm. The requirement of the co-localization of the two Abs for the formation of this new epitope is revealed in the minimum distance matrices for all the Cα atoms for the receptor and trastuzumab in the simulations of the binary and ternary complexes (see Figure 1C upper and lower diagonal of the matrix, respectively). This can also be observed from the reported distances between two residues in this novel epitope and residues at the interface in trastuzumab (Additional file 6). For this epitope to emerge in the Her2:trastuzumab interface, the presence of pertuzumab is strictly required. This novel association is mainly driven by long-range interactions, and could be the origin of the observed synergism.

Quantification of synergism between pertuzumab and trastuzumab

The MM-GBSA (molecular mechanics Generalized Born surface area) [25] is regularly employed to obtain estimates of the free binding energy of large systems and these are largely reasonable [3032]. These only refer to the enthalpic component of the overall free energies of binding because computations of entropies of such large systems are prohibitively expensive. The derived association energetics of trastuzumab and pertuzumab to Her2 are shown in Table 1; the computed values are very large as is the case with computations involving large systems; indeed in comparison the binding of trastuzumab to the extracellular domain of Her2 and to Her2 positive cells ranges from approximately -12.5 to -14 kcal/mol [33]. However, despite being large, qualitative analyses of the differences in the computed binding energies are very useful in ranking interactions. Our findings suggest (the more negative an energy, the more stable the associated system):

Table 1 In silico binding free energy for the association of Her2 ectodomain and the two antibodies.

- trastuzumab has a higher affinity for apo Her2 than pertuzumab (Figure 2A); this arises from stronger electrostatic and van der Waals interactions between trastuzumab and Her2 (the net charge on Her2 is -15, while on trastuzumab and pertuzumab it is +6 and +5 respectively); the buried area at the interface is approximately 1,675 and approximately 1,530 Å2 for trastuzumab and pertuzumab, respectively

Figure 2

Thermodynamic cycles to study the cooperative binding of Her2 to trastuzumab and pertuzumab. (A) The cycle with the energy estimations calculated from the independent association simulations. (B) The cycle analysed using the conformations of receptor and binary complexes from the co-localized interaction simulations. The energetics for the binary complexes have been obtained from the ternary complex of Her2:pertuzumab:trastuzumab complex, in an attempt to account for the influence of the other antibody when both are co-localized. (C) The cycle representing the gain in energy between the binding of the two antibodies to each other in an isolated system or when co-localized onto the structure of the receptor. The ectodomain of Her2 has been coloured in green, while trastuzumab and pertuzumab are represented in black and orange, respectively. The differences in free energy binding estimated for every system are shown. The computations were for the following schemes: A) her2-trastuzumab < -her2+ trastuzumab; her2-trast-pert < -her2-trast + pertuzumab her2-trast + pertuzumab < -her2+ trastuzumab + pertuzumab her2-pertuzumab < -her2+ pertuzumab; her2-trast-pert < -her2-pert + trastuzumab. B) her2-trastuzumab (extracted from her2:pert:trast) < -her2+ trastuzumab; her2-trast-pert < -her2-trast (extracted from her2-pert:trast) + pertuzumab her2-pertuzumab (extracted from her2:pert:trast) < -her2+ pertuzumab; her2-trast-pert < -her2-pert (extracted from her2-pert:trast) + trastuzumab. C) Abs (extracted from her2:pert:trast) < -pertuzumab + trastuzumab; Abs < -pertuzumab + trastuzumab.

- a marked gain in affinity when both Abs are bound simultaneously to the receptor (synergism). This effect arises mainly from the interactions between these two antibodies upon co-localization, with a difference of approximately 27 kcal/mol between isolated and simultaneously bound antibodies (Figure 2C).

The computed energetics within these thermodynamic cycles suggest a sequential mechanism, where the initial binding involves trastuzumab (Figure 2A) followed by the binding of pertuzumab (Figure 2B). This seems to be partly driven by the conformational changes triggered upon the binding of trastuzumab that enhance pertuzumab binding (Figure 3) and the new additional contacts formed between trastuzumab and Her2. The enhanced interactions that accrue when the two antibodies are co-localized on Her2 leading to the cooperative effects, may provide a partial explanation for the clinically observed synergy. In order to support the effect of this newly revealed pertuzumab-induced trastuzumab epitope, we have constructed an in silico truncated form of Her2, where this β-hairpin has been deleted (T323C-AEDG-T328C), and we find that the binding energetics undergo a marked destabilization (Table 1).

Figure 3

Structural changes on Her2 ectodomain upon binding of trastuzumab and/or pertuzumab. To understand the structural changes that occur as a result of the binding of the two antibodies, we analyse the conformational transition/populations using cross 2D-RMSD plots (a plot of RMSD vs RMSD) of the ectodomain of Her2 in the trajectory of the ternary complex of Her2:pertuzumab:trastuzumab and the different apo and binary trajectories. The latter were fitted onto the full receptor structure, except the pertuzumab epitope region, and the rmsd has been calculated for the Cα atoms using the initial structure of the ternary complex simulation as reference. These maps reveal the differences of the various conformations sampled during the simulation with respect to the ternary complex. Different populations of structures are observed in the plots, with structures closest to the conformation adopted for the receptor at the beginning of the ternary simulations shown in blue, to those that are most altered in red (the largest variations are up to approximately 10 Å RMSD). The binding of trastuzumab displaces the conformation of the receptor towards that found in the ternary complex (see plot at the left and density distribution for such a complex). These conformational changes enhance the binding of pertuzumab and thus the affinity of the two antibodies for the receptor is increased, as shown in Table 1, leading to cooperativity.

Abs-induced bending of the Her2 ectodomain

Principal Component Analysis (Additional file 2) combined with dynamical domain motion analysis using DynDom suggests that domains I to III (residues 1 to 480) generally behave as a rigid body, with hinge regions varying in their location in domain IV, depending on the binding partner of the receptor (Additional files 7 and 8). Two different motions are detected in the apo state of Her2; the first one involves a shearing motion between the stalk portion of the receptor and the large rigid-body unit; whereas the second one corresponds to a hinge motion within domain IV. The overall change in conformation involves a large conformational space scanned by the dimerization arm with respect to the membrane surface (with an angle of approximately 50° for the parallel and approximately 78° for the motion perpendicular to the membrane), which represents periscope-like movements that this receptor might execute in search of binding partners. In the case of the ternary complex, the binding of trastuzumab causes a shift in the hinge region along domain IV. This is associated with small swivelling motions across an axis orthogonal to the membrane; if the receptor lies parallel to the membrane [34] then this would represent "out of the membrane" motions. The dimerization arm spans approximately 16 Å in comparison to the apo form which covers approximately 33 Å (measured from the Cα of the atom (Pro252) at the base of this motif). In contrast, in the binary states, the pattern of motion differs: domain IV behaves as a rigid body and only one hinge region is apparent. When trastuzumab is bound, a hinge motion is observed (18°) between the bulky domain I-to-III and domain IV (crystallographic data show approximately 7° in the angle subtended between domains III and IV); the Her2:pertuzumab complex undergoes more complex motions. The deletion in the truncated mutant results in attenuation of two-body rotations along the parallel axis (from 48° in the "full-length" receptor to approximately 28° in the truncated form). In contrast, the bound form of this truncated receptor undergoes a hinged motion spanning approximately 20°.

The overall comparison of the conformations that Her2 ectodomain adopts in the different conditions analysed here indicates that domains I, II and III adopt a relatively fixed interdomain relationship, but their relative orientations with respect to domain IV varies considerably (up to approximately 80°), and appear to be pivoted across the region 504 to 532 that lies between domains III and IV. This is indicative of the plasticity inherent in these receptors and suggests that the extracellular domain can swivel quite considerably across domain IV that itself is constrained by being embedded partly in the membrane.


In summary, MD simulations suggest that the clinical synergism observed for the two antibodies arises partly from an enhanced affinity that originates in cooperative interactions between these two antibodies when they are co-localized on Her2 and "clamp" it; this may inhibit dimerization and possibly higher oligomerizations with neighbouring receptors. In the presence of trastuzumab, the receptor becomes highly plastic, especially domains I-III, and this appears to promote increased association with pertuzumab. While combining antibodies is a feature that increasingly appears to have major effects in therapy (see, for example, the work of Chao and co-workers [35], which combines two different antibodies against two different cell surface molecules for treatment of non-Hodgkin's lymphoma), a model incorporating the idea that two antibodies bind to the same receptor, albeit at distinct epitopes, is relatively new; indeed there has been only one report of such a feature as highlighted for targeting the IGF-1 receptor [13]. These findings could be tested, for example, in 2D and 3D preclinical cell culture models to explore trastuzumab efficacy [36, 37]. While we await experimental investigations of this novel mechanism in Her2, if validated, it offers routes towards the design of new antibodies whereby combinations of trastuzumab and pertuzumab (or parts thereof) may be constructed that have the potential to generate therapeutics with much higher affinities and hopefully efficacies. The technology of generating hybrid antibodies or bi-specific antibodies [38] is known; however, this technique has been used for linking two antibodies that target different molecules. In contrast our hypotheses open the avenue for linking two antibodies that target distinct epitopes on the same molecule.



binding free energy


electrostatics energy


non-polar energy


van der Waals energy




complementary determining regions


epidermal growth factor receptor


fragment antigen binding




human epidermal growth factor 2


insulin-like growth factor


linear combination of pairwise overlaps


molecular dynamics


molecular mechanics/Generalized-Born surface Area


Principal component analysis


particle mesh Ewald


receptor tyrosine kinase


solvent-accessible surface area


  1. 1.

    Yarden Y, Sliwkowski MX: Untangling the ErbB signalling network. Nat Rev Mol Cell Biol. 2001, 2: 127-137. 10.1038/35052073.

    CAS  Article  Google Scholar 

  2. 2.

    Holbro T, Hynes NE: ErbB receptors: directing key signaling networks throughout life. Annu Rev Pharmacol Toxicol. 2004, 44: 195-217. 10.1146/annurev.pharmtox.44.101802.121440.

    CAS  Article  Google Scholar 

  3. 3.

    Cho HS, Mason K, Ramyar KX, Stanley AM, Gabelli SB, Denney DW, Leahy DJ: Structure of the extracellular region of Her2 alone and in complex with the Herceptin Fab. Nature. 2003, 421: 756-60. 10.1038/nature01392.

    CAS  Article  Google Scholar 

  4. 4.

    Garrett TP, McKern NM, Lou M, Elleman TC, Adams TE, Lovrecz GO, Kofler M, Jorissen RN, Nice EC, Burgess AW, Ward CW: The crystal structure of a truncated ErbB2 ectodomain reveals an active conformation, poised to interact with other ErbB receptors. Mol Cell. 2003, 11: 495-505. 10.1016/S1097-2765(03)00048-0.

    CAS  Article  Google Scholar 

  5. 5.

    Jorissen RN, Walker F, Pouliot N, Garrett TP, Ward CW, Burgess AW: Epidermal growth factor receptor: mechanisms of activation and signalling. Exp Cell Res. 2003, 284: 31-53. 10.1016/S0014-4827(02)00098-8.

    CAS  Article  Google Scholar 

  6. 6.

    Schlessinger J: Cell signaling by receptor tyrosine kinases. Cell. 2000, 103: 211-225. 10.1016/S0092-8674(00)00114-8.

    CAS  Article  Google Scholar 

  7. 7.

    Slamon DJ, Godolphin W, Jones LA, Holt JA, Wong SG, Keith DE, Levin WJ, Stuart SG, Udove J, Ullrich A: Studies of the HER-2/neu proto-oncogene in human breast and ovarian cancer. Science. 1989, 244: 707-712. 10.1126/science.2470152.

    CAS  Article  Google Scholar 

  8. 8.

    Slamon DJ: Proto-oncogenes and human cancers. N Engl J Med. 1987, 317: 955-957. 10.1056/NEJM198710083171509.

    CAS  Article  Google Scholar 

  9. 9.

    Baselga J, Swain SM: Novel anticancer targets: revisiting ERBB2 and discovering ERBB3. Nat Rev Cancer. 2009, 9: 463-475. 10.1038/nrc2656.

    CAS  Article  Google Scholar 

  10. 10.

    Harari D, Yarden Y: Molecular mechanisms underlying ErbB2/Her2 action in breast cancer. Oncogene. 2000, 19: 6102-6114. 10.1038/sj.onc.1203973.

    CAS  Article  Google Scholar 

  11. 11.

    Sliwkowski MX, Lofgren JA, Lewis GD, Hotaling TE, Fendly BM, Fox JA: Nonclinical studies addressing the mechanism of action of trastuzumab (Herceptin). Semin Oncol. 1999, 26: 60-70.

    CAS  PubMed  Google Scholar 

  12. 12.

    Baselga J, Gelmon KA, Verma S, Wardley A, Conte P, Miles D, Bianchi G, Cortes J, McNally VA, Ross GA, Fumoleau P, Gianni L: Phase II trial of pertuzumab and trastuzumab in patients with human epidermal growth factor receptor 2-positive metastatic breast cancer that progressed during prior trastuzumab therapy. J Clin Oncol. 2010, 28: 1138-1144. 10.1200/JCO.2009.24.2024.

    CAS  Article  Google Scholar 

  13. 13.

    Dong J, Demarest SJ, Sereno A, Tamraz S, Langley E, Doern A, Snipas T, Perron K, Joseph I, Glaser SM, Ho SN, Reff ME, Hariharan K: Combination of two insulin-like growth factor-I receptor inhibitory antibodies targeting distinct epitopes leads to an enhanced antitumor response. Mol Cancer Ther. 2010, 9: 2593-2604. 10.1158/1535-7163.MCT-09-1018.

    CAS  Article  Google Scholar 

  14. 14.

    Franklin MC, Carey KD, Vajdos FF, Leahy DJ, de Vos AM, Sliwkowski MX: Insights into ErbB signaling from the structure of the ErbB2-pertuzumab complex. Cancer Cell. 2004, 5: 317-328. 10.1016/S1535-6108(04)00083-2.

    CAS  Article  Google Scholar 

  15. 15.

    Case DA, Darden TA, Cheatham TE, Simmerling CL, Wang J, Duke RE, Luo R, Merz KM, Pearlman DA, Crowley M, Walker RC, Zhang W, Wang B, Havik S, Riotberg A, Seabra G, Wong KF, Paesani F, Wu X, Brozell S, Tsui V, Gohlke H, Yang L, Tan C, Mongan J, Hornak V, Cui G, Beroza P, Mathews DH, Schafmeister C, et al: AMBER 9. 2006, University of California, San Francisco

    Google Scholar 

  16. 16.

    Wang J, Cieplak P, Kollman PA: How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules?. J Comput Chem. 2000, 21: 1049-1074. 10.1002/1096-987X(200009)21:12<1049::AID-JCC3>3.0.CO;2-F.

    CAS  Article  Google Scholar 

  17. 17.

    Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA: PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 2004, 32: W665-10.1093/nar/gkh381.

    CAS  Article  Google Scholar 

  18. 18.

    Izaguirre JA, Catarello DP, Wozniak JM, Skeel RD: Langevin stabilization of molecular dynamics. J Chem Phys. 2001, 114: 2090-10.1063/1.1332996.

    CAS  Article  Google Scholar 

  19. 19.

    Berendsen HJC, Postma JPM, Van Gunsteren WF, DiNola A, Haak JR: Molecular dynamics with coupling to an external bath. J Chem Phys. 1984, 81: 3684-10.1063/1.448118.

    CAS  Article  Google Scholar 

  20. 20.

    Darden T, York D, Pedersen L: Particle mesh Ewald: An N log(N) method for Ewald sums in large systems. J Chem Phys. 1993, 98: 10089-10092. 10.1063/1.464397.

    CAS  Article  Google Scholar 

  21. 21.

    Ryckaert JP, Ciccotti G, Berendsen HJC: Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys. 1977, 22: 327-341.

    Article  Google Scholar 

  22. 22.

    Grant BJ, Rodrigues APC, ElSawy KM, McCammon JA, Caves LSD: Bio 3 d: an R package for the comparative analysis of protein structures. Bioinformatics. 2006, 22: 2695-10.1093/bioinformatics/btl461.

    CAS  Article  Google Scholar 

  23. 23.

    Hayward S, Berendsen HJC: Systematic analysis of domain motions in proteins from conformational change: new results on citrate synthase and T4 lysozyme. Proteins. 1998, 30: 144-154. 10.1002/(SICI)1097-0134(19980201)30:2<144::AID-PROT4>3.0.CO;2-N.

    CAS  Article  Google Scholar 

  24. 24.

    PyMOL. DeLano Scientific, Palo Alto, CA, []

  25. 25.

    Kollman PA, Massova I, Reyes C, Kuhn B, Huo S, Chong L, Lee M, Lee T, Duan Y, Wang W, Donini O, Cieplak P, Srinivasan J, Case DA, Cheatham TE: Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Acc Chem Res. 2000, 33: 889-897. 10.1021/ar000033j.

    CAS  Article  Google Scholar 

  26. 26.

    Weiser J, Shenkin PS, Still WC: Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO). J Comput Chem. 1999, 20: 217-230. 10.1002/(SICI)1096-987X(19990130)20:2<217::AID-JCC4>3.0.CO;2-A.

    CAS  Article  Google Scholar 

  27. 27.

    Amadei A, Linssen AB, Berendsen HJ: Essential dynamics of proteins. Proteins. 1993, 17: 412-425. 10.1002/prot.340170408.

    CAS  Article  Google Scholar 

  28. 28.

    Levy RM, Srinivasan AR, Olson WK, McCammon JA: Quasi-harmonic method for studying very low frequency modes in proteins. Biopolymers. 1984, 23: 1099-1112. 10.1002/bip.360230610.

    CAS  Article  Google Scholar 

  29. 29.

    Stepanova M: Dynamics of essential collective motions in proteins: theory. Phys Rev E. 2007, 76: 51918-

    Article  Google Scholar 

  30. 30.

    Gohlke H, Kiel C, Case DA: Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. J Mol Biol. 2003, 330: 891-913. 10.1016/S0022-2836(03)00610-7.

    CAS  Article  Google Scholar 

  31. 31.

    Balius TE, Rizzo RC: Quantitative Prediction of Fold Resistance for Inhibitors of EGFR. Biochemistry. 2009, 48: 8435-8448. 10.1021/bi900729a.

    CAS  Article  Google Scholar 

  32. 32.

    McGillick BE, Balius TE, Mukherjee S, Rizzo RC: Origins of Resistance to the HIVgp41 Viral Entry Inhibitor T20. Biochemistry. 2010, 49: 3575-3579. 10.1021/bi901915g.

    CAS  Article  Google Scholar 

  33. 33.

    Troise F, Cafaro V, Giancola C, D'Alessio G, De Lorenzo C: Differential binding of human immunoagents and Herceptin to the ErbB2 receptor. FEBS J. 2008, 275: 4967-4979. 10.1111/j.1742-4658.2008.06625.x.

    CAS  Article  Google Scholar 

  34. 34.

    Kästner J, Loeffler HH, Roberts SK, Martin-Fernandez ML, Winn MD: Ectodomain orientation, conformational plasticity and oligomerization of ErbB1 receptors investigated by molecular dynamics. J Struct Biol. 2009, 167: 117-128. 10.1016/j.jsb.2009.04.007.

    Article  Google Scholar 

  35. 35.

    Chao MP, Alizadeh AA, Tang C, Myklebust JH, Varghese B, Gill S, Jan M, Cha AC, Chan CK, Tan BT, Park CY, Zhao F, Kohrt HE, Malumbres R, Briones J, Gascoyne RD, Lossos IS, Levy R, Weissman IL, Majeti R: Anti-CD47 antibody synergizes with rituximab to promote phagocytosis and eradicate non-Hodgkin lymphoma. Cell. 2010, 142: 699-713. 10.1016/j.cell.2010.07.044.

    CAS  Article  Google Scholar 

  36. 36.

    Pickl M, Ries CH: Comparison of 3D and 2D tumor models reveals enhanced HER2 activation in 3D associated with an increased response to trastuzumab. Oncogene. 2009, 28: 461-468. 10.1038/onc.2008.394.

    CAS  Article  Google Scholar 

  37. 37.

    Weigelt B, Lo AT, Park CC, Gray JW, Bissell MJ: HER2 signaling pathway activation and response of breast cancer cells to HER2-targeting agents is dependent strongly on the 3D microenvironment. Breast Cancer Res Treat. 2010, 122: 35-43. 10.1007/s10549-009-0502-2.

    CAS  Article  Google Scholar 

  38. 38.

    Thakur A, Lum LG: Cancer therapy with bispecific antibodies: Clinical experience. Curr Opin Mol Ther. 2010, 12: 340-349.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references


CSV is adjunct at DBS (NUS) and SBS (NTU); Grant from European Research Council (AdG09 250244) is acknowledged.

Author information



Corresponding authors

Correspondence to Gloria Fuentes or Chandra S Verma.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

GF conceived of the study, participated in study design and data interpretation, and drafted the manuscript. MS and JB conceived of the study and revised the manuscript critically for important intellectual content. CV conceived of the study, participated in the analysis and interpretation, and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Supplemental Figure S1

Additional file 1: . Cartoon representation of all the systems involved in the study; A) Apo Her2; B) binary complex between Her2:pertuzumab; C) binary complex between Her2:trastuzumab; D) ternary complex between Her2:pertuzumab:trastuzumab; E) pertuzumab; F) trastuzumab; G) "binary" interaction between trastuzumab:pertuzumab. Trastuzumab has been coloured in black; while pertuzumab is shown in orange. (PDF 528 KB)

Supplemental Figure S2

Additional file 2: . The characterization of the collective protein dynamics have been analyzed using Principal Component Analysis (PCA). PCA is a statistical method that highlight similarities and differences in a complex data set by reducing the number of dimensions without losing much of the information. In MD simulations, the application of PCA can separate the configurational space into two subspaces; the essential one containing the functionally important correlated motions (comprising only few degrees of freedom or eigenvectors) and the irrelevant subspace with the independent Gaussian fluctuations with little functional relevance. In this study, PCA has been performed on an ensemble of conformations derived from the different MD simulations. Results of PCA obtained from the diagonalization of the atomic displacement correlation matrix of Cα atoms of the ectodomain from six of the MD simulations carried out in this study (using standard Euclidean distances) are shown here. The resulting Principal Components (PCs) are sorted according to their contribution to the total fluctuation along the ensemble of conformations, and only a small subset of these PCs are necessary to describe the great majority of the total atomic displacements which capture the essential dynamics of the conformational changes occurring in the system (the above so-referred essential subspace). Every figure contains the percentage of the cumulative eigenvalues as a function of the number of eigenvectors considered. The inset includes the projection on the first two eigenvectors of the coordinates extracted from the time interval of 8 to 16 ns. Every figure represents the coordinates for the ectodomain of Her2 from the following simulations: (A) apo Her2; (B) Her2:pertuzumab:trastuzumab; (C) Her2:pertuzumab; (D) Her2:trastuzumab; (E) apo truncated Her2; (F) truncated Her2:pertuzumab:trastuzumab. The application of PCA on data obtained from MD simulations provides a useful way of detecting global, correlated motions of the system that are separate from the tightly-constrained harmonic motions of more constrained atoms. The structures projected on the first eigenvectors have facilitated our understanding of the dynamics of the six systems analysed here (see Supplementary Figure S6). (PDF 746 KB)

Supplemental Figure S3

Additional file 3: . Cα RMS deviation with respect to the starting frame of the 2-ns equilibration phase versus time. (A) apo Her2 simulation; (B) ternary complex Her2:trastuzumab:pertuzumab; (C) binary complex Her2:trastuzumab; (D) binary complex Her2:pertuzumab. The top panel for all the plots represent the RMSD for the whole system; this is broken into the different components in the association in the bottom plots. (PDF 101 KB)

Supplemental Figure S4

Additional file 4: . RSM fluctuation plots and cartoon figures for the different systems of Her2 in the apo-and holo-states with pertuzumab and trastuzumab. A) RMS fluctuations for trastuzumab-containing systems (top); RMS fluctuations for pertuzumab-containing systems (middle panel) and RMS fluctuations for the ectodomain of the receptor in the four different MD simulations; B) putty/sausage cartoon representation for apo-receptor (right), Her2 in complex with the two antibodies, with trastuzumab and with pertuzumab (from top to bottom) and colored accordingly to the fluctuation values per residue. The complementarity determining regions (CDR) of both antibodies have been marked in the plot with blue boxes. The trastuzumab and pertuzumab epitopes have been highlighted in shaded boxes in grey and cream and black, respectively. (PDF 488 KB)

Supplemental Table S1

Additional file 5: . Cα-Cα distances for the residues at the receptor:Ab interfaces detected along the simulations between Her2 and the two Abs: pertuzumab and trastuzumab. The different shading boxes indicate the CDR of the antibodies. (DOC 342 KB)

Supplemental Figure S5

Additional file 6: . Interactions between the two antibodies and new epitope revealed in the MD simulations. (A) Density distributions for the intermolecular Cα-Cα distance between Glu294 and Asp295 (labeled in the Figure C) in the pertuzumab-induced trastuzumab epitope of the receptor and their putative interacting residues at the interface in trastuzumab. The dashed lines represent the population in the Her2:pertuzumab:trastuzumab simulation, while the straight lines show the distribution found in the Her2:trastuzumab simulation. (B) Density distribution for the intermolecular Cα-Cα distance between residues in the different Abs. In horizontal lines the average distance for the ternary (red) and binary (blue) simulations are shown. (C) Schematic representation of the residues for which the distances have been measured. The color coding is the same as in Figure S1. (PDF 318 KB)

Supplemental Figure S6

Additional file 7: . Cartoon representation of the maximum and minimum conformations of Her2 along the first eigenvector extracted from the PCA analysis and colored according to the Dyndom output (see Additional file 8) for the following simulations: (A) apo Her2; (B) Her2:pertuzumab:trastuzumab; (C) Her2:pertuzumab; (D) Her2:trastuzumab; (E) apo truncated Her2; (F) truncated Her2:pertuzumab:trastuzumab. All the models were fit onto the last 20 residues of the ectodomain. (PDF 515 KB)

Supplemental Table S2

Additional file 8: . Dynamical domains and hinge bending residues determined by the program DYNDOM for the superposition of the Her2 receptor of the trajectories of this domain in the apo, dimeric and trimeric forms with trastuzumab and pertuzumab. (DOC 58 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Cite this article

Fuentes, G., Scaltriti, M., Baselga, J. et al. Synergy between trastuzumab and pertuzumab for human epidermal growth factor 2 (Her2) from colocalization: an in silicobased mechanism. Breast Cancer Res 13, R54 (2011).

Download citation


  • Trastuzumab
  • Molecular Dynamic Simulation
  • Lapatinib
  • Ternary Complex
  • Pertuzumab