High-resolution solution structure of the beryllofluoride-activated NtrC receiver domain

Bacterial receiver domains mediate the cellular response to environmental changes through conformational changes induced by phosphorylation of a conserved aspartate residue. While the structures of several activated receiver domains have recently been determined, there is substantial variation in the conformational changes occurring upon activation. Here we present the high-resolution structure of the activated NtrC receiver domain (BeF(3)(-)-NtrC(r) complex) determined using NMR data, including residual dipolar couplings, yielding a family of structures with a backbone rmsd of 0.57 +/- 0.08 A, which is compared with the previous lower-resolution structure of the phosphorylated protein. Both phosphorylation and beryllofluoride addition induce a shift in register and an axial rotation of alpha-helix 4. In this high-resolution structure, we are able to observe a concerted change in the positions of Thr82 and Tyr101; this correlated change in two conserved residues (termed Y-T coupling) has been considered a general feature of the conformational change in receiver domains upon activation. In NtrC, this correlated side chain shift, leading to the helix reorientation, is distinctly different from the smaller reorganization seen in other activated receiver domains, and involves numerous other residues which do not participate in conformational changes seen in the other systems. Titration of the activated receiver domain with peptides from the NtrC ATPase domain provides direct evidence for interactions on the rearranged face of the receiver domain, which are likely to be responsible for enabling assembly into the active aggregate. Analysis of the active structure also suggests that His84 may play a role in controlling the phosphate hydrolysis rate.     

Molecular dynamics simulations of Trp side-chain conformational flexibility in the gramicidin A channel

Gramicidin A (gA) is prototypical peptide antibiotic and a model ion channel former. Configured in the solid-state NMR beta(6.5)-helix channel conformation, gA was subjected to 1-ns molecular dynamics (MD) gas phase simulations using the all-atom charmm22 force field to ascertain the conformational stability of the Trp side chains as governed by backbone and neighboring side-chain contacts. Three microcanonical trajectories were computed using different initial atomic velocities for each of twenty different initial structures. For each set, one of the four Trp side chains in each monomer was initially positioned in one of the five non-native conformations (A. E. Dorigo et al., Biophysical Journal, 1999, Vol. 76, 1897-1908), the other Trps being positioned in the native state, o1. In three additional control simulations, all Trps were initiated in the native conformation. After equilibration, constraints were removed and subsequent conformational changes of the initially constrained Trp were measured. The chi(1) was more flexible than chi(2.1). The energetically optimal orientation, o1 (Dorigo et al., 1999), was the most stable in all four Trp positions (9, 11, 13, 15) and remained unchanged for the entire 1 ns simulation in 19 of 24 trials. Changes in chi(1) from each of the 5 suboptimal states occur readily. Two of the non-native conformations reverted readily to o1, whereas the other three converted to an intermediate state, i2. There were frequent interconversions between i2 and o1. We speculate that experimentally observed Trp stability is caused by interactions with the lipid-water interface, and that stabilization of one of the suboptimal conformations in gA, such as i2, by lipid headgroups could produce a secondary, metastable conformational state. This could explain recent experimental studies of differences in the channel conductance dispersity between gA and a Trp-to-Phe gA analog, gramicidin M (gM, J. C. Markham et al., Biochimica et Biophysica Acta, 2001, Vol. 1513, 185-192).     

Disaccharide conformational flexibility. II. Molecular dynamics simulations of sucrose

Molecular dynamics simulations have been used to study the motions in vacuum of the disaccharide sucrose. Ensembles of trajectories were calculated for each of the five local minimum energy conformations identified in the adiabatic conformational energy mapping of this molecule. The model sucrose molecules were found to exhibit a variety of motions, although the global minimum energy conformation was found to be dynamically stable, and no transitions away from this structure were observed to occur spontaneously. In all but one of these vacuum trajectories, the intramolecular hydrogen bond between residues was maintained, in accord with recent nmr studies of this molecule in aqueous solution. Considerable flexibility of the furanoid ring was found in the trajectories. No "flips" to the opposite puckering for this ring were found in the simulations starting from the global minimum, although such a transition was observed for a trajectory initiated with one of the higher local minimum energy conformations. Overall, the observed structural fluctuations were consistent with the experimental picture of sucrose as a relatively rigid molecule.     

