Peak refreezing in the Greenland firn layer under future warming scenarios.

Firn (compressed snow) covers approximately 90[Formula: see text] of the Greenland ice sheet (GrIS) and currently retains about half of rain and meltwater through refreezing, reducing runoff and subsequent mass loss. The loss of firn could mark a tipping point for sustained GrIS mass loss, since decades to centuries of cold summers would be required to rebuild the firn buffer. Here we estimate the warming required for GrIS firn to reach peak refreezing, using 51 climate simulations statistically downscaled to 1 km resolution, that project the long-term firn layer evolution under multiple emission scenarios (1850-2300). We predict that refreezing stabilises under low warming scenarios, whereas under extreme warming, refreezing could peak and permanently decline starting in southwest Greenland by 2100, and further expanding GrIS-wide in the early 22[Formula: see text] century. After passing this peak, the GrIS contribution to global sea level rise would increase over twenty-fold compared to the last three decades.

Description and Demonstration of the Coupled Community Earth System Model v2 - Community Ice Sheet Model v2 (CESM2-CISM2).

Earth system/ice-sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea-level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy-balance-based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice-sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO2 scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions.

Projected land ice contributions to twenty-first-century sea level rise.

The land ice contribution to global mean sea level rise has not yet been predicted1 using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models2-8, but primarily used previous-generation scenarios9 and climate models10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

An Efficient Ice Sheet/Earth System Model Spin-up Procedure for CESM2-CISM2: Description, Evaluation, and Broader Applicability.

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin-up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet-enabled simulations with CESM2-CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6).

Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison.

Earlier large-scale Greenland ice sheet sea-level projections (e.g., those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of the initMIP-Greenland intercomparison exercise is to compare, evaluate and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of 1) the initial present-day state of the ice sheet and 2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly), and should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap, but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.

Rising Oceans Guaranteed: Arctic Land Ice Loss and Sea Level Rise.

PURPOSE OF REVIEW: This paper reviews sea level contributions from land ice across the Arctic, including Greenland. We summarize ice loss measurement methods, ice loss mechanisms, and recent observations and projections, and highlight research advances over the last 3-5 years and remaining scientific challenges. RECENT FINDINGS: Mass loss across the Arctic began to accelerate during the late twentieth century, with projections of continued loss across all future greenhouse gas emission scenarios. Recent research has improved knowledge of ice hydrology and surface processes, influences of atmospheric and oceanic changes on land ice, and boundary conditions such as subglacial topography. New computer models can also more accurately simulate glacier and ice sheet evolution. SUMMARY: Rapid Arctic ice loss is underway, and future ice loss and sea level rise are guaranteed. Research continues to better understand and model physical processes and to improve projections of ice loss rates, especially after 2050.

An ice sheet model validation framework for the Greenland ice sheet.

