Improvement of diffraction quality upon rehydration of dehydrated icosahedral Enterococcus faecalis pyruvate dehydrogenase core crystals.

Members of the family of 2-oxoacid dehydrogenase multienzyme complexes catalyze the oxidative decarboxylation of alpha-keto acids and are among the most remarkable enzymatic machineries in the living cell. These multienzyme complexes combine a highly symmetric (cubic or icosahedral) core with a dynamic and flexible arrangement of numerous subunits and domains surrounding the core. The center of the complex is formed by either 24 or 60 copies of dihydrolipoamide acetyltransferase (E2)-a multidomain enzyme. The hollow icosahedral cores are composed of 60 identical subunits of the catalytic domain of E2 with a molecular weight of about 1.8 million Da. Bipyramidal crystals suitable for X-ray diffraction of the icosahedral core of the pyruvate dehydrogenase multienzyme complex from Enterococcus faecalis were grown up to 0.7 mm in each dimension. The crystals belong to space group R32 with a = b = 244.3 A (hexagonal setting), and have a solvent content of 73%. The asymmetric unit contains one-third of the molecule, i.e., 20 of the 60 subunits. Initial X-ray crystallographic data to 7 A resolution were collected at cryotemperatures at synchrotron facilities. Interestingly, the diffraction was improved significantly upon rehydrating dehydrated crystals and extended to 4.2 A.     

Exploration of the Dynamic Properties of Protein Complexes Predicted from Spatially Constrained Protein-Protein Interaction Networks

Protein complexes are not static, but rather highly dynamic with subunits that undergo 1-dimensional diffusion with respect to each other. Interactions within protein complexes are modulated through regulatory inputs that alter interactions and introduce new components and deplete existing components through exchange. While it is clear that the structure and function of any given protein complex is coupled to its dynamical properties, it remains a challenge to predict the possible conformations that complexes can adopt. Protein-fragment Complementation Assays detect physical interactions between protein pairs constrained to ≤8 nm from each other in living cells. This method has been used to build networks composed of 1000s of pair-wise interactions. Significantly, these networks contain a wealth of dynamic information, as the assay is fully reversible and the proteins are expressed in their natural context. In this study, we describe a method that extracts this valuable information in the form of predicted conformations, allowing the user to explore the conformational landscape, to search for structures that correlate with an activity state, and estimate the abundance of conformations in the living cell. The generator is based on a Markov Chain Monte Carlo simulation that uses the interaction dataset as input and is constrained by the physical resolution of the assay. We applied this method to an 18-member protein complex composed of the seven core proteins of the budding yeast Arp2/3 complex and 11 associated regulators and effector proteins. We generated 20,480 output structures and identified conformational states using principle component analysis. We interrogated the conformation landscape and found evidence of symmetry breaking, a mixture of likely active and inactive conformational states and dynamic exchange of the core protein Arc15 between core and regulatory components. Our method provides a novel tool for prediction and visualization of the hidden dynamics within protein interaction networks.     

Dynamic and coordinating domain motions in the active subunits of the F1-ATPase molecular motor

F1-ATPase is a rotary molecular motor crucial for various cellular functions. In F1-ATPase, the rotation of the gammadeltaepsilon subunits against the hexameric alpha(3)beta(3) subunits is highly coordinative, driven by ATP hydrolysis and structural changes at three beta subunits. However, the dynamical and coordinating structural transitions in the beta subunits are not fully understood at the molecular level. Here we examine structural transitions and domain motions in the active subunits of F1-ATPase via dynamical domain analysis of the alpha(3)beta(3)gammadeltaepsilon complex. The domain movement and hinge axes and bending residues have been identified and determined for various conformational changes of the beta-subunits. P-loop and the ATP-binding pocket are for the first time found to play essential mechanical functions additional to the catalytic roles. The cooperative conformational changes pertaining to the rotary mechanism of F1-ATPase appears to be more complex than Boyer's 'bi-site' activity. These findings provide unique molecular insights into dynamic and cooperative domain motions in F1-ATPase.     

