Biochemical and mutational analysis of EcoRII functional domains reveals evolutionary links between restriction enzymes

The archetypal Type IIE restriction endonuclease EcoRII is a dimer that has a modular structure. DNA binding studies indicate that the isolated C-terminal domain dimer has an interface that binds a single cognate DNA molecule whereas the N-terminal domain is a monomer that also binds a single copy of cognate DNA. Hence, the full-length EcoRII contains three putative DNA binding interfaces: one at the C-terminal domain dimer and two at each of the N-terminal domains. Mutational analysis indicates that the C-terminal domain shares conserved active site architecture and DNA binding elements with the tetrameric restriction enzyme NgoMIV. Data provided here suggest possible evolutionary relationships between different subfamilies of restriction enzymes.     

C.EcoO109I,a regulatory protein for production of EcoO109I restriction endonuclease,specifically binds to and bends DNA upstream of its translational start site

The EcoO109I restriction-modification system, which recognizes 5′-(A/G)GGNCC(C/T)-3′, has been cloned, and contains convergently transcribed endonuclease and methylase. The role and action mechanism of the gene product, C.EcoO109I, of a small open reading frame located upstream of ecoO109IR were investigated in vivo and in vitro. The results of deletion analysis suggested that C.EcoO109I acts as a positive regulator of ecoO109IR expression but has little effect on ecoO109IM expression. Assaying of promoter activity showed that the expression of ecoO109IC was regulated by its own gene product, C.EcoO109I. C.EcoO109I was overproduced as a His-tag fusion protein in recombinant Escherichia coli HB101 and purified to homogeneity. C.EcoO109I exists as a homodimer, and recognizes and binds to the DNA sequence 5′-CTAAG(N)₅CTTAG-3′ upstream of the ecoO109IC translational start site. It was also shown that C.EcoO109I bent the target DNA by 54 ± 4°.     

Dynamic and Structural Changes in the Minimally Restructuring EcoRI Bound to a Minimally Mutated DNA Chain

Abstract

The dynamics of a protein plays an important role in protein functionality. Here, we examine the differences in the dynamics of a minimally restructuring protein, EcoRI, when it is bound to its cognate DNA and to a noncognate sequence which differs by just a single basepair. Molecular dynamics simulations of the complexes and essential dynamics analyses reveal that the overall dynamics of the protein subunits change from a coordinated motion in the cognate complex to a scrambled motion in the noncognate complex. This dynamical difference extends to the protein-DNA interface where EcoRI tries to constrict the DNA in the cognate complex. In the noncognate complex, absence of the constricting motion of interfacial residues, overall change in backbone dynamics and structural relaxation of the arms enfolding the DNA leave the DNA less-kinked relative to the situation in the cognate complex, thus indicating that the protein is poised for linear diffusion along the DNA rather than for catalytic action. In a larger context, the results imply that the DNA sequences dictate protein dynamics and that when a protein chances upon the recognition sequence some of the key domains of the protein undergo dynamical changes that prepare the protein for eventual catalytic action.     

Dynamic and Structural Changes in the Minimally Restructuring EcoRI Bound to a Minimally Mutated DNA Chain

The dynamics of a protein plays an important role in protein functionality. Here, we examine the differences in the dynamics of a minimally restructuring protein, EcoRI, when it is bound to its cognate DNA and to a noncognate sequence which differs by just a single basepair. Molecular dynamics simulations of the complexes and essential dynamics analyses reveal that the overall dynamics of the protein subunits change from a coordinated motion in the cognate complex to a scrambled motion in the noncognate complex. This dynamical difference extends to the protein-DNA interface where EcoRI tries to constrict the DNA in the cognate complex. In the noncognate complex, absence of the constricting motion of interfacial residues, overall change in backbone dynamics and structural relaxation of the arms enfolding the DNA leave the DNA less-kinked relative to the situation in the cognate complex, thus indicating that the protein is poised for linear diffusion along the DNA rather than for catalytic action. In a larger context, the results imply that the DNA sequences dictate protein dynamics and that when a protein chances upon the recognition sequence some of the key domains of the protein undergo dynamical changes that prepare the protein for eventual catalytic action.     

