Bond-optimized ring closure for proline: comparison of conformations and semiempirical energies with small molecule X-ray structures

A method is described for generating proline ring structures by successive addition of atoms, wherein ring closure is achieved by optimizing the fit to known ring bond-angles and one closing bond-length ("bond-optimized ring closure"). Two ring torsion angles are fixed independently within broad, allowed ranges, and the remaining torsion angles are determined uniquely in most cases. The independent torsion angles are chosen as phi and chi 2, and ring closure is achieved without prohibitive strain through most of the ranges -130 degrees less than phi less than -20 degrees and -60 degrees less than chi 2 less than 60 degrees. Comparisons of predicted ring structures to 191 X-ray diffraction structures from the literature, starting with the known values of phi and chi 2, yielded root-mean-square deviations of 4.8 degrees in chi 1, 4.7 degrees in chi 3, 8.3 degrees in chi 4, and 0.3-2% in the ring bond angles and the N-C delta distance. Semiempirical energies were calculated for the optimized structures using three sets of energy parameters from the literature. The energy surfaces show broad minima coinciding with the torsion angle regions in which the highest concentrations of observed structures are found. Two of the sets of energy parameters produce double minima corresponding to the "up" and "down" puckered conformations.     

Dynamic conformations of the CD38-mediated NAD cyclization captured in a single crystal

The extracellular domain of human CD38 is a multifunctional enzyme involved in the metabolism of two Ca²⁺ messengers: cyclic ADP-ribose and nicotinic acid adenine dinucleotide phosphate. When NAD is used as substrate, CD38 predominantly hydrolyzes it to ADP-ribose, with a trace amount of cyclic ADP-ribose produced through cyclization of the substrate. However, mutation of a key residue at the active site, E146, inhibits the hydrolysis activity of CD38 but greatly increases its cyclization activity. To understand the role of the residue E146 in the catalytic process, we determined the crystal structure of the E146A mutant protein with a substrate analogue, arabinosyl-2′-fluoro-deoxy-nicotinamide adenine dinucleotide. The structure captured the enzymatic reaction intermediates in six different conformations in a crystallographic asymmetric unit. The structural results indicate a folding-back process for the adenine ring of the substrate and provide the first multiple snapshots of the process. Our approach of utilizing multiple molecules in the crystallographic asymmetric unit should be generally applicable for capturing the dynamic nature of enzymatic catalysis.     

Synergy within structural biology of single crystal optical spectroscopy and X-ray crystallography

Advances in the adaptation of optical spectroscopy to monitor photo-induced or enzyme-catalyzed reactions in the crystalline state have enabled X-ray crystal structures to be accurately linked with spectroscopically defined intermediates. This, in turn, has led to a deeper understanding of the role protein structural changes play in function. The integration of optical spectroscopy with X-ray crystallography is growing and now extends beyond linking crystal structure to reaction intermediate. Recent examples of this synergy include applications in protein crystallization, X-ray data acquisition, radiation damage, and acquisition of phase information important for structure determination.     

Structure of tumour necrosis factor by X-ray solution scattering and preliminary studies by single crystal X-ray diffraction

The structure of tumour necrosis factor has been investigated by X-ray small-angle scattering and X-ray diffraction using synchrotron radiation. The overall radius of gyration is 25.5 A. A plausible model accounting for the scattering curves consists of an elongated trimer with an axial ratio of 3 to 4 and a maximal chord with a lower limit of 80 A. Tumour necrosis factor has been crystallized in a trigonal space group. Our results are in favour of a single trimer in the asymmetric unit. The diffraction extends to 3.5 A.     

Hominid and gelada baboon evolution: Agreement between molecular and fossil time scales

The time of origin of the hominid lineage has long been debated. Macromolecular studies have consistently shown genetic distances between living humans and African apes to be quite small. The molecular clock hypothesis proposes that the time of separation of these lineages is relatively recent (in the range of 4–8 million years ago) and not 15 million years or more ago as usually suggested. Three independent molecular comparisons yield a mean estimate of 4.6 million years for the hominid-African pongid divergence. The relationship of Theropithecusand Papiois a parallel case within Primates of two taxa which are quite similar at the molecular level, but which are usually thought to have separated relatively long ago. The two cases of seeming discordance between different lines of evidence are analogous. Each involves a speciation event which eventually resulted in one substantially derived lineage and one or more relatively unchanged lineages. In each case, claims of the antiquity of the divergence event extend to at least twice the age of the first certain appearance of the more derived lineage in the fossil record. Finally, in each case, the molecular clock model suggests a range of possible divergence times that overlaps with the first appearances of undoubted hominids and Theropithecusin the fossil record. This test involving paleontological evidence supports the molecular clock hypothesis.     