Structure of the pyruvate dehydrogenase multienzyme complex E1 component from Escherichia coli at 1.85 A resolution

The crystal structure of the recombinant thiamin diphosphate-dependent E1 component from the Escherichia coli pyruvate dehydrogenase multienzyme complex (PDHc) has been determined at a resolution of 1.85 A. The E. coli PDHc E1 component E1p is a homodimeric enzyme and crystallizes with an intact dimer in an asymmetric unit. Each E1p subunit consists of three domains: N-terminal, middle, and C-terminal, with all having alpha/beta folds. The functional dimer contains two catalytic centers located at the interface between subunits. The ThDP cofactors are bound in the "V" conformation in clefts between the two subunits (binding involves the N-terminal and middle domains), and there is a common ThDP binding fold. The cofactors are completely buried, as only the C2 atoms are accessible from solution through the active site clefts. Significant structural differences are observed between individual domains of E1p relative to heterotetrameric multienzyme complex E1 components operating on branched chain substrates. These differences may be responsible for reported alternative E1p binding modes to E2 components within the respective complexes. This paper represents the first structural example of a functional pyruvate dehydrogenase E1p component from any species. It also provides the first representative example for the entire family of homodimeric (alpha2) E1 multienzyme complex components, and should serve as a model for this class of enzymes.     

Computer simulations of cyclic enkephalin analogues

Molecular dynamics simulations and energy minimization studies of cyclic enkephalin analogues incorporating retro-inverso modifications have been carried out. The dynamic trajectories are analyzed in terms of the relative mobility of the 14-membered rings, conformational transitions among equilibrium states, and hydrogen-bonding patterns. The cyclization of the molecules reduces the motion of the ring structures substantially. Time-correlated conformational transitions resulting in the reorientation of peptide units are observed. Hydrogen bonds form principally C7 structures. Because of the incorporation of retro-inverso residues, C6 and C8 structures are also formed. Starting conformations for energy minimizations were obtained from the molecular dynamics simulations and from a systematic search of the conformational space available to the molecules. Several minimum energy backbone and side-chain conformations were found for each analogue. The effect of retro-inverso residues on hydrogen-bonding patterns and backbone conformations is discussed.     

Conformational sampling and nucleotide-dependent transitions of the GroEL subunit probed by unbiased molecular dynamics simulations

GroEL is an ATP dependent molecular chaperone that promotes the folding of a large number of substrate proteins in E. coli. Large-scale conformational transitions occurring during the reaction cycle have been characterized from extensive crystallographic studies. However, the link between the observed conformations and the mechanisms involved in the allosteric response to ATP and the nucleotide-driven reaction cycle are not completely established. Here we describe extensive (in total long) unbiased molecular dynamics (MD) simulations that probe the response of GroEL subunits to ATP binding. We observe nucleotide dependent conformational transitions, and show with multiple 100 ns long simulations that the ligand-induced shift in the conformational populations are intrinsically coded in the structure-dynamics relationship of the protein subunit. Thus, these simulations reveal a stabilization of the equatorial domain upon nucleotide binding and a concomitant "opening" of the subunit, which reaches a conformation close to that observed in the crystal structure of the subunits within the ADP-bound oligomer. Moreover, we identify changes in a set of unique intrasubunit interactions potentially important for the conformational transition.     

Reconciling Mediating and Slaving Roles of Water in Protein Conformational Dynamics

Proteins accomplish their physiological functions with remarkably organized dynamic transitions among a hierarchical network of conformational substates. Despite the essential contribution of water molecules in shaping functionally important protein dynamics, their exact role is still controversial. Water molecules were reported either as mediators that facilitate or as masters that slave protein dynamics. Since dynamic behaviour of a given protein is ultimately determined by the underlying energy landscape, we systematically analysed protein self energies and protein-water interaction energies obtained from extensive molecular dynamics simulation trajectories of barstar. We found that protein-water interaction energy plays the dominant role when compared with protein self energy, and these two energy terms on average have negative correlation that increases with increasingly longer time scales ranging from 10 femtoseconds to 100 nanoseconds. Water molecules effectively roughen potential energy surface of proteins in the majority part of observed conformational space and smooth in the remaining part. These findings support a scenario wherein water on average slave protein conformational dynamics but facilitate a fraction of transitions among different conformational substates, and reconcile the controversy on the facilitating and slaving roles of water molecules in protein conformational dynamics.     