Conformational Dynamics of Thermus aquaticus DNA Polymerase I during Catalysis

Despite the fact that DNA polymerases have been investigated for many years and are commonly used as tools in a number of molecular biology assays, many details of the kinetic mechanism they use to catalyze DNA synthesis remain unclear. Structural and kinetic studies have characterized a rapid, pre-catalytic open-to-close conformational change of the Finger domain during nucleotide binding for many DNA polymerases including Thermus aquaticus DNA polymerase I (Taq Pol), a thermostable enzyme commonly used for DNA amplification in PCR. However, little has been performed to characterize the motions of other structural domains of Taq Pol or any other DNA polymerase during catalysis. Here, we used stopped-flow Förster resonance energy transfer to investigate the conformational dynamics of all five structural domains of the full-length Taq Pol relative to the DNA substrate during nucleotide binding and incorporation. Our study provides evidence for a rapid conformational change step induced by dNTP binding and a subsequent global conformational transition involving all domains of Taq Pol during catalysis. Additionally, our study shows that the rate of the global transition was greatly increased with the truncated form of Taq Pol lacking the N-terminal domain. Finally, we utilized a mutant of Taq Pol containing a de novo disulfide bond to demonstrate that limiting protein conformational flexibility greatly reduced the polymerization activity of Taq Pol.     

Coupling between Catalytic Loop Motions and Enzyme Global Dynamics

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.     

Structural basis of cellulosome efficiency explored by small angle X-ray scattering

Cellulose, the main structural component of plant cell walls, is the most abundant carbohydrate polymer in nature. To break down plant cell walls, anaerobic microorganisms have evolved a large extracellular enzyme complex termed cellulosome. This megadalton catalytic machinery organizes an enzymatic assembly, tenaciously bound to a scaffolding protein via specialized intermodular "cohesin-dockerin" interactions that serve to enhance synergistic activity among the different catalytic subunits. Here, we report the solution structure properties of cellulosome-like assemblies analyzed by small angle x-ray scattering and molecular dynamics. The atomic models, generated by our strategy for the free chimeric scaffoldin and for binary and ternary complexes, reveal the existence of various conformations due to intrinsic structural flexibility with no, or only coincidental, inter-cohesin interactions. These results provide primary evidence concerning the mechanisms by which these protein assemblies attain their remarkable synergy. The data suggest that the motional freedom of the scaffoldin allows precise positioning of the complexed enzymes according to the topography of the substrate, whereas short-scale motions permitted by residual flexibility of the enzyme linkers allow "fine-tuning" of individual catalytic domains.     

Structure and Function in Homodimeric Enzymes: Simulations of Cooperative and Independent Functional Motions

Large-scale conformational change is a common feature in the catalytic cycles of enzymes. Many enzymes function as homodimers with active sites that contain elements from both chains. Symmetric and anti-symmetric cooperative motions in homodimers can potentially lead to correlated active site opening and/or closure, likely to be important for ligand binding and release. Here, we examine such motions in two different domain-swapped homodimeric enzymes: the DcpS scavenger decapping enzyme and citrate synthase. We use and compare two types of all-atom simulations: conventional molecular dynamics simulations to identify physically meaningful conformational ensembles, and rapid geometric simulations of flexible motion, biased along normal mode directions, to identify relevant motions encoded in the protein structure. The results indicate that the opening/closure motions are intrinsic features of both unliganded enzymes. In DcpS, conformational change is dominated by an anti-symmetric cooperative motion, causing one active site to close as the other opens; however a symmetric motion is also significant. In CS, we identify that both symmetric (suggested by crystallography) and asymmetric motions are features of the protein structure, and as a result the behaviour in solution is largely non-cooperative. The agreement between two modelling approaches using very different levels of theory indicates that the behaviours are indeed intrinsic to the protein structures. Geometric simulations correctly identify and explore large amplitudes of motion, while molecular dynamics simulations indicate the ranges of motion that are energetically feasible. Together, the simulation approaches are able to reveal unexpected functionally relevant motions, and highlight differences between enzymes.     

