Search

A fundamental appreciation for how biological macromolecules work requires knowledge of structure and dynamics. Molecular dynamics simulations provide powerful tools for the exploration of the conformational energy landscape accessible to these molecules, and the rapid increase in computational power coupled with improvements in methodology makes this an exciting time for the application of simulation to structural biology. In this Perspective we survey two areas, protein folding and enzymatic catalysis, in which simulations have contributed to a general understanding of mechanism. We also describe results for the [F.sub.1] ATPase molecular motor and the Src family of signaling proteins as examples of applications of simulations to specific biological systems.

A fundamental appreciation for how biological macromolecules work requires knowledge of structure and dynamics. Molecular dynamics simulations provide powerful tools for the exploration of the conformational energy landscape accessible to these molecules, and the rapid increase in computational power coupled with improvements in methodology makes this an exciting time for the application of simulation to structural biology. In this Perspective we survey two areas, protein folding and enzymatic catalysis, in which simulations have contributed to a general understanding of mechanism. We also describe results for the Fi ATPase molecular motor and the Src family of signaling proteins as examples of applications of simulations to specific biological systems. [ABSTRACT FROM AUTHOR]

Variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) challenge currently available coronavirus disease 2019 vaccines and monoclonal antibody therapies through epitope change on the receptor binding domain of the viral spike glycoprotein. Hence, there is a specific urgent need for alternative antivirals that target processes less likely to be affected by mutation, such as the membrane fusion step of viral entry into the host cell. One such antiviral class includes peptide inhibitors, which block formation of the so-called heptad repeat 1 and 2 (HR1HR2) six-helix bundle of the SARS-CoV-2 spike (S) protein and thus interfere with viral membrane fusion. We performed structural studies of the HR1HR2 bundle, revealing an extended, well-folded N-terminal region of HR2 that interacts with the HR1 triple helix. Based on this structure, we designed an extended HR2 peptide that achieves single-digit nanomolar inhibition of SARS-CoV-2 in cell-based and virus-based assays without the need for modifications such as lipidation or chemical stapling. The peptide also strongly inhibits all major SARS-CoV-2 variants to date. This extended peptide is ∼100-fold more potent than all previously published short, unmodified HR2 peptides, and it has a very long inhibition lifetime after washout in virus infection assays, suggesting that it targets a prehairpin intermediate of the SARS-CoV-2 S protein. Together, these results suggest that regions outside the HR2 helical region may offer new opportunities for potent peptide-derived therapeutics for SARS-CoV-2 and its variants, and even more distantly related viruses, and provide further support for the prehairpin intermediate of the S protein.

We reconstitute a phosphotyrosine-mediated protein condensation phase transition of the ∼200 residue cytoplasmic tail of the epidermal growth factor receptor (EGFR) and the adaptor protein, Grb2, on a membrane surface. The phase transition depends on phosphorylation of the EGFR tail, which recruits Grb2, and crosslinking through a Grb2-Grb2 binding interface. The Grb2 Y160 residue plays a structurally critical role in the Grb2-Grb2 interaction, and phosphorylation or mutation of Y160 prevents EGFR:Grb2 condensation. By extending the reconstitution experiment to include the guanine nucleotide exchange factor, SOS, and its substrate Ras, we further find that the condensation state of the EGFR tail controls the ability of SOS, recruited via Grb2, to activate Ras. These results identify an EGFR:Grb2 protein condensation phase transition as a regulator of signal propagation from EGFR to the MAPK pathway.

Son of Sevenless (SOS) is a Ras guanine nucleotide exchange factor (GEF) that plays a central role in numerous cellular signaling pathways. Like many other signaling molecules, SOS is autoinhibited in the cytosol and activates only after recruitment to the membrane. The mean activation time of individual SOS molecules has recently been measured to be ∼60 s, which is unexpectedly long and seemingly contradictory with cellular signaling timescales, which have been measured to be as fast as several seconds. Here, we rectify this discrepancy using a first-passage time analysis to reconstruct the effective signaling timescale of multiple SOS molecules from their single-molecule activation kinetics. Along with corresponding experimental measurements, this analysis reveals how the functional response time, comprised of many slowly activating molecules, can become substantially faster than the average molecular kinetics. This consequence stems from the enzymatic processivity of SOS in a highly out-of-equilibrium reaction cycle during receptor triggering. Ultimately, rare, early activation events dominate the macroscopic reaction dynamics., Competing Interests: The authors declare no competing interest.

