Analysis of correlated motion in antibody combining sites from molecular dynamics simulations

The crystal structures of the NC6.8-antisweet taste ligand complex and the uncomplexed antibody structures display significant differences in the conformations of residues in the combining site. A molecular dynamics method was employed to understand the flexibility and correlated motion of key combining site residues in the uncomplexed antibody. The simulations reveal that residues that show conformational differences between the complex and uncomplexed structures display strong dynamical correlations. Extensive analysis of the dynamics trajectory using time correlation methods is presented.     

Structural predictions of the binding site architecture for monoclonal antibody NC6.8 using computer-aided molecular modeling, ligand binding, and spectroscopy.

Monoclonal antibody NC6.8 binds the superpotent sweetener ligand N-(p-cyanophenyl)-N'-(diphenylmethyl) guanidineacetic acid with high affinity (Kd = 53 nM). Using computer-aided molecular modeling and several experimental techniques, such as competitive ligand binding, absorbance spectroscopy, and fluorescence spectroscopy, we have predicted the structure of the variable domain fragment (Fv) and identified the key residues in the combining site of the antibody. We have identified nine specific amino acids as being involved in ligand recognition and complexation. Most notable are H:33W, which is responsible for ligand-induced tryptophan fluorescence quenching, H:56R, which forms a salt bridge with the carboxylate moiety of the ligand, and L:34H, which, deep in the binding site, interacts with the cyanophenyl portion of the ligand. Two residues located deep in the putative binding pocket, H:35E and H:50E, provide the negatively charged potential for interaction with the protonated aryl nitrogen and the positive guanidinium group. These modeling predictions were made before the solution of high-resolution structures of the native Fab (2.6 A) and the Fab-ligand complex (2.2 A). Comparisons between the theoretical model and experimental native and liganded Fab structures are made.     

Crystal structure of an anti-interleukin-2 monoclonal antibody Fab complexed with an antigenic nonapeptide

The three-dimensional structure of the Fab fragment of a monoclonal antibody (LNKB-2) to human interleukin-2 (IL-2) complexed with a synthetic antigenic nonapeptide, Ac-Lys-Pro-Leu-Glu-Glu-Val-Leu-Asn-Leu-OMe, has been determined at 3.0 A resolution. In the structure, four out of the six hypervariable loops of the Fab (complementarity determining regions [CDRs] L1, H1, H2, and H3) are involved in peptide association through hydrogen bonding, salt bridge formation, and hydrophobic interactions. The Tyr residues in the Fab antigen binding site play a major role in antigen-antibody recognition. The structures of the complexed and uncomplexed Fab were compared. In the antigen binding site the CDR-L1 loop of the antibody shows the largest structural changes upon peptide binding. The peptide adopts a mostly alpha-helical conformation similar to that in the epitope fragment 64-72 of the IL-2 antigen. The side chains of residues Leu 66, Val 69, and Leu 70, which are shielded internally in the IL-2 structure, are involved in interactions with the Fab in the complex studied. This indicates that antibody-antigen complexation involves a significant rearrangement of the epitope-containing region of the IL-2 with retention of the alpha-helical character of the epitope fragment.     

Alternative oligomeric states of the yeast Rvb1/Rvb2 complex induced by histidine tags

The crystal structure of the cdk5/p25 complex has provided information on possible molecular mechanisms of the ligand binding, specificity, and regulation of the kinase. Comparative molecular dynamics simulations are reported here for physiological conditions. This study provides new insight on the mechanisms that modulate such processes, which may be exploited to control pathological activation by p25. The structural changes observed in the kinase are stabilized by a network of interactions involving highly conserved residues within the cyclin-dependent kinase (cdk) family. Collective motions of the proteins (cdk5, p25, and CIP) and their complexes are identified by principal component analysis, revealing two conformational states of the activation loop upon p25 complexation, which are absent in the uncomplexed kinase and not apparent from the crystal. Simulations of the uncomplexed inhibitor CIP show structural rearrangements and increased flexibility of the interfacial loop containing the critical residue E240, which becomes fully hydrated and available for interactions with one of several positively charged residues in the kinase. These changes provide a rationale for the observed high affinity and enhanced inhibitory action of CIP when compared to either p25 or the physiological activators of cdk5.     

Construction and characterization of two anti‐sweetener single chain antibodies using radioligand binding,fluorescence and circular dichroism spectroscopy