Molecular motions and conformational changes of HPPK

6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) belongs to a class of catalytic enzymes involved in phosphoryl transfer and is a new target for the development of novel antimicrobial agents. In the present study, the fundamental consideration is to view the overall structure of HPPK as a network of interacting residues and to extract the most cooperative collective motions that define its global dynamics. A coarse-grained model, harmonically constrained according to HPPK's crystal structure is used. Four crystal structures of HPPK (one apo and three holo forms with different nucleotide and pterin analogs) are studied with the goal of providing insights about the function-dynamic correlation and ligand induced conformational changes. The dynamic differences are examined between HPPK's apo- and holo-forms, because they are involved in the catalytic reaction steps. Our results indicate that the palm-like structure of HPPK is nearly rigid, whereas the two flexible loops: L2 (residues 43-53) and L3 (residues 82-92) exhibit the most concerted motions for ligand recognition and presumably, catalysis. These two flexible loops are involved in the recognition of HPPKs nucleotide and pterin ligands, whereas the rigid palm region is associated with binding of these cognate ligands. Six domains of collective motions are identified, comprised of structurally close but not necessarily sequential residues. Two of these domains correspond to the flexible loops (L2 and L3), whereas the remaining domains correspond to the rigid part of the molecule.     

Conformational Dynamics of Damage Processing by Human DNA Glycosylase NEIL1

Endonuclease VIII-like protein 1 (NEIL1) is a DNA repair enzyme found in higher eukaryotes, including humans. It belongs to the helix–two turn–helix (H2TH) structural superfamily together with Escherichia coli formamidopyrimidine–DNA glycosylase (Fpg) and endonuclease VIII (Nei), and removes a variety of oxidized purine and pyrimidine bases from DNA. Structural, modeling and kinetic studies have established that the bacterial H2TH superfamily enzymes proceed through several conformational intermediates while recognizing and removing their cognate lesions. Here we apply stopped-flow kinetics with detection of intrinsic Trp fluorescence and Förster resonance energy transfer fluorescence to follow the conformational dynamics of human NEIL1 and DNA when the enzyme interacts with undamaged DNA, or DNA containing cleavable or non-cleavable abasic sites, or dihydrouracil lesions. NEIL1 processed a natural abasic site and a damaged base in DNA equally well but showed an additional fluorescently discernible step when DHU was present, likely reflecting additional rearrangements during base eversion into the enzyme's active site. With undamaged DNA and DNA containing a non-cleavable abasic site analog, (3-hydroxytetrahydrofuran-2-yl)methyl phosphate, NEIL1 was diverted to a non-productive DNA conformation early in the reaction. Our results support the view of NEIL1 as an enzyme that actively destabilizes damaged DNA and uses multiple checkpoints along the reaction coordinate to drive substrate lesions into the active site while rejecting normal bases and non-substrate lesions.     

A Physical Picture of Protein Dynamics and Conformational Changes

A physical model is reviewed which explains different aspects of protein dynamics consistently. At low temperatures, the molecules are frozen in conformational substates. Their average energy is 3/2RT. Solid-state vibrations occur on a time scale of femtoseconds to nanoseconds. Above a characteristic temperature, often called the dynamical transition temperature, slow modes of motions can be observed occurring on a time scale between about 140 and 1 ns. These motions are overdamped, quasidiffusive, and involve collective motions of segments of the size of an α-helix. Molecules performing these types of motion are in the “flexible state”. This state is reached by thermal activation. It is shown that these motions are essential for conformational relaxation. Based on this picture, a new approach is proposed to understand conformational changes. It connects structural fluctuations and conformational transitions.     

Combination of Neutron Scattering and Molecular Dynamics to Determine Internal Motions in Biomolecules

Abstract

