Multiscale modeling of the plastic behaviour in single crystal tungsten: from atomistic to crystal plasticity simulations

En una planta de fusión, los materiales en contacto con el plasma así como los materiales de primera pared experimentan condiciones particularmente hostiles al estar expuestos a altos flujos de partículas, neutrones y grandes cargas térmicas. Como consecuencia de estas diferentes y complejas condiciones de trabajo, el estudio, desarrollo y diseño de estos materiales es uno de los más importantes retos que ha surgido en los últimos años para la comunidad científica en el campo de los materiales y la energía. Debido a su baja tasa de erosión, alta resistencia al sputtering, alta conductividad térmica, muy alto punto de fusión y baja retención de tritio, el tungsteno (wolframio) es un importante candidato como material de primera pared y como posible material estructural avanzado en fusión por confinamiento magnético e inercial. Sin embargo, el tiempo de vida del tungsteno viene controlado por diversos factores como son su respuesta termo-mecánica en la superficie, la posibilidad de fusión y el fallo por acumulación de helio. Es por ello que el tiempo de vida limitado por la respuesta mecánica del tungsteno (W), y en particular su fragilidad, sean dos importantes aspectos que tienes que ser investigados. El comportamiento plástico en materiales refractarios con estructura cristalina cúbica centrada en las caras (bcc) como el tungsteno está gobernado por las dislocaciones de tipo tornillo a escala atómica y por conjuntos e interacciones de dislocaciones a escalas más grandes. El modelado de este complejo comportamiento requiere la aplicación de métodos capaces de resolver de forma rigurosa cada una de las escalas. El trabajo que se presenta en esta tesis propone un modelado multiescala que es capaz de dar respuestas ingenieriles a las solicitudes técnicas del tungsteno, y que a su vez está apoyado por la rigurosa física subyacente a extensas simulaciones atomísticas. En primer lugar, las propiedades estáticas y dinámicas de las dislocaciones de tipo tornillo en cinco potenciales interatómicos de tungsteno son comparadas, determinando cuáles de ellos garantizan una mayor fidelidad física y eficiencia computacional. Las grandes tasas de deformación asociadas a las técnicas de dinámica molecular hacen que las funciones de movilidad de las dislocaciones obtenidas no puedan ser utilizadas en los siguientes pasos del modelado multiescala. En este trabajo, proponemos dos métodos alternativos para obtener las funciones de movilidad de las dislocaciones: un modelo Monte Cario cinético y expresiones analíticas. El conjunto de parámetros necesarios para formular el modelo de Monte Cario cinético y la ley de movilidad analítica son calculados atomísticamente. Estos parámetros incluyen, pero no se limitan a: la determinación de las entalpias y energías de formación de las parejas de escalones que forman las dislocaciones, la parametrización de los efectos de no Schmid característicos en materiales bcc,etc. Conociendo la ley de movilidad de las dislocaciones en función del esfuerzo aplicado y la temperatura, se introduce esta relación como ecuación de flujo dentro de un modelo de plasticidad cristalina. La predicción del modelo sobre la dependencia del límite de fluencia con la temperatura es validada experimentalmente con ensayos uniaxiales en tungsteno monocristalino. A continuación, se calcula el límite de fluencia al aplicar ensayos uniaxiales de tensión para un conjunto de orientaciones cristalográticas dentro del triángulo estándar variando la tasa de deformación y la temperatura de los ensayos. Finalmente, y con el objetivo de ser capaces de predecir una respuesta más dúctil del tungsteno para una variedad de estados de carga, se realizan ensayos biaxiales de tensión sobre algunas de las orientaciones cristalográficas ya estudiadas en función de la temperatura.-------------------------------------------------------------------------ABSTRACT ----------------------------------------------------------Tungsten and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield and low longterm disposal radioactive footprint. However, tungsten also presents a very low fracture toughness, mostly associated with inter-granular failure and bulk plasticity, that limits its applications. As a result of these various and complex conditions of work, the study, development and design of these materials is one of the most important challenges that have emerged in recent years to the scientific community in the field of materials for energy applications. The plastic behavior of body-centered cubic (bcc) refractory metals like tungsten is governed by the kink-pair mediated thermally activated motion of h¿ (\1 11)i screw dislocations on the atomistic scale and by ensembles and interactions of dislocations at larger scales. Modeling this complex behavior requires the application of methods capable of resolving rigorously each relevant scale. The work presented in this thesis proposes a multiscale model approach that gives engineering-level responses to the technical specifications required for the use of tungsten in fusion energy reactors, and it is also supported by the rigorous underlying physics of extensive atomistic simulations. First, the static and dynamic properties of screw dislocations in five interatomic potentials for tungsten are compared, determining which of these ensure greater physical fidelity and computational efficiency. The large strain rates associated with molecular dynamics techniques make the dislocation mobility functions obtained not suitable to be used in the next steps of the multiscale model. Therefore, it is necessary to employ mobility laws obtained from a different method. In this work, we suggest two alternative methods to get the dislocation mobility functions: a kinetic Monte Carlo model and analytical expressions. The set of parameters needed to formulate the kinetic Monte Carlo model and the analytical mobility law are calculated atomistically. These parameters include, but are not limited to: enthalpy and energy barriers of kink-pairs as a function of the stress, width of the kink-pairs, non-Schmid effects ( both twinning-antitwinning asymmetry and non-glide stresses), etc. The function relating dislocation velocity with applied stress and temperature is used as the main source of constitutive information into a dislocation-based crystal plasticity framework. We validate the dependence of the yield strength with the temperature predicted by the model against existing experimental data of tensile tests in singlecrystal tungsten, with excellent agreement between the simulations and the measured data. We then extend the model to a number of crystallographic orientations uniformly distributed in the standard triangle and study the effects of temperature and strain rate. Finally, we perform biaxial tensile tests and provide the yield surface as a function of the temperature for some of the crystallographic orientations explored in the uniaxial tensile tests.

