Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

Vascular endothelial growth factor: Crystal structure and functional mapping of the kinase domain receptor binding site

Vascular endothelial growth factor (VEGF) is a homodimeric member of the cystine knot family of growth factors, with limited sequence homology to platelet-derived growth factor (PDGF) and transforming growth factor β2 (TGF-β). We have determined its crystal structure at a resolution of 2.5 Å, and identified its kinase domain receptor (KDR) binding site using mutational analysis. Overall, the VEGF monomer resembles that of PDGF, but its N-terminal segment is helical rather than extended. The dimerization mode of VEGF is similar to that of PDGF and very different from that of TGF-β. Mutational analysis of VEGF reveals that symmetrical binding sites for KDR are located at each pole of the VEGF homodimer. Each site contains two functional “hot spots” composed of binding determinants presented across the subunit interface. The two most important determinants are located within the largest hot spot on a short, three-stranded sheet that is conserved in PDGF and TGF-β. Functional analysis of the binding epitopes for two receptor-blocking antibodies reveal different binding determinants near each of the KDR binding hot spots.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Excited-state dynamics in photosystem II: Insights from the x-ray crystal structure

The heart of oxygenic photosynthesis is photosystem II (PSII), a multisubunit protein complex that uses solar energy to drive the splitting of water and production of molecular oxygen. The effectiveness of the photochemical reaction center of PSII depends on the efficient transfer of excitation energy from the surrounding antenna chlorophylls. A kinetic model for PSII, based on the x-ray crystal structure coordinates of 37 antenna and reaction center pigment molecules, allows us to map the major energy transfer routes from the antenna chlorophylls to the reaction center chromophores. The model shows that energy transfer to the reaction center is slow compared with the rate of primary electron transport and depends on a few bridging chlorophyll molecules. This unexpected energetic isolation of the reaction center in PSII is similar to that found in the bacterial photosystem, conflicts with the established view of the photophysics of PSII, and may be a functional requirement for primary photochemistry in photosynthesis. In addition, the model predicts a value for the intrinsic photochemical rate constant that is 4 times that found in bacterial reaction centers.

A cooperative model for receptor recognition and cell adhesion: evidence from the molecular packing in the 1.6-A crystal structure of the pheromone Er-1 from the ciliated protozoan Euplotes raikovi.

The crystal structure of the pheromone Er-1 from the unicellular eukaryotic organism Euplotes raikovi was determined at 1.6 A resolution and refined to a crystallographic R factor of 19.9%. In the tightly packed crystal, two extensive intermolecular helix-helix interactions arrange the Er-1 molecules into layers. Since the putative receptor of the pheromone is a membrane-bound protein, whose extracellular C-terminal domain is identical in amino acid sequence to the soluble pheromone, the interactions found in the crystal may mimic the pheromone-receptor interactions as they occur on a cell surface. Based on this, we propose a model for the interaction between soluble pheromone molecules and their receptors. In this model, strong pheromone-receptor binding emerges as a consequence of the cooperative utilization of several weak interactions. The model offers an explanation for the results of binding studies and may also explain the adhesion between cells that occurs during mating.

Crystal and molecular structure of paclitaxel (taxol).

Paclitaxel (formerly called taxol), an important anticancer drug, inhibits cell replication by binding to and stabilizing microtubule polymers. As drug-receptor interactions are governed by the three-dimensional stereochemistries of both participants, we have determined the crystal structure of paclitaxel to identify its conformational preferences that may be related to biological activity. The monoclinic crystals contain two independent paclitaxel molecules in the asymmetric unit plus several water and dioxane solvent molecules. Taxane ring conformation is very similar in both paclitaxel molecules and is similar to the taxane ring conformation found in the crystal structure of the paclitaxel analogue docetaxel (formerly called taxotere). The two paclitaxel molecules have carbon-13 side-chain conformations that differ from each other and from that of the corresponding side chain in the docetaxel crystal structure. The carbon-13 side-chain conformation of one paclitaxel molecule is similar to what was proposed from NMR studies done in polar solvents, while that of the other paclitaxel molecule is different and hitherto unobserved. The paclitaxel molecules interact with each other and with solvent atoms through an extensive network of hydrogen bonds. Analysis of the hydrogen-bonding network together with structure-activity studies may suggest which atoms of paclitaxel are important for binding to microtubule receptors.

