Immunologic memory to phosphocholine. IV. Hybridomas representative of Group I (T15-like) and Group II (non-T15-like) antibodies utilize distinct VH genes

The anti-phosphocholine (PC) memory response elicited in BALB/c mice by phosphocholine-keyhole limpet hemocyanin (PC-KLH) contains two groups of antibodies distinguished by their fine specificity for PC and p-nitrophenylphosphocholine (NPPC). Group I antibodies are inhibited by both PC and NPPC, while Group II antibodies are inhibited appreciably only by NPPC; only Group I antibodies are dominated by the T15 idiotype. Anti-PC hybridomas representative of the memory response to PC-KLH were produced to examine the variable region genes expressed by memory B cells. Two IgM hybridomas were of the Group I type, because they were inhibited by both PC and NPPC and they bound to the pneumococcus R36A. However, only one of these antibodies (PCM-2) expressed a T15 idiotope, while the other (PCM-1) did not express any of three T15 idiotopes. Despite its negative T15 idiotype profile, N-terminal amino acid sequencing of PCM-1 purified heavy chain and Southern blots of the hybridoma DNA indicated that it utilizes the T15 VH and JH1 genes. Three hybridomas, IgG1, IgM, and IgE, typical of Group II antibodies, were examined; these were negative for three T15 idiotopes and displayed measurable avidity only for NPPC in a PC-protein binding inhibition assay. These three hybridoma antibodies, like serum Group II IgG1, did not measurably bind to the bacterium R36A. The heavy chain amino termini of all three of these antibodies were inaccessible for Edman degradation. Southern blots of DNA from the IgG1 hybridoma revealed the T15 VH gene to be in the germ line configuration only and unassociated with any JH segment, indicating that this Group II antibody utilizes a VH gene different from the T15 family. These results signify that, whereas some diversity of the (anti-PC) memory response may be generated by somatic diversification of variable regions important in the primary response, a significant contribution to the overall heterogeneity of memory antibodies originates in the expression of additional variable region genes.     

Expression, purification and the 1.8 angstroms resolution crystal structure of human neuron specific enolase

Human neuron-specific enolase (NSE) or isozyme gamma has been expressed with a C-terminal His-tag in Escherichia coli. The enzyme has been purified, crystallized and its crystal structure determined. In the crystals the enzyme forms the asymmetric complex NSE x Mg2 x SO4/NSE x Mg x Cl, where "/" separates the dimer subunits. The subunit that contains the sulfate (or phosphate) ion and two magnesium ions is in the closed conformation observed in enolase complexes with the substrate or its analogues; the other subunit is in the open conformation observed in enolase subunits without bound substrate or analogues. This indicates negative cooperativity for ligand binding between subunits. Electrostatic charge differences between isozymes alpha and gamma, -19 at physiological pH, are concentrated in the regions of the molecular surface that are negatively charged in alpha, i.e. surface areas negatively charged in alpha are more negatively charged in gamma, while areas that are neutral or positively charged tend to be charge-conserved.     

Binding interactions between the core central domain of 16S rRNA and the ribosomal protein S15 determined by molecular dynamics simulations

The goal of the current study is to utilize molecular dynamic (MD) simulations to investigate the dynamic behavior of 16S rRNA in the presence and absence of S15 and to identify the binding interactions between these two molecules. The simulations show that: (i) 16S rRNA remains in a highly folded structure when it is bound to S15; (ii) in the absence of S15, 16S rRNA significantly alters its conformation and transiently forms conformations that are similar to the bound structure that make it available for binding with S15; (iii) the unbound rRNA spends the majority of its time in extended conformations. The formation of the extended conformations is a result of the molecule reaching a lower electrostatic energy and the formation of the highly folded, crystal-like conformation is a result of achieving a lower solvation energy. In addition, our MD simulations show that 16S rRNA and S15 bind across the major groove of helix 22 (H22) via electrostatic interactions. The negatively charged phosphate groups of G658, U740, G741 and G742 bind to the positively charged S15 residues Lys7, Arg34 and Arg37. The current study provides a dynamic view of the binding of 16S rRNA with S15.     

Domain:domain interactions within Hop, the Hsp70/Hsp90 organizing protein, are required for protein stability and structure