Crystal structure of human p32, a doughnut-shaped acidic mitochondrial matrix protein

Human p32 (also known as SF2-associated p32, p32/TAP, and gC1qR) is a conserved eukaryotic protein that localizes predominantly in the mitochondrial matrix. It is thought to be involved in mitochondrial oxidative phosphorylation and in nucleus–mitochondrion interactions. We report the crystal structure of p32 determined at 2.25 Å resolution. The structure reveals that p32 adopts a novel fold with seven consecutive antiparallel β-strands flanked by one N-terminal and two C-terminal α-helices. Three monomers form a doughnut-shaped quaternary structure with an unusually asymmetric charge distribution on the surface. The implications of the structure on previously proposed functions of p32 are discussed and new specific functional properties are suggested.

Crystal structure of a thermostable type B DNA polymerase from Thermococcus gorgonarius

Most known archaeal DNA polymerases belong to the type B family, which also includes the DNA replication polymerases of eukaryotes, but maintain high fidelity at extreme conditions. We describe here the 2.5 Å resolution crystal structure of a DNA polymerase from the Archaea Thermococcus gorgonarius and identify structural features of the fold and the active site that are likely responsible for its thermostable function. Comparison with the mesophilic B type DNA polymerase gp43 of the bacteriophage RB69 highlights thermophilic adaptations, which include the presence of two disulfide bonds and an enhanced electrostatic complementarity at the DNA–protein interface. In contrast to gp43, several loops in the exonuclease and thumb domains are more closely packed; this apparently blocks primer binding to the exonuclease active site. A physiological role of this “closed” conformation is unknown but may represent a polymerase mode, in contrast to an editing mode with an open exonuclease site. This archaeal B DNA polymerase structure provides a starting point for structure-based design of polymerases or ligands with applications in biotechnology and the development of antiviral or anticancer agents.

Crystal structure of the HLA-Cw3 allotype-specific killer cell inhibitory receptor KIR2DL2