Two single‐chain antibodies (scFv) that bind the superpotent sweetener ligand, NC‐174, were generated from mouse monoclonal antibodies (mAb) NC6.8 (IgG, κ) and NC10.14 (IgG, λ). These scFv were constructed by cloning the variable region sequences of the mAb, connecting them in tandem with a 25‐amino‐acid polypeptide linker, and expressing them in E. coli using the pET‐11a system. The recombinant proteins were purified using Ni²⁺–NTA–agarose by virtue of a hexahistidine sequence introduced to the C‐terminus of the heavy chain variable region during the cloning process. The secondary structure and ligand binding properties of the two scFv, the parent mAbs and proteolytically derived Fab fragments were examined using radioligand binding, circular dichroism (CD) and fluorescence spectroscopy. The far‐UV CD spectra of both scFv possessed predominantly β character, as did those of the Fab, and the near‐UV CD spectral data for scFvNC10.14, NC6.8 and NC10.14 Fab indicated that chromophore perturbation occurred upon ligand binding. The affinity constants determined for the two scFv, Fab and mAb were nearly equivalent. Copyright © 1999 John Wiley & Sons, Ltd.     

Mutual conformational adaptations in antigen and antibody upon complex formation between an Fab and HIV-1 capsid protein p24

BACKGROUND: Elucidating the structural basis of antigen-antibody recognition ideally requires a structural comparison of free and complexed components. To this end we have studied a mouse monoclonal antibody, denoted 13B5, raised against p24, the capsid protein of HIV-1. We have previously described the first crystal structure of intact p24 as visualized in the Fab13B5-p24 complex. Here we report the structure of the uncomplexed Fab13B5 at 1.8 A resolution and analyze the Fab-p24 interface and the conformational changes occurring upon complex formation. RESULTS: Fab13B5 recognizes a nearly continuous epitope comprising a helix-turn-helix motif in the C-terminal domain of p24. Only 4 complementarity-determining regions (CDRs) are in contact with p24 with most interactions being by the heavy chain. Comparison of the free and complexed Fab reveals that structural changes upon binding are localized to a few side chains of CDR-H1 and -H2 but involve a larger, concerted displacement of CDR-H3. Antigen binding is also associated with an 8 degrees relative rotation of the heavy and light chain variable regions. In p24, small conformational changes localized to the turn between the two helices comprising the epitope result from Fab binding. CONCLUSIONS: The relatively small area of contact between Fab13B5 and p24 may be related to the fact that the epitope is a continuous peptide rather than a more complex protein surface and correlates with a relatively low affinity of antigen and antibody. Despite this, a significant quaternary structural change occurs in the Fab upon complex formation, with additional smaller adaptations of both antigen and antibody.     

The three-dimensional structure of a complex of a murine Fab (NC10. 14) with a potent sweetener (NC174): an illustration of structural diversity in antigen recognition by immunoglobulins

The three-dimensional structure of a complex of an Fab from a murine IgG2b(lambda) antibody (NC10.14) with a high potency sweet tasting hap- ten, N-(p-cyanophenyl)-N'-(diphenylmethyl)-N"-(carboxymethyl)guan idine (NC174), has been determined to 2.6 A resolution by X-ray crystallography. This complex crystallized in the triclinic space group P1, with two molecules in the asymmetric unit. In contrast to a companion monoclonal antibody (NC6.8) with a kappa-type light chain and similar high affinity for the NC174 ligand, the NC10.14 antibody possessed a large and deep antigen combining site bounded primarily by the third complementarity-determining regions (CDR3s) of the light and heavy chains. CDR3 of the heavy chain dominated the site and its crown protruded into the external solvent as a type 1' beta-turn. NC174 was nested against HCDR3 and was held in place by two tryptophan side-chains (L91 and L96) from LCDR3. The diphenyl rings were accommodated on an upper tier of the binding pocket that is largely hydrophobic. At the floor of the site, a positively charged arginine side-chain (H95) stabilized the orientation of the electronegative cyano group of the hapten. The negative charge on the acetate group was partially neutralized by a hydrogen bond with the phenolic hydroxyl group of tyrosine H58. Comparisons of the modes of binding of NC174 to the NC6.8 and NC10.14 antibodies illustrate the enormous structural and mechanistic diversity manifest by immune responses.     