Development of a single detector ring micro crystal element scanner: QuickPET II

This article describes a single ring version of the micro crystal element scanner (MiCES) and investigation of its spatial resolution imaging characteristics for mouse positron emission tomography (PET) imaging. This single ring version of the MiCES system, referred to as QuickPET II, consists of 18 MiCE detector modules mounted as a single ring in a vertical gantry. The system has a 5.76-cm transverse field of view and a 1.98-cm axial field of view. In addition to the scanner and data acquisition system, we have developed an iterative reconstruction that includes a model of the system's detector response function. Evaluation images of line sources and mice have been acquired. Using filtered backprojection, the resolution for a reconstructed line source has been measured at 1.2 mm full width at half maximum. F-18-2-fluoro-2-deoxyglucose mouse PET images are provided. The result shows that QuickPET II has the imaging characteristics to support high-resolution, static mouse PET studies using 18-F labeled compounds.     

Correlation between glycosyl torsion angle and sugar ring pucker does not always exist.

A survey of the conformational parameters of deoxypyrimidine nucleosides and nucleotides shows that the correlation between glycosyl torsion angle and sugar pucker, which has often been considered to be well-established, does not always exist. This may be of significance when interactions between DNA and other molecules are considered.     

Quartz crystal biosensor for real-time monitoring of molecular recognition between protein and small molecular medicinal agents

A quartz crystal microbalance (QCM) biosensor integrated into a flow injection analysis (FIA) system was used for the real-time investigation of molecular recognition between a protein and small molecular medicinal agents. Two sulfa-drugs, sulfamethazine (SMZ) and sulfamethoxazole (SMO), were, respectively, immobilized on the gold electrodes of the piezoelectric crystals using appropriate procedures based on self-assembly of the dithiothreitol (DTT). The binding interactions of the two immobilized drug ligands, with various proteins in solution, were followed as changes in the resonant frequency of the modified crystals. Results obtained from this rapid screen analysis clearly indicated that the two drug ligands appeared quite different in this molecular recognition procedure although their structures were similar. SMZ-immobilized sensor showed specific interaction only with IgG, while SMO-immobilized sensor showed negligible specific binding with IgG, but binding with trypsin and chymotrypsin. Further studies on the specific interaction between immobilized SMZ and three different species of IgG--human IgG, goat IgG and mouse IgG were carried out and the marked species-dependent difference was observed. The resultant sensorgrams were rapidly analyzed by using an in-house kinetic analysis software based on genetic algorithm (GA) to derive both the kinetic rate constants (kass and kdiss) and equilibrium association constants (KA) for IgG-SMZ interactions. For the interactions, KA were 5.48 x 10(5), 2.75 x 10(5) and 1.86 x 10(5) M(-1) for human IgG, goat IgG and mouse IgG, respectively. The kinetic data provided further insight into the structural/functional relationships of different IgG on a molecular level.     

Cesium cholate: determination of X-ray crystal structure indicates participation of the ring hydroxyl groups in metal binding

