Large-scale Domain Motions Research Articles

Conformational Selection Governs Carrier Domain Positioning in Staphylococcus aureus Pyruvate Carboxylase.

Biotin-dependent enzymes employ a carrier domain to efficiently transport reaction intermediates between distant active sites. The translocation of this carrier domain is critical to the interpretation of kinetic and structural studies, but there have been few direct attempts to investigate the dynamic interplay between ligand binding and carrier domain positioning in biotin-dependent enzymes. Pyruvate carboxylase (PC) catalyzes the MgATP-dependent carboxylation of pyruvate where the biotinylated carrier domain must translocate ∼70 Å from the biotin carboxylase domain to the carboxyltransferase domain. Many prior studies have assumed that carrier domain movement is governed by ligand-induced conformational changes, but the mechanism underlying this movement has not been confirmed. Here, we have developed a system to directly observe PC carrier domain positioning in both the presence and absence of ligands, independent of catalytic turnover. We have incorporated a cross-linking trap that reports on the interdomain conformation of the carrier domain when it is positioned in proximity to a neighboring carboxyltransferase domain. Cross-linking was monitored by gel electrophoresis, inactivation kinetics, and intrinsic tryptophan fluorescence. We demonstrate that the carrier domain positioning equilibrium is sensitive to substrate analogues and the allosteric activator acetyl-CoA. Notably, saturating concentrations of biotin carboxylase ligands do not prevent carrier domain trapping proximal to the neighboring carboxyltransferase domain, demonstrating that carrier domain positioning is governed by conformational selection. This model of carrier domain translocation in PC can be applied to other multi-domain enzymes that employ large-scale domain motions to transfer intermediates during catalysis.

Revealing the activation mechanism of autoinhibited RalF by integrated simulation and experimental approaches

RalF is an Arf GEF from Legionella pneumophilia, the bacterium that causes severe pneumonia. In its crystal structure, RalF is in the autoinhibited form. A large-scale domain motion is expected to lift the autoinhibition, the mechanism of which is still unknown. Since RalF is activated in the presence of the membrane, its active structure and the structure of the RalF-Arf1 complex could not have been determined experimentally. On the simulation side, it has been proven that classical Molecular Dynamics (MD) alone is not efficient enough to map motions of such amplitude and determine the active conformation of RalF. In this article, using Molecular Dynamics with excited Normal Modes (MDeNM) combined with previous experimental findings we were able to determine the active RalF structure and the structure of the RalF-Arf1 complex in the presence of the membrane, bridging the gap between experiments and simulation.

Conformational Dynamics in NF-κB Transcriptional Regulation

The NF-κB family of transcription factor is a central mediator of immune and imflammatory responses. NF-κB functions as a dimer and recognize DNA through the two N-terminal domains (NTDs). Previous hydrogen-deuterium exchange mass spectrometry, stopped-flow experiments and molecular dynamics simulations suggested the two NTDs can undergo relative motions which plays an important role in the kinetics of DNA association and dissociation. Here we showed by single-molecule FRET the direct evidence of relative domain motions in the NF-κB RelA-p50 dimer. Surprisingly, the two NTDs underwent slow heterogeneous motions on the timescale of seconds to minutes. These motions are scaled down when NF-κB is bound to DNA and allosterically altered by the inhibitor protein IκBα, which is known to accelerate DNA dissociation and prevent DNA binding. A new DNA-bound conformation unexpected from crystal structure was also observed. Together, our results provide a direct visualization and quantitative characterization of large-scale domain motions of NF-κB and shed light on the understanding of the role of protein dynamics in protein-protein interactions.

Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis.

Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.

Protein Conformational Change Is Essential for Reductive Activation of Lytic Polysaccharide Monooxygenase by Cellobiose Dehydrogenase.