Molecular dynamic simulations of cisplatin- and oxaliplatin-d(GG) intrastand cross-links reveal differences in their conformational dynamics

Mismatch repair proteins, DNA damage-recognition proteins and translesion DNA polymerases discriminate between Pt-GG adducts containing cis-diammine ligands (formed by cisplatin (CP) and carboplatin) and trans-RR-diaminocyclohexane ligands (formed by oxaliplatin (OX)) and this discrimination is thought to be important in determining differences in the efficacy, toxicity and mutagenicity of these platinum anticancer agents. We have postulated that these proteins recognize differences in conformation and/or conformational dynamics of the DNA containing the adducts. We have previously determined the NMR solution structure of OX-DNA, CP-DNA and undamaged duplex DNA in the 5'-d(CCTCAGGCCTCC)-3' sequence context and have shown the existence of several conformational differences in the vicinity of the Pt-GG adduct. Here we have used molecular dynamics simulations to explore differences in the conformational dynamics between OX-DNA, CP-DNA and undamaged DNA in the same sequence context. Twenty-five 10 ns unrestrained fully solvated molecular dynamics simulations were performed starting from two different DNA conformations using AMBER v8.0. All 25 simulations reached equilibrium within 4 ns, were independent of the starting structure and were in close agreement with previous crystal and NMR structures. Our data show that the cis-diammine (CP) ligand preferentially forms hydrogen bonds on the 5' side of the Pt-GG adduct, while the trans-RR-diaminocyclohexane (OX) ligand preferentially forms hydrogen bonds on the 3' side of the adduct. In addition, our data show that these differences in hydrogen bond formation are strongly correlated with differences in conformational dynamics, specifically the fraction of time spent in different DNA conformations in the vicinity of the adduct, for CP- and OX-DNA adducts. We postulate that differential recognition of CP- and OX-GG adducts by mismatch repair proteins, DNA damage-recognition proteins and DNA polymerases may be due, in part, to differences in the fraction of time that the adducts spend in a conformation favorable for protein binding.     

Characterization of the conformational state and flexibility of HIV-1 glycoprotein gp120 core domain

gp120 is key to the human immunodeficiency virus type 1 viral cell entry. Knowledge of the detailed conformational states of gp120 is crucial to intervention, yet the unbound form is still resistant to structural characterization probably because of its flexibility. Toward this goal, we performed molecular dynamics simulations on the wild type gp120 core domain extracted from its ternary crystal structure and on a modeled mutant, S375W, that experimentally has a significantly different phenotype from the wild type. Although the mutant retained a bound-like conformation, the wild type drifted to a different conformational state. The wild type beta strands 2 and 3 of the bridging sheet were very mobile and partially unfolded, and the organization among the inner and outer domains and beta strands 20 and 21 of the bridging sheet, near the mutation site, was more open than in the bound form, although the overall structure was maintained. These differences were apparently a result of the strengthening of the hydrophobic core in the mutant. This stabilization further explains the experimentally significantly different thermodynamic properties between the wild type and the mutant. Taken together, our results suggest that the free form, although different from the bound state, shares many of the bound structural features. The observed loss of freedom near the binding site, rather than the previously hypothesized more dramatic conformational transition from the unbound to the bound state, appears to be the major contributor to the large entropy cost for the CD4 binding to the wild type.     

Membrane-bound 3D structures reveal the intrinsic flexibility of annexin VI

Several quasi-ordered arrays and three two-dimensional crystal forms of annexin VI were obtained on artificial lipid monolayers. Three-dimensional reconstructions of the crystal forms exhibit marked differences in the orientations of the two lobes, revealing flexibility of the linker between the two lobes of annexin VI. Evidence is presented that the lobes may bind the monolayer in a parallel orientation, or an antiparallel orientation, in which the second lobe is turned away from the monolayer. It is hypothesized that annexin VI may also adopt several conformations in vivo, underlying different functional roles.     

Crystal structures of histone and p53 methyltransferase SmyD2 reveal a conformational flexibility of the autoinhibitory C-terminal domain

