Ultrafast diffraction and structural dynamics: The nature of complex molecules far from equilibrium

Studies of molecular structures at or near their equilibrium configurations have long provided information on their geometry in terms of bond distances and angles. Far-from-equilibrium structures are relatively unknown—especially for complex systems—and generally, neither their dynamics nor their average geometries can be extrapolated from equilibrium values. For such nonequilibrium structures, vibrational amplitudes and bond distances play a central role in phenomena such as energy redistribution and chemical reactivity. Ultrafast electron diffraction, which was developed to study transient molecular structures, provides a direct method for probing the nature of complex molecules far from equilibrium. Here we present our ultrafast electron diffraction observations of transient structures for two cyclic hydrocarbons. At high internal energies of ≈4 eV, these molecules display markedly different behavior. For 1,3,5-cycloheptatriene, excitation results in the formation of hot ground-state structures with bond distances similar to those of the initial structure, but with nearly three times the average vibrational amplitude. Energy is redistributed within 5 ps, but with a negative temperature characterizing the nonequilibrium population. In contrast, the ring-opening reaction of 1,3-cyclohexadiene is shown to result in hot structures with a C—C bond distance of over 1.7 Å, which is 0.2 Å away from any expected equilibrium value. Even up to 400 ps, energy remains trapped in large-amplitude motions comprised of torsion and asymmetric stretching. These studies promise a new direction for studying structural dynamics in nonequilibrium complex systems.

Antiviral agent blocks breathing of the common cold virus

A dynamic capsid is critical to the events that shape the viral life cycle; events such as cell attachment, cell entry, and nucleic acid release demand a highly mobile viral surface. Protein mass mapping of the common cold virus, human rhinovirus 14 (HRV14), revealed both viral structural dynamics and the inhibition of such dynamics with an antiviral agent, WIN 52084. Viral capsid digestion fragments resulting from proteolytic time-course experiments provided structural information in good agreement with the HRV14 three-dimensional crystal structure. As expected, initial digestion fragments included peptides from the capsid protein VP1. This observation was expected because VP1 is the most external viral protein. Initial digestion fragments also included peptides belonging to VP4, the most internal capsid protein. The mass spectral results together with x-ray crystallography data provide information consistent with a “breathing” model of the viral capsid. Whereas the crystal structure of HRV14 shows VP4 to be the most internal capsid protein, mass spectral results show VP4 fragments to be among the first digestion fragments observed. Taken together this information demonstrates that VP4 is transiently exposed to the viral surface via viral breathing. Comparative digests of HRV14 in the presence and absence of WIN 52084 revealed a dramatic inhibition of digestion. These results indicate that the binding of the antiviral agent not only causes local conformational changes in the drug binding pocket but actually stabilizes the entire viral capsid against enzymatic degradation. Viral capsid mass mapping provides a fast and sensitive method for probing viral structural dynamics as well as providing a means for investigating antiviral drug efficacy.

Specific interaction of mutant p53 with regions of matrix attachment region DNA elements (MARs) with a high potential for base-unpairing

Mutant, but not wild-type p53 binds with high affinity to a variety of MAR-DNA elements (MARs), suggesting that MAR-binding of mutant p53 relates to the dominant-oncogenic activities proposed for mutant p53. MARs recognized by mutant p53 share AT richness and contain variations of an AATATATTT “DNA-unwinding motif,” which enhances the structural dynamics of chromatin and promotes regional DNA base-unpairing. Mutant p53 specifically interacted with MAR-derived oligonucleotides carrying such unwinding motifs, catalyzing DNA strand separation when this motif was located within a structurally labile sequence environment. Addition of GC-clamps to the respective MAR-oligonucleotides or introducing mutations into the unwinding motif strongly reduced DNA strand separation, but supported the formation of tight complexes between mutant p53 and such oligonucleotides. We conclude that the specific interaction of mutant p53 with regions of MAR-DNA with a high potential for base-unpairing provides the basis for the high-affinity binding of mutant p53 to MAR-DNA.

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.

Identification of two distinct human SMC protein complexes involved in mitotic chromosome dynamics

The structural maintenance of chromosomes (SMC) family member proteins previously were shown to play a critical role in mitotic chromosome condensation and segregation in yeast and Xenopus. Other family members were demonstrated to be required for DNA repair in yeast and mammals. Although several different SMC proteins were identified in different organisms, little is known about the SMC proteins in humans. Here, we report the identification of four human SMC proteins that form two distinct heterodimeric complexes in the cell, the human chromosome-associated protein (hCAP)-C and hCAP-E protein complex (hCAP-C/hCAP-E), and the human SMC1 (hSMC1) and hSMC3 protein complex (hSMC1/hSMC3). The hCAP-C/hCAP-E complex is the human ortholog of the Xenopus chromosome-associated protein (XCAP)-C/XCAP-E complex required for mitotic chromosome condensation. We found that a second complex, hSMC1/hSMC3, is required for metaphase progression in mitotic cells. Punctate vs. diffuse distribution patterns of the hCAP-C/hCAP-E and hSMC1/hSMC3 complexes in the interphase nucleus indicate independent behaviors of the two complexes during the cell cycle. These results suggest that two distinct classes of SMC protein complexes are involved in different aspects of mitotic chromosome organization in human cells.