The forms and frequencies of atomic dynamics on the pico- and nanosecond timescales are accessible experimentally using incoherent neutron scattering. Molecular dynamics simulations cover the same space and time domains and neutron scattering intensities can be calculated from the simulations for direct comparison with experiment. To illustrate the complementarity of neutron scattering and molecular dynamics we examine measured and simulation-derived elastic incoherent scattering profiles from myoglobin and from the crystalline alanine dipeptide. Elastic incoherent scattering gives information on the geometry of the volume accessible to the atoms in the samples. The simulation-derived dipeptide elastic scattering profiles are in reasonable accord with experiment, deviations being due to the sampling limitations in the simulations and experimental detector normalisation procedures. The simulated dynamics is decomposed, revealing characteristic profiles due to rotational diffusional and translational vibrational motions of the methyl groups. In myoglobin, for which the timescale of the simulation matches more closely that accessible to the experiment, good agreement is seen for the elastic incoherent structure factor. This indicates that the space sampled by the hydrogen atoms in the protein on the timescale <100 ps is well represented by the simulation. Part of the helix atom fluctuations can be described in terms of rigid helix motions.     

Simultaneous Measurements of Solvent Dynamics and Functional Kinetics in a Light-Activated Enzyme

Solvent fluctuations play a key role in controlling protein motions and biological function. Here, we have studied how individual steps of the reaction catalyzed by the light-activated enzyme protochlorophyllide oxidoreductase (POR) couple with solvent dynamics. To simultaneously monitor the catalytic cycle of the enzyme and the dynamical behavior of the solvent, we designed temperature-dependent UV-visible microspectrophotometry experiments, using flash-cooled nanodroplets of POR to which an exogenous soluble fluorophore was added. The formation and decay of the first two intermediates in the POR-catalyzed reaction were measured, together with the solvent glass transition and the buildup of crystalline ice at cryogenic temperatures. We find that formation of the first intermediate occurs below the glass transition temperature (T_g), and is not affected by changes in solvent dynamics induced by modifying the glycerol content. In contrast, formation of the second intermediate occurs above T_g and is influenced by changes in glycerol concentration in a manner remarkably similar to the buildup of crystalline ice. These results suggest that internal, nonslaved protein motions drive the first step of the POR-catalyzed reaction whereas solvent-slaved motions control the second step. We propose that the concept of solvent slaving applies to complex enzymes such as POR.     

The elastic network model reveals a consistent picture on intrinsic functional dynamics of type II restriction endonucleases

The vibrational dynamics of various type II restriction endonucleases, in complex with cognate/non-cognate DNA and in the apo form, are investigated with the elastic network model in order to reveal common functional mechanisms in this enzyme family. Scissor-like and tong-like motions observed in the slowest modes of all enzymes and their complexes point to common DNA recognition and cleavage mechanisms. Normal mode analysis further points out that the scissor-like motion has an important role in differentiating between cognate and non-cognate sequences at the recognition site, thus implying its catalytic relevance. Flexible regions observed around the DNA-binding site of the enzyme usually concentrate on the highly conserved β-strands, especially after DNA binding. These β-strands may have a structurally stabilizing role in functional dynamics for target site recognition and cleavage. In addition, hot spot residues based on high-frequency modes reveal possible communication pathways between the two distant cleavage sites in the enzyme family. Some of these hot spots also exist on the shortest path between the catalytic sites and are highly conserved.     

Conformational Dynamics and Ligand Binding in the Multi-Domain Protein PDC109

PDC109 is a modular multi-domain protein with two fibronectin type II (Fn2) repeats joined by a linker. It plays a major role in bull sperm binding to the oviductal epithelium through its interactions with phosphorylcholines (PhCs), a head group of sperm cell membrane lipids. The crystal structure of the PDC109-PhC complex shows that each PhC binds to the corresponding Fn2 domain, while the two domains are on the same face of the protein. Long timescale explicit solvent molecular dynamics (MD) simulations of PDC109, in the presence and absence of PhC, suggest that PhC binding strongly correlates with the relative orientation of choline-phospholipid binding sites of the two Fn2 domains; unless the two domains tightly bind PhCs, they tend to change their relative orientation by deforming the flexible linker. The effective PDC109-PhC association constant of 28 M, estimated from their potential of mean force is consistent with the experimental result. Principal component analysis of the long timescale MD simulations was compared to the significantly less expensive normal mode analysis of minimized structures. The comparison indicates that difference between relative domain motions of PDC109 with bound and unbound PhC is captured by the first principal component in the principal component analysis as well as the three lowest normal modes in the normal mode analysis. The present study illustrates the use of detailed MD simulations to clarify the energetics of specific ligand-domain interactions revealed by a static crystallographic model, as well as their influence on relative domain motions in a multi-domain protein.     