The crystal structure of cesium cholate, C(24)H(36)(OH)(3) COOCs has been determined with three-dimensional X-ray diffractometer data. It crystallized in the monoclinic space group P2(1) with unit-cell dimensions a = 11.543(5) A, b = 8.614(3) A, and c = 12.662(5) A, beta(deg) = 107.95(2), V = 1197.7 A(3) and Z = 2. The atomic parameters were refined to a final r = 0.0269 and R(omega) = 0.0280 for 2342 observed reflections. Each Cs(+) is coordinated to 7 oxygen atoms from 5 different cholate anions with Cs-O distances ranging from 2.957(4) A to 3.678(5) A. In this crystal, 5 cholates are coordinated with 1 Cs(+), and 5 Cs(+) are coordinated with 1 cholate anion. Carboxyl and all the 3 ring hydroxyl groups of cholate anion participate in binding to Cs(+) simultaneously, and there is no water molecule coordinated with the Cs(+). The pattern of successive rows arranged with polar (p) and non-polar (n) faces in apposition leads to the formation of a sandwich sheet structure with polar and non-polar channels. The Cs ions lie within the polar interior of the sandwich. The H-bond network is reorganized in forming cesium cholate from cholic acid. All the oxygen atoms in cholate anion are involved in H-bonding reciprocally or with water molecules to form an extensive 3-dimensional network of H-bonds. Compared with cholic acid and other similar type of steroids, the coordination structure and H-bonding of Cs cholate crystal are distinct.     

Characterization of D-glucaric acid using NMR, X-ray crystal structure, and MM3 molecular modeling analyses

D-Glucaric acid was characterized in solution by comparing NMR spectra from the isotopically unlabeled molecule with those from D-glucaric acid labeled with deuterium or carbon-13 atoms. The NMR studies provided unequivocal assignments for all carbon atoms and non-hydroxyl protons of the molecule. The crystal structure of D-glucaric acid was obtained by X-ray diffraction techniques and the structure was a close match to the low energy conformation generated from a Monte-Carlo-based searching protocol employing the MM3 molecular mechanics program. The molecule adopts a bent structure in both the crystalline and computationally generated lowest-energy structure, a conformation that is devoid of destabilizing eclipsed 1,3-hydroxyl interactions.     

Structural features of human initiation factor 4E,studied by X-ray crystal analyses and molecular dynamics simulations

The structural features of human eIF4E were investigated by X-ray crystal analyses of its cap analog (m(7)GTP and m(7)GpppA) complexes and molecular dynamics (MD) simulations of cap-free and cap-bound eIF4Es, as well as the cap-bound Ser209-phosphorylated eIF4E. Crystal structure analyses at 2.0 A resolution revealed that the molecule forms a temple-bell-shaped surface of eight antiparallel beta-structures, three alpha-helices and ten loop structures, where the N-terminal region corresponds to the handle of the bell. This concave backbone provides a scaffold for the mRNA cap-recognition pocket consisting of three receiving parts for the 5'-terminal m(7)G base, the triphosphate, and the second nucleotide. The m(7)G base is sandwiched between the two aromatic side-chains of Trp102 and Trp56. The two (m(7)G)NH-O (Glu103 carboxy group) hydrogen bonds stabilize the stacking interaction. The basic residues of Arg157 and Lys162 and water molecules construct a binding pocket for the triphosphate moiety, where a universal hydrogen-bonding network is formed. The flexible C-terminal loop region unobserved in the m(7)GTP complex was clearly observed in the m(7)GpppA complex, as a result of the fixation of this loop by the interaction with the adenosine moiety, indicating the function of this loop as a receiving pocket for the second nucleotide. On the other hand, MD simulation in an aqueous solution system revealed that the cap-binding pocket, especially its C-terminal loop structure, is flexible in the cap-free eIF4E, and the entrance of the cap-binding pocket becomes narrow, although the depth is relatively unchanged. SDS-PAGE analyses showed that this structural instability is highly related to the fast degradation of cap-free eIF4E, compared with cap-bound or 4E-BP/cap-bound eIF4E, indicating the conferment of structural stability of eIF4E by the binary or ternary complex formation. MD simulation of m(7)GpppA-bound Ser209-phosphorylated eIF4E showed that the size of the cap-binding entrance is dependent on the ionization state in the Ser209 phosphorylation, which is associated with the regulatory function through the switching on/off of eIF4E phosphorylation.     

X-ray crystal structure and molecular dynamics simulation of bovine pancreas phospholipase A2–n-dodecylphosphorylcholine complex