SmyD2 belongs to a new class of chromatin regulators that control gene expression in heart development and tumorigenesis. Besides methylation of histone H3 K4, SmyD2 can methylate non-histone targets including p53 and the retinoblastoma tumor suppressor. The methyltransferase activity of SmyD proteins has been proposed to be regulated by autoinhibition via the intra- and interdomain bending of the conserved C-terminal domain (CTD). However, there has been no direct evidence of a conformational change in the CTD. Here, we report two crystal structures of SmyD2 bound either to the cofactor product S-adenosylhomocysteine or to the inhibitor sinefungin. SmyD2 has a two-lobed structure with the active site located at the bottom of a deep crevice formed between the CTD and the catalytic domain. By extensive engagement with the methyltransferase domain, the CTD stabilizes the autoinhibited conformation of SmyD2 and restricts access to the catalytic site. Unexpectedly, despite that the two SmyD2 structures are highly superimposable, significant differences are observed in the first two helices of the CTDs: the two helices bend outwards and move away from the catalytic domain to generate a less closed conformation in the sinefungin-bound structure. Although the overall fold of the individual domains is structurally conserved among SmyD proteins, SmyD2 appear to be a conformational "intermediate" between a close form of SmyD3 and an open form of SmyD1. In addition, the structures reveal that the CTD is structurally similar to tetratricopeptide repeats (TPR), a motif through which many cochaperones bind to the heat shock protein Hsp90. Our results thus provide the first evidence for the intradomain flexibility of the TPR-like CTD, which may be important for the activation of SmyD proteins by Hsp90.     

Molecular dynamics simulations on the free and complexed N-terminal SH2 domain of SHP-2

Signal transduction events are often mediated by small protein domains such as SH2 (Src homology 2) domains that recognize phosphotyrosines (pY) and flanking sequences. In case of the SHP-2 receptor tyrosine phosphatase an N-terminal SH2 domain binds and inactivates the phosphatase (PTP) domain. The pY-peptide-binding site on the N-terminal SH2 domain does not overlap with the PTP binding region. Nevertheless, pY-peptide binding causes domain dissociation and phosphatase activation. Comparative multi-nanosecond molecular dynamics simulations on the N-SH2 domain in ligand-bound and free states have been performed to study the allosteric mechanism that leads to domain dissociation upon pY-peptide binding. Significant ligand-dependent differences in the conformational flexibility of regions that are involved in SH2-PTP domain association have been observed. The results support a mechanism of signal transduction where SH2-peptide binding modulates the domain flexibility and reduces its capacity to fit into the entrance of the PTP catalytic domain of SHP-2.     

Coarse-grained dynamics of the receiver domain of NtrC: fluctuations, correlations and implications for allosteric cooperativity

Receiver domains are key molecular switches in bacterial signaling. Structural studies have shown that the receiver domain of the nitrogen regulatory protein C (NtrC) exists in a conformational equilibrium encompassing both inactive and active states, with phosphorylation of Asp54 allosterically shifting the equilibrium towards the active state. To analyze dynamical fluctuations and correlations in NtrC as it undergoes activation, we have applied a coarse-grained dynamics algorithm using elastic network models. Normal mode analysis reveals possible dynamical pathways for the transition of NtrC from the inactive state to the active state. The diagonalized correlation between the inactive and the active (phosphorylated) state shows that most correlated motions occur around the active site of Asp54 and in the region Thr82 to Tyr101. This indicates a coupled correlation of dynamics in the "Thr82-Tyr101" motion. With phosphorylation inducing significant flexibility changes around the active site and alpha3 and alpha4 helices, we find that this activation makes the active-site region and the loops of alpha3/beta4 and alpha4/beta5 more stable. This means that phosphorylation entropically favors the receiver domain in its active state, and the induced conformational changes occur in an allosteric manner. Analyses of the local flexibility and long-range correlated motion also suggest a dynamics criterion for determining the allosteric cooperativity of NtrC, and may be applicable to other proteins.     

Multiple distinct assemblies reveal conformational flexibility in the small heat shock protein Hsp26