Purification and characterization of a cellulolytic multienzyme complex produced by Neocallimastix patriciarum J11

Understanding the roles of the components of the multienzyme complex of the anaerobial cellulase system, acting on complex substrates, is crucial to the development of efficient cellulase systems for industrial applications such as converting lignocellulose to sugars for bioethanol production. In this study, we purified the multienzyme complex of Neocallimastix patriciarum J11 from a broth through cellulose affinity purification. The multienzyme complex is composed of at least 12 comprised proteins, based on sodium dodecyl sulfate polyacrylamide gel electrophoresis. Eight of these constituents have demonstrated β-glucanase activity on zymogram analysis. The multienzyme complex contained scaffoldings that respond to the gathering of the cellulolytic components. The levels and subunit ratio of the multienzyme complex from N. patriciarum J11 might have been affected by their utilized carbon sources, whereas the components of the complexes were consistent. The trypsin-digested peptides of six proteins were matched to the sequences of cellulases originating from rumen fungi, based on identification through liquid chromatography/mass spectrometry, revealing that at least three types of cellulase, including one endoglucanase and two exoglucanases, could be found in the multienzyme complex of N. patriciarum J11. The cellulolytic subunits could hydrolyze synergistically on both the internal bonds and the reducing and nonreducing ends of cellulose. Based on our research, our findings are the first to depict the composition of the multienzyme complex produced by N. patriciarum J11, and this complex is composed of scaffoldin and three types of cellulase.     

X-ray small-angle studies of the pyruvate dehydrogenase core complex from Escherichia coli K-12

The pyruvate dehydrogenase core complex from E. coli K-12, defined as the multienzyme complex which can be obtained with a unique polypeptide chain composition, has been investigated in solution with the X-ray small-angle technique. The molecular mass of the core complex of 3.78·10⁶ daltons verifies the ratio of polypeptide chains of 161616 of the three enzyme components, pyruvate dehydrogenase, dihydrolipoamide transacetylase, and dihydrolipoamide dehydrogenase, present in the complex. In connection with the values obtained for the radius of gyration (156.5 å), volume (1.07·10⁷ å³) and amount of solvent associated with the complex (1.03 g/g) a loose packing of subunits in the complex has to be assumed. The maximum diameter of the core complex of 433 å, as determined from the correlation function, corroborates the large extension of the complex. The comparison of experimental and theoretical scattering curves reveals a relatively isometric overall shape of the core complex.Enzymes: Pyruvate dehydrogenase complex = pyruvate dehydrogenase (EC 1.2.4.1) plus dihydrolipoamide transacetylase (EC 2.3.1.12) plus dihydrolipoamide dehydrogenase (EC 1.6.4.3).     

Conformational Ensemble Modulates Cooperativity in the Rate-determining Catalytic Step in the E1 Component of the Escherichia coli Pyruvate Dehydrogenase Multienzyme Complex

Cooperativity is extensively used by enzymes, particularly those acting at key metabolic branch points, to “fine tune” catalysis. Thus, cooperativity and enzyme catalysis are intimately linked, yet their linkage is poorly understood. Here we show that negative cooperativity in the rate-determining step in the E1 component of the Escherichia coli pyruvate dehydrogenase multienzyme complex is an outcome of redistribution of a “rate-promoting” conformational pre-equilibrium. An array of biophysical and biochemical studies indicates that non-catalytic but conserved residues directly regulate the redistribution. Furthermore, factors such as ligands and temperature, individually or in concert, also strongly influence the redistribution. As a consequence, these factors also exert their influence on catalysis by profoundly influencing the pre-equilibrium facilitated dynamics of communication between multienzyme components. Our observations suggest a mode of cooperativity in the E1 component that is consistent with the dynamical hypothesis shown to satisfactorily explain cooperativity in many well studied enzymes. The results point to the likely existence of multiple modes of communication between subunits when the entire class of thiamin diphosphate-dependent enzymes is considered.     

