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.

Molecular basis for the recognition of two structurally different major histocompatibility complex/peptide complexes by a single T-cell receptor

2C is a typical alloreactive cytotoxic T lymphocyte clone that recognizes two different ligands. These ligands are adducts of the allo-major histocompatibility complex (MHC) molecule H-2Ld and an endogenous octapeptide, and of the self-MHC molecule H-2Kb and another peptide. MHC-binding and T-cell assays with synthetic peptides in combination with molecular modeling studies were employed to analyze the structural basis for this crossreactivity. The molecular surfaces of the two complexes differ greatly in densities and distributions of positive and negative charges. However, modifications of the peptides that increase similarity decrease the capacities of the resulting MHC peptide complexes to induce T-cell responses. Moreover, the roles of the peptides in ligand recognition are different for self- and allo-MHC-restricted T-cell responses. The self-MHC-restricted T-cell responses were finely tuned to recognition of the peptide. The allo-MHC-restricted responses, on the other hand, largely ignore modifications of the peptide. The results strongly suggest that adaptation of the T-cell receptor to the different ligand structures, rather than molecular mimicry by the ligands, is the basis for the crossreactivity of 2C. This conclusion has important implications for T-cell immunology and for the understanding of immunological disorders.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Attenuated T2 relaxation by mutual cancellation of dipole–dipole coupling and chemical shift anisotropy indicates an avenue to NMR structures of very large biological macromolecules in solution

Fast transverse relaxation of 1H, 15N, and 13C by dipole-dipole coupling (DD) and chemical shift anisotropy (CSA) modulated by rotational molecular motions has a dominant impact on the size limit for biomacromolecular structures that can be studied by NMR spectroscopy in solution. Transverse relaxation-optimized spectroscopy (TROSY) is an approach for suppression of transverse relaxation in multidimensional NMR experiments, which is based on constructive use of interference between DD coupling and CSA. For example, a TROSY-type two-dimensional 1H,15N-correlation experiment with a uniformly 15N-labeled protein in a DNA complex of molecular mass 17 kDa at a 1H frequency of 750 MHz showed that 15N relaxation during 15N chemical shift evolution and 1HN relaxation during signal acquisition both are significantly reduced by mutual compensation of the DD and CSA interactions. The reduction of the linewidths when compared with a conventional two-dimensional 1H,15N-correlation experiment was 60% and 40%, respectively, and the residual linewidths were 5 Hz for 15N and 15 Hz for 1HN at 4°C. Because the ratio of the DD and CSA relaxation rates is nearly independent of the molecular size, a similar percentagewise reduction of the overall transverse relaxation rates is expected for larger proteins. For a 15N-labeled protein of 150 kDa at 750 MHz and 20°C one predicts residual linewidths of 10 Hz for 15N and 45 Hz for 1HN, and for the corresponding uniformly 15N,2H-labeled protein the residual linewidths are predicted to be smaller than 5 Hz and 15 Hz, respectively. The TROSY principle should benefit a variety of multidimensional solution NMR experiments, especially with future use of yet somewhat higher polarizing magnetic fields than are presently available, and thus largely eliminate one of the key factors that limit work with larger molecules.

Molecular Cloning and Characterization of a Radial Spoke Head Protein of Sea Urchin Sperm Axonemes: Involvement of the Protein in the Regulation of Sperm Motility

Monoclonal antibodies raised against axonemal proteins of sea urchin spermatozoa have been used to study regulatory mechanisms involved in flagellar motility. Here, we report that one of these antibodies, monoclonal antibody D-316, has an unusual perturbating effect on the motility of sea urchin sperm models; it does not affect the beat frequency, the amplitude of beating or the percentage of motile sperm models, but instead promotes a marked transformation of the flagellar beating pattern which changes from a two-dimensional to a three-dimensional type of movement. On immunoblots of axonemal proteins separated by SDS-PAGE, D-316 recognized a single polypeptide of 90 kDa. This protein was purified following its extraction by exposure of axonemes to a brief heat treatment at 40°C. The protein copurified and coimmunoprecipitated with proteins of 43 and 34 kDa, suggesting that it exists as a complex in its native form. Using D-316 as a probe, a full-length cDNA clone encoding the 90-kDa protein was obtained from a sea urchin cDNA library. The sequence predicts a highly acidic (pI = 4.0) protein of 552 amino acids with a mass of 62,720 Da (p63). Comparison with protein sequences in databases indicated that the protein is related to radial spoke proteins 4 and 6 (RSP4 and RSP6) of Chlamydomonas reinhardtii, which share 37% and 25% similarity, respectively, with p63. However, the sea urchin protein possesses structural features distinct from RSP4 and RSP6, such as the presence of three major acidic stretches which contains 25, 17, and 12 aspartate and glutamate residues of 34-, 22-, and 14-amino acid long stretches, respectively, that are predicted to form α-helical coiled-coil secondary structures. These results suggest a major role for p63 in the maintenance of a planar form of sperm flagellar beating and provide new tools to study the function of radial spoke heads in more evolved species.

Kinetic equilibrium of forces and molecular events in muscle contraction