Large-scale protein domain dynamics and electron transfer are often associated. However, as protein motions span a broad range of time and length scales, it is often challenging to identify and thus link functionally relevant dynamic changes to electron transfer in proteins. It is hypothesized that large-scale domain motions direct electrons through a FAD and a heme b cofactor of the fungal cellobiose dehydrogenase (CDH) enzymes to the type-II copper center (T2Cu) of the polysaccharide-degrading lytic polysaccharide monooxygenases (LPMOs). However, as of yet, domain motions in CDH have not been linked formally to enzyme-catalyzed electron transfer reactions. The detailed structural features of CDH, which govern the functional conformational landscapes of the enzyme, have only been partially resolved. Here, we use a combination of pressure, viscosity, ionic strength, and temperature perturbation stopped-flow studies to probe the conformational landscape associated with the electron transfer reactions of CDH. Through the use of molecular dynamics simulations, potentiometry, and stopped-flow spectroscopy, we investigated how a conserved Tyr99 residue plays a key role in shaping the conformational landscapes for both the interdomain electron transfer reactions of CDH (from FAD to heme) and the delivery of electrons from the reduced heme cofactor to the LPMO T2Cu. Our studies show how motions gate the electron transfer within CDH and from CDH to LPMO and illustrate the conformational landscape for interdomain and interprotein electron transfer in this extracellular fungal electron transfer chain.

A Multi-State Coarse-Grained Simulation Model Captures Conformational Cycling in P-Type ATPases

P-type ATPases are a large family of ion pumps that precisely move ions into and out of cell membrane against their concentration gradients by consuming ATP. All the family members share a common basic pump cycle (also known as the Post-Albers cycle) consisting of two main conformational states (termed E1 and E2) with distinct conformations and affinities for the primary transported ions. Over 120 high-resolution structures of P-type ATPases have been accumulated in the Protein Data Bank (PDB), with most of the available structures being of the sarcoplasmic reticulum Ca2+-ATPase (SERCA, ∼60%) and the Na+/K+-ATPase (NKA, ∼30%). While these experimentally-derived structures capture the ATPases in a variety of conformational states, and thus provide crucial insight into transport mechanism, they do not directly reveal how the conformational changes occur, and how signals are transferred between the transmembrane regions and the cytoplasmic domains consisting of the phosphorylation (P) domain, the nucleotide binding (N) domain and the actuator (A) domain. To bridge the missing dynamics information between static structures during the Post-Albers cycle, we have developed a coarse-grained model by combining a multi-basin structure-based potential with a modified implicit solvent MARTINI model and enhance sampling using metadynamics simulations. The model allows access the millisecond dynamics and large-scale domain motions of such large membrane proteins with more than a thousand of residues. In this way we are able to capture conformational changes between the E1-E2 states of SERCA and NKA and obtain more detailed information about the mechanisms underlying transport function and regulation.

Multiscale enhanced sampling of glucokinase: Regulation of the enzymatic reaction via a large scale domain motion

Enhanced sampling yields a comprehensive structural ensemble or a free energy landscape, which is beyond the capability of a conventional molecular dynamics simulation. Our recently developed multiscale enhanced sampling (MSES) method employs a coarse-grained model coupled with the target physical system for the efficient acceleration of the dynamics. MSES has demonstrated applicability to large protein systems in solution, such as intrinsically disordered proteins and protein-protein and protein-ligand interactions. Here, we applied the MSES simulation to an important drug discovery target, glucokinase (GCK), to elucidate the structural basis of the positive cooperativity of the enzymatic reaction at an atomistic resolution. MSES enabled us to compare two sets of the free energy landscapes of GCK, for the glucose-bound and glucose-unbound forms, and thus demonstrated the drastic change of the free energy surface depending on the glucose concentration. In the glucose-bound form, we found two distinct basins separated by a high energy barrier originating from the domain motion and the folding/unfolding of the α13 helix. By contrast, in the glucose-unbound form, a single flat basin extended to the open and super-open states. These features illustrated the two distinct phases achieving the cooperativity, the fast reaction cycle staying in the closed state at a high glucose concentration and the slow cycle primarily in the open/super-open state at a low concentration. The weighted ensemble simulations revealed the kinetics of the structural changes in GCK with the synergetic use of the MSES results; the rate constant of the transition between the closed state and the open/super-open states, kC/O = 1.1 ms-1, is on the same order as the experimental catalytic rate, kcat = 0.22 ms-1. Finally, we discuss the pharmacological activities of GCK activators (small molecular drugs modulating the GCK activity) in terms of the slight changes in the domain motion, depending on their chemical structures as regulators. The present study demonstrated the capability of the enhanced sampling and the associated kinetic calculations for understanding the atomistic structural dynamics of protein systems in physiological environments.