Killer cell inhibitory receptors (KIR) protect class I HLAs expressing target cells from natural killer (NK) cell-mediated lysis. To understand the molecular basis of this receptor-ligand recognition, we have crystallized the extracellular ligand-binding domains of KIR2DL2, a member of the Ig superfamily receptors that recognize HLA-Cw1, 3, 7, and 8 allotypes. The structure was determined in two different crystal forms, an orthorhombic P212121 and a trigonal P3221 space group, to resolutions of 3.0 and 2.9 Å, respectively. The overall fold of this structure, like KIR2DL1, exhibits K-type Ig topology with cis-proline residues in both domains that define β-strand switching, which sets KIR apart from the C2-type hematopoietic growth hormone receptor fold. The hinge angle of KIR2DL2 is approximately 80°, 14° larger than that observed in KIR2DL1 despite the existence of conserved hydrophobic residues near the hinge region. There is also a 5° difference in the observed hinge angles in two crystal forms of 2DL2, suggesting that the interdomain hinge angle is not fixed. The putative ligand-binding site is formed by residues from several variable loops with charge distribution apparently complementary to that of HLA-C. The packing of the receptors in the orthorhombic crystal form offers an intriguing model for receptor aggregation on the cell surface.

A dimeric crystal structure for the N-terminal two domains of intercellular adhesion molecule-1

The 3.0-Å structure of a 190-residue fragment of intercellular adhesion molecule-1 (ICAM-1, CD54) reveals two tandem Ig-superfamily (IgSF) domains. Each of two independent molecules dimerizes identically with a symmetry-related molecule over a hydrophobic interface on the BED sheet of domain 1, in agreement with dimerization of ICAM-1 on the cell surface. The residues that bind to the integrin LFA-1 are well oriented for bivalent binding in the dimer, with the critical Glu-34 residues pointing away from each other on the periphery. Residues that bind to rhinovirus are in the flexible BC and FG loops at the tip of domain 1, and these and the upper half of domain 1 are well exposed in the dimer for docking to virus. By contrast, a residue important for binding to Plasmodium falciparum-infected erythrocytes is in the dimer interface. The presence of A′ strands in both domains 1 and 2, conserved hydrogen bonds at domain junctions, and elaborate hydrogen bond networks around the key integrin binding residues in domain 1 make these domains suited to resist tensile forces during adhesive interactions. A subdivision of the intermediate (I) set of IgSF domains is proposed in which domain 1 of ICAM-1 and previously described I set domains belong to the I1 set and domain 2 of ICAM-1, ICAM-2, and vascular cell adhesion molecule-1 belong to the I2 set.

Crystal structure of chemically synthesized [N33A] stromal cell-derived factor 1α, a potent ligand for the HIV-1 “fusin” coreceptor

Stromal cell-derived factor-1α (SDF-1α ) is a member of the chemokine superfamily and functions as a growth factor and chemoattractant through activation of CXCR4/LESTR/Fusin, a G protein-coupled receptor. This receptor also functions as a coreceptor for T-tropic syncytium-inducing strains of HIV-1. SDF-1α antagonizes infectivity of these strains by competing with gp120 for binding to the receptor. The crystal structure of a variant SDF-1α ([N33A]SDF-1α ) prepared by total chemical synthesis has been refined to 2.2-Å resolution. Although SDF-1α adopts a typical chemokine β-β-β-α topology, the packing of the α-helix against the β-sheet is strikingly different. Comparison of SDF-1α with other chemokine structures confirms the hypothesis that SDF-1α may be either an ancestral protein from which all other chemokines evolved or the chemokine that is the least divergent from a primordial chemokine. The structure of SDF-1α reveals a positively charged surface ideal for binding to the negatively charged extracellular loops of the CXCR4 HIV-1 coreceptor. This ionic complementarity is likely to promote the interaction of the mobile N-terminal segment of SDF-1α with interhelical sites of the receptor, resulting in a biological response.

Crystal structure of the BTB domain from PLZF