We propose a new ice sheet model validation framework - the Cryospheric Model Comparison Tool (CmCt) - that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013 using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quantitative metrics for use in evaluating the different model simulations against the observations. We find that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin- and whole-ice-sheet scale metrics, we find that simulations using both idealized conceptual models and dynamic, numerical models provide an equally reasonable representation of the ice sheet surface (mean elevation differences of <1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CmCt, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate predictive skill with respect to observed dynamic changes occurring on Greenland over the past few decades. An extensible design will allow for continued use of the CmCt as future altimetry, gravimetry, and other remotely sensed data become available for use in ice sheet model validation.

Ice Sheet Model Intercomparison Project (ISMIP6) contribution to CMIP6.

Reducing the uncertainty in the past, present and future contribution of ice sheets to sea-level change requires a coordinated effort between the climate and glaciology communities. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the Greenland and Antarctic Ice Sheets. In this paper, we describe the framework for ISMIP6 and its relationship to other activities within CMIP6. The ISMIP6 experimental design relies on CMIP6 climate models and includes, for the first time within CMIP, coupled ice sheet - climate models as well as standalone ice sheet models. To facilitate analysis of the multi-model ensemble and to generate a set of standard climate inputs for standalone ice sheet models, ISMIP6 defines a protocol for all variables related to ice sheets. ISMIP6 will provide a basis for investigating the feedbacks, impacts, and sea-level changes associated with dynamic ice sheets and for quantifying the uncertainty in ice-sheet-sourced global sea-level change.

Structure and mechanisms of Escherichia coli aspartate transcarbamoylase.

Enzymes catalyze a particular reaction in cells, but only a few control the rate of this reaction and the metabolic pathway that follows. One specific mechanism for such enzymatic control of a metabolic pathway involves molecular feedback, whereby a metabolite further down the pathway acts at a unique site on the control enzyme to alter its activity allosterically. This regulation may be positive or negative (or both), depending upon the particular system. Another method of enzymatic control involves the cooperative binding of the substrate, which allows a large change in enzyme activity to emanate from only a small change in substrate concentration. Allosteric regulation and homotropic cooperativity are often known to involve significant conformational changes in the structure of the protein. Escherichia coli aspartate transcarbamoylase (ATCase) is the textbook example of an enzyme that regulates a metabolic pathway, namely, pyrimidine nucleotide biosynthesis, by feedback control and by the cooperative binding of the substrate, L-aspartate. The catalytic and regulatory mechanisms of this enzyme have been extensively studied. A series of X-ray crystal structures of the enzyme in the presence and absence of substrates, products, and analogues have provided details, at the molecular level, of the conformational changes that the enzyme undergoes as it shifts between its low-activity, low-affinity form (T state) to its high-activity, high-affinity form (R state). These structural data provide insights into not only how this enzyme catalyzes the reaction between l-aspartate and carbamoyl phosphate to form N-carbamoyl-L-aspartate and inorganic phosphate, but also how the allosteric effectors modulate this activity. In this Account, we summarize studies on the structure of the enzyme and describe how these structural data provide insights into the catalytic and regulatory mechanisms of the enzyme. The ATCase-catalyzed reaction is regulated by nucleotide binding some 60 Å from the active site, inducing structural alterations that modulate catalytic activity. The delineation of the structure and function in this particular model system will help in understanding the molecular basis of cooperativity and allosteric regulation in other systems as well.

Boron atoms as spin carriers in two- and three-dimensional systems.

Paramagnetic compounds with at least partially boron-centered electron spin can be constructed using either the prototypically electron-accepting boron atoms bridged by planar pi-conjugated organic systems, or by taking advantage of the three-dimensional delocalized bonding in oligonuclear borane, haloborane, or carborane clusters. The concept of mixed valency can thus be transferred from organic and transition-metal chemistry to main-group-element molecules, and density functional theory is capable of reproducing the variable spin distribution.

Structure-guided design of AMP mimics that inhibit fructose-1,6-bisphosphatase with high affinity and specificity.

AMP binding sites are commonly used by nature for allosteric regulation of enzymes controlling the production and metabolism of carbohydrates and lipids. Since many of these enzymes represent potential drug targets for metabolic diseases, efforts were initiated to discover AMP mimics that bind to AMP-binding sites with high affinity and high enzyme specificity. Herein we report the structure-guided design of potent fructose 1,6-bisphosphatase (FBPase) inhibitors that interact with the AMP binding site on FBPase despite their structural dissimilarity to AMP. Molecular modeling, free-energy perturbation calculations, X-ray crystallography, and enzyme kinetic data guided our redesign of AMP, which began by replacing the 5'-phosphate with a phosphonic acid attached to C8 of the adenine base via a 3-atom spacer. Additional binding affinity was gained by replacing the ribose with an alkyl group that formed van der Waals interactions with a hydrophobic region within the AMP binding site and by replacing the purine nitrogens N1 and N3 with carbons to minimize desolvation energy expenditures. The resulting benzimidazole phosphonic acid, 16, inhibited human FBPase (IC50 = 90 nM) 11-fold more potently than AMP and exhibited high specificity for the AMP binding site on FBPase. 16 also inhibited FBPase in primary rat hepatocytes and correspondingly resulted in concentration-dependent inhibition of the gluconeogenesis pathway. Accordingly, these results suggest that the AMP site of FBPase may represent a potential drug target for reducing the excessive glucose produced by the gluconeogenesis pathway in patients with type 2 diabetes.

Discovery of potent and specific fructose-1,6-bisphosphatase inhibitors and a series of orally-bioavailable phosphoramidase-sensitive prodrugs for the treatment of type 2 diabetes.

Excessive glucose production by the liver coupled with decreased glucose uptake and metabolism by muscle, fat, and liver results in chronically elevated blood glucose levels in patients with type 2 diabetes. Efforts to treat diabetes by reducing glucose production have largely focused on the gluconeogenesis pathway and rate-limiting enzymes within this pathway such as fructose-1,6-bisphosphatase (FBPase). The first potent FBPase inhibitors were identified using a structure-guided drug design strategy (Erion, M. D.; et al. J. Am. Chem. Soc. 2007, 129, 15480-15490) but proved difficult to deliver orally. Herein, we report the synthesis and characterization of a series of orally bioavailable FBPase inhibitors identified following the combined discoveries of a low molecular weight inhibitor series with increased potency and a phosphonate prodrug class suitable for their oral delivery. The lead inhibitor, 10A, was designed with the aid of X-ray crystallography and molecular modeling to bind to the allosteric AMP binding site of FBPase. High potency (IC50 = 16 nM) and FBPase specificity were achieved by linking a 2-aminothiazole with a phosphonic acid. Free-energy perturbation calculations provided insight into the factors that contributed to the high binding affinity. 10A and standard phosphonate prodrugs of 10A exhibited poor oral bioavailability (0.2-11%). Improved oral bioavailability (22-47%) was achieved using phosphonate diamides that convert to the corresponding phosphonic acid by sequential action of an esterase and a phosphoramidase. Oral administration of the lead prodrug, MB06322 (30, CS-917), to Zucker Diabetic Fatty rats led to dose-dependent inhibition of gluconeogenesis and endogenous glucose production and consequently to significant blood glucose reduction.

An eight-subunit COP9 signalosome with an intact JAMM motif is required for fungal fruit body formation.

Fruit body formation in filamentous fungi is a complex and yet hardly understood process. We show here that protein turnover control is crucial for Aspergillus nidulans development. Deletion of genes encoding COP9 signalosome (CSN) subunits 1, 2, 4, or 5 resulted in identical blocks in fruit body formation. The CSN multiprotein complex controls ubiquitin-dependent protein degradation in eukaryotes. Six CSN subunits interacted in a yeast two-hybrid analysis, and the complete eight-subunit CSN was recruited by a functional tandem affinity purification tag fusion of subunit 5 (CsnE). The tagged CsnE was unable to recruit any CSN subunit in a strain deleted for subunit 1 or subunit 4. Mutations in the JAMM metalloprotease core of CsnE resulted in mutant phenotypes identical to those of csn deletion strains. We propose that a correctly assembled CSN including a functional JAMM links protein turnover to fungal sexual development.

Nature of "hydrogen bond" in the diborane-benzene complex: covalent, electrostatic, or dispersive?

Motivated by the recent discovery of unusual "hydrogen bonding"-like interaction between a borane system and benzene molecules in a molecular crystal, we carried out quantum mechanical calculations on a model complex, diborane-benzene cluster. The aim is to understand the nature of this unique interaction, which is expected to play an essential role in this novel class of molecular crystals. As analyzed in the present study, the interaction between diborane and benzene is special in the following aspects: (1) this interaction is mostly dispersive; (2) the observed pseudodirectionality with one of the diborane bridge hydrogen directed toward the benzene centroid minimizes the van der Waals contact; and (3) in the "hydrogen bond" map, this interaction is located in a unique region, which is presently populated by a few known molecular complexes with very different chemical characteristics. It is anticipated that the results from the present analysis will provide meaningful guidance for molecular engineering with diborane-benzene as a building block and for stabilization of this and possible other hydrogen bonds by dispersive contributions.

T-state active site of aspartate transcarbamylase: crystal structure of the carbamyl phosphate and L-alanosine ligated enzyme.

An X-ray diffraction study to 2.0 A resolution shows that this enzyme, ATCase, is in the T-state (the c3 to c3 distance is 45.2 A) when ATCase is bound to carbamyl phosphate (CP) and to L-alanosine (an analogue of aspartate). This result strongly supports the kinetic results that alanosine did not inhibit the carbamylation of aspartate in the normal reaction of native ATCase plus CP and aspartate [Baillon, J., Tauc, P., and Hervé, G. (1985) Biochemistry 24, 7182-7187]. The structure further reveals that the phosphate of CP is 4 A away from its known position in the R-state and is in the T-state position of P(i) in a recent study of ATCase complexed with products, phosphate (P(i)) and N-carbamyl-L-aspartate [Huang, J., and Lipscomb, W. N. (2004) Biochemistry 43, 6422-6426]. Moreover, the alanosine position in this T-state is somewhat displaced from that expected for its analogue, aspartate, from the R-state position. The relations of these structural aspects to the kinetics are presented.

The allosteric mechanism of yeast chorismate mutase: a dynamic analysis.

The effector-regulated allosteric mechanism of yeast chorismate mutase (YCM) was studied by normal mode analysis and targeted molecular dynamics. The normal mode analysis shows that the conformational change between YCM in the R state and in the T state can be represented by a relatively small number of low-frequency modes. This suggests that the transition is coded in the structure and is likely to have a low energetic barrier. Quantitative comparisons (i.e. frequencies) between the low-frequency modes of YCM with and without effectors (modeled structures) reveal that the binding of Trp increases the global flexibility, whereas Tyr decreases global flexibility. The targeted molecular dynamics simulation of substrate analog release from the YCM active site suggests that a series of residues are critical for orienting and "recruiting" the substrate. The simulation led to the switching of a series of substrate-release-coupled salt-bridge partners in the ligand-binding domain; similar changes occur in the transition between YCM R-state and T-state crystal structures. Thus, the normal mode analysis and targeted molecular dynamics results provide evidence that the effectors regulate YCM activity by influencing the global flexibility. The change in flexibility is coupled to the binding of substrate to the T state and release of the product from the R state, respectively.

Evolution of 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase-encoding genes in the yeast Saccharomyces cerevisiae.

The shikimate pathway resulting in three aromatic amino acids is initiated in different organisms by two and three 3-deoxy-d-arabino-heptulosonate-7-phosphate synthases, respectively. Aro3p and Aro4p are the yeast enzymes feedback-inhibited by phenylalanine and tyrosine, respectively. A yeast strain deficient in the general control transcriptional regulatory system of amino acid biosynthesis is unable to live in the presence of high amounts of phenylalanine and tyrosine. Here, we show that this yeast strain can be rescued by the expression of aroH from Escherichia coli encoding the tryptophan-regulated AroH as third isoenzyme. Yeast carrying Ec AroH as the only enzyme for the initial step of the shikimate pathway can grow in the absence of tryptophan. Without aromatic amino acids, this yeast strain survives only when the yeast ARO3 promoter instead of the ARO4 promoter drives E. coli aroH. The detailed analysis of Aro3p and Aro4p revealed a triple feedback control by tyrosine/phenylalanine and tryptophan. Dissecting this control allowed engineering of Aro4p S195A as an enzyme, which is inhibited like AroH only by tryptophan. In addition, Aro4p variants were constructed that show an equally strong inhibition by tyrosine and tryptophan (Aro4p P165G Q302R) and in which the regulation by tyrosine and tryptophan was reversed (Aro4p P165G). Our data suggest that yeast possesses only two instead of three isogenes encoding 3-deoxy-D-arabino-heptulosonate-7-phosphate synthases because both isoenzymes can be fine tuned by tryptophan as additional effector and because transcriptional regulation by the general control system can be induced as backup when aromatic amino acids in the environment are imbalanced.

MB06322 (CS-917): A potent and selective inhibitor of fructose 1,6-bisphosphatase for controlling gluconeogenesis in type 2 diabetes.

In type 2 diabetes, the liver produces excessive amounts of glucose through the gluconeogenesis (GNG) pathway and consequently is partly responsible for the elevated glucose levels characteristic of the disease. In an effort to find safe and efficacious GNG inhibitors, we targeted the AMP binding site of fructose 1,6-bisphosphatase (FBPase). The hydrophilic nature of AMP binding sites and their widespread use for allosteric regulation of enzymes in metabolic pathways has historically made discovery of AMP mimetics suitable for drug development difficult. By using a structure-based drug design strategy, we discovered a series of compounds that mimic AMP but bear little structural resemblance. The lead compound, MB05032, exhibited high potency and specificity for human FBPase. Oral delivery of MB05032 was achieved by using the bisamidate prodrug MB06322 (CS-917), which is converted to MB05032 in two steps through the action of an esterase and a phosphoramidase. MB06322 inhibited glucose production from a variety of GNG substrates in rat hepatocytes and from bicarbonate in male Zucker diabetic fatty rats. Analysis of liver GNG pathway intermediates confirmed FBPase as the site of action. Oral administration of MB06322 to Zucker diabetic fatty rats led to a dose-dependent decrease in plasma glucose levels independent of insulin levels and nutritional status. Glucose lowering occurred without signs of hypoglycemia or significant elevations in plasma lactate or triglyceride levels. The findings suggest that potent and specific FBPase inhibitors represent a drug class with potential to treat type 2 diabetes through inhibition of GNG.

Aspartate transcarbamylase (ATCase) of Escherichia coli: a new crystalline R-state bound to PALA, or to product analogues citrate and phosphate.

Structures of the R-state of Escherichia coli ATCase maintained with carbamyl phosphate and succinate, phosphonoacetamide and malonate, or N-phosphonacetyl-l-aspartate (PALA) have previously been made in the space group P321, in which the two independent r (regulatory) and two independent c (catalytic) chains are repeated by crystallographic symmetry to yield the holoenzyme c(6)r(6), ((c(3))(2)(r(2))(3)). The exploration of a new crystalline R-state P2(1)2(1)2(1) was undertaken to examine the c(3).c(3) expansion of 11 A in the T-to-R transition, and to further test whether intermolecular contacts influence the binding of PALA. The results show that the expansion along the 3-fold axis is 10 A, and that the binding modes of the six crystallographic independent PALA molecules are virtually identical to one another, and to modes described previously. As further test, the PALA, a bisubstrate analogue, was displaced by citrate and phosphate, where citrate is an analogue of product carbamylaspartate. The results support the conclusions about the binding of the three previously studied analogues, and further support, within about 0.5 A, the structure proposed for the transition state [Gouaux, J. E., Krause, K. L., and Lipscomb, W. N. (1987) Biochem. Biophys. Res. Commun. 142, 893-897; Jin, L., Stec, B., Lipscomb, W. N., and Kantrowitz, E. R. (1999) Proteins: Struct., Funct., Genet. 37, 729-742].

Products in the T-state of aspartate transcarbamylase: crystal structure of the phosphate and N-carbamyl-L-aspartate ligated enzyme.

The structure of aspartate transcarbamylase of Escherichia coli ligated to products (phosphate and N-carbamyl-l-aspartate) has been determined at 2.37 A resolution (R-factor = 0.23, R(free) = 0.27). Results might indicate a product release mode, rather than close analogues to the transition state like those found in our earlier studies of other ligands (N-phosphonacetyl-L-aspartate, carbamyl phosphate plus malonate, phosphonoacetamide plus malonate, or citrate plus phosphate). Ordered product release, first carbamylaspartate (CLA) and then phosphate, might be facilitated by a 4 A movement of phosphate from the substrate-analogue position to the product (phosphate) binding position, and by a somewhat similar release movement of the other product (CLA) relative to its analogue (citrate). This movement is consistent with earlier studies of binding of either pyrophosphate or phosphate alone [Honzatko, R. B., and Lipscomb, W. N. (1982) J. Mol. Biol. 160, 265-286].