Ca 2+ /calmodulin-dependent protein kinase II alpha subunit (CaMKIIα) is a key neuronal signaling protein and an emerging drug target. The central hub domain regulates the activity of CaMKIIα by organizing the holoenzyme complex into functional oligomers, yet pharmacological modulation of the hub domain has never been demonstrated. Here, using a combination of photoaffinity labeling and chemical proteomics, we show that compounds related to the natural substance γ-hydroxybutyrate (GHB) bind selectively to CaMKIIα. By means of a 2.2-Å x-ray crystal structure of ligand-bound CaMKIIα hub, we reveal the molecular details of the binding site deep within the hub. Furthermore, we show that binding of GHB and related analogs to this site promotes concentration-dependent increases in hub thermal stability believed to alter holoenzyme functionality. Selectively under states of pathological CaMKIIα activation, hub ligands provide a significant and sustained neuroprotection, which is both time and dose dependent. This is demonstrated in neurons exposed to excitotoxicity and in a mouse model of cerebral ischemia with the selective GHB analog, HOCPCA (3-hydroxycyclopent-1-enecarboxylic acid). Together, our results indicate a hitherto unknown mechanism for neuroprotection by a highly specific and unforeseen interaction between the CaMKIIα hub domain and small molecule brain-penetrant GHB analogs. This establishes GHB analogs as powerful tools for investigating CaMKII neuropharmacology in general and as potential therapeutic compounds for cerebral ischemia in particular., Competing Interests: Competing interest statement: The University of Copenhagen and Otago Innovation Ltd. have licensed the patent rights for GHB derivatives and their uses (WO/2019/149329) to Ceremedy Ltd., of which B.F., B.R.K., and P.W. are cofounders., (Copyright © 2021 the Author(s). Published by PNAS.)

The transformation of molecular binding events into cellular decisions is the basis of most biological signal transduction. A fundamental challenge faced by these systems is that reliance on protein-ligand chemical affinities alone generally results in poor sensitivity to ligand concentration, endangering the system to error. Here, we examine the lipid-binding pleckstrin homology and Tec homology (PH-TH) module of Bruton's tyrosine kinase (Btk). Using fluorescence correlation spectroscopy (FCS) and membrane-binding kinetic measurements, we identify a phosphatidylinositol (3-5)-trisphosphate (PIP 3 ) sensing mechanism that achieves switch-like sensitivity to PIP 3 levels, surpassing the intrinsic affinity discrimination of PIP 3 :PH binding. This mechanism employs multiple PIP 3 binding as well as dimerization of Btk on the membrane surface. Studies in live cells confirm that mutations at the dimer interface and peripheral site produce effects comparable to that of the kinase-dead Btk in vivo. These results demonstrate how a single protein module can institute an allosteric counting mechanism to achieve high-precision discrimination of ligand concentration. Furthermore, this activation mechanism distinguishes Btk from other Tec family member kinases, Tec and Itk, which we show are not capable of dimerization through their PH-TH modules. This suggests that Btk plays a critical role in the stringency of the B cell response, whereas T cells rely on other mechanisms to achieve stringency., Competing Interests: The authors declare no competing interests.

Upon activation, the epidermal growth factor receptor (EGFR) phosphorylates tyrosine residues in its cytoplasmic tail, which triggers the binding of Src homology 2 (SH2) and phosphotyrosine-binding (PTB) domains and initiates downstream signaling. The sequences flanking the tyrosine residues (referred to as "phosphosites") must be compatible with phosphorylation by the EGFR kinase domain and the recruitment of adapter proteins, while minimizing phosphorylation that would reduce the fidelity of signal transmission. To understand how phosphosite sequences encode these functions within a small set of residues, we carried out high-throughput mutational analysis of three phosphosite sequences in the EGFR tail. We used bacterial surface display of peptides coupled with deep sequencing to monitor phosphorylation efficiency and the binding of the SH2 and PTB domains of the adapter proteins Grb2 and Shc1, respectively. We found that the sequences of phosphosites in the EGFR tail are restricted to a subset of the range of sequences that can be phosphorylated efficiently by EGFR. Although efficient phosphorylation by EGFR can occur with either acidic or large hydrophobic residues at the -1 position with respect to the tyrosine, hydrophobic residues are generally excluded from this position in tail sequences. The mutational data suggest that this restriction results in weaker binding to adapter proteins but also disfavors phosphorylation by the cytoplasmic tyrosine kinases c-Src and c-Abl. Our results show how EGFR-family phosphosites achieve a trade-off between minimizing off-pathway phosphorylation and maintaining the ability to recruit the diverse complement of effectors required for downstream pathway activation., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