In biomolecular systems, the mechanical transfer of free energy occurs with both high efficiency and high speed. It is shown here that such a transfer can be achieved only if the participating free-energy-storing elements exhibit opposing relationships between their content of free energy and the force they exert in the transfer direction. A kinetic equilibrium of forces (KEF) results, in which the transfer of free energy is mediated essentially by thermal molecular motion. On the basis of present evidence, KEF is used as a guiding principle in developing a mechanical model of the crossbridge cycle in muscle contraction. The model allows the basic features of molecular events to be visualized in terms of plausible structures. Real understanding of the process will require identification of the elements that perform the functions described here. Besides chemomechanical energy transduction, KEF may have a role in other biomolecular processes in which free energy is transferred mechanically over large distances.

NMR structures of three single-residue variants of the human prion protein

The NMR structures of three single-amino acid variants of the C-terminal domain of the human prion protein, hPrP(121–230), are presented. In hPrP(M166V) and hPrP(R220K) the substitution is with the corresponding residue in murine PrP, and in hPrP(S170N) it is with the corresponding Syrian hamster residue. All three substitutions are in the surface region of the structure of the cellular form of PrP (PrPC) that is formed by the C-terminal part of helix 3, with residues 218–230, and a loop of residues 166–172. This molecular region shows high species variability and has been implicated in specific interactions with a so far not further characterized “protein X,” and it is related to the species barrier for transmission of prion diseases. As expected, the three variant hPrP(121–230) structures have the same global architecture as the previously determined wild-type bovine, human, murine, and Syrian hamster prion proteins, but with the present study two localized “conformational markers” could be related with single amino acid exchanges. These are the length and quality of definition of helix 3, and the NMR-observability of the residues in the loop 166–172. Poor definition of the C-terminal part of helix 3 is characteristic for murine PrP and has now been observed also for hPrP(R220K), and NMR observation of the complete loop 166–172 has so far been unique for Syrian hamster PrP and is now also documented for hPrP(S170N).

Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.

Crystal structures of eukaryotic translation initiation factor 5A from Methanococcus jannaschii at 1.8 Å resolution

Eukaryotic translation initiation factor 5A (eIF-5A) is a ubiquitous protein found in all eukaryotic cells. The protein is closely associated with cell proliferation in the G1–S stage of the cell cycle. Recent findings show that the eIF-5A proteins are highly expressed in tumor cells and act as a cofactor of the Rev protein in HIV-1-infected cells. The mature eIF is the only protein known to have the unusual amino acid hypusine, a post-translationally modified lysine. The crystal structure of eIF-5A from Methanococcus jannaschii (MJ eIF-5A) has been determined at 1.9 Å and 1.8 Å resolution in two crystal forms by using the multiple isomorphous replacement method and the multiwavelength anomalous diffraction method for the first crystal form and the molecular replacement method for the second crystal form. The structure consists of two folding domains, one of which is similar to the oligonucleotide-binding domain found in the prokaryotic cold shock protein and the translation initiation factor IF1 despite the absence of any significant sequence similarities. The 12 highly conserved amino acid residues found among eIF-5As include the hypusine site and form a long protruding loop at one end of the elongated molecule.

Synthetic zeolites and other microporous oxide molecular sieves

Use of synthetic zeolites and other microporous oxides since 1950 has improved insulated windows, automobile air-conditioning, refrigerators, air brakes on trucks, laundry detergents, etc. Their large internal pore volumes, molecular-size pores, regularity of crystal structures, and the diverse framework chemical compositions allow “tailoring” of structure and properties. Thus, highly active and selective catalysts as well as adsorbents and ion exchangers with high capacities and selectivities were developed. In the petroleum refining and petrochemical industries, zeolites have made possible cheaper and lead-free gasoline, higher performance and lower-cost synthetic fibers and plastics, and many improvements in process efficiency and quality and in performance. Zeolites also help protect the environment by improving energy efficiency, reducing automobile exhaust and other emissions, cleaning up hazardous wastes (including the Three Mile Island nuclear power plant and other radioactive wastes), and, as specially tailored desiccants, facilitating the substitution of new refrigerants for the ozone-depleting chlorofluorocarbons banned by the Montreal Protocol.

An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

A Chloroplast DNA Helicase II from Pea That Prefers Fork-Like Replication Structures

A DNA helicase, called chloroplast DNA (ctDNA) helicase II, was purified to apparent homogeneity from pea (Pisum sativum). The enzyme contained intrinsic, single-stranded, DNA-dependent ATPase activity and an apparent molecular mass of 78 kD on sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The DNA helicase was markedly stimulated by DNA substrates with fork-like replication structures. A 5′-tailed fork was more active than the 3′-tailed fork, which itself was more active than substrates without a fork. The direction of unwinding was 3′ to 5′ along the bound strand, and it failed to unwind blunt-ended duplex DNA. DNA helicase activity required only ATP or dATP hydrolysis. The enzyme also required a divalent cation (Mg2+>Mn2+>Ca2+) for its unwinding activity and was inhibited at 200 mm KCl or NaCl. This enzyme could be involved in the replication of ctDNA. The DNA major groove-intercalating ligands nogalamycin and daunorubicin were inhibitory to unwinding (Ki approximately 0.85 μm and 2.2 μm, respectively) and ATPase (Ki approximately 1.3 μm and 3.0 μm, respectively) activities of pea ctDNA helicase II, whereas ellipticine, etoposide (VP-16), and camptothecin had no effect on the enzyme activity. These ligands may be useful in further studies of the mechanisms of chloroplast helicase activities.

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

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.