Modeling of the two wave response of ATP receptors to jumps of the agonist concentration in pheochromocytoma cells

A model for the kinetics of conformational transitions of ionotropic ATP receptors in pheochromocytoma cells was elaborated. The contribution of the states of ionotropic receptors (upon the blockage of the "open" channel state) to the kinetics of postsynaptic currents was estimated at mediator concentrations studied. The model enables one to determine the contribution of various conformational states of the receptor, in particular in the "closed" state, to the dynamics of ionic current that is registered upon stimulation of ATP receptors. It is shown that after the cessation of the agonist application, a secondary current wave can arise. The rate constants for conformational transitions of ATP receptors were determined.     

Combined use of molecular dynamics simulations and NMR to explore peptide bond isomerization and multiple intramolecular hydrogen-bonding possibilities in a cyclic pentapeptide, cyclo(Gly-Pro-D-Phe-Gly-Val).

The conformational behavior of a model cyclic pentapeptide--cyclo(Gly-L-Pro-D-Phe-Gly-L-Val)--has been explored through the combined use of in vacuo molecular dynamics simulations and a range of nmr experiments (preceding paper). The molecular dynamics analysis suggests that, despite the conformational constraints imposed by formation of the pentapeptide cycle, this pentapeptide undergoes conformational transitions between various hydrogen-bonded conformations, characterized by low energy barriers. An inverse gamma turn with Pro in position i + 1 and a gamma turn with D-Phe in position i + 1 are two alternatives occurring frequently. Like other DLDDL cyclic pentapeptides, cyclo(Gly-Pro-D-Phe-Gly-Val) is also stabilized by an inverse gamma-turn structure with the beta-branched Val residue in position i + 1, and this hydrogen bond is retained in the different conformational families. The gamma-turn around D-Phe3 and the inverse gamma turn around Val5 are consistent with the nmr observations. 3JNH-CH alpha coupling constants of the all-trans forms were calculated from one of the molecular dynamics trajectories and are comparable to nmr experimental data, suggesting that the conformational states visited during the simulation are representative of the conformational distribution in solution. In addition to the equilibrium among various hydrogen-bonded all-trans conformers, the observation in nmr spectra of two sets of resonances for all peptide protons indicated a slow conformational interconversion of the Gly-Pro peptide bond between trans and cis isomers. The activation energy between these two conformers was determined experimentally by magnetization transfer and was calculated by high temperature constrained molecular dynamics simulation. Both methods yield a free energy of activation of ca. 20 kcal/mol. Furthermore, the free energy of activation is dependent on the direction of rotation of the Gly-Pro peptide bond.     

Structure of the bc1 complex from Seculamonas ecuadoriensis,a jakobid flagellate with an ancestral mitochondrial genome

In eubacteria, the respiratory bc(1) complex (complex III) consists of three or four different subunits, whereas that of mitochondria, which have descended from an alpha-proteobacterial endosymbiont, contains about seven additional subunits. To understand better how mitochondrial protein complexes evolved from their simpler bacterial predecessors, we purified complex III of Seculamonas ecuadoriensis, a member of the jakobid protists, which possess the most bacteria-like mitochondrial genomes known. The S. ecuadoriensis complex III has an apparent molecular mass of 460 kDa and exhibits antimycin-sensitive quinol:cytochrome c oxidoreductase activity. It is composed of at least eight subunits between 6 and 46 kDa in size, including two large "core" subunits and the three "respiratory" subunits. The molecular mass of the S. ecuadoriensis bc(1) complex is slightly lower than that reported for other eukaryotes, but about 2x as large as complex III in bacteria. This indicates that the departure from the small bacteria-like complex III took place at an early stage in mitochondrial evolution, prior to the divergence of jakobids. We posit that the recruitment of additional subunits in mitochondrial respiratory complexes is a consequence of the migration of originally alpha-proteobacterial genes to the nucleus.     