Lipidated small GTPases and their regulators need to bind to membranes to propagate actions in the cell, but an integrated understanding of how the lipid bilayer exerts its effect has remained elusive. Here we focused on ADP ribosylation factor (Arf) GTPases, which orchestrate a variety of regulatory functions in lipid and membrane trafficking, and their activation by the guanine-nucleotide exchange factor (GEF) Brag2, which controls integrin endocytosis and cell adhesion and is impaired in cancer and developmental diseases. Biochemical and structural data are available that showed the exceptional efficiency of Arf activation by Brag2 on membranes. We determined the high-resolution crystal structure of unbound Brag2 containing the GEF (Sec7) and membrane-binding (pleckstrin homology) domains, revealing that it has a constitutively active conformation. We used this structure to analyze the interaction of uncomplexed Brag2 and of the myristoylated Arf1/Brag2 complex with a phosphatidylinositol bisphosphate (PIP 2 ) -containing lipid bilayer, using coarse-grained molecular dynamics. These simulations revealed that the system forms a close-packed, oriented interaction with the membrane, in which multiple PIP 2 lipids bind the canonical lipid-binding site and unique peripheral sites of the PH domain, the Arf GTPase and, unexpectedly, the Sec7 domain. We cross-validated these predictions by reconstituting the binding and kinetics of Arf and Brag2 in artificial membranes. Our coarse-grained structural model thus suggests that the high efficiency of Brag2 requires interaction with multiple lipids and a well-defined orientation on the membrane, resulting in a local PIP 2 enrichment, which has the potential to signal toward the Arf pathway., Competing Interests: The authors declare no conflict of interest.

The assembly of cell surface receptors with downstream signaling molecules is a commonly occurring theme in multiple signaling systems. However, little is known about how these assemblies modulate reaction kinetics and the ultimate propagation of signals. Here, we reconstitute phosphotyrosine-mediated assembly of extended linker for the activation of T cells (LAT):growth factor receptor-bound protein 2 (Grb2):Son of Sevenless (SOS) networks, derived from the T-cell receptor signaling system, on supported membranes. Single-molecule dwell time distributions reveal two, well-differentiated kinetic species for both Grb2 and SOS on the LAT assemblies. The majority fraction of membrane-recruited Grb2 and SOS both exhibit fast kinetics and single exponential dwell time distributions, with average dwell times of hundreds of milliseconds. The minor fraction exhibits much slower kinetics, extending the dwell times to tens of seconds. Considering this result in the context of the multistep process by which the Ras GEF (guanine nucleotide exchange factor) activity of SOS is activated indicates that kinetic stabilization from the LAT assembly may be important. This kinetic proofreading effect would additionally serve as a stochastic noise filter by reducing the relative probability of spontaneous SOS activation in the absence of receptor triggering. The generality of receptor-mediated assembly suggests that such effects may play a role in multiple receptor proximal signaling processes.

Heme-nitric oxide/oxygen (H-NOX) binding domains are a recently discovered family of heme-based gas sensor proteins that are conserved across eukaryotes and bacteria. Nitric oxide (NO) binding to the heme cofactor of H-NOX proteins has been implicated as a regulatory mechanism for processes ranging from vasodilation in mammals to communal behavior in bacteria. A key molecular event during NO-dependent activation of H-NOX proteins is rupture of the heme-histidine bond and formation of a five-coordinate nitrosyl complex. Although extensive biochemical studies have provided insight into the NO activation mechanism, precise molecular-level details have remained elusive. In the present study, high-resolution crystal structures of the H-NOX protein from Shewanella oneidensis in the unligated, intermediate six-coordinate and activated five-coordinate, NO-bound states are reported. From these structures, it is evident that several structural features in the heme pocket of the unligated protein function to maintain the heme distorted from planarity. NO-induced scission of the iron-histidine bond triggers structural rearrangements in the heme pocket that permit the heme to relax toward planarity, yielding the signaling-competent NO-bound conformation. Here, we also provide characterization of a nonheme metal coordination site occupied by zinc in an H-NOX protein.