Essential role of the flexible linker on the conformational equilibrium of bacterial peroxiredoxin reductase for effective regeneration of peroxiredoxin

Reactive oxygen species (ROS) can damage DNA, proteins, and lipids, so cells have antioxidant systems that regulate ROS. In many bacteria, a dedicated peroxiredoxin reductase, alkyl hydroperoxide reductase subunit F (AhpF), catalyzes the rapid reduction of the redox-active disulfide center of the antioxidant protein peroxiredoxin (AhpC) to detoxify ROS such as hydrogen peroxide, organic hydroperoxide, and peroxynitrite. AhpF is a flexible multidomain protein that enables a series of electron transfers among the redox centers by accepting reducing equivalents from NADH. A flexible linker connecting the N-terminal domain (NTD) and C-terminal domain (CTD) of AhpF suggests that the enzyme adopts a large-scale domain motion that alternates between the closed and open states to shuttle electrons from the CTD via the NTD to AhpC. Here, we conducted comprehensive mutational, biochemical, and biophysical analyses to gain insights into the role of the flexible linker and the residues critical for the domain motions of Escherichia coli AhpF (EcAhpF) during electron transfer. Small-angle X-ray scattering studies of linker mutants revealed that a group of charged residues, 200EKR202, is crucial for the swiveling motion of the NTD. Moreover, NADH binding significantly affected EcAhpF flexibility and the movement of the NTD relative to the CTD. The mutants also exhibited a decrease in H2O2 reduction by the AhpF-AhpC ensemble. We propose that a concerted movement involving the NTD, C-terminal NADH, and FAD domains, and the flexible linker between them is essential for optimal intra-domain cross-talk and for efficient electron transfer to the redox partner AhpC required for peroxidation.

High Affinity vs. Native Fibronectin in the Modulation of αvβ3 Integrin Conformational Dynamics: Insights from Computational Analyses and Implications for Molecular Design.

Understanding how binding events modulate functional motions of multidomain proteins is a major issue in chemical biology. We address several aspects of this problem by analyzing the differential dynamics of αvβ3 integrin bound to wild type (wtFN10, agonist) or high affinity (hFN10, antagonist) mutants of fibronectin. We compare the dynamics of complexes from large-scale domain motions to inter-residue coordinated fluctuations to characterize the distinctive traits of conformational evolution and shed light on the determinants of differential αvβ3 activation induced by different FN sequences. We propose an allosteric model for ligand-based integrin modulation: the conserved integrin binding pocket anchors the ligand, while different residues on the two FN10’s act as the drivers that reorganize relevant interaction networks, guiding the shift towards inactive (hFN10-bound) or active states (wtFN10-bound). We discuss the implications of results for the design of integrin inhibitors.

Bridging Enzymatic Structure Function via Mechanics: A Coarse-Grain Approach.

Flexibility is a central aspect of protein function, and ligand binding in enzymes involves a wide range of structural changes, ranging from large-scale domain movements to small loop or side-chain rearrangements. In order to understand how the mechanical properties of enzymes, and the mechanical variations that are induced by ligand binding, relate to enzymatic activity, we carried out coarse-grain Brownian dynamics simulations on a set of enzymes whose structures in the unbound and ligand-bound forms are available in the Protein Data Bank. Our results show that enzymes are remarkably heterogeneous objects from a mechanical point of view and that the local rigidity of individual residues is tightly connected to their part in the protein's overall structure and function. The systematic comparison of the rigidity of enzymes in their unbound and bound forms highlights the fact that small conformational changes can induce large mechanical effects, leading to either more or less flexibility depending on the enzyme's architecture and the location of its ligand-biding site. These mechanical variations target a limited number of specific residues that occupy key locations for enzymatic activity, and our approach thus offers a mean to detect perturbation-sensitive sites in enzymes, where the addition or removal of a few interactions will lead to important changes in the proteins internal dynamics.