Small heat shock proteins are a superfamily of molecular chaperones that suppress protein aggregation and provide protection from cell stress. A key issue for understanding their action is to define the interactions of subunit domains in these oligomeric assemblies. Cryo-electron microscopy of yeast Hsp26 reveals two distinct forms, each comprising 24 subunits arranged in a porous shell with tetrahedral symmetry. The subunits form elongated, asymmetric dimers that assemble via trimeric contacts. Modifications of both termini cause rearrangements that yield a further four assemblies. Each subunit contains an N-terminal region, a globular middle domain, the alpha-crystallin domain, and a C-terminal tail. Twelve of the C termini form 3-fold assembly contacts which are inserted into the interior of the shell, while the other 12 C termini form contacts on the surface. Hinge points between the domains allow a variety of assembly contacts, providing the flexibility required for formation of supercomplexes with non-native proteins.     

Copper-induced conformational changes in the N-terminal domain of the Wilson disease copper-transporting ATPase

The Wilson disease copper-transporting ATPase plays a critical role in the intracellular trafficking of copper. Mutations in this protein lead to the accumulation of a toxic level of copper in the liver, kidney, and brain followed by extensive tissue damage and death. The ATPase has a novel amino-terminal domain ( approximately 70 kDa) which contains six repeats of the copper binding motif GMTCXXC. We have expressed and characterized this domain with respect to the copper binding sites and the conformational consequences of copper binding. A detailed analysis of this domain by X-ray absorption spectroscopy (XAS) has revealed that each binding site ligates copper in the +1 oxidation state using two cysteine side chains with distorted linear geometry. Analysis of copper-induced conformational changes in the amino-terminal domain indicates that both secondary and tertiary structure changes take place upon copper binding. These copper-induced conformational changes could play an important role in the function and regulation of the ATPase in vivo. In addition to providing important insights on copper binding to the protein, these results suggest a possible mechanism of copper trafficking by the Wilson disease ATPase.     

Solution structure of Atg8 reveals conformational polymorphism of the N-terminal domain

During autophagy a crescent shaped like membrane is formed, which engulfs the material that is to be degraded. This membrane grows further until its edges fuse to form the double membrane covered autophagosome. Atg8 is a protein, which is required for this initial step of autophagy. Therefore, a multistage conjugation process of newly synthesized Atg8 to phosphatidylethanolamine is of critical importance. Here we present the high resolution structure of unprocessed Atg8 determined by nuclear magnetic resonance spectroscopy. Its C-terminal subdomain shows a well-defined ubiquitin-like fold with slightly elevated mobility in the pico- to nanosecond timescale as determined by heteronuclear NOE data. In comparison to unprocessed Atg8, cleaved Atg8^G116 shows a decreased mobility behaviour. The N-terminal domain adopts different conformations within the micro- to millisecond timescale. The possible biological relevance of the differences in dynamic behaviours between both subdomains as well as between the cleaved and uncleaved forms is discussed.     

Molecular dynamics simulations of aptamer-binding reveal generalized allostery in thrombin

Thrombin is an attractive target for antithrombotic therapy due to its central role in thrombosis and hemostasis as well as its role in inducing tumor growth, metastasis, and tumor invasion. The thrombin-binding DNA aptamer (TBA), is under investigation for anticoagulant drugs. Although aptamer binding experiments have been revealed various effects on thrombin’s enzymatic activities, the detailed picture of the thrombin’s allostery from TBA binding is still unclear. To investigate thrombin’s response to the aptamer-binding at the molecular level, we compare the mechanical properties and free energy landscapes of the free and aptamer-bound thrombin using microsecond-scale all-atom GPU-based molecular dynamics simulations. Our calculations on residue fluctuations and coupling illustrate the allosteric effects of aptamer-binding at the atomic level, highlighting the exosite II, 60s, γ and the sodium loops, and the alpha helix region in the light chains involved in the allosteric changes. This level of details clarifies the mechanisms of previous experimentally demonstrated phenomena, and provides a prediction of the reduced autolysis rate after aptamer-binding. The shifts in thrombin’s ensemble of conformations and free energy surfaces after aptamer-binding demonstrate that the presence of bound-aptamer restricts the conformational freedom of thrombin suggesting that conformational selection, i.e. generalized allostery, is the dominant mechanism of thrombin-aptamer binding. The profound perturbation on thrombin’s mechanical and thermodynamic properties due to the aptamer-binding, which was revealed comprehensively as a generalized allostery in this work, may be exploited in further drug discovery and development.     