The major heat shock protein (Hsp) chaperones Hsp70 and Hsp90 both bind the co-chaperone Hop (Hsp70/Hsp90 organizing protein), which coordinates Hsp actions in folding protein substrates. Hop contains three tetratricopeptide repeat (TPR) domains that have binding sites for the conserved EEVD C termini of Hsp70 and Hsp90. Crystallographic studies have shown that EEVD interacts with positively charged amino acids in Hop TPR-binding pockets (called carboxylate clamps), and point mutations of these carboxylate clamp positions can disrupt Hsp binding. In this report, we use circular dichroism to assess the effects of point mutations and Hsp70/Hsp90 peptide binding on Hop conformation. Our results show that Hop global conformation is destabilized by single point mutations in carboxylate clamp positions at pH 5, while the structure of individual TPR domains is unaffected. Binding of peptides corresponding to the C termini of Hsp70 and Hsp90 alters the global conformation of wild-type Hop, whereas peptide binding does not alter conformation of individual TPR domains. These results provide biophysical evidence that Hop-binding pockets are directly involved with domain:domain interactions, both influencing Hop global conformation and Hsp binding, and contributing to proper coordination of Hsp70 and Hsp90 interactions with protein substrates.     

Structural studies of the streptavidin binding loop.

The streptavidin-biotin complex provides the basis for many important biotechnological applications and is an interesting model system for studying high-affinity protein-ligand interactions. We report here crystallographic studies elucidating the conformation of the flexible binding loop of streptavidin (residues 45 to 52) in the unbound and bound forms. The crystal structures of unbound streptavidin have been determined in two monoclinic crystal forms. The binding loop generally adopts an open conformation in the unbound species. In one subunit of one crystal form, the flexible loop adopts the closed conformation and an analysis of packing interactions suggests that protein-protein contacts stabilize the closed loop conformation. In the other crystal form all loops adopt an open conformation. Co-crystallization of streptavidin and biotin resulted in two additional, different crystal forms, with ligand bound in all four binding sites of the first crystal form and biotin bound in only two subunits in a second. The major change associated with binding of biotin is the closure of the surface loop incorporating residues 45 to 52. Residues 49 to 52 display a 3(10) helical conformation in unbound subunits of our structures as opposed to the disordered loops observed in other structure determinations of streptavidin. In addition, the open conformation is stabilized by a beta-sheet hydrogen bond between residues 45 and 52, which cannot occur in the closed conformation. The 3(10) helix is observed in nearly all unbound subunits of both the co-crystallized and ligand-free structures. An analysis of the temperature factors of the binding loop regions suggests that the mobility of the closed loops in the complexed structures is lower than in the open loops of the ligand-free structures. The two biotin bound subunits in the tetramer found in the MONO-b1 crystal form are those that contribute Trp 120 across their respective binding pockets, suggesting a structural link between these binding sites in the tetramer. However, there are no obvious signatures of binding site communication observed upon ligand binding, such as quaternary structure changes or shifts in the region of Trp 120. These studies demonstrate that while crystallographic packing interactions can stabilize both the open and closed forms of the flexible loop, in their absence the loop is open in the unbound state and closed in the presence of biotin. If present in solution, the helical structure in the open loop conformation could moderate the entropic penalty associated with biotin binding by contributing an order-to-disorder component to the loop closure.     

A syn-anti conformational difference allows SRSF2 to recognize guanines and cytosines equally well

SRSF2 (SC35) is a key player in the regulation of alternative splicing events and binds degenerated RNA sequences with similar affinity in nanomolar range. We have determined the solution structure of the SRSF2 RRM bound to the 5'-UCCAGU-3' and 5'-UGGAGU-3' RNA, both identified as SRSF2 binding sites in the HIV-1 tat exon 2. RNA recognition is achieved through a novel sandwich-like structure with both termini forming a positively charged cavity to accommodate the four central nucleotides. To bind both RNA sequences equally well, SRSF2 forms a nearly identical network of intermolecular interactions by simply flipping the bases of the two consecutive C or G nucleotides into either anti or syn conformation. We validate this unusual mode of RNA recognition functionally by in-vitro and in-vivo splicing assays and propose a 5'-SSNG-3' (S=C/G) high-affinity binding consensus sequence for SRSF2. In conclusion, in addition to describe for the first time the RNA recognition mode of SRSF2, we provide the precise consensus sequence to identify new putative binding sites for this splicing factor.     

Selected-fit versus induced-fit protein binding: kinetic differences and mutational analysis