Probing flexibility and "induced-fit" phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations

Aldose reductase is a promising target for the treatment of diabetic complications, and as such, has become the focus of various drug design projects. As revealed by a survey of available crystal structures, the protein shows pronounced induced-fit effects upon ligand binding. Although helping to explain the enzyme's substrate promiscuity, phenomena of this kind are still responsible for significant complications in structure-based design efforts directed to aldose reductase. Accordingly, a deeper understanding of the principles governing conformational alterations in this enzyme would be of utmost practical importance. As a first step in addressing this issue, molecular dynamics (MD) simulations have been carried out. The ultrahigh resolution crystal structure of aldose reductase complexed with inhibitor IDD594 served as ideal starting point for a set of different simulations of nanosecond time scale: the native complexed state with bound inhibitor, the uncomplexed state (after removal of the inhibitor) at standard temperature, and the uncomplexed state at elevated temperature. The reference simulation of the complex exhibits extraordinary stability of the overall fold, whereas two distinct conformational substates are found for the binding-site region. In contrast, already at standard temperature pronounced changes are observed in the binding region during the simulation of the uncomplexed state. Leu300, for example, closes the access to the pocket opened by IDD594. On the other hand, conformations around the catalytic site are highly conserved, with the His110-Tyr48-NADP+ orientation being stabilized by a water molecule. Detailed analysis of the trajectories allows to reveal a set of distinct conformational substates that may prove useful as alternative structural templates in virtual screening for new aldose reductase inhibitors.     

Computer-aided molecular modeling of the binding site architecture for eight monoclonal antibodies that bind a high potency guanidinium sweetener

