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.

Three-dimensional cryoelectron microscopy of dimeric kinesin and ncd motor domains on microtubules.

Kinesin and ncd motor proteins are homologous in sequence yet move in opposite directions along microtubules. We have previously shown that monomeric kinesin and ncd bind in the same orientation on equivalent sites relative to the ends of tubulin sheets of known polarity. We now report cryoelectron microscope images of 16-protofilament microtubules decorated with both single- and double-headed kinesin and double-headed ncd. Three-dimensional density maps and difference maps show that, in adenosine 5'-[beta,gamma-imido]triphosphate, both dimeric motors bind tightly to microtubules via one head, leaving the other free, though apparently in a fixed position. The attached heads of dimers bind to tubulin in the same way as single kinesin heads. The second heads are connected to the tops of the first but, whereas the second kinesin head is closely associated with the first, pairs of ncd heads are splayed apart. There is also a distinct difference in orientation: the second kinesin head is tilted toward the microtubule plus end, while the second head of ncd points toward the minus end.

Structure of the ripple phase in lecithin bilayers.

The phases of the x-ray form factors are derived for the ripple (Pbeta') thermodynamic phase in the lecithin bilayer system. By combining these phases with experimental intensity data, the electron density map of the ripple phase of dimyristoyl-phosphatidylcholine is constructed. The phases are derived by fitting the intensity data to two-dimensional electron density models, which are created by convolving an asymmetric triangular ripple profile with a transbilayer electron density profile. The robustness of the model method is indicated by the result that many different models of the transbilayer profile yield essentially the same phases, except for the weaker, purely ripple (0,k) peaks. Even with this residual ambiguity, the ripple profile is well determined, resulting in 19 angstroms for the ripple amplitude and 10 degrees and 26 degrees for the slopes of the major and the minor sides, respectively. Estimates for the bilayer head-head spacings show that the major side of the ripple is consistent with gel-like structure, and the minor side appears to be thinner with lower electron density.

X-ray structure determination of a metastable state of carbonmonoxy myoglobin after photodissociation.

The x-ray structure of carbon monoxide (CO)-ligated myoglobin illuminated during data collection by a laser diode at the wavelength lambda = 690 nm has been determined to a resolution of 1.7 A at T = 36 K. For comparison, we also measured data sets of deoxymyoglobin and CO-ligated myoglobin. In the photon-induced structure the electron density associated with the CO ligand can be described by a tube extending from the iron into the heme pocket over more than 4 A. This density can be interpreted by two discrete positions of the CO molecule. One is close to the heme iron and can be identified to be bound CO. In the second, the CO is dissociated from the heme iron and lies on top of pyrrole ring C. At our experimental conditions the overall structure of myoglobin in the metastable state is close to the structure of a CO-ligated molecule. However, the iron has essentially relaxed into the position of deoxymyoglobin. We compare our results with those of Schlichting el al. [Schlichting, I., Berendzen, J., Phillips, G. N., Jr., & Sweet, R. M. (1994) Nature 317, 808-812], who worked with the myoglobin mutant (D122N) that crystallizes in the space group P6 and Teng et al. [Teng, T. Y., Srajer, V. & Moffat, K. (1994) Nat. Struct. Biol. 1, 701-705], who used native myoglobin crystals of the space group P2(1). Possible reasons for the structural differences are discussed.

Two-dimensional protein crystallization via metal-ion coordination by naturally occurring surface histidines.

A powerful and potentially general approach to the targeting and crystallization of proteins on lipid interfaces through coordination of surface histidine residues to lipid-chelated divalent metal ions is presented. This approach, which should be applicable to the crystallization of a wide range of naturally occurring or engineered proteins, is illustrated here by the crystallization of streptavidin on a monolayer of an iminodiacetate-Cu(II) lipid spread at the air-water interface. This method allows control of the protein orientation at interfaces, which is significant for the facile production of highly ordered protein arrays and for electron density mapping in structural analysis of two-dimensional crystals. Binding of native streptavidin to the iminodiacetate-Cu lipids occurs via His-87, located on the protein surface near the biotin binding pocket. The two-dimensional streptavidin crystals show a previously undescribed microscopic shape that differs from that of crystals formed beneath biotinylated lipids.

Crystallographic analysis of endogenous peptides associated with HLA-DR1 suggests a common, polyproline II-like conformation for bound peptides.

The structure of the human major histocompatibility complex (MHC) class II molecule HLA-DR1 derived from the human lymphoblastoid cell line LG-2 has been determined in a complex with the Staphylococcus aureus enterotoxin B superantigen. The HLA-DR1 molecule contains a mixture of endogenous peptides derived from cellular or serum proteins bound in the antigen-binding site, which copurify with the class II molecule. Continuous electron density for 13 amino acid residues is observed in the MHC peptide-binding site, suggesting that this is the core length of peptide that forms common interactions with the MHC molecule. Electron density is also observed for side chains of the endogenous peptides. The electron density corresponding to peptide side chains that interact with the DR1-binding site is more clearly defined than the electron density that extends out of the binding site. The regions of the endogenous peptides that interact with DRI are therefore either more restricted in conformation or sequence than the peptide side chains or amino acids that project out of the peptide-binding site. The hydrogen-bond interactions and conformation of a peptide model built into the electron density are similar to other HLA-DR-peptide structures. The bound peptides assume a regular conformation that is similar to a polyproline type II helix. The side-chain pockets and conserved asparagine residues of the DR1 molecule are well-positioned to interact with peptides in the polyproline type II conformation and may restrict the range of acceptable peptide conformations.

Crystal structure of Pseudomonas aeruginosa catabolic ornithine transcarbamoylase at 3.0-A resolution: a different oligomeric organization in the transcarbamoylase family.

The crystal structure of the Glu-105-->Gly mutant of catabolic ornithine transcarbamoylase (OTCase; carbamoyl phosphate + L-ornithine = orthophosphate + L-citrulline, EC 2.1.3.3) from Pseudomonas aeruginosa has been determined at 3.0-A resolution. This mutant is blocked in the active R (relaxed) state. The structure was solved by the molecular replacement method, starting from a crude molecular model built from a trimer of the catalytic subunit of another transcarbamoylase, the extensively studied aspartate transcarbamoylase (ATCase) from Escherichia coli. This model was used to generate initial low-resolution phases at 8-A resolution, which were extended to 3-A by noncrystallographic symmetry averaging. Four phase extensions were required to obtain an electron density map of very high quality from which the final model was built. The structure, including 4020 residues, has been refined to 3-A, and the current crystallographic R value is 0.216. No solvent molecules have been added to the model. The catabolic OTCase is a dodecamer composed of four trimers organized in a tetrahedral manner. Each monomer is composed of two domains. The carbamoyl phosphate binding domain shows a strong structural homology with the equivalent ATCase part. In contrast, the other domain, mainly implicated in the binding of the second substrate (ornithine for OTCase and aspartate for ATCase) is poorly conserved. The quaternary structures of these two allosteric transcarbamoylases are quite divergent: the E. coli ATCase has pseudo-32 point-group symmetry, with six catalytic and six regulatory chains; the catabolic OTCase has 23 point-group symmetry and only catalytic chains. However, both enzymes display homotropic and heterotropic cooperativity.

X-ray structure of clotting factor IXa: active site and module structure related to Xase activity and hemophilia B.

Hereditary deficiency of factor IXa (fIXa), a key enzyme in blood coagulation, causes hemophilia B, a severe X chromosome-linked bleeding disorder afflicting 1 in 30,000 males; clinical studies have identified nearly 500 deleterious variants. The x-ray structure of porcine fIXa described here shows the atomic origins of the disease, while the spatial distribution of mutation sites suggests a structural model for factor X activation by phospholipid-bound fIXa and cofactor VIIIa. The 3.0-A-resolution diffraction data clearly show the structures of the serine proteinase module and the two preceding epidermal growth factor (EGF)-like modules; the N-terminal Gla module is partially disordered. The catalytic module, with covalent inhibitor D-Phe-1I-Pro-2I-Arg-3I chloromethyl ketone, most closely resembles fXa but differs significantly at several positions. Particularly noteworthy is the strained conformation of Glu-388, a residue strictly conserved in known fIXa sequences but conserved as Gly among other trypsin-like serine proteinases. Flexibility apparent in electron density together with modeling studies suggests that this may cause incomplete active site formation, even after zymogen, and hence the low catalytic activity of fIXa. The principal axes of the oblong EGF-like domains define an angle of 110 degrees, stabilized by a strictly conserved and fIX-specific interdomain salt bridge. The disorder of the Gla module, whose hydrophobic helix is apparent in electron density, can be attributed to the absence of calcium in the crystals; we have modeled the Gla module in its calcium form by using prothrombin fragment 1. The arched module arrangement agrees with fluorescence energy transfer experiments. Most hemophilic mutation sites of surface fIX residues occur on the concave surface of the bent molecule and suggest a plausible model for the membrane-bound ternary fIXa-FVIIIa-fX complex structure: fIXa and an equivalently arranged fX arch across an underlying fVIIIa subdomain from opposite sides; the stabilizing fVIIIa interactions force the catalytic modules together, completing fIXa active site formation and catalytic enhancement.

The crystal structure of Pseudomonas aeruginosa exotoxin domain III with nicotinamide and AMP: conformational differences with the intact exotoxin.

Domain III of Pseudomonas aeruginosa exotoxin A catalyses the transfer of ADP-ribose from NAD to a modified histidine residue of elongation factor 2 in eukaryotic cells, thus inactivating elongation factor 2. This domain III is inactive in the intact toxin but is active in the isolated form. We report here the 2.5-A crystal structure of this isolated domain crystallized in the presence of NAD and compare it with the corresponding structure in the intact Pseudomonas aeruginosa exotoxin A. We observe a significant conformational difference in the active site region from Arg-458 to Asp-463. Contacts with part of domain II in the intact toxin prevent the adoption of the isolated domain conformation and provide a structural explanation for the observed inactivity. Additional electron density in the active site region corresponds to separate AMP and nicotinamide and indicates that the NAD has been hydrolyzed. The structure has been compared with the catalytic domain of the diphtheria toxin, which was crystallized with ApUp.

Simulação e modelagem computacional de dados de espalhamento à baixos ângulos - enfoque em estruturas de alta simetria

Esta tese apresenta uma abordagem para a criação rápida de modelos em diferentes geometrias (complexas ou de alta simetria) com objetivo de calcular a correspondente intensidade espalhada, podendo esta ser utilizada na descrição de experimentos de es- palhamento à baixos ângulos. A modelagem pode ser realizada com mais de 100 geome- trias catalogadas em um Banco de Dados, além da possibilidade de construir estruturas a partir de posições aleatórias distribuídas na superfície de uma esfera. Em todos os casos os modelos são gerados por meio do método de elementos finitos compondo uma única geometria, ou ainda, compondo diferentes geometrias, combinadas entre si a partir de um número baixo de parâmetros. Para realizar essa tarefa foi desenvolvido um programa em Fortran, chamado de Polygen, que permite modelar geometrias convexas em diferentes formas, como sólidos, cascas, ou ainda com esferas ou estruturas do tipo DNA nas arestas, além de usar esses modelos para simular a curva de intensidade espalhada para sistemas orientados e aleatoriamente orientados. A curva de intensidade de espalhamento é calculada por meio da equação de Debye e os parâmetros que compõe cada um dos modelos, podem ser otimizados pelo ajuste contra dados experimentais, por meio de métodos de minimização baseados em simulated annealing, Levenberg-Marquardt e algorítmicos genéticos. A minimização permite ajustar os parâmetros do modelo (ou composição de modelos) como tamanho, densidade eletrônica, raio das subunidades, entre outros, contribuindo para fornecer uma nova ferramenta para modelagem e análise de dados de espalhamento. Em outra etapa desta tese, é apresentado o design de modelos atomísticos e a sua respectiva simulação por Dinâmica Molecular. A geometria de dois sistemas auto-organizado de DNA na forma de octaedro truncado, um com linkers de 7 Adeninas e outro com linkers de ATATATA, foram escolhidas para realizar a modelagem atomística e a simulação por Dinâmica Molecular. Para este sistema são apresentados os resultados de Root Mean Square Deviations (RMSD), Root Mean Square Fluctuations (RMSF), raio de giro, torção das hélices duplas de DNA além da avaliação das ligações de Hidrogênio, todos obtidos por meio da análise de uma trajetória de 50 ns.

Estudo de efeitos da polarização eletrostática periférica no Tokamak TCABR

Efeitos da polarização eletrostática de eletrodos na periferia de tokamaks têm sido investigados em pequenos tokamaks e mesmo em alguns tokamaks de grande porte. Em geral as experiências são realizadas em condições em que bifurcação do campo elétrico radial é obtida, processo este identificado como modo H de polarização. No Tokamak TCABR, as experiências indicam que o confinamento aumenta para tensões aplicadas até +300 volts, atingindo um máximo de duas vezes o tempo de confinamento do modo L, mas sem bifurcação. Indícios de bifurcação foram notados com +400 V de polarização, mas a descarga termina devido à excitação da atividade MHD, ainda sob investigação. No presente trabalho, a pesquisa é aprofundada com a utilização de uma sonda de Langmuir com 18 pinos dispostos em duas fileiras sob a forma de um ancinho (rake probe) o que permite a medição da temperatura, densidade e flutuação de potencial ao longo do raio menor na periferia do Tokamak. A resolução temporal desse sistema é de cerca de 0,5 ms, para a temperatura, e 5 microssegundos para densidade e potencial flutuante do plasma. Outra sonda eletrostática com 5-pinos na mesma posição radial, mas em diferentes posições poloidal e toroidal foi usada para medições de turbulência e transporte de partículas. Os efeitos da polarização foram investigados e indicam que os níveis de turbulência e transporte começam a diminuir entre +150 e +200 V e para +300 V chegam a atingir uma quase supressão. Nesse mesmo intervalo de tensão a densidade começa a aumentar e para +300 V chega a ser um fator de aproximadamente 2. Quanto ao perfil de temperatura a variação é pouco significativa, mas as incertezas das medidas são maiores. Esses dados são compatíveis com a criação de uma barreira de transporte na região entre o eletrodo em r = 17 cm e o limitador em a = 18 cm. Além disso, o campo elétrico radial mostra forte cisalhamento nessa região. Tomando o início da subida do potencial flutuante como origem de uma escala de tempo, o atraso temporal do início da subida da densidade de elétrons e o atraso do início do decréscimo do transporte de partículas foram medidos. Os resultados são 50 microssegundos para a densidade de elétrons e 60 microssegundos para o transporte de partículas. A questão dos limiares de potência é discutida no texto. Os dados desta experiência indicam que o campo elétrico radial desempenha o papel principal para a melhoria do confinamento.

Enhanced spin injection efficiency in ferromagnet/semiconductor tunnel junctions

Within the ballistic transport picture, we have investigated the spin-polarized transport properties of a ferromagnetic metal/two-dimensional semiconductor (FM/SM) hybrid junction and an FM/FM/SM structure using quantum tunnelling theory. Our calculations indicate explicitly that the low spin injection efficiency (SIE) from an FM into an SM, compared with a ferromagnet/normal metal junction, originates from the mismatch of electron densities in the FM and SM. To enhance the SIE from an FM into an SM, we introduce another FM film between them to form FM/FM/SM double tunnel junctions, in which the quantum interference effect will lead to the current polarization exhibiting periodically oscillating behaviour, with a variation according to the thickness of the middle FM film and/or its exchange energy strength. Our results show that, for some suitable values of these parameters, the SIE can reach a very high level, which can also be affected by the electron density in the SM electrode.

The electric field effect on spin injection in tunneling regime

Using the quantum tunneling theory, we investigate the spin-dependent transport properties of the ferromagnetic metal/Schottky barrier/semiconductor heterojunction under the influence of an external electric field. It is shown that increasing the electric field, similar to increasing the electron density in semiconductor, will result in a slight enhancement of spin injection in tunneling regime, and this enhancement is significantly weakened when the tunneling Schottky barrier becomes stronger. Temperature effect on spin injection is also discussed. (C) 2003 Elsevier B.V. All rights reserved.

Crystallization and preliminary X-ray diffraction analysis of importin-alpha acomplexed with NLS peptidomimetics

Importin-alpha is the nuclear import receptor that recognizes cargo proteins with nuclear localization sequences (NLSs). Tile study of NLS peptidomimetics can provide a better understanding of the requirements for the molecular recognition of cargo proteins by importin-alpha, and potentially engender a large number of applications in medicine. Importin-a was crystallized with a set of six NLS peptidomimetics, and X-ray diffraction data were collected in the range 2.1-2.5 angstrom resolution. Preliminary electron density calculations show that the ligands are present in the crystals. (c) 2005 Elsevier B.V All rights reserved.

Enhancement of hydrogen physisorption on single-walled carbon nanotubes resulting from defects created by carbon bombardment

The defect effect on hydrogen adsorption on single-walled carbon nanotubes (SWNTs) has been studied by using extensive molecular dynamics simulations and density functional theory (DFT) calculations. It indicates that the defects created on the exterior wall of the SWNTs by bombarding the tube wall with carbon atoms and C-2 dimers at a collision energy of 20 eV can enhance the hydrogen adsorption potential of the SWNTs substantially. The average adsorption energy for a H-2 molecule adsorbed on the exterior wall of a defected (10,10) SWNT is similar to 150 meV, while that for a H-2 molecule adsorbed on the exterior wall of a perfect (10,10) SWNT is similar to 104 meV. The H-2 sticking coefficient is very sensitive to temperature, and has a maximum value around 70 to 90 K. The electron density contours, the local density of states, and the electron transfers obtained from the DFT calculations clearly indicate that the H-2 molecules are all physisorbed on the SWNTs. At temperatures above 200 K, most of the H-2 molecules adsorbed on the perfect SWNT are soon desorbed, but the H-2 molecules can still remain on the defected SWNTs at 300 K. The detailed processes of H-2 molecules adsorbing on and desorbing from the (10,10) SWNTs are demonstrated.