The binding of a ligand molecule to a protein is often accompanied by conformational changes of the protein. A central question is whether the ligand induces the conformational change (induced-fit), or rather selects and stabilizes a complementary conformation from a pre-existing equilibrium of ground and excited states of the protein (selected-fit). We consider here the binding kinetics in a simple four-state model of ligand-protein binding. In this model, the protein has two conformations, which can both bind the ligand. The first conformation is the ground state of the protein when the ligand is off, and the second conformation is the ground state when the ligand is bound. The induced-fit mechanism corresponds to ligand binding in the unbound ground state, and the selected-fit mechanism to ligand binding in the excited state. We find a simple, characteristic difference between the on- and off-rates in the two mechanisms if the conformational relaxation into the ground states is fast. In the case of selected-fit binding, the on-rate depends on the conformational equilibrium constant, whereas the off-rate is independent. In the case of induced-fit binding, in contrast, the off-rate depends on the conformational equilibrium, while the on-rate is independent. Whether a protein binds a ligand via selected-fit or induced-fit thus may be revealed by mutations far from the protein's binding pocket, or other "perturbations" that only affect the conformational equilibrium. In the case of selected-fit, such mutations will only change the on-rate, and in the case of induced-fit, only the off-rate.     

Minimal ensembles of side chain conformers for modeling protein–protein interactions

The goal of this article is to reduce the complexity of the side chain search within docking problems. We apply six methods of generating side chain conformers to unbound protein structures and determine their ability of obtaining the bound conformation in small ensembles of conformers. Methods are evaluated in terms of the positions of side chain end groups. Results for 68 protein complexes yield two important observations. First, the end‐group positions change less than 1 Å on association for over 60% of interface side chains. Thus, the unbound protein structure carries substantial information about the side chains in the bound state, and the inclusion of the unbound conformation into the ensemble of conformers is very beneficial. Second, considering each surface side chain separately in its protein environment, small ensembles of low‐energy states include the bound conformation for a large fraction of side chains. In particular, the ensemble consisting of the unbound conformation and the two highest probability predicted conformers includes the bound conformer with an accuracy of 1 Å for 78% of interface side chains. As more than 60% of the interface side chains have only one conformer and many others only a few, these ensembles of low‐energy states substantially reduce the complexity of side chain search in docking problems. This approach was already used for finding pockets in protein–protein interfaces that can bind small molecules to potentially disrupt protein–protein interactions. Side‐chain search with the reduced search space will also be incorporated into protein docking algorithms. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Escherichia coli peptide binding protein OppA has a preference for positively charged peptides

The Escherichia coli peptide binding protein OppA is an essential component of the oligopeptide transporter Opp. Based on studies on its orthologue from Salmonella typhimurium, it has been proposed that OppA binds peptides between two and five amino acids long, with no apparent sequence selectivity. Here, we studied peptide binding to E. coli OppA directly and show that the protein has an unexpected preference for basic peptides. OppA was expressed in the periplasm, where it bound to available peptides. The protein was purified in complex with tightly bound peptides. The crystal structure (up to 2.0 Å) of OppA liganded with the peptides indicated that the protein has a preference for peptides containing a lysine. Mass spectrometry analysis of the bound peptides showed that peptides between two and five amino acids long bind to the protein and indeed hinted at a preference for positively charged peptides. The preference of OppA for peptides with basic residues, in particular lysines, was corroborated by binding studies with peptides of defined sequence using isothermal titration calorimetry and intrinsic protein fluorescence titration. The protein bound tripeptides and tetrapeptides containing positively charged residues with high affinity, whereas related peptides without lysines/arginines were bound with low affinity. A structure of OppA in an open conformation in the absence of ligands was also determined to 2.0 Å, revealing that the initial binding site displays a negative surface charge, consistent with the observed preference for positively charged peptides. Taken together, E. coli OppA appears to have a preference for basic peptides.     

Structure-activity relationship and site of binding of polyamine derivatives at the nicotinic acetylcholine receptor.