Dissociation of ribulose-1,5-bisphosphate carboxylase/oxygenase from spinach by urea

The dissociation of D-ribulose-1,5-bisphosphate carboxylase/oxygenase from spinach, which consists of eight large subunits (L, 53 kDa) and eight small subunits (S, 14 kDa) and thus has a quarternary structure L8S8, has been investigated using a variety of physical techniques. Gel chromatography using Sephadex G-100 indicates the quantitative dissociation of the small subunit S from the complex at 3-4 M urea (50 mM Tris/Cl pH 8.0, 0.5 mM EDTA, 1 mM dithiothreitol and 5 mM 2-mercaptoethanol). The dissociated S is monomeric. Analytical ultracentrifuge studies show that the core of large subunits, L, remaining at 3-4 M urea sediments with S20, w = 15.0 S, whereas the intact enzyme (L8S8) sediments with S20, w = 17.7S. The observed value is consistent with a quarternary structure L8. The dissociation reaction in 3-4 M urea can thus be represented by L8S8----L8 + 8S. At urea concentrations c greater than 5 M the L8 core dissociates into monomeric, unfolded large subunits. A large decrease in fluorescence emission intensity accompanies the dissociation of the small subunit S. This change is completed at 4 M urea. No changes are observed upon dissociating the L8 core. The kinetics of dissociation of the small subunit, as monitored by fluorescence spectroscopy, closely follow the kinetics of loss of carboxylase activity of the enzyme. Studies of the circular dichroism of D-ribulose-1,5-bisphosphate carboxylase in the wavelength region 200-260 nm indicate two conformational transitions. The first one ([0]220 from -8000 to -3500 deg cm2 dmol-1) is completed at 4 M urea and corresponds to the dissociation of the small subunit and coupled conformational changes. The second one ([0]220 from -3500 to -1200 deg cm2 dmol-1) is completed at 6 M urea and reflects the dissociation and unfolding of large subunits from the core. The effect of activation of the enzyme by addition of MgCl2 (10 mM) and NaHCO3 (10 mM) on these conformational transitions was investigated. The first conformational transition is then shifted to higher urea concentrations: a single transition ([0]220 from -8000 to -1200 deg cm2 dmol-1) is observed for the activated enzyme. From the urea dissociation experiments we conclude that both large (L) and small (S) subunits are important for carboxylase activity of spinach D-ribulose-1,5-bisphosphate carboxylase: the L-S subunit interactions tighten upon activation and dissociation of S leads to a coupled, proportional loss of enzyme activity.     

Conformational relaxation following hydride transfer plays a limiting role in dihydrofolate reductase catalysis

The catalytic cycle of an enzyme is frequently associated with conformational changes that may limit maximum catalytic throughput. In Escherichia coli dihydrofolate reductase, release of the tetrahydrofolate (THF) product is the rate-determining step under physiological conditions and is associated with an "occluded" to "closed" conformational change. In this study, we demonstrate that in dihydrofolate reductase the closed to occluded conformational change in the product ternary complex (E.THF.NADP (+)) also gates progression through the catalytic cycle. Using NMR relaxation dispersion, we have measured the temperature and pH dependence of microsecond to millisecond time scale backbone dynamics of the occluded E.THF.NADP (+) complex. Our studies indicate the presence of three independent dynamic regions, associated with the active-site loops, the cofactor binding cleft, and the C-terminus and an adjacent loop, which fluctuate into discrete conformational substates with different kinetic and thermodynamic parameters. The dynamics of the C-terminally associated region is pH-dependent (p K a < 6), but the dynamics of the active-site loops and cofactor binding cleft are pH-independent. The active-site loop dynamics access a closed conformation, and the accompanying closed to occluded rate constant is comparable to the maximum pH-independent hydride transfer rate constant. Together, these results strongly suggest that the closed to occluded conformational transition in the product ternary complex is a prerequisite for progression through the catalytic cycle and that the rate of this process places an effective limit on the maximum rate of the hydride transfer step.     