Computer-aided molecular modeling of the antibody binding site of eight different monoclonal antibodies (mAb) that bind the intense sweetener ligand (N-(p-cyanophenyl)-N'-diphenylmethyl) guanidine acetic acid was completed using canonical loop structures and framework regions from known immuno globulins as “parent structures” for the molecular scaffoldings. The models of the fragment variable (Fv) region of the mAb were analyzed/or the presence and location of residues predicted to be involved in ligand binding. Several binding site tryptophan residues in these models were located in positions that support previous flurospectroscopic observations of the mAb-ligand complexation. Computer-aided renderings of the electrostatic potential at the van der Waals surface of the Fv region were compared and found to be consistent with the ligand binding specificity profiles for the different mAb. The Fv model of mAb NC6.8 was consistent with the binding site features determined in the Fab structure recently solved by X-ray diffraction techniques. These Fv models should provide an adequate basis for site-directed mutagenesis experiments in order to characterize interactive motifs in the mAb binding site.     

Models for binding cooperativities of inhibitors with transthyretin

Here, molecular dynamics (MD) simulations are performed to study the differences of binding channel shapes of TTR with two inhibitors, flufenamic acid (FLU) and one kind of N-phenyl phenoxazine (BPD). The asymmetries of global structure including the central binding channel are found to be intrinsic. Moreover, the conformational changes of the binding channel are responsible for negative cooperativity (NC) or independent cooperativity (IC) of ligands. The results suggested a possible binding mechanism addressing NC of FLU and IC of BPD. For FLU, when the first ligand binds with TTR, it leads to expansion of the second binding site which may weaken the interaction of the second FLU with TTR. But for BPD, the first ligand's binding changes the second site's shape slightly, the second ligand has similar binding ability with TTR in the second site like the first binding event.     

Coarse‐grained and all‐atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen‐binding protein composed of four subunits (α1, α2, β1, and β2), is a well‐known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse‐grained (CG) modeling, all‐atom simulation, and structural data from X‐ray crystallography and wide‐angle X‐ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot‐spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all‐atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all‐atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Three-dimensional structures of influenza virus neuraminidase-antibody complexes

X-ray diffraction analysis of crystals of a monoclonal Fab fragment NC41 bound to a viral antigen, influenza virus neuraminidase, shows an epitope involving five surface loops of the antigen. In addition it reveals an unusual pairing pattern between the domains of light and heavy chains in the variable module of the antibody. We interpret this result to imply that association with antigen can induce changes in the quaternary structure of the Fab, through a sliding of domains at the variable light/variable heavy chains (VL-VH) interface. In addition, Fab binding has altered the conformation of some of the surface loops of the antigen. The structure of the NC10 Fab-neuraminidase complex has now also been solved. It binds an epitope that overlaps the NC41 epitope. In this structure, there is no electron density for the C-module of the Fab fragment, implying it is disordered in the crystal lattice. The implications of these, and other antibody-antigen structures, for immune recognition are discussed.     

Cocrystal structures of NC6.8 Fab identify key interactions for high potency sweetener recognition: implications for the design of synthetic sweeteners

The crystal structures of the murine monoclonal IgG2b(kappa) antibody NC6.8 Fab fragment complexed with high-potency sweetener compound SC45647 and nontasting high-affinity antagonist TES have been determined. The crystal structures show how sweetener potency is fine-tuned by multiple interactions between specific receptor residues and the functionally different groups of the sweeteners. Comparative analysis with the structure of NC6.8 complexed with the super-potency sweetener NC174 reveals that although the same residues in the antigen binding pocket of NC6.8 interact with the zwitterionic, trisubstituted guanidinium sweeteners as well as TES, specific differences exist and provide guidance for the design of new artificial sweeteners. In case of the nonsweetener TES, the interactions with the receptor are indirectly mediated through a hydrogen bonded water network, while the sweeteners bind with high affinity directly to the receptor. The presence of a hydrophobic group interacting with multiple receptor residues as a major determinant for sweet taste has been confirmed. The nature of the hydrophobic group is likely a discriminator for super- versus high-potency sweeteners, which can be exploited in the design of new, highly potent sweetener compounds. Overall similarities and partial conservation of interactions indicate that the NC6.8 Fab surrogate is representing crucial features of the T1R2 taste receptor VFTM binding site.     

Crystal structures of two mutant neuraminidase-antibody complexes with amino acid substitutions in the interface.

The site on influenza virus N9 neuraminidase recognized by NC41 monoclonal antibody comprises 19 amino acid residues that are in direct contact with 17 residues on the antibody. Single sequence changes in some of the neuraminidase residues in the site markedly reduce antibody binding. However, two mutants have been found within the site, Ile368 to Arg and Asn329 to Asp selected by antibodies other than NC41, and these mutants bind NC41 antibody with only slightly reduced affinity. The three-dimensional structures of the two mutant N9-NC41 antibody complexes as derived from the wild-type complex are presented. Both structures show that some amino acid substitutions can be accommodated within an antigen-antibody interface by local structural rearrangements around the mutation site. In the Ile368 to Arg mutant complex, the side-chain of Arg368 is shifted by 2.9 A from its position in the uncomplexed mutant and a shift of 1.3 A in the position of the light chain residue HisL55 with respect to the wild-type complex is also observed. In the other mutant, the side-chain of Asp329 appears rotated by 150 degrees around C alpha-C beta with respect to the uncomplexed mutant, so that the carboxylate group is moved to the periphery of the antigen-antibody interface. The results provide a basis for understanding some of the potential structural effects of somatic hypermutation on antigen-antibody binding in those cases where the mutation in the antibody occurs at antigen-contacting residues, and demonstrate again the importance of structural context in evaluating the effect of amino acid substitutions on protein structure and function.     

Molecular dynamics simulations of the ligand-binding domain of the ionotropic glutamate receptor GluR2.

Ionotropic glutamate receptors are essential for fast synaptic nerve transmission. Recent x-ray structures for the ligand-binding (S1S2) region of the GluR2 alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA)-sensitive receptor have suggested how differences in protein/ligand interactions may determine whether a ligand will behave as a full agonist. We have used multiple molecular dynamics simulations of 2-5 ns duration to explore the structural dynamics of GluR2 S1S2 in the presence and absence of glutamate and in a complex with kainate. Our studies indicate that not only is the degree of domain closure dependent upon interactions with the ligand, but also that protein/ligand interactions influence the motion of the S2 domain with respect to S1. Differences in domain mobility between the three states (apo-S1S2, glutamate-bound, and kainate-bound) are surprisingly clear-cut. We discuss how these changes in dynamics may provide an explanation relating the mechanism of transmission of the agonist-binding event to channel opening. We also show here how the glutamate may adopt an alternative mode of binding not seen in the x-ray structure, which involves a key threonine (T480) side chain flipping into a new conformation. This new conformation results in an altered pattern of hydrogen bonding at the agonist-binding site.     

Induced‐fit or preexisting equilibrium dynamics? Lessons from protein crystallography and MD simulations on acetylcholinesterase and implications for structure‐based drug design