Tyrosine kinase inhibitors (TKIs) represent transformative therapies for several malignancies. Two critical features necessary for maximizing TKI tolerability and response duration are kinase selectivity and invulnerability to resistance-conferring kinase domain (KD) mutations in the intended target. No prior TKI has demonstrated both of these properties. Aiming to maximize selectivity, medicinal chemists have largely sought to create TKIs that bind to an inactive (type II) kinase conformation. Here we demonstrate that the investigational type I TKI crenolanib is a potent inhibitor of Fms tyrosine kinase-3 (FLT3) internal tandem duplication, a validated therapeutic target in human acute myeloid leukemia (AML), as well as all secondary KD mutants previously shown to confer resistance to the first highly active FLT3 TKI quizartinib. Moreover, crenolanib is highly selective for FLT3 relative to the closely related protein tyrosine kinase KIT, demonstrating that simultaneous FLT3/KIT inhibition, a prominent feature of other clinically active FLT3 TKIs, is not required for AML cell cytotoxicity in vitro and may contribute to undesirable toxicity in patients. A saturation mutagenesis screen of FLT3-internal tandem duplication failed to recover any resistant colonies in the presence of a crenolanib concentration well below what has been safely achieved in humans, suggesting that crenolanib has the potential to suppress KD mutation-mediated clinical resistance. Crenolanib represents the first TKI to exhibit both kinase selectivity and invulnerability to resistance-conferring KD mutations, which is unexpected of a type I inhibitor. Crenolanib has significant promise for achieving deep and durable responses in FLT3-mutant AML, and may have a profound impact upon future medicinal chemistry efforts in oncology.

Interior topological features, such as pockets and channels, have evolved in proteins to regulate biological functions by facilitating the diffusion of biomolecules. Decades of research using the globins as model heme proteins have clearly highlighted the importance of gas pockets around the heme in controlling the capture and release of O(2). However, much less is known about how ligand migration contributes to the diverse functions of other heme protein scaffolds. Heme nitric oxide/oxygen binding (H-NOX) domains are a conserved family of gas-sensing heme proteins with a divergent fold that are critical to numerous signaling pathways. Utilizing X-ray crystallography with xenon, a tunnel network has been shown to serve as a molecular pathway for ligand diffusion. Structure-guided mutagenesis results show that the tunnels have unexpected effects on gas-sensing properties in H-NOX domains. The findings provide insights on how the flux of biomolecules through protein scaffolds modulates protein chemistry.

Membrane-bound Ras is activated by translocation of the Son of Sevenless (SOS) protein to the plasma membrane. SOS is inactive unless Ras is bound to an allosteric site on SOS, and the Dbl homology (DH) and Pleckstrin homology (PH) domains of SOS (the DH-PH unit) block allosteric Ras binding. We showed previously that the activity of SOS at the membrane increases with the density of PIP(2) and the local concentration of Ras-GTP, which synergize to release the DH-PH unit. Here we present a new crystal structure of SOS that contains the N-terminal histone domain in addition to the DH-PH unit and the catalytic unit (SOS(HDFC), residues 1-1049). The structure reveals that the histone domain plays a dual role in occluding the allosteric site and in stabilizing the autoinhibitory conformation of the DH-PH unit. Additional insight is provided by kinetic analysis of the activation of membrane-bound Ras by mutant forms of SOS that contain mutations in the histone and the PH domains (E108K, C441Y, and E433K) that are associated with Noonan syndrome, a disease caused by hyperactive Ras signaling. Our results indicate that the histone domain and the DH-PH unit are conformationally coupled, and that the simultaneous engagement of the membrane by a PH domain PIP(2)-binding interaction and electrostatic interactions between a conserved positively charged patch on the histone domain and the negatively charged membrane coincides with a productive reorientation of SOS at the membrane and increased accessibility of both Ras binding sites on SOS.