Mitochondrial Disease-related Mutation G167P in Cytochrome b of Rhodobacter capsulatus Cytochrome bc1 (S151P in Human) Affects the Equilibrium Distribution of [2Fe-2S] Cluster and Generation of Superoxide

Cytochrome bc1 is one of the key enzymes of many bioenergetic systems. Its operation involves a large scale movement of a head domain of iron-sulfur protein (ISP-HD), which functionally connects the catalytic quinol oxidation Qo site in cytochrome b with cytochrome c1. The Qo site under certain conditions can generate reactive oxygen species in the reaction scheme depending on the actual position of ISP-HD in respect to the Qo site. Here, using a bacterial system, we show that mutation G167P in cytochrome b shifts the equilibrium distribution of ISP-HD toward positions remote from the Qo site. This renders cytochrome bc1 non-functional in vivo. This effect is remediated by addition of alanine insertions (1Ala and 2Ala) in the neck region of the ISP subunit. These insertions, which on their own shift the equilibrium distribution of ISP-HD in the opposite direction (i.e. toward the Qo site), also act in this manner in the presence of G167P. Changes in the equilibrium distribution of ISP-HD in G167P lead to an increased propensity of cytochrome bc1 to generate superoxide, which becomes evident when the concentration of quinone increases. This result corroborates the recently proposed model in which "semireverse" electron transfer back to the Qo site, occurring when ISP-HD is remote from the site, favors reactive oxygen species production. G167P suggests possible molecular effects of S151P (corresponding in sequence to G167P) identified as a mitochondrial disease-related mutation in human cytochrome b. These effects may be valid for other human mutations that change the equilibrium distribution of ISP-HD in a manner similar to G167P.

Fast internal dynamics in alcohol dehydrogenase.

Large-scale domain motions in alcohol dehydrogenase (ADH) have been observed previously by neutron spin-echo spectroscopy (NSE). We have extended the investigation on the dynamics of ADH in solution by using high-resolution neutron time-of-flight (TOF) and neutron backscattering (BS) spectroscopy in the incoherent scattering range. The observed hydrogen dynamics were interpreted in terms of three mobility classes, which allowed a simultaneous description of the measured TOF and BS spectra. In addition to the slow global protein diffusion and domain motions observed by NSE, a fast internal process could be identified. Around one third of the protons in ADH participate in the fast localized diffusive motion. The diffusion coefficient of the fast internal motions is around two third of the value of the surrounding D2O solvent. It is tempting to associate the fast internal process with solvent exposed amino acid residues with dangling side chains.

Large-Scale Domain Motions and Pyridoxal-5'-Phosphate Assisted Radical Catalysis in Coenzyme B12-Dependent Aminomutases

Lysine 5,6-aminomutase (5,6-LAM) and ornithine 4,5-aminomutase (4,5-OAM) are two of the rare enzymes that use assistance of two vitamins as cofactors. These enzymes employ radical generating capability of coenzyme B12 (5′-deoxyadenosylcobalamin, dAdoCbl) and ability of pyridoxal-5′-phosphate (PLP, vitamin B6) to stabilize high-energy intermediates for performing challenging 1,2-amino rearrangements between adjacent carbons. A large-scale domain movement is required for interconversion between the catalytically inactive open form and the catalytically active closed form. In spite of all the similarities, these enzymes differ in substrate specificities. 4,5-OAM is highly specific for d-ornithine as a substrate while 5,6-LAM can accept d-lysine and l-β-lysine. This review focuses on recent computational, spectroscopic and structural studies of these enzymes and their implications on the related enzymes. Additionally, we also discuss the potential biosynthetic application of 5,6-LAM.