The BTB domain (also known as the POZ domain) is an evolutionarily conserved protein–protein interaction motif found at the N terminus of 5–10% of C2H2-type zinc-finger transcription factors, as well as in some actin-associated proteins bearing the kelch motif. Many BTB proteins are transcriptional regulators that mediate gene expression through the control of chromatin conformation. In the human promyelocytic leukemia zinc finger (PLZF) protein, the BTB domain has transcriptional repression activity, directs the protein to a nuclear punctate pattern, and interacts with components of the histone deacetylase complex. The association of the PLZF BTB domain with the histone deacetylase complex provides a mechanism of linking the transcription factor with enzymatic activities that regulate chromatin conformation. The crystal structure of the BTB domain of PLZF was determined at 1.9 Å resolution and reveals a tightly intertwined dimer with an extensive hydrophobic interface. Approximately one-quarter of the monomer surface area is involved in the dimer intermolecular contact. These features are typical of obligate homodimers, and we expect the full-length PLZF protein to exist as a branched transcription factor with two C-terminal DNA-binding regions. A surface-exposed groove lined with conserved amino acids is formed at the dimer interface, suggestive of a peptide-binding site. This groove may represent the site of interaction of the PLZF BTB domain with nuclear corepressors or other nuclear proteins.

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.

Crystal structure of a multifunctional 2-Cys peroxiredoxin heme-binding protein 23 kDa/proliferation-associated gene product

Heme-binding protein 23 kDa (HBP23), a rat isoform of human proliferation-associated gene product (PAG), is a member of the peroxiredoxin family of peroxidases, having two conserved cysteine residues. Recent biochemical studies have shown that HBP23/PAG is an oxidative stress-induced and proliferation-coupled multifunctional protein that exhibits specific bindings to c-Abl protein tyrosine kinase and heme, as well as a peroxidase activity. A 2.6-Å resolution crystal structure of rat HBP23 in oxidized form revealed an unusual dimer structure in which the active residue Cys-52 forms a disulfide bond with conserved Cys-173 from another subunit by C-terminal tail swapping. The active site is largely hydrophobic with partially exposed Cys-173, suggesting a reduction mechanism of oxidized HBP23 by thioredoxin. Thus, the unusual cysteine disulfide bond is involved in peroxidation catalysis by using thioredoxin as the source of reducing equivalents. The structure also provides a clue to possible interaction surfaces for c-Abl and heme. Several significant structural differences have been found from a 1-Cys peroxiredoxin, ORF6, which lacks the C-terminal conserved cysteine corresponding to Cys-173 of HBP23.

The 2.0-Å resolution crystal structure of a trimeric antibody fragment with noncognate VH–VL domain pairs shows a rearrangement of VH CDR3

The 2.0-Å resolution x-ray crystal structure of a novel trimeric antibody fragment, a “triabody,” has been determined. The trimer is made up of polypeptides constructed in a manner identical to that previously described for some “diabodies”: a VL domain directly fused to the C terminus of a VH domain—i.e., without any linker sequence. The trimer has three Fv heads with the polypeptides arranged in a cyclic, head-to-tail fashion. For the particular structure reported here, the polypeptide was constructed with a VH domain from one antibody fused to the VL domain from an unrelated antibody giving rise to “combinatorial” Fvs upon formation of the trimer. The structure shows that the exchange of the VL domain from antibody B1-8, a Vλ domain, with the VL domain from antibody NQ11, a Vκ domain, leads to a dramatic conformational change in the VH CDR3 loop of antibody B1-8. The magnitude of this change is similar to the largest of the conformational changes observed in antibody fragments in response to antigen binding. Combinatorial pairing of VH and VL domains constitutes a major component of antibody diversity. Conformationally flexible antigen-binding sites capable of adapting to the specific CDR3 loop context created upon VH–VL pairing may be employed by the immune system to maximize the structural diversity of the immune response.

The crystal structure of human interferon β at 2.2-Å resolution