The kinase domain of human epidermal growth factor receptor (HER) 3/ErbB3, a member of the EGF receptor (EGFR) family, lacks several residues that are critical for catalysis. Because catalytic activity in EGFR family members is switched on by an allosteric interaction between kinase domains in an asymmetric kinase domain dimer, HER3 might be specialized to serve as an activator of other EGFR family members. We have determined the crystal structure of the HER3 kinase domain and show that it appears to be locked into an inactive conformation that resembles that of EGFR and HER4. Although the crystal structure shows that the HER3 kinase domain binds ATP, we confirm that it is catalytically inactive but can serve as an activator of the EGFR kinase domain. The HER3 kinase domain forms a dimer in the crystal, mediated by hydrophobic contacts between the N-terminal lobes of the kinase domains. This N-lobe dimer closely resembles a dimer formed by inactive HER4 kinase domains in crystal structures determined previously, and molecular dynamics simulations suggest that the HER3 and HER4 N-lobe dimers are stable. The kinase domains of HER3 and HER4 form similar chains in their respective crystal lattices, in which N-lobe dimers are linked together by reciprocal exchange of C-terminal tails. The conservation of this tiling pattern in HER3 and HER4, which is the closest evolutionary homolog of HER3, might represent a general mechanism by which this branch of the HER receptors restricts ligand-independent formation of active heterodimers with other members of the EGFR family.

The delivery of signals from the activated T cell antigen receptor (TCR) inside the cell relies on the protein tyrosine kinase ZAP-70 (zeta-associated protein of 70 kDa). A recent crystal structure of inactive full-length ZAP-70 suggests that a central interface formed by the docking of the two SH2 domains of ZAP-70 onto the kinase domain is crucial for suppressing catalytic activity. Here we validate the significance of this autoinhibitory interface for the regulation of ZAP-70 catalytic activity and the T cell response. For this purpose, we perform in vitro catalytic activity assays and binding experiments using ZAP-70 proteins purified from insect cells to examine activation of ZAP-70. Furthermore, we use cell lines stably expressing wild-type or mutant ZAP-70 to monitor proximal events in T cell signaling, including TCR-induced phosphorylation of ZAP-70 substrates, activation of the MAP kinase pathway, and intracellular Ca(2+) levels. Taken together, our results directly correlate the stability of the autoinhibitory interface with the activation of these key events in the T cell response.

In many protein kinases, a characteristic conformational change (the "DFG flip") connects catalytically active and inactive conformations. Many kinase inhibitors--including the cancer drug imatinib--selectively target a specific DFG conformation, but the function and mechanism of the flip remain unclear. Using long molecular dynamics simulations of the Abl kinase, we visualized the DFG flip in atomic-level detail and formulated an energetic model predicting that protonation of the DFG aspartate controls the flip. Consistent with our model's predictions, we demonstrated experimentally that the kinetics of imatinib binding to Abl kinase have a pH dependence that disappears when the DFG aspartate is mutated. Our model suggests a possible explanation for the high degree of conservation of the DFG motif: that the flip, modulated by electrostatic changes inherent to the catalytic cycle, allows the kinase to access flexible conformations facilitating nucleotide binding and release.

DNA polymerases attach to the DNA sliding clamp through a common overlapping binding site. We identify a small-molecule compound that binds the protein-binding site in the Escherichia coli beta-clamp and differentially affects the activity of DNA polymerases II, III, and IV. To understand the molecular basis of this discrimination, the cocrystal structure of the chemical inhibitor is solved in complex with beta and is compared with the structures of Pol II, Pol III, and Pol IV peptides bound to beta. The analysis reveals that the small molecule localizes in a region of the clamp to which the DNA polymerases attach in different ways. The results suggest that the small molecule may be useful in the future to probe polymerase function with beta, and that the beta-clamp may represent an antibiotic target.

The Ras-specific guanine nucleotide-exchange factors Son of sevenless (Sos) and Ras guanine nucleotide-releasing factor 1 (RasGRF1) transduce extracellular stimuli into Ras activation by catalyzing the exchange of Ras-bound GDP for GTP. A truncated form of RasGRF1 containing only the core catalytic Cdc25 domain is sufficient for stimulating Ras nucleotide exchange, whereas the isolated Cdc25 domain of Sos is inactive. At a site distal to the catalytic site, nucleotide-bound Ras binds to Sos, making contacts with the Cdc25 domain and with a Ras exchanger motif (Rem) domain. This allosteric Ras binding stimulates nucleotide exchange by Sos, but the mechanism by which this stimulation occurs has not been defined. We present a crystal structure of the Rem and Cdc25 domains of Sos determined at 2.0-A resolution in the absence of Ras. Differences between this structure and that of Sos bound to two Ras molecules show that allosteric activation of Sos by Ras occurs through a rotation of the Rem domain that is coupled to a rotation of a helical hairpin at the Sos catalytic site. This motion relieves steric occlusion of the catalytic site, allowing substrate Ras binding and nucleotide exchange. A structure of the isolated RasGRF1 Cdc25 domain determined at 2.2-A resolution, combined with computational analyses, suggests that the Cdc25 domain of RasGRF1 is able to maintain an active conformation in isolation because the helical hairpin has strengthened interactions with the Cdc25 domain core. These results indicate that RasGRF1 lacks the allosteric activation switch that is crucial for Sos activity.