Molecular dynamics simulations of a guaiacyl beta-O-4 lignin model compound: examination of intramolecular hydrogen bonding and conformational flexibility

The dynamical conformational behavior of a guaiacyl beta-O-4 lignin model compound has been investigated by molecular simulations. The potential energy surface of the molecule in vacuum has been examined by means of an adiabatic map, showing a large accessible conformational space with multiple energy minima separated by low barriers. Molecular dynamics simulations have been performed in vacuum and with explicit solvent molecules for 10 and 2.1 ns, respectively. Molecular dynamics trajectories recorded in vacuum have shown the molecule to be flexible and to visit a large number of conformations. Many intramolecular H-bonds have been observed, existing for more than 90% of the total simulation time. The presence of explicit solvent molecules induces a significant broadening of some regions of the accessible conformational space and also largely reduces the statistical significance of intramolecular H-bonding. Intramolecular H-bonds observed in vacuum do not persist significantly and are preferentially exchanged with intermolecular H-bonds to the surrounding solvent molecules. The theoretical results are in good agreement with experimental NMR data that do not support the existence of strong and persistent intramolecular H-bonds in solution but instead indicate that H-bonds to solvent predominate. Finally, both molecular modeling and NMR approaches predict the guaiacyl beta-O-4 structure to be flexible and indicate that intramolecular H-bonds are not strong and persistent enough to confer rigidity to the molecule in solution.     

Solution NMR studies reveal no global flexibility in the catalytic domain of CDC25B

The CDC25B phosphatase is a critical regulator of the cell cycle and has been validated as an important therapeutic target in cancer. Previous studies using molecular dynamics simulations have concluded that the catalytic domain of CDC25B may experience a significant degree of dynamics or be partially disordered in solution, a finding that has a pronounced impact on the structure‐based development of CDC25B inhibitors. We have probed the backbone dynamics of the CDC25B catalytic domain in solution using NMR relaxation experiments and found that the core of the protein is relatively rigid and does not experience any large‐scale dynamics over a broad range of time scales. Furthermore, based on residual dipolar coupling measurements we have concluded that the conformation in solution is very similar to that observed in the crystal form. Importantly, these findings rationalize the application of the CDC25B crystal structure in structure‐based drug development. Proteins 2014; 82:2889–2895. © 2014 Wiley Periodicals, Inc.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

The conformational flexibility of the helicase-like domain from Thermotoga maritima reverse gyrase is restricted by the topoisomerase domain

Reverse gyrase is the only enzyme known to introduce positive supercoils into DNA. Positive supercoiling is achieved by the functional cooperation of a helicase-like and a topoisomerase domain. The isolated helicase-like domain is a DNA-stimulated ATPase, and the isolated topoisomerase domain can relax supercoiled DNA. In the context of reverse gyrase, these individual activities are suppressed or attenuated. The helicase-like domain of Thermotoga maritima reverse gyrase is a nucleotide-dependent conformational switch that binds DNA and ATP cooperatively. It provides a nucleotide-dependent DNA-binding site to reverse gyrase and thus serves as a valuable model for the investigation of the effect of nucleotides on DNA processing by reverse gyrase that is key to its supercoiling activity. To improve our understanding of the structural basis for the functional cooperation of a helicase domain with a DNA topoisomerase, we have determined the structures of the isolated helicase-like domain of T. maritima reverse gyrase in five different conformations. Comparison of these structures reveals extensive domain flexibility in the absence of conformational restrictions by the topoisomerase that is consistent with single-molecule Fo?rster resonance energy transfer experiments presented here. The structure of the first ADP-bound form provides novel details about nucleotide binding to reverse gyrase. It demonstrates that reverse gyrases use the canonical nucleotide binding mode common to superfamily 2 helicases despite large deviations in the conserved motifs. A characteristic insert region adopts drastically different structures in different reverse gyrases. Counterparts of this insert region are located at very different positions in other DNA-processing enzymes but may point toward a general role in DNA strand separation.