Type I interferons (IFNs) are helical cytokines that have diverse biological activities despite the fact that they appear to interact with the same receptor system. To achieve a better understanding of the structural basis for the different activities of α and β IFNs, we have determined the crystal structure of glycosylated human IFN-β at 2.2-Å resolution by molecular replacement. The molecule adopts a fold similar to that of the previously determined structures of murine IFN-β and human IFN-α2b but displays several distinct structural features. Like human IFN-α2b, human IFN-β contains a zinc-binding site at the interface of the two molecules in the asymmetric unit, raising the question of functional relevance for IFN-β dimers. However, unlike the human IFN-α2b dimer, in which homologous surfaces form the interface, human IFN-β dimerizes with contact surfaces from opposite sides of the molecule. The relevance of the structure to the effects of point mutations in IFN-β at specific exposed residues is discussed. A potential role of ligand–ligand interactions in the conformational assembly of IFN receptor components is discussed.

Crystal structure of the RNA-dependent RNA polymerase of hepatitis C virus

We report the crystal structure of the RNA-dependent RNA polymerase of hepatitis C virus, a major human pathogen, to 2.8-Å resolution. This enzyme is a key target for developing specific antiviral therapy. The structure of the catalytic domain contains 531 residues folded in the characteristic fingers, palm, and thumb subdomains. The fingers subdomain contains a region, the “fingertips,” that shares the same fold with reverse transcriptases. Superposition to the available structures of the latter shows that residues from the palm and fingertips are structurally equivalent. In addition, it shows that the hepatitis C virus polymerase was crystallized in a closed fingers conformation, similar to HIV-1 reverse transcriptase in ternary complex with DNA and dTTP [Huang H., Chopra, R., Verdine, G. L. & Harrison, S. C. (1998) Science 282, 1669–1675]. This superposition reveals the majority of the amino acid residues of the hepatitis C virus enzyme that are likely to be implicated in binding to the replicating RNA molecule and to the incoming NTP. It also suggests a rearrangement of the thumb domain as well as a possible concerted movement of thumb and fingertips during translocation of the RNA template-primer in successive polymerization rounds.

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.

Crystal structure of SQD1, an enzyme involved in the biosynthesis of the plant sulfolipid headgroup donor UDP-sulfoquinovose

The SQD1 enzyme is believed to be involved in the biosynthesis of the sulfoquinovosyl headgroup of plant sulfolipids, catalyzing the transfer of SO3− to UDP-glucose. We have determined the structure of the complex of SQD1 from Arabidopsis thaliana with NAD+ and the putative substrate UDP-glucose at 1.6-Å resolution. Both bound ligands are completely buried within the binding cleft, along with an internal solvent cavity which is the likely binding site for the, as yet, unidentified sulfur-donor substrate. SQD1 is a member of the short-chain dehydrogenase/reductase (SDR) family of enzymes, and its structure shows a conservation of the SDR catalytic residues. Among several highly conserved catalytic residues, Thr-145 forms unusually short hydrogen bonds with both susceptible hydroxyls of UDP-glucose. A His side chain may also be catalytically important in the sulfonation.

Crystal structure of the Sec18p N-terminal domain

Yeast Sec18p and its mammalian orthologue N-ethylmaleimide-sensitive fusion protein (NSF) are hexameric ATPases with a central role in vesicle trafficking. Aided by soluble adapter factors (SNAPs), Sec18p/NSF induces ATP-dependent disassembly of a complex of integral membrane proteins from the vesicle and target membranes (SNAP receptors). During the ATP hydrolysis cycle, the Sec18p/NSF homohexamer undergoes a large-scale conformational change involving repositioning of the most N terminal of the three domains of each protomer, a domain that is required for SNAP-mediated interaction with SNAP receptors. Whether an internal conformational change in the N-terminal domains accompanies their reorientation with respect to the rest of the hexamer remains to be addressed. We have determined the structure of the N-terminal domain from Sec18p by x-ray crystallography. The Sec18p N-terminal domain consists of two β-sheet-rich subdomains connected by a short linker. A conserved basic cleft opposite the linker may constitute a SNAP-binding site. Despite structural variability in the linker region and in an adjacent loop, all three independent molecules in the crystal asymmetric unit have the identical subdomain interface, supporting the notion that this interface is a preferred packing arrangement. However, the linker flexibility allows for the possibility that other subdomain orientations may be sampled.