The Ras-specific nucleotide exchange factor son of sevenless (SOS) is a large, multidomain protein with complex regulation, including a Ras-dependent allosteric mechanism. The N-terminal segment of SOS, the histone domain, contains two histone folds, which is highly unusual for a cytoplasmic protein. Using a combination of computational docking, small-angle x-ray scattering, mutagenesis, and calorimetry, we show that the histone domain folds into the rest of SOS and docks onto a helical linker that connects the pleckstrin-homology (PH) and Dbl-homology (DH) domains of SOS to the catalytic domain. In this model, a positively charged surface region on the histone domain is positioned so as to provide a fourth potential anchorage site on the membrane for SOS in addition to the PH domain, the allosteric Ras molecule, and the C-terminal adapter-binding site. The histone domain in SOS interacts with the helical linker, using a region of the surface that in nucleosomes is involved in histone tetramerization. Adjacent surface elements on the histone domain that correspond to the DNA-binding surface of nucleosomes form the predicted interaction site with the membrane. The orientation and position of the histone domain in the SOS model implicates it as a potential mediator of membrane-dependent activation signals.

Sliding clamps are ring-like multimeric proteins that encircle duplex DNA and serve as mobile DNA-bound platforms that are essential for efficient DNA replication and repair. Sliding clamps are placed on DNA by clamp loader complexes, in which the clamp-interacting elements are organized in a right-handed spiral assembly. To understand how the flat, ring-like clamps might interact with the spiral interaction surface of the clamp loader complex, we have performed molecular dynamics simulations of sliding clamps (proliferating cell nuclear antigen from the budding yeast, humans, and an archaeal species) in which we have removed one of the three subunits so as to release the constraint of ring closure. The simulations reveal significant structural fluctuations corresponding to lateral opening and out-of-plane distortions of the clamp, which result principally from bending and twisting of the beta-sheets that span the intermolecular interfaces, with smaller but similar contributions from beta-sheets that span the intramolecular interfaces within each subunit. With the integrity of these beta-sheets intact, the predominant fluctuations seen in the simulations are oscillations between lateral openings and right-handed spirals. The tendency for clamps to adopt a right-handed spiral conformation implies that once opened, the conformation of the clamp can easily match the spiraling of clamp loader subunits, a feature that is intrinsic to the recognition of DNA and subsequent hydrolysis of ATP by the clamp-bound clamp loader complex.

Clamp-loader complexes are heteropentameric AAA+ ATPases that load sliding clamps onto DNA. The structure of the nucleotide-free Escherichia coli clamp loader had been determined previously and led to the proposal that the clamp-loader cycles between an inactive state, in which the ATPase domains form a closed ring, and an active state that opens up to form a "C" shape. The crystal structure was interpreted as being closer to the active state than the inactive state. The crystal structure of a nucleotide-bound eukaryotic clamp loader [replication factor C (RFC)] revealed a different and more tightly packed spiral organization of the ATPase domains, raising questions about the significance of the conformation seen earlier for the bacterial clamp loader. We describe crystal structures of the E. coli clamp-loader complex bound to the ATP analog ATPgammaS (at a resolution of 3.5 A) and ADP (at a resolution of 4.1 A). These structures are similar to that of the nucleotide-free clamp-loader complex. Only two of the three functional ATP-binding sites are occupied by ATPgammaS or ADP in these structures, and the bound nucleotides make no interfacial contacts in the complex. These results, along with data from isothermal titration calorimetry, molecular dynamics simulations, and comparison with the RFC structure, suggest that the more open form of the E. coli clamp loader described earlier and in the present work corresponds to a stable inactive state of the clamp loader in which the ATPase domains are prevented from engaging the clamp in the highly cooperative manner seen in the fully ATP-loaded RFC-clamp structure.

