Leonard M. Thomas

Crystal structure of HLA-A2 bound to LIR-1, a host and viral major histocompatibility complex receptor.

Leukocyte immunoglobulin-like receptor 1 (LIR-1), an inhibitory receptor expressed on monocytes, dendritic cells and lymphocytes, regulates cellular function by binding a broad range of classical and nonclassical major histocompatibility complex (MHC) class I molecules, and the human cytomegalovirus MHC class I homolog UL18. Here we describe the 3.4-Å crystal structure of a complex between the LIR-1 D1D2 domains and the MHC class I molecule HLA-A2. LIR-1 contacts the mostly conserved β2-microglobulin and α3 domains of HLA-A2. The LIR-1 binding site comprises residues at the interdomain hinge, and a patch at the D1 tip. The structure shows how LIR-1 recognizes UL18 and diverse MHC class I molecules, and indicates that a similar mode of MHC class I recognition is used by other LIR family members.

Iterative approach to computational enzyme design.

A general approach for the computational design of enzymes to catalyze arbitrary reactions is a goal at the forefront of the field of protein design. Recently, computationally designed enzymes have been produced for three chemical reactions through the synthesis and screening of a large number of variants. Here, we present an iterative approach that has led to the development of the most catalytically efficient computationally designed enzyme for the Kemp elimination to date. Previously established computational techniques were used to generate an initial design, HG-1, which was catalytically inactive. Analysis of HG-1 with molecular dynamics simulations (MD) and X-ray crystallography indicated that the inactivity might be due to bound waters and high flexibility of residues within the active site. This analysis guided changes to our design procedure, moved the design deeper into the interior of the protein, and resulted in an active Kemp eliminase, HG-2. The cocrystal structure of this enzyme with a transition state analog (TSA) revealed that the TSA was bound in the active site, interacted with the intended catalytic base in a catalytically relevant manner, but was flipped relative to the design model. MD analysis of HG-2 led to an additional point mutation, HG-3, that produced a further threefold improvement in activity. This iterative approach to computational enzyme design, including detailed MD and structural analysis of both active and inactive designs, promises a more complete understanding of the underlying principles of enzymatic catalysis and furthers progress toward reliably producing active enzymes.

Crystal structure of prostate-specific membrane antigen, a tumor marker and peptidase

Prostate-specific membrane antigen (PSMA) is highly expressed in prostate cancer cells and nonprostatic solid tumor neovasculature and is a target for anticancer imaging and therapeutic agents. PSMA acts as a glutamate carboxypeptidase (GCPII) on small molecule substrates, including folate, the anticancer drug methotrexate, and the neuropeptide N-acetyl-l-aspartyl-l-glutamate. Here we present the 3.5-Å crystal structure of the PSMA ectodomain, which reveals a homodimer with structural similarity to transferrin receptor, a receptor for iron-loaded transferrin that lacks protease activity. Unlike transferrin receptor, the protease domain of PSMA contains a binuclear zinc site, catalytic residues, and a proposed substrate-binding arginine patch. Elucidation of the PSMA structure combined with docking studies and a proposed catalytic mechanism provides insight into the recognition of inhibitors and the natural substrate N-acetyl-l-aspartyl-l-glutamate. The PSMA structure will facilitate development of chemotherapeutics, cancer-imaging agents, and agents for treatment of neurological disorders.

Crystal structure of LIR-2 (ILT4) at 1.8 Å: differences from LIR-1 (ILT2) in regions implicated in the binding of the Human Cytomegalovirus class I MHC homolog UL18