Several wasp venoms contain philanthotoxins (PhTXs), natural polyamine amides, which act as noncompetitive inhibitors (NCIs) on the nicotinic acetylcholine receptor (nAChR). Effects of varying the structure of PhTXs and poly(methylene tetramine)s on the binding affinity have been investigated. Using the fluorescent NCI ethidium in a displacement assay Kapp values of these compounds have been determined. We found that an increase in size of the PhTX's hydrophobic head group significantly increased the binding affinity, while inserting positive charge almost completely destroyed it. Elongating the PhTX polyamine chain by introducing an additional aminomethylene group decreased the binding affinity, whereas a terminal lysine improved it. In general, poly(methylene tetramine)s showed higher binding affinities than PhTX analogues. The stoichiometry of PhTX binding was determined to be two PhTX molecules per receptor monomer. PhTXs appeared to bind to a single class of nonallosterically interacting binding sites and bound PhTX was found to be completely displaced by well-characterized luminal NCIs. To elucidate the site of PhTX binding, a photolabile, radioactive PhTX derivative was photocross-linked to the nAChR in its closed channel conformation resulting in labeling yields for the two alpha and the beta, gamma and delta subunits of 10.4, 11.1, 4.0 and 7.4%, respectively. Based on these findings we suggest that PhTXs and poly(methylene tetramine)s enter the receptor's ionic channel from the extracellular side. The hydrophobic head groups most likely bind to the high-affinity NCI site, while the positively charged polyamine chains presumably interact with the negatively charged selectivity filter located deep in the channel lumen.     

Design, synthesis, and conformational study of biologically active photolabeled analogues of the main immunogenic region of the acetylcholine receptor

Photoaffinity labeling is a powerful tool for the characterization of the molecular basis of ligand binding to acceptor molecules, which provides important insights for mapping the bimolecular interfaces. The autoimmune disease myasthenia gravis is caused by autoantibodies against the acetylcholine receptor (AChR). The majority of the anti-AChR antibodies bind to the "main immunogenic region" (MIR) of the AChR. To identify the contact points between the complementarity determining regions of the anti-MIR antibodies that recognize the MIR contact sites of the AChR, we present here three photoreactive dodecapeptide MIR analogues containing the photolabel p-benzoyl-L-phenylalanine (Bpa) moiety, either in position 1 or 11. The structure of the produced 12-mers was analyzed using two-dimensional (1)H-NMR spectroscopy, whereas their binding to anti-MIR monoclonal antibodies (mAbs) was determined by immunochemical assays. In all cases the modifications resulted in conservation of the beta-turn conformation of the N-terminus, which has been proved essential for antibody recognition and increased anti-MIR binding relative to the MIR decapeptide.     

Essential motions and energetic contributions of individual residues in a peptide bound to an SH3 domain

We have studied protein-ligand interactions by molecular dynamics simulations using software designed to exploit parallel computing architectures. The trajectories were analyzed to extract the essential motions and to estimate the individual contributions of fragments of the ligand to overall binding enthalpy. Two forms of the bound ligand are compared, one with the termini blocked by covalent derivatization, and one in the underivatized, zwitterionic form. The ends of the peptide tend to bind more loosely in the capped form. We can observe significant motions in the bound ligand and distinguish between motions of the peptide backbone and of the side chains. This could be useful in designing ligands, which fit optimally to the binding protein. We show that it is possible to determine the different contributions of each residue in a peptide to the enthalpy of binding. Proline is a major net contributor to binding enthalpy, in keeping with the known propensity for this family of proteins to bind proline-rich peptides.     

Insights into binding of fatty acids by fatty acid binding proteins

Members of the phylogenetically related intracellular lipid binding protein (iLBP) are characterized by a highly conserved tertiary structure, but reveal distinct binding preferences with regard to ligand structure and conformation, when binding is assessed by the Lipidex method (removal of unbound ligand by hydrophobic polymer) or by isothermal titration calorimetry, a true equilibrium method. Subfamily proteins bind retinoids, subfamily II proteins bind bulky ligands, examples are intestinal bile acid binding protein (I-BABP) and liver fatty acid binding protein (L-FABP) which binds 2 ligand molecules, preferably monounsaturated and n-3 fatty acids. Subfamily III intestinal fatty acid binding protein (I-FABP) binds fatty acid in a bent conformation. The fatty acid bound by subfamily IV FABPs has a U-shaped conformation; here heart (H-) FABP preferably binds n-6, brain (B-) FABP n-3 fatty acids. The ADIFAB-method is a fluorescent test for fatty acid in equilibrium with iLBP and reveals some correlation of binding affinity to fatty acid solubility in the aqueous phase; these data are often at variance with those obtained by the other methods. Thus, in this review published binding data are critically discussed, taking into account on the one hand binding increments calculated for fatty acid double bonds on the basis of the solubility hypothesis, on the other hand the interpretation of calorimetric data on the basis of crystallographic and solution structures of iLBPs.     

A preformed binding interface in the unbound ensemble of an intrinsically disordered protein: evidence from molecular simulations