Structural basis of resistance to anti-cytochrome bc₁ complex inhibitors: implication for drug improvement.

The emergence of drug resistance has devastating economic and social consequences, a testimonial of which is the rise and fall of inhibitors against the respiratory component cytochrome bc₁ complex, a time tested and highly effective target for disease control. Unfortunately, the mechanism of resistance is a multivariate problem, including primarily mutations in the gene of the cytochrome b subunit but also activation of alternative pathways of ubiquinol oxidation and pharmacokinetic effects. There is a considerable interest in designing new bc₁ inhibitors with novel modes of binding and lower propensity to induce the development of resistance. The accumulation of crystallographic data of bc₁ complexes with and without inhibitors bound provides the structural basis for rational drug design. In particular, the cytochrome b subunit offers two distinct active sites that can be targeted for inhibition - the quinol oxidation site and the quinone reduction site. This review brings together available structural information of inhibited bc₁ by various quinol oxidation- and reductionsite inhibitors, the inhibitor binding modes, conformational changes upon inhibitor binding of side chains in the active site and large scale domain movements of the iron-sulfur protein subunit. Structural data analysis provides a clear understanding of where and why existing inhibitors fail and points towards promising alternatives.

Motions and mechanics: investigating conformational transitions in multi-domain proteins with coarse-grain simulations

Protein flexibility is an essential aspect of protein function, and proteic activity involves a wide range of structural changes, varying from small side-chain movements to large-scale domain motion. In order to understand how these large-scale rearrangements impact on proteins internal dynamics and mechanics, we carried out coarse-grain simulations on a set of proteins presenting conformational changes due to domain–domain motions, and investigated the resulting variations of their mechanical properties. These changes are highly heterogeneous along the protein sequence, and our results show that the residues undergoing important variations of their force constant occupy key positions for protein function, as they are mostly located in the ligand-binding site or on the domain–domain interface.

Thermodynamic dissection of large‐scale domain motions coupled with ligand binding of enzyme I

Domain motions are central to the biological functions of many proteins. The energetics of the motions, however, is often difficult to characterize when motions are coupled with the ligand binding. Here, we determined the thermodynamic parameters of individual domain motions and ligand binding of enzyme I (EI) using strategic domain-deletion mutants that selectively removed particular motions. Upon ligand binding, EI employs two large-scale domain motions, the hinge motion and the swivel motion, to switch between conformational states of distinct domain-domain orientations. Calorimetric analysis of the EI mutants separated the free energy changes of the binding and motions, demonstrating that the unfavorable hinge motion (ΔG = 1.5 kcal mol(-1)) was driven by the favorable swivel motion (ΔG = -5.2 kcal mol(-1)). The large free energy differences could be explained by the physicochemical nature of the domain interfaces associated with the motions; the hinge motion employed much narrower interface than the swivel motion without any hydrogen bonds or salt bridges. The small heat capacity further suggested that the packing of the domain interfaces associated with the hinge motion was less compact than that commonly observed in proteins. Lastly, thermodynamic analysis of phosphorylated EI suggests that the domain motions are regulated by the ligand binding and the phosphorylation states. Taken together, the thermodynamic dissection approach illustrates how multiple motions and ligand binding are energetically connected during the functional cycle of EI.

Redox-Linked Domain Movements in the Catalytic Cycle of Cytochrome P450 Reductase

SummaryNADPH-cytochrome P450 reductase is a key component of the P450 mono-oxygenase drug-metabolizing system. There is evidence for a conformational equilibrium involving large-scale domain motions in this enzyme. We now show, using small-angle X-ray scattering (SAXS) and small-angle neutron scattering, that delivery of two electrons to cytochrome P450 reductase leads to a shift in this equilibrium from a compact form, similar to the crystal structure, toward an extended form, while coenzyme binding favors the compact form. We present a model for the extended form of the enzyme based on nuclear magnetic resonance and SAXS data. Using the effects of changes in solution conditions and of site-directed mutagenesis, we demonstrate that the conversion to the extended form leads to an enhanced ability to transfer electrons to cytochrome c. This structural evidence shows that domain motion is linked closely to the individual steps of the catalytic cycle of cytochrome P450 reductase, and we propose a mechanism for this.