Molecular Dynamics of HIV1-Integrase in Complex with 93del-A Structural Perspective on the Mechanism of Inhibition

HIV1 integrase is an important target for the antiviral therapy. Guanine-rich quadruplex, such as 93del, have been shown to be potent inhibitors of this enzyme and thus representing a new class of antiviral agents. Although X-ray and NMR structures of HIV1 integrase and 93del have been reported, there is no structural information of the complex and the mechanism of inhibition still remains unexplored. A number of computational methods including automated protein-DNA docking and molecular dynamics simulation in explicit solvent were used to model the binding of 93del to HIV1 integrase. Analysis of the dynamic behaviour of the complex using principal components analysis and elastic network modelling techniques allow us to understand how the binding of 93del aptamer and its interactions with key residues affect the intrinsic motions of the catalytic loops by stabilising them in catalytically inactive conformations. Such insights into the structural mechanism of inhibition can aid in improving the design of anti-HIV aptamers.     

Cooperative dynamics of a DNA polymerase replicating complex

Engineered DNA polymerases continue to be the workhorses of many applications in biotechnology, medicine and nanotechnology. However, the dynamic interplay between the enzyme and the DNA remains unclear. In this study, we performed an extensive replica exchange with flexible tempering (REFT) molecular dynamics simulation of the ternary replicating complex of the archaeal family B DNA polymerase from the thermophile Thermococcus gorgonarius, right before the chemical step. The convoluted dynamics of the enzyme are reducible to rigid-body motions of six subdomains. Upon binding to the enzyme, the DNA double helix conformation changes from a twisted state to a partially untwisted state. The twisted state displays strong bending motion, whereby the DNA oscillates between a straight and a bent conformation. The dynamics of double-stranded DNA are strongly correlated with rotations of the thumb toward the palm, which suggests an assisting role of the enzyme during DNA translocation. In the complex, the primer–template duplex displays increased preference for the B-DNA conformation at the n ? 2 and n ? 3 dinucleotide steps. Interactions at the primer 3′ end indicate that Thr541 and Asp540 are the acceptors of the first proton transfer in the chemical step, whereas in the translocation step both residues hold the primer 3′ terminus in the vicinity of the priming site, which is crucial for high processivity.     

Promotion of Enzyme Flexibility by Dephosphorylation and Coupling to the Catalytic Mechanism of a Phosphohexomutase

The enzyme phosphomannomutase/phosphoglucomutase (PMM/PGM) from Pseudomonas aeruginosa catalyzes an intramolecular phosphoryl transfer across its phosphosugar substrates, which are precursors in the synthesis of exoproducts involved in bacterial virulence. Previous structural studies of PMM/PGM have established a key role for conformational change in its multistep reaction, which requires a dramatic 180° reorientation of the intermediate within the active site. Here hydrogen-deuterium exchange by mass spectrometry and small angle x-ray scattering were used to probe the conformational flexibility of different forms of PMM/PGM in solution, including its active, phosphorylated state and the unphosphorylated state that occurs transiently during the catalytic cycle. In addition, the effects of ligand binding were assessed through use of a substrate analog. We found that both phosphorylation and binding of ligand produce significant effects on deuterium incorporation. Phosphorylation of the conserved catalytic serine has broad effects on residues in multiple domains and is supported by small angle x-ray scattering data showing that the unphosphorylated enzyme is less compact in solution. The effects of ligand binding are generally manifested near the active site cleft and at a domain interface that is a site of conformational change. These results suggest that dephosphorylation of the enzyme may play two critical functional roles: a direct role in the chemical step of phosphoryl transfer and secondly through propagation of structural flexibility. We propose a model whereby increased enzyme flexibility facilitates the reorientation of the reaction intermediate, coupling changes in structural dynamics with the unique catalytic mechanism of this enzyme.