The crystal structure of n-dodecylphosphorylcholine (n-C₁₂PC)–bovine pancreas phospholipase A₂ (PLA₂) complex provided the following structural.characteristics: (1) the dodecyl chain of n-C₁₂PC was located at the PLA₂ N -terminal helical region by hydrophobic interactions, which corresponds to the binding pocket of 2-acyl fatty acid chain (β-chain) of the substrate phospholipid, (2) the region from Lys-53 to Lys-56 creates a cholinereceiving pocket of n-C₁₂PC and (3) the N-termillal group of Ala-1 shifts significantly toward the Tyr-52 OH group by the binding of the n-C1₂PC inhibitor. Since the accuracy of the X-ray analysis (R = 0.275 at 2.3 Å resolution) was insufficient to establish these important X-ray insights, the complex structure was further investigated through the molecular dynamics (M D) simulation, assuming a system in aqueous solution at 310K. The M D simulation covering 176 ps showed that the structural characteristics observed by X-ray analysis are intrinsic and also stable in the dynamic state. Furthermore, the M D simulation made clear that the PLA₂ binding pocket is large enough to permit the conformational fluctuation of the n-C₁₂PC hydrocarbon chain. © 1994 Wiley-Liss, Inc. © 1994 Wiley-Liss, Inc.     

3 Nsec molecular dynamics simulation of the protein ubiquitin and comparison with X-ray crystal and solution NMR structures.

Mainly due to computational limitations, past protein molecular dynamics simulations have rarely been extended to 300 psec; we are not aware of any published results beyond 350 psec. The present work compares a 3000 psec simulation of the protein ubiquitin with the available x-ray crystallographic and solution NMR structures. Aside from experimental structure availability, ubiquitin was studied because of its relatively small size (76 amino acids) and lack of disulfide bridges. An implicit solvent model was used except for explicit treatment of waters of crystallization. We found that the simulated average structure retains most of the character of the starting x-ray crystal structure. In two highly surface accessible regions, the simulation was not in agreement with the x-ray structure. In addition, there are six backbone-backbone hydrogen bonds that are in conflict between the solution NMR and x-ray crystallographic structures; two are bonds that the NMR does not locate, and four are ones that the two methods disagree upon the donor. Concerning these six backbone-backbone hydrogen bonds, the present simulation agrees with the solution NMR structure in five out-of-the six cases, in that if a hydrogen bond is present in the x-ray structure and not in the NMR structure, the bond breaks within 700 psec. Of the two hydrogen bonds that are found in the NMR structure and not in the x-ray structure, one forms at 1400 psec and the other forms rarely. The present results suggest that relatively long molecular dynamics simulations, that use protein x-ray crystal coordinates for the starting structure and a computationally efficient solvent representation, may be used to gain an understanding of conformational and dynamic differences between the solid-crystal and dilute-solution states.     

Nuclear magnetic resonance study of ring conformations of cyclic tetradepsipeptide AM-toxins and related peptides

Cyclic tetradepsipeptides, AM-toxin I and II, are the host-specific phytotoxins of Alternaria mali. In order to elucidate conformation-toxicity relationships, we analyzed the 270-MHz proton nmr spectra of AM-toxins and hydrogenated analogs, (D -Ala²)AM-toxin I (toxic) and (L -Ala²)AM-toxin I (not toxic), in (C²H₃)₂SO. These cyclic tetradepsipeptides do not contain N-substituted amino acid residues, and all the peptide and ester groups have been found to be transoid. Two conformers with very unequal populations have been found for AM-toxin I and II; the C^β?C^α? C?O conformations of the Dha² residues are nonplanar S-trans in the major conformer and nonplanar S-cis in the minor conformer. Only one ring conformation has been found for each of (L -Ala²) and (D -Ala²)AM-toxin I. (L -Ala²)AM-toxin I takes a C₄-type ring conformation; all the C?O groups and C^α-H bonds are oriented to the same side of the ring. (D -Ala²)AM-toxin I takes a new ring conformation; the side chain and C?O group of the L -Amp¹ residue are oriented to the same side of the ring. This new conformation is also found for the major conformers of AM-toxin I and II and thus appears to be required for the toxicity. The ring conformations of Tyr(OCH₃)¹-bearing analog tetradepsipeptides have been found to be much the same as those of Amp¹-bearing depsipeptides. Furthermore, on the basis of the two distinct conformations of (D -Ala²) and (L -Ala²)AM-toxin I, an empirical rule is proposed for the stable ring conformations of cyclic tetra-D ,L -peptides, not containing N-substituted amino acid residues.