Entropy and barrier-controlled fluctuations determine conformational viscoelasticity of single biomolecules

Biological macromolecules have complex and nontrivial energy landscapes, endowing them with a unique conformational adaptability and diversity in function. Hence, understanding the processes of elasticity and dissipation at the nanoscale is important to molecular biology and emerging fields such as nanotechnology. Here we analyze single molecule fluctuations in an atomic force microscope, using a generic model of biopolymer viscoelasticity that includes local "internal" conformational dissipation. Comparing two biopolymers, dextran and cellulose (polysaccharides with and without local bistable transitions), demonstrates that signatures of simple conformational change are minima in both the elastic and internal friction constants around a characteristic force. A novel analysis of dynamics on a bistable energy landscape provides a simple explanation: an elasticity driven by the entropy, and friction by a barrier-controlled hopping time of populations between states, which is surprisingly distinct to the well-known relaxation time. This nonequilibrium microscopic analysis thus provides a means of quantifying new dynamical features of the energy landscape of the glucopyranose ring, revealing an unexpected underlying roughness and information on the shape of the barrier of the chair-boat transition in dextran. The results presented herein provide a basis toward probing the viscoelasticity of macromolecular conformational transitions on more complex energy landscapes, such as during protein folding.     

Three-dimensional structure of the truncated core of the Saccharomyces cerevisiae pyruvate dehydrogenase complex determined from negative stain and cryoelectron microscopy images.

Dihydrolipoamide acyltransferase (E2), a catalytic and structural component of the three functional classes of multienzyme complexes that catalyze the oxidative decarboxylation of alpha-keto acids, forms the central core to which the other components are attached. We have imaged by negative stain and cryoelectron microscopy the truncated dihydrolipoamide acetyltransferase core (60 subunits; M(r) = 2.7 x 10(6)) of the Saccharomyces cerevisiae pyruvate dehydrogenase complex. Using icosahedral particle reconstruction techniques, we determined its structure to 25 A resolution. Although the model derived from the negative stain reconstruction was approximately 20% smaller than the model derived from the frozen-hydrated data, when corrected for the effects of the electron microscope contrast transfer functions, the reconstructions showed excellent correspondence. The pentagonal dodecahedron-shaped macromolecule has a maximum diameter, as measured along the 3-fold axis, of approximately 226 A (frozen-hydrated value), and 12 large openings (approximately 63 A in diameter) on the 5-fold axes that lead into a large solvent-accessible cavity (approximately 76-140 A diameter). The 20 vertices consist of cone-shaped trimers, each with a flattened base on the outside of the structure and an apex directed toward the center. The trimers are interconnected by 20 A thick "bridges" on the 2-fold axes. These studies also show that the highest resolution features apparent in the frozen-hydrated reconstruction are revealed in a filtered reconstruction of the stained molecule.     

Microtubule dynamics and microtubule caps: a time-resolved cryo- electron microscopy study

Microtubules display the unique property of dynamic instability characterized by phase changes between growth and shrinkage, even in constant environmental conditions. The phases can be synchronized, leading to bulk oscillations of microtubules. To study the structural basis of dynamic instability we have examined growing, shrinking, and oscillating microtubules by time-resolved cryo-EM. In particular we have addressed three questions which are currently a matter of debate: (a) What is the relationship between microtubules, tubulin subunits, and tubulin oligomers in microtubule dynamics?; (b) How do microtubules shrink? By release of subunits or via oligomers?; and (c) Is there a conformational change at microtubule ends during the transitions from growth to shrinkage and vice versa? The results show that (a) oscillating microtubules coexist with a substantial fraction of oligomers, even at a maximum of microtubule assembly; (b) microtubules disassemble primarily into oligomers; and (c) the ends of growing microtubules have straight protofilaments, shrinking microtubules have protofilaments coiled inside out. This is interpreted as a transition from a tense to a relaxed conformation which could be used to perform work, as suggested by some models of poleward chromosome movement during anaphase.