Intrinsically disordered proteins play an important role in cellular signalling, mediated by their interactions with other biomolecules. A key question concerns the nature of their binding mechanism, and whether the bound structure is induced only by proximity to the binding partner. This is difficult to answer through experiment alone because of the very heterogeneous nature of the unbound ensemble, and the probable rapid interconversion of the various unbound structures. Here we report the most extensive set of simulations on NCBD to date: we use large-scale replica exchange molecular dynamics to explore the unbound state. An important feature of the study is the use of an atomistic force field that has been parametrised against experimental data for weakly structured peptides, together with an accurate explicit water model. Neither the force field nor the starting conformations are biased towards a particular structure. The regions of NCBD that have high helical propensity in the simulations correspond closely to helices in the 'core' unbound conformation determined by NMR, although no single member of the simulated unbound ensemble closely resembles the core conformation, or either of the two known bound conformations. We have validated the results against NMR spectroscopy and SAXS measurements, obtaining reasonable agreement. The two helices which most stabilise the binding of NCBD with ACTR are formed readily; the third helix, which is less important for binding but is involved in most of the intraprotein contacts of NCBD in the bound conformation, is formed more rarely, and tends not to coexist with the other helices. These results support a mechanism by which NCBD gains the advantages of disorder, while forming binding-competent structures in the unbound state. We obtain support for this mechanism from coarse-grained simulations of NCBD with, and without, its binding partner.     

Atomic-level insights into metabolite recognition and specificity of the SAM-II riboswitch

Although S-adenosylhomocysteine (SAH), a metabolic by-product of S-adenosylmethionine (SAM), differs from SAM only by a single methyl group and an overall positive charge, SAH binds the SAM-II riboswitch with more than 1000-fold less affinity than SAM. Using atomistic molecular dynamics simulations, we investigated the molecular basis of such high selectivity in ligand recognition by SAM-II riboswitch. The biosynthesis of SAM exclusively generates the (S,S) stereoisomer, and (S,S)-SAM can spontaneously convert to the (R,S) form. We, therefore, also examined the effects of (R,S)-SAM binding to SAM-II and its potential biological function. We find that the unfavorable loss in entropy in SAM-II binding is greater for (S,S)- and (R,S)-SAM than SAH, which is compensated by stabilizing electrostatic interactions with the riboswitch. The positively charged sulfonium moiety on SAM acts as the crucial anchor point responsible for the formation of key ionic interactions as it fits favorably in the negatively charged binding pocket. In contrast, SAH, with its lone pair of electrons on the sulfur, experiences repulsion in the binding pocket of SAM-II and is enthalpically destabilized. In the presence of SAH, similar to the unbound riboswitch, the pseudoknot structure of SAM-II is not completely formed, thus exposing the Shine-Dalgarno sequence. Unlike SAM, this may further facilitate ribosomal assembly and translation initiation. Our analysis of the conformational ensemble sampled by SAM-II in the absence of ligands and when bound to SAM or SAH reveals that ligand binding follows a combination of conformational selection and induced-fit mechanisms.     

Stabilization of antibody structure upon association to a human carbonic anhydrase IX epitope studied by X-ray crystallography, microcalorimetry, and molecular dynamics simulations

Specific antibodies interfere with the function of human tumor-associated carbonic anhydrase IX (CA IX), and show potential as tools for anticancer interventions. In this work, a correlation between structural elements and thermodynamic parameters of the association of antibody fragment Fab M75 to a peptide corresponding to its epitope in the proteoglycan-like domain of CA IX, is presented. Comparisons of the crystal structures of free Fab M75 and its complex with the epitope peptide reveal major readjustments of CDR-H1 and CDR-H3. In contrast, the overall conformations and positions of CDR-H2 and CDR-L2 remain unaltered, and their positively charged residues may thus present a fixed frame for epitope recognition. Adoption of the altered CDR-H3 conformation in the structure of the complex is accompanied by an apparent local stabilization. Analysis of domain mobility with translation-libration-screw (TLS) method shows that librations of the entire heavy chain variable domain (V(H)) decrease and reorient in the complex, which correlates well with participation of the heavy chain in ligand binding. Isothermal titration microcalorimetry (ITC) experiments revealed a highly unfavorable entropy term, which can be attributed mainly to the decrease in the degrees of freedom of the system, the loss of conformational freedom of peptide and partially to a local stabilization of CDR-H3. Moreover, it was observed that one proton is transferred from the environment to the protein-ligand complex upon binding. Molecular dynamics simulations followed by molecular mechanics/generalized Born surface area (MM-GBSA) calculations of the ligand (epitope peptide) binding energy yielded energy values that were in agreement with the ITC measurements and indicated that the charged residues play crucial role in the epitope binding. Theoretical arguments presented in this work indicate that two adjacent arginine residues (ArgH50 and ArgH52) are responsible for the observed proton transfer.     