The Increasingly Complex Mechanism of HMG-CoA Reductase

HMG-CoA reductase (HMGR) is the target of statins, cholesterol-lowering drugs prescribed to millions of patients worldwide. More recent research indicates that HMGR could be a useful target in the development of antimicrobial agents. Over the last seven decades, researchers have proposed a series of increasingly complex reaction mechanisms for this biomedically important enzyme. The maturation of the mechanistic proposals for HMGR have paralleled advances in a diverse set of research areas, such as molecular biology and computational chemistry. Thus, the development of the HMGR mechanism provides a useful case study for following the advances in state-of-the-art methods in enzyme mechanism research. Similarly, the questions raised by these mechanism proposals reflect the limitations of the methods used to develop them. The mechanism of HMGR, a four-electron oxidoreductase, is unique and far more complex than originally thought. The reaction contains multiple chemical steps, coupled to large-scale domain motions of the homodimeric enzyme. The first proposals for the HMGR mechanism were based on kinetic and labeling experiments, drawing analogies to the mechanism of known dehydrogenases. Advances in molecular biology and bioinformatics enabled researchers to use site-directed mutagenesis experiments and protein sequencing to identify catalytically important glutamate, aspartate, and histidine residues. These studies, in turn, have generated new and more complicated mechanistic proposals. With the development of protein crystallography, researchers solved HMGR crystal structures to reveal an unexpected lysine residue at the center of the active site. The many crystal structures of HMGR led to increasingly complex mechanistic proposals, but the inherent limitations of the protein crystallography left a number of questions unresolved. For example, the protonation state of the glutamate residue within the active site cannot be clearly determined from the crystal structure. The differing protonation state of this residue leads to different proposed mechanisms for the enzyme. As computational analysis of large biomolecules has become more feasible, the application of methods such as hybrid quantum mechanics/molecular mechanics (QM/MM) calculations to the HMGR mechanism have led to the most detailed mechanistic proposal yet. As these methodologies continue to improve, they prove to be very powerful for the study of enzyme mechanisms in conjunction with protein crystallography. Nevertheless, even the most current mechanistic proposal for HMGR remains incomplete due to limitations of the current computational methodologies. Thus, HMGR serves as a model for how the combination of increasingly sophisticated experimental and computational methods can elucidate very complex enzyme mechanisms.

Molecular dynamics studies on the conformational transitions of adenylate kinase: a computational evidence for the conformational selection mechanism.

Escherichia coli adenylate kinase (ADK) is a monomeric phosphotransferase enzyme that catalyzes reversible transfer of phosphoryl group from ATP to AMP with a large-scale domain motion. The detailed mechanism for this conformational transition remains unknown. In the current study, we performed long time-scale molecular dynamics simulations on both open and closed states of ADK. Based on the structural analyses of the simulation trajectories, we detected over 20 times conformational transitions between the open and closed states of ADK and identified two novel conformations as intermediate states in the catalytic processes. With these findings, we proposed a possible mechanism for the large-scale domain motion of Escherichia coli ADK and its catalytic process: (1) the substrate free ADK adopted an open conformation; (2) ATP bound with LID domain closure; (3) AMP bound with NMP domain closure; (4) phosphoryl transfer occurred with ATP, and AMP converted into two ADPs, and no conformational transition was detected in the enzyme; (5) LID domain opened with one ADP released; (6) another ADP released with NMP domain open. As both open and closed states sampled a wide range of conformation transitions, our simulation strongly supported the conformational selection mechanism for Escherichia coli ADK.

Application of Systems Biology and Bioinformatics Methods in Biochemistry and Biomedicine

Application of Systems Biology and Bioinformatics Methods in Biochemistry and Biomedicine