Dynamics and orientation of N+(CD3)3-bromoacetylcholine bound to its binding site on the nicotinic acetylcholine receptor

Dynamic and structural information has been obtained for an analogue of acetylcholine while bound to the agonist binding site on the nicotinic acetylcholine receptor (nAcChoR), using wide-line deuterium solid-state NMR. Analysis of the deuterium lineshape obtained at various temperatures from unoriented nAcChoR membranes labeled with deuterated bromoacetylcholine (BAC) showed that the quaternary ammonium group of the ligand is well constrained within the agonist binding site when compared with the dynamics observed in the crystalline solids. This motional restriction would suggest that a high degree of complementarity exists between the quaternary ammonium group of the ligand and the protein within the agonist binding site. nAcChoR membranes were uniaxially oriented by isopotential centrifugation as determined by phosphorous NMR of the membrane phospholipids. Analysis of the deuterium NMR lineshape of these oriented membranes enriched with the nAcChoR labeled with N+(CD3)3-BAC has enabled us to determine that the angle formed between the quaternary ammonium group of the BAC and the membrane normal is 42° in the desensitized form of the receptor. This measurement allows us to orient in part the bound ligand within the proposed receptor binding site.

The role of structure, energy landscape, dynamics, and allostery in the enzymatic function of myoglobin

The grail of protein science is the connection between structure and function. For myoglobin (Mb) this goal is close. Described as only a passive dioxygen storage protein in texts, we argue here that Mb is actually an allosteric enzyme that can catalyze reactions among small molecules. Studies of the structural, spectroscopic, and kinetic properties of Mb lead to a model that relates structure, energy landscape, dynamics, and function. Mb functions as a miniature chemical reactor, concentrating and orienting diatomic molecules such as NO, CO, O2, and H2O2 in highly conserved internal cavities. Reactions can be controlled because Mb exists in distinct taxonomic substates with different catalytic properties and connectivities of internal cavities.

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.

Mantle dynamics and seismic tomography

Three-dimensional imaging of the Earth's interior, called seismic tomography, has achieved breakthrough advances in the last two decades, revealing fundamental geodynamical processes throughout the Earth's mantle and core. Convective circulation of the entire mantle is taking place, with subducted oceanic lithosphere sinking into the lower mantle, overcoming the resistance to penetration provided by the phase boundary near 650-km depth that separates the upper and lower mantle. The boundary layer at the base of the mantle has been revealed to have complex structure, involving local stratification, extensive structural anisotropy, and massive regions of partial melt. The Earth's high Rayleigh number convective regime now is recognized to be much more interesting and complex than suggested by textbook cartoons, and continued advances in seismic tomography, geodynamical modeling, and high-pressure–high-temperature mineral physics will be needed to fully quantify the complex dynamics of our planet's interior.

Symmetry, stability, and dynamics of multidomain and multicomponent protein systems

Symmetry is commonly observed in many biological systems. Here we discuss representative examples of the role of symmetry in structural molecular biology. Point group symmetries are observed in many protein oligomers whose three-dimensional atomic structures have been elucidated by x-ray crystallography. Approximate symmetry also occurs in multidomain proteins. Symmetry often confers stability on the molecular system and results in economical usage of basic components to build the macromolecular structure. Symmetry is also associated with cooperativity. Mild perturbation from perfect symmetry may be essential in some systems for dynamic functions.

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.

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.

A structural basis for a phosphoramide mustard-induced DNA interstrand cross-link at 5'-d(GAC).

Phosphoramide mustard-induced DNA interstrand cross-links were studied both in vitro and by computer simulation. The local determinants for the formation of phosphoramide mustard-induced DNA interstrand cross-links were defined by using different pairs of synthetic oligonucleotide duplexes, each of which contained a single potentially cross-linkable site. Phosphoramide mustard was found to cross-link dG to dG at a 5'-d(GAC)-3'. The structural basis for the formation of this 1,3 cross-link was studied by molecular dynamics and quantum chemistry. Molecular dynamics indicated that the geometrical proximity of the binding sites also favored a 1,3 dG-to-dG linkage over a 1,2 dG-to-dG linkage in a 5'-d(GCC)-3' sequence. While the enthalpies of 1,2 and 1,3 mustard cross-linked DNA were found to be very close, a 1,3 structure was more flexible and may therefore be in a considerably higher entropic state.