Soluble guanylate cyclases are nitric oxide-responsive signaling proteins in which the nitric oxide sensor is a heme-binding domain of unknown structure that we have termed the heme-NO and oxygen binding (H-NOX) domain. H-NOX domains are also found in bacteria, either as isolated domains, or are fused through a membrane-spanning region to methyl-accepting chemotaxis proteins. We have determined the crystal structure of an oxygen-binding H-NOX domain of one such signaling protein from the obligate anaerobe Thermoanaerobacter tengcongensis at 1.77-angstroms resolution, revealing a protein fold unrelated to known structures. Particularly striking is the structure of the protoporphyrin IX group, which is distorted from planarity to an extent not seen before in protein-bound heme groups. Comparison of the structure of the H-NOX domain in two different crystal forms suggests a mechanism whereby alteration in the degree of distortion of the heme group is coupled to changes on the molecular surface of the H-NOX domain and potentially to changes in intermolecular interactions., (Copyright 2004 The National Academy of Sciencs of the USA)

Glutathione reductase (EC 1.6.4.2; CAS registry number 9001-48-3) and trypanothione reductase (CAS registry number 102210-35-5), which are related flavoprotein disulfide oxidoreductases, have marked specificities for glutathione and trypanothione, respectively. A combination of primary sequence alignments and molecular modeling, together with the high-resolution crystal structure of human glutathione reductase, identified certain residues as potentially being responsible for substrate discrimination. Site-directed mutagenesis of Escherichia coli glutathione reductase was used to test these predictions. The mutation of Asn-21 to Arg demonstrated that this single change was insufficient to generate the greater discrimination against trypanothione shown by human glutathione reductase compared with the E. coli enzyme. However, the mutation of Ala-18, Asn-21, and Arg-22 to the amino acid residues (Glu, Trp, and Asn, respectively) in corresponding positions in Trypanosoma congolense trypanothione reductase confirmed that this region of polypeptide chain is intimately involved in substrate recognition. It led to a mutant form of E. coli glutathione reductase that possessed essentially no activity with glutathione but that was able to catalyze trypanothione reduction with a kcat/Km value that was 10% of that measured for natural trypanothione reductases. These results should be of considerable importance in the design of trypanocidal drugs targeted at the differences between glutathione and trypanothione metabolism in trypanosomatids and their hosts.

Trypanosomes and related protozoan parasites lack glutathione reductase and possess instead a closely related enzyme that serves as the reductant of a bis(glutathione)-spermidine conjugate, trypanothione. The human and parasite enzymes have mutually exclusive substrate specificities, providing a route for the design of therapeutic agents by specific inhibition of the parasite enzyme. We report here the three-dimensional structure of trypanothione reductase from Crithidia fasciculata and show that it closely resembles the structure of human glutathione reductase. In particular, the core structure surrounding the catalytic machinery is almost identical in the two enzymes. However, significant differences are found at the substrate binding sites. A cluster of basic residues in glutathione reductase is replaced by neutral, hydrophobic, or acidic residues in trypanothione reductase, consistent with the nature of the spermidine linkage and the change in overall charge of the substrate from -2 to +1, respectively. The binding site is more open in trypanothione reductase due to rotations of about 4 degrees in the domains that form the site, with relative shifts of as much as 2-3 A in residue positions. These results provide a detailed view of the residues that can interact with potential inhibitors and complement previous modeling and mutagenesis studies on the two enzymes.

The extent to which the librations of rigid molecules can model the crystallographic temperature factor profiles of proteins has been examined. For all proteins considered, including influenza virus hemagglutinin, glutathione reductase, myohemerythrin, myoglobin, and streptavidin, a simple 10-parameter model [V. Schomaker and K. N. Trueblood (1968) Acta Crystallogr. Sect. B 24, 63-76] is found to reproduce qualitatively the patterns of maxima and minima in the isotropic backbone meansquare displacements. Large deviations between the rigid molecule and individual atomic temperature factors are found to be correlated with a region in hemagglutinin for which the refined structural model is unsatisfactory and with errors in the structure in a partially incorrect model of myohemerythrin. For the high-resolution glutathione reductase structure, better results are obtained on treating each of the compact domains in the structure as independent rigid bodies. The method allows for the refinement of reliable temperature factors with the introduction of minimal parameters and may prove useful for the evaluation of models in the early stages of x-ray structure refinement. While these results by themselves do not establish the nature of the underlying displacements, the success of the rigid protein model in reproducing qualitative features of temperature factor profiles suggests that rigid body refinement results should be considered in any interpretation of crystallographic thermal parameters.