BackgroundLeukocyte Immunoglobulin-like Receptor-1 (LIR-1) and LIR-2 (also known as ILT2 and ILT4 respectively) are highly related cell surface receptors that bind a broad range of class I MHC molecules with low (μM) affinities. Expressed on monocytic cells and macrophages, both molecules transmit inhibitory signals after binding ligands. In addition to binding host class I MHC, the LIR-1 molecule, which is also expressed on lymphoid tissues, binds with a high (nM) affinity to UL18, a class I MHC homolog encoded by Human Cytomegalovirus (HCMV). In comparison, LIR-2 binds UL18 only weakly (μM KD). To understand how HCMV preferentially targets the more broadly expressed LIR-1 molecule, we determined the crystal structure of a ligand-binding fragment of LIR-2, and compared this to the existing high-resolution crystal structure of LIR-1.ResultsRecombinant LIR-2 (domains 1 and 2) was produced in E. coli and crystallized using streak seeding to optimize the crystal morphology. A data set complete to 1.8 Å was collected at 100 K from a single crystal in the P41212 spacegroup. The structure was solved by molecular replacement, using a search model based on the LIR-1 structure.ConclusionsThe overall structure of LIR-2 D1D2 resembles both LIR-1, and Killer Inhibitory Receptors, in that the A strand in each domain forms hydrogen bonds to both β sheets, and there is a sharp angle between the two immunoglobulin-like domains. However, differences from LIR-1 are observed in each domain, with two key changes apparent in the ligand-binding domain, D1. The region corresponding to the residue 44–57 helix of LIR-1 adopts a topology distinct from that of both LIR-1 and the KIR structures, involving a shortened 310 helix. Secondly, the predicted UL18 binding region of LIR-1 is altered substantially in LIR-2: the 76–84 loop mainchain is displaced 11 Å with respect to LIR-1, and Tyrosine 38 adopts an alternative rotamer conformation. In summary, the structure of LIR-2 has revealed significant differences to LIR-1, including ones that may help to explain the >1000-fold lower affinity of LIR-2 for UL18.

Crystallographic Trapping of Heme Loss Intermediates during the Nitrite-Induced Degradation of Human Hemoglobin.

Heme is an important cofactor in a large number of essential proteins and is often involved in small molecule binding and activation. Loss of heme from proteins thus negatively affects the function of these proteins but is also an important component of iron recycling. The characterization of intermediates that form during the loss of heme from proteins has been problematic, in a large part, because of the instability of such intermediates. We have characterized, by X-ray crystallography, three compounds that form during the nitrite-induced degradation of human α(2)β(2) hemoglobin (Hb). The first is an unprecedented complex that exhibits a large β heme displacement of 4.8 Å toward the protein exterior; the heme displacement is stabilized by the binding of the distal His residue to the heme Fe, which in turn allows for the unusual binding of an exogenous ligand on the proximal face of the heme. We have also structurally characterized complexes that display regiospecific nitration of the heme at the 2-vinyl position; we show that heme nitration is not a prerequisite for heme loss. Our results provide structural insight into a possible pathway for nitrite-induced loss of heme from human Hb.

Structure of human R-state aquomethemoglobin at 2.0 Å resolution.

The crystal structure of tetrameric (αβ)(2) R-state human adult aquomethemoglobin is reported at 2.0 Å resolution. The asymmetric unit contained one αβ subunit pair. The R-state crystal belonged to space group P4(1)2(1)2, with unit-cell parameters a = b = 53.6, c = 192.8 Å. An Fe-bound water molecule was modeled into the heme distal pockets of each of the α and β subunits. In the α subunit, a highly ordered liganded water was modeled with an Fe-O(water) distance of 2.2 Å and appears to be protected against escape from the distal pocket by the conformation of the heme propionate groups, which point upwards towards the distal His58 residue aided by a hydrogen-bonding network involving the solvent. In the β subunit, the liganded water exhibited greater motion and was modeled with a longer Fe-O(water) distance of 2.5 Å; in this subunit both propionate groups point downwards away from the distal His63 residue, presumably allowing greater motion of the liganded water in and out of the distal pocket.

Structure of Escherichia coli AdhP (ethanol-inducible dehydrogenase) with bound NAD.