Structural homologies with ATP- and folate-binding enzymes in the crystal structure of folylpolyglutamate synthetase

Folylpolyglutamate synthetase, which is responsible for the addition of a polyglutamate tail to folate and folate derivatives, is an ATP-dependent enzyme isolated from eukaryotic and bacterial sources, where it plays a key role in the retention of the intracellular folate pool. Here, we report the 2.4-Å resolution crystal structure of the MgATP complex of the enzyme from Lactobacillus casei. The structural analysis reveals that folylpolyglutamate synthetase is a modular protein consisting of two domains, one with a typical mononucleotide-binding fold and the other strikingly similar to the folate-binding enzyme dihydrofolate reductase. We have located the active site of the enzyme in a large interdomain cleft adjacent to an ATP-binding P-loop motif. Opposite this site, in the C domain, a cavity likely to be the folate binding site has been identified, and inspection of this cavity and the surrounding protein structure suggests that the glutamate tail of the substrate may project into the active site. A further feature of the structure is a well defined Ω loop, which contributes both to the active site and to interdomain interactions. The determination of the structure of this enzyme represents the first step toward the elucidation of the molecular mechanism of polyglutamylation of folates and antifolates.

Crystal structure of Taq DNA polymerase in complex with an inhibitory Fab: The Fab is directed against an intermediate in the helix-coil dynamics of the enzyme

We report the crystal structure of Thermus aquaticus DNA polymerase I in complex with an inhibitory Fab, TP7, directed against the native enzyme. Some of the residues present in a helical conformation in the native enzyme have adopted a γ turn conformation in the complex. Taken together, structural information that describes alteration of helical structure and solution studies that demonstrate the ability of TP7 to inhibit 100% of the polymerase activity of the enzyme suggest that the change in conformation is probably caused by trapping of an intermediate in the helix-coil dynamics of this helix by the Fab. Antibodies directed against modified helices in proteins have long been anticipated. The present structure provides direct crystallographic evidence. The Fab binds within the DNA binding cleft of the polymerase domain, interacting with several residues that are used by the enzyme in binding the primer:template complex. This result unequivocally corroborates inferences drawn from binding experiments and modeling calculations that the inhibitory activity of this Fab is directly attributable to its interference with DNA binding by the polymerase domain of the enzyme. The combination of interactions made by the Fab residues in both the polymerase and the vestigial editing nuclease domain of the enzyme reveal the structural basis of its preference for binding to DNA polymerases of the Thermus species. The orientation of the structure-specific nuclease domain with respect to the polymerase domain is significantly different from that seen in other structures of this polymerase. This reorientation does not appear to be antibody-induced and implies remarkably high relative mobility between these two domains.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

Crystal structure of Bacillus subtilis YabJ, a purine regulatory protein and member of the highly conserved YjgF family

The yabJ gene in Bacillus subtilis is required for adenine-mediated repression of purine biosynthetic genes in vivo and codes for an acid-soluble, 14-kDa protein. The molecular mechanism of YabJ is unknown. YabJ is a member of a large, widely distributed family of proteins of unknown biochemical function. The 1.7-Å crystal structure of YabJ reveals a trimeric organization with extensive buried hydrophobic surface and an internal water-filled cavity. The most important finding in the structure is a deep, narrow cleft between subunits lined with nine side chains that are invariant among the 25 most similar homologs. This conserved site is proposed to be a binding or catalytic site for a ligand or substrate that is common to YabJ and other members of the YER057c/YjgF/UK114 family of proteins.

The crystal structure of the Rev binding element of HIV-1 reveals novel base pairing and conformational variability