Crystal structures of acetylcholinesterase complexed with ligands are compared with side-chain conformations accessed by native acetylcholinesterase in molecular dynamics (MD) simulations. Several crystallographic conformations of a key residue in a specific binding site are accessed in a simulation of native acetylcholinesterase, although not seen in rotomer plots. Conformational changes upon ligand binding thus involve preexisting equilibrium dynamics. Consequently, rational drug design could benefit significantly from conformations monitored by MD simulations of native targets.     

Molecular dynamics simulation studies of induced fit and conformational capture in U1A–RNA binding: Do molecular substates code for specificity?

Molecular dynamics (MD) simulations on stem loop 2 of U1 small nuclear RNA and a construct of the U1A protein were carried out to obtain predictions of the structures for the unbound forms in solution and to elucidate dynamical aspects of induced fit upon binding. A crystal structure of the complex between the U1A protein and stem loop 2 RNA and an NMR structure for the uncomplexed form of the U1A protein are available from Oubridge et al. (Nature, 1994, Vol. 372, pp. 432-438) and Avis et al. (Journal of Molecular Biology, 1996, Vol. 257, pp. 398-411), respectively. As a consequence, U1A-RNA binding is a particularly attractive case for investigations of induced fit in protein-nucleic acid complexation. When combined with the available structural data, the results from simulations indicate that structural adaptation of U1A protein and RNA define distinct mechanisms for induced fit. For the protein, the calculations indicate that induced fit upon binding involves a non-native thermodynamic substate in which the structure is preorganized for binding. In contrast, induced fit of the RNA involves a distortion of the native structure in solution to an unstable form. However, the RNA solution structures predicted from simulation show evidence that structures in which groups of bases are favorably oriented for binding the U1A protein are thermally accessible. These results, which quantify with computational modeling recent proposals on induced fit and conformational capture by Leuillot and Varani (Biochemistry, 2001, Vol. 40, pp. 7947-7956) and by Williamson (Nature Structural Biology, 2000, Vol. 7, pp. 834-837) suggest an important role for intrinsic molecular architecture and substates other than the native form in the specificity of protein-RNA interactions.     

Molecular Dynamics Simulations of Calcium-Free Calmodulin in Solution

Abstract

A 4-ns molecular dynamics simulation of calcium-free calmodulin in solution has been performed, using Ewald summation to treat electrostatic interactions. Our simulation results were mostly consistent with solution experimental studies, including NMR, fluorescence and x-ray scattering. The secondary structures within the N- and C-terminal domains were conserved in the simulation, with trajectory structures similar to the NMR-derived model structure 1CFD. However, the relative orientations of the domains, for which there are no NMR restraints, differed in details between the simulation and the 1CFD model. The most interesting information provided by the simulations is that the dynamics of calcium-free calmod- ulin in solution is dominated by slow rigid body reorientations of the domains. The interdomain distance fluctuated between 29 and 39 Å, and interdomain orientation angle, defined as the pseudo-dihedral formed by the four calcium binding sites, varied between ?2° and 108°. Similarly, the domain linker region also exhibited significant fluctuations, with its length varying in the 34–45 Å range and its bend angle in the 10–100° range. The simulations are in accord with fluorescence results suggesting that calcium-free calmodulin is more compact and more flexible than the calcium activated form. Surprisingly, quite similar solvent accessibilities of the hydrophobic patches were seen in the calcium-free trajectory described in this work and previously generated calcium-loaded calmodulin simulations. Thus, our simulations suggest a reexamination of the standard model of the structural change of calmodulin upon calcium binding, involving exposure of the hydrophobic patches to solvent.     

Retinal Flip in Rhodopsin Activation?

Rhodopsin is a well-characterized structural model of a G protein-coupled receptor. Photoisomerization of the covalently bound retinal triggers activation. Surprisingly, the x-ray crystal structure of the active Meta-II state has a 180° rotation about the long-axis of the retinal polyene chain. Unbiased microsecond-timescale all-atom molecular dynamics simulations show that the retinal cofactor can flip back to the orientation observed in the inactive state of rhodopsin under conditions favoring the Meta-I state. Our results provide, to our knowledge, the first evidence from molecular dynamics simulations showing how rotation of the retinal ligand within its binding pocket can occur in the activation mechanism of rhodopsin.