High-resolution crystal structure of aldose reductase complexed with the novel sulfonyl-pyridazinone inhibitor exhibiting an alternative active site anchoring group

The crystal structure of a novel sulfonyl-pyridazinone inhibitor in complex with aldose reductase, the first enzyme of the polyol pathway, has been determined to 1.43 angstroms and 0.95 angstroms resolution. The ternary complex of inhibitor, cofactor and enzyme has been obtained by soaking of preformed crystals. Supposedly due to low solubility in the crystallisation buffer, in both structures the inhibitor shows reduced occupancy of 74% and 46% population, respectively. The pyridazinone head group of the inhibitor occupies the catalytic site, whereas the chloro-benzofuran moiety penetrates into the opened specificity pocket. The high-resolution structure provides some evidence that the pyridazinone group binds in a negatively charged deprotonated state, whereas the neighbouring His110 residue most likely adopts a neutral uncharged status. Since the latter structure is populated by the ligand to only 46%, a second conformation of the C-terminal ligand-binding region can be detected. This conformation corresponds to the closed state of the specificity pocket when no or only small ligands are bound to aldose reductase. The two conformational states are in good agreement with frames observed along a molecular dynamics trajectory describing the transition from closed to open situation. Accordingly, both geometries, superimposed in the averaged crystal structure, correspond to snapshots of the ligand-bound and the unbound state. Isothermal titration calorimetry has been applied to determine the binding constants of the investigated pyridazinone in comparison to the hydantoin sorbinil and the carboxylate-type inhibitors IDD 594 and tolrestat. The pyridazinone exhibits a binding affinity similar to those of tolrestat and sorbinil, and shows slightly reduced affinity compared to IDD 594. These studies elucidating the binding mode and providing information about protonation states of protein side-chains involved in binding of this novel class of inhibitors establish the platform for further structure-based drug design.     

Docking unbound proteins using shape complementarity, desolvation, and electrostatics

A comprehensive docking study was performed on 27 distinct protein-protein complexes. For 13 test systems, docking was performed with the unbound X-ray structures of both the receptor and the ligand. For the remaining systems, the unbound X-ray structure of only molecule was available; therefore the bound structure for the other molecule was used. Our method optimizes desolvation, shape complementarity, and electrostatics using a Fast Fourier Transform algorithm. A global search in the rotational and translational space without any knowledge of the binding sites was performed for all proteins except nine antibodies recognizing antigens. For these antibodies, we docked their well-characterized binding site-the complementarity-determining region defined without information of the antigen-to the entire surface of the antigen. For 24 systems, we were able to find near-native ligand orientations (interface C(alpha) root mean square deviation less than 2.5 A from the crystal complex) among the top 2,000 choices. For three systems, our algorithm could identify the correct complex structure unambiguously. For 13 other complexes, we either ranked a near-native structure in the top 20 or obtained 20 or more near-native structures in the top 2,000 or both. The key feature of our algorithm is the use of target functions that are highly tolerant to conformational changes upon binding. If combined with a post-processing method, our algorithm may provide a general solution to the unbound docking problem. Our program, called ZDOCK, is freely available to academic users (http://zlab.bu.edu/~rong/dock/).     

Flexible ligand docking using conformational ensembles.

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.     

Monoclonal antibodies to cruciform DNA structures

Two monoclonal antibodies, 2D3 and 4B4, have been raised against a cruciform structure in a heteroduplex DNA molecule. Antibody binding to DNA fragments was determined by a radioimmunoassay in which DNA--antibody complexes were separated from unbound DNA by acrylamide gel electrophoresis. These antibodies seem to recognize conformational determinants specific to cruciform structures. 2D3 and 4B4 antibodies do not bind to linear double-stranded homoduplex DNA fragments, linear single-stranded DNA or single-stranded simian virus 40 DNA containing a stem--loop structure, but do bind to the original cruciform and to a different cruciform with one shortened arm. 2D3 also bound to a T-shaped double-stranded DNA molecule, while 4B4 binding to this structure was weak. The monoclonal antibodies 2D3 and 4B4 were found to be immunoglobulin G1 and immunoglobulin M, respectively.