The crystal and molecular structure of an RNA duplex corresponding to the high affinity Rev protein binding element (RBE) has been determined at 2.1-Å resolution. Four unique duplexes are present in the crystal, comprising two structural variants. In each duplex, the RNA double helix consists of an annealed 12-mer and 14-mer that form an asymmetric internal loop consisting of G-G and G-A noncanonical base pairs and a flipped-out uridine. The 12-mer strand has an A-form conformation, whereas the 14-mer strand is distorted to accommodate the bulges and noncanonical base pairing. In contrast to the NMR model of the unbound RBE, an asymmetric G-G pair with N2-N7 and N1-O6 hydrogen bonding, is formed in each helix. The G-A base pairing agrees with the NMR structure in one structural variant, but forms a novel water-mediated pair in the other. A backbone flip and reorientation of the G-G base pair is required to assume the RBE conformation present in the NMR model of the complex between the RBE and the Rev peptide.

Crystal structure of a Staphylococcus aureus protein A domain complexed with the Fab fragment of a human IgM antibody: Structural basis for recognition of B-cell receptors and superantigen activity

Staphylococcus aureus produces a virulence factor, protein A (SpA), that contains five homologous Ig-binding domains. The interactions of SpA with the Fab region of membrane-anchored Igs can stimulate a large fraction of B cells, contributing to lymphocyte clonal selection. To understand the molecular basis for this activity, we have solved the crystal structure of the complex between domain D of SpA and the Fab fragment of a human IgM antibody to 2.7-Å resolution. In the complex, helices II and III of domain D interact with the variable region of the Fab heavy chain (VH) through framework residues, without the involvement of the hypervariable regions implicated in antigen recognition. The contact residues are highly conserved in human VH3 antibodies but not in other families. The contact residues from domain D also are conserved among all SpA Ig-binding domains, suggesting that each could bind in a similar manner. Features of this interaction parallel those reported for staphylococcal enterotoxins that are superantigens for many T cells. The structural homology between Ig VH regions and the T-cell receptor Vβ regions facilitates their comparison, and both types of interactions involve lymphocyte receptor surface remote from the antigen binding site. However, T-cell superantigens reportedly interact through hydrogen bonds with T-cell receptor Vβ backbone atoms in a primary sequence-independent manner, whereas SpA relies on a sequence-restricted conformational binding with residue side chains, suggesting that this common bacterial pathogen has adopted distinct molecular recognition strategies for affecting large sets of B and T lymphocytes.

Crystal structure of an anti-carbohydrate antibody directed against Vibrio cholerae O1 in complex with antigen: Molecular basis for serotype specificity

The crystal structure of the murine Fab S-20-4 from a protective anti-cholera Ab specific for the lipopolysaccharide Ag of the Ogawa serotype has been determined in its unliganded form and in complex with synthetic fragments of the Ogawa O-specific polysaccharide (O-SP). The upstream terminal O-SP monosaccharide is shown to be the primary antigenic determinant. Additional perosamine residues protrude outwards from the Ab surface and contribute only marginally to the binding affinity and specificity. A complementary water-excluding hydrophobic interface and five Ab–Ag hydrogen bonds are crucial for carbohydrate recognition. The structure reported here explains the serotype specificity of anti-Ogawa Abs and provides a rational basis toward the development of a synthetic carbohydrate-based anti-cholera vaccine.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Crystal structure of the T state of allosteric yeast chorismate mutase and comparison with the R state.

The crystal structure of the tyrosine-bound T state of allosteric yeast Saccharomyces cerevisiae chorismate mutase was solved by molecular replacement at a resolution of 2.8 angstroms using a monomer of the R-state structure as the search model. The allosteric inhibitor tyrosine was found to bind in the T state at the same binding site as the allosteric activator tryptophan binds in the R state, thus defining one regulatory binding site for each monomer. Activation by tryptophan is caused by the larger steric size of its side chain, thereby pushing apart the allosteric domain of one monomer and helix H8 of the catalytic domain of the other monomer. Inhibition is caused by polar contacts of tyrosine with Arg-75 and Arg-76 of one monomer and with Gly-141, Ser-142, and Thr-145 of the other monomer, thereby bringing the allosteric and catalytic domains closer together. The allosteric transition includes an 8 degree rotation of each of the two catalytic domains relative to the allosteric domains of each monomer (domain closure). Alternatively, this transition can be described as a 15 degree rotation of the catalytic domains of the dimer relative to each other.