The crystal structure of AdhP, a recombinantly expressed alcohol dehydrogenase from Escherichia coli K-12 (substrain MG1655), was determined to 2.01 Å resolution. The structure, which was solved using molecular replacement, also included the structural and catalytic zinc ions and the cofactor nicotinamide adenine dinucleotide (NAD). The crystals belonged to space group P21, with unit-cell parameters a = 68.18, b = 118.92, c = 97.87 Å, β = 106.41°. The final R factor and Rfree were 0.138 and 0.184, respectively. The structure of the active site of AdhP suggested a number of residues that may participate in a proton relay, and the overall structure of AdhP, including the coordination to structural and active-site zinc ions, is similar to those of other tetrameric alcohol dehydrogenase enzymes.

Distal pocket control of nitrite binding in myoglobin

An important function of heme protein distal pockets is to orient substrates to allow selective reactions involving heme iron. The mammalian muscle protein myoglobin (Mb) has been utilized successfully as a model for distal pocket control of heme protein structure–function relationships. Herein, we unambiguously demonstrate that the Mb active site pocket alters the preferred binding mode of a ligand from what is normally observed in its bioinorganic non-protein model. The interaction of the nitrite anion (NO2 ) with Mb has a long and fascinating history. The curing of meat with nitrite, a century-old process, results primarily from the reaction of Mb with nitrite to generate the red-pink pigment derived from MbNO that is appealing in color to the meat purchaser. 4] The ability of Mb to reduce nitrite to bioactive nitric oxide, a vasodilator for blood pressure control, received renewed attention when it was shown that nitrite was protective against heart attacks in Mb mice, but not in Mb / knockout mice, implicating this heme protein as a genuine nitrite reductase (NiR) protein under hypoxic conditions. Although it had been assumed for decades that nitrite binds to Mb through the N-binding mode, the X-ray crystal structure of the ferric Mb–nitrite adduct revealed an unexpected O-binding mode (Mb(ONO)) that was, at the time, the sole example of nitrite O-binding to any heme protein.

Dioxane contributes to the altered conformation and oligomerization state of a designed engrailed homeodomain variant

Our goal was to compute a stable, full‐sequence design of the Drosophila melanogaster engrailed homeodomain. Thermal and chemical denaturation data indicated the design was significantly more stable than was the wild‐type protein. The data were also nearly identical to those for a similar, later full‐sequence design, which was shown by NMR to adopt the homeodomain fold: a three‐helix, globular monomer. However, a 1.65 Å crystal structure of the design described here turned out to be of a completely different fold: a four‐helix, rodlike tetramer. The crystallization conditions included ∼25% dioxane, and subsequent experiments by circular dichroism and sedimentation velocity analytical ultracentrifugation indicated that dioxane increases the helicity and oligomerization state of the designed protein. We attribute at least part of the discrepancy between the target fold and the crystal structure to the presence of a high concentration of dioxane.

Using Molecular Dynamics Simulations as an Aid in the Prediction of Domain Swapping of Computationally Designed Protein Variants

In standard implementations of computational protein design, a positive-design approach is used to predict sequences that will be stable on a given backbone structure. Possible competing states are typically not considered, primarily because appropriate structural models are not available. One potential competing state, the domain-swapped dimer, is especially compelling because it is often nearly identical with its monomeric counterpart, differing by just a few mutations in a hinge region. Molecular dynamics (MD) simulations provide a computational method to sample different conformational states of a structure. Here, we tested whether MD simulations could be used as a post-design screening tool to identify sequence mutations leading to domain-swapped dimers. We hypothesized that a successful computationally designed sequence would have backbone structure and dynamics characteristics similar to that of the input structure and that, in contrast, domain-swapped dimers would exhibit increased backbone flexibility and/or altered structure in the hinge-loop region to accommodate the large conformational change required for domain swapping. While attempting to engineer a homodimer from a 51-amino-acid fragment of the monomeric protein engrailed homeodomain (ENH), we had instead generated a domain-swapped dimer (ENH_DsD). MD simulations on these proteins showed increased B-factors derived from MD simulation in the hinge loop of the ENH_DsD domain-swapped dimer relative to monomeric ENH. Two point mutants of ENH_DsD designed to recover the monomeric fold were then tested with an MD simulation protocol. The MD simulations suggested that one of these mutants would adopt the target monomeric structure, which was subsequently confirmed by X-ray crystallography.