Ligand binding and protein dynamics in lactate dehydrogenase

Recent experimental studies suggest that lactate dehydrogenase (LDH) binds its substrate via the formation of a LDH/NADH.substrate encounter complex through a select-fit mechanism, whereby only a minority population of LDH/NADH is binding-competent. In this study, we perform molecular dynamics calculations to explore the variations in structure accessible to the binary complex with a focus on identifying structures that seem likely to be binding-competent and which are in accord with the known experimental characterization of forming binding-competent species. We find that LDH/NADH samples quite a range of protein conformations within our 2.148 ns calculations, some of which yield quite facile access of solvent to the active site. The results suggest that the mobile loop of LDH is perhaps just partially open in these conformations and that multiple open conformations, yielding multiple binding pathways, are likely. These open conformations do not require large-scale unfolding/melting of the binary complex. Rather, open versus closed conformations are due to subtle protein and water rearrangements. Nevertheless, the large heat capacity change observed between binding-competent and binding-incompetent can be explained by changes in solvation and an internal rearrangement of hydrogen bonds. We speculate that such a strategy for binding may be necessary to get a ligand efficiently to a binding pocket that is located fairly deep within the protein's interior.     

Interdomain interactions in hinge-bending transitions.

BACKGROUND: The mechanisms that allow or constrain protein movement have not been understood. Here we study interdomain interactions in proteins to investigate hinge-bending motions. RESULTS: We find a limited number of salt bridges and hydrogen bonds at the interdomain interface, in both the "closed" and the "open" conformations. Consistently, analysis of 222 salt bridges in an independently selected database indicates that most salt bridges form within rather than between independently folding hydrophobic units. Calculations show that these interdomain salt bridges either destabilize or only marginally stabilize the closed conformation in most proteins. In contrast, the nonpolar buried surface area between the moving parts can be extensive in the closed conformations. However, when the nonpolar buried surface area is large, we find that at the interdomain interface in the open conformation it may be as large or larger than in the closed conformation. Hence, the energetic penalty of opening the closed conformation is overcome. Consistently, a large nonpolar surface area buried in the closed interdomain interface accompanies limited opening of the domains, yielding a larger interface. CONCLUSIONS: Short-range electrostatic interactions are largely absent between moving domains. Interdomain nonpolar buried surface area may be large in the closed conformation, but it is largely offset by the area buried in the open conformation. In such cases the opening of the domains appears to be relatively small. This may allow prediction of the extent of domain opening. Such predictions may have implications for the shape and size of the binding pockets in drug/protein design.     

Use of normal modes for structural modeling of proteins: the case study of rat heme oxygenase 1

We present an original approach based on full-atom normal mode analysis (NMA) aimed to expand the general framework of homology modeling. Using the rat heme-free oxygenase 1 as a case system, we show how NMA can be used to model different physiologically relevant conformations of the same protein. Starting from a unique heme-bound X-ray structure, and using two structural templates corresponding to a human and an incomplete rat heme-free structures, we generate models of the rat unbound species with open and closed conformations. Less than 100 lowest frequency modes of the target were sufficient to obtain the heme-free conformations, the closest to the templates. The rat HO-1 model built for the open form shows features similar to the open form of the human heme-free oxygenase, and the one built for the closed form was similar to the incompletely resolved X-ray structure of the same protein available in the Protein DataBank. In the latter case, the use of NMA was particularly useful since it allowed to build a complete structure and therefore to discuss on the reason of the structural differences between open and closed forms. This study shows that the amount of main chain flexibility provided by the normal modes can lead to major improvements in homology modeling approaches. Such applications will allow the characterization of alternative conformations of a target protein with respect to the templates and/or the construction of good quality 3D models based on existing templates with unresolved parts in their tertiary structure.     

Structural bioinformatics study of EPSP synthase from Mycobacterium tuberculosis

The shikimate pathway is an attractive target for herbicides and antimicrobial agent development because it is essential in algae, higher plants, bacteria, and fungi, but absent from mammals. Homologues to enzymes in the shikimate pathway have been identified in the genome sequence of Mycobacterium tuberculosis. Among them, the EPSP synthase was proposed to be present by sequence homology. Accordingly, in order to pave the way for structural and functional efforts towards anti-mycobacterial agent development, here we describe the molecular modeling of 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase isolated from M. tuberculosis that should provide a structural framework on which the design of specific inhibitors may be based on. Significant differences in the relative orientation of the domains in the two models result in "open" and "closed" conformations. The possible relevance of this structural transition in the ligand biding is discussed.     

Insights into the Structure of Invisible Conformations of Large Methyl Group Labeled Molecular Machines from High Pressure NMR

Most proteins are highly flexible and can adopt conformations that deviate from the energetically most favorable ground state. Structural information on these lowly populated, alternative conformations is often lacking, despite the functional importance of these states. Here, we study the pathway by which the Dcp1:Dcp2 mRNA decapping complex exchanges between an autoinhibited closed and an open conformation. We make use of methyl Carr–Purcell–Meiboom–Gill (CPMG) NMR relaxation dispersion (RD) experiments that report on the population of the sparsely populated open conformation as well as on the exchange rate between the two conformations. To obtain volumetric information on the open conformation as well as on the transition state structure we made use of RD measurements at elevated pressures. We found that the open Dcp1:Dcp2 conformation has a lower molecular volume than the closed conformation and that the transition state is close in volume to the closed state. In the presence of ATP the volume change upon opening of the complex increases and the volume of the transition state lies in-between the volumes of the closed and open state. These findings show that ATP has an effect on the volume changes that are associated with the opening-closing pathway of the complex. Our results highlight the strength of pressure dependent NMR methods to obtain insights into structural features of protein conformations that are not directly observable. As our work makes use of methyl groups as NMR probes we conclude that the applied methodology is also applicable to high molecular weight complexes.     

Exploring the role of large conformational changes in the fidelity of DNA polymerase beta

The relationships between the conformational landscape, nucleotide insertion catalysis and fidelity of DNA polymerase beta are explored by means of computational simulations. The simulations indicate that the transition states for incorporation of right (R) and wrong (W) nucleotides reside in substantially different protein conformations. The protein conformational changes that reproduce the experimentally observed fidelity are significantly larger than the small rearrangements that usually accompany motions from the reactant state to the transition state in common enzymatic reactions. Once substrate binding has occurred, different constraints imposed on the transition states for insertion of R and W nucleotides render it highly unlikely that both transition states can occur in the same closed structure, because the predicted fidelity would then be many orders of magnitude too large. Since the conformational changes reduce the transition state energy of W incorporation drastically they decrease fidelity rather than increase it. Overall, a better agreement with experimental data is attained when the R is incorporated through a transition state in a closed conformation and W is incorporated through a transition state in one or perhaps several partially open conformations. The generation of free energy surfaces for R and W also allow us to analyze proposals about the relationship between induced fit and fidelity.     

Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.     

Large amplitude conformational change in proteins explored with a plastic network model: adenylate kinase

The plastic network model (PNM) is used to generate a conformational change pathway for Escherichia coli adenylate kinase based on two crystal structures, namely that of an open and a closed conformer. In this model, the energy basins corresponding to known conformers are connected at their lowest common energies. The results are used to evaluate and analyze the minimal energy pathways between these basins. The open to closed transition analysis provides an identification of hinges that is in agreement with the existing definitions based on the available X-ray structures. The elastic energy distribution and the C(alpha) pseudo-dihedral variation provide similar information on these hinges. The ensemble of the 45 published structures for this protein and closely related proteins is shown to always be within 3.0 A of the pathway, which corresponds to a conformational change between two end structures that differ by a C(alpha)-atom root-mean-squared deviation of 7.1A.     

Conformational flexibility of the leucine binding protein examined by protein domain coarse-grained molecular dynamics

Periplasmic binding proteins are the initial receptors for the transport of various substrates over the inner membrane of gram-negative bacteria. The binding proteins are composed of two domains, and the substrate is entrapped between these domains. For several of the binding proteins it has been established that a closed-up conformation exists even without substrate present, suggesting a highly flexible apo-structure which would compete with the ligand-bound protein for the transporter interaction. For the leucine binding protein (LBP), structures of both open and closed conformations are known, but no closed-up structure without substrate has been reported. Here we present molecular dynamics simulations exploring the conformational flexibility of LBP. Coarse grained models based on the MARTINI force field are used to access the microsecond timescale. We show that a standard MARTINI model cannot maintain the structural stability of the protein whereas the ELNEDIN extension to MARTINI enables simulations showing a stable protein structure and nanosecond dynamics comparable to atomistic simulations, but does not allow the simulation of conformational flexibility. A modification to the MARTINI-ELNEDIN setup, referred to as domELNEDIN, is therefore presented. The domELNEDIN setup allows the protein domains to move independently and thus allows for the simulation of conformational changes. Microsecond domELNEDIN simulations starting from either the open or the closed conformations consistently show that also for LBP, the apo-structure is flexible and can exist in a closed form.

Direct Observation of T4 Lysozyme Hinge-Bending Motion by Fluorescence Correlation Spectroscopy

Bacteriophage T4 Lysozyme (T4L) catalyzes the hydrolysis of the peptidoglycan layer of the bacterial cell wall late in the infection cycle. It has long been postulated that equilibrium dynamics enable substrate access to the active site located at the interface between the N- and C-terminal domains. Crystal structures of WT-T4L and point mutants captured a range of conformations that differ by the hinge-bending angle between the two domains. Evidence of equilibrium between open and closed conformations in solution was gleaned from distance measurements between the two domains but the nature of the equilibrium and the timescale of the underlying motion have not been investigated. Here, we used fluorescence fluctuation spectroscopy to directly detect T4L equilibrium conformational fluctuations in solution. For this purpose, Tetramethylrhodamine probes were introduced at pairs of cysteines in regions of the molecule that undergo relative displacement upon transition from open to closed conformations. Correlation analysis of Tetramethylrhodamine intensity fluctuations reveals hinge-bending motion that changes the relative distance and orientation of the N- and C-terminal domains with ≅15 μs relaxation time. That this motion involves interconversion between open and closed conformations was further confirmed by the dampening of its amplitude upon covalent substrate trapping. In contrast to the prevalent two-state model of T4L equilibrium, molecular brightness and number of particles obtained from cumulant analysis suggest that T4L populates multiple intermediate states, consistent with the wide range of hinge-bending angles trapped in the crystal structure of T4L mutants.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

Human P-glycoprotein is active when the two halves are clamped together in the closed conformation

The P-glycoprotein (P-gp, ABCB1) drug pump protects us from toxic compounds and confers multidrug resistance. Each of the two homologous halves of P-gp is composed of a transmembrane domain (TMD) with six TM segments followed by a nucleotide-binding domain (NBD). The drug- and ATP-binding sites reside at the interface between the TMDs and NBDs, respectively. Crystal structures show drug pumps in the open and closed conformations, where the drug-binding pocket and NBDs are open or closed at the cytoplasmic side, respectively. Although it has been postulated that drug substrates enter the drug-binding pocket in the open conformation, it is unknown if they can enter in the closed conformation. To determine this, we introduced cysteines into regions of TM3 (residues 175-178) and TM9 (residues 820-822) that extend into the cytoplasm and are 4 Å and 20 Å apart in the closed and open conformations, respectively. The 12 double cysteine mutants were then cross-linked with a short cross-linker, M1M (4 Å) at 0 °C to reduce thermal motion in the protein. Only mutant L175C/N820C was cross-linked. Cross-linking was not increased in the presence of ATP or drug substrates. Cross-linking increased its basal ATPase activity about 3-fold. Activity could be increased further by drug substrates such as verapamil and rhodamine B. These results suggest that P-gp in the membrane is in the closed conformation that has a high affinity for drug substrates.     

Exploring conformational states and helical packings in the P2X receptor transmembrane domain by molecular dynamics simulation

The P2X receptor is a trimeric transmembrane protein that acts as an ATP-gated ion channel. Its transmembrane domain (TMD) contains only six helices and three of them, the M2 helices, line the ion conduction pathway. Here, using molecular dynamics simulation, I identify four conformational states of the TMD that are associated with four types of packing between M2 helices. Packing in the extracellular half of the M2 helix produces closed conformations, while packing in the intracellular half produces both open and closed conformations. State transition is observed and supports a mechanism where iris-like twisting of the M2 helices switches the location of helical packing between the extracellular and the intracellular halves of the helices. In addition, this twisting motion alters the position and orientation of residue side-chains relative to the pore and therefore influences the pore geometry and possibly ion permeation. Helical packing, on the other hand, may restrict the twisting motion and generate discrete conformational states.     

Rapid Sampling of Molecular Motions with Prior Information Constraints

Proteins are active, flexible machines that perform a range of different functions. Innovative experimental approaches may now provide limited partial information about conformational changes along motion pathways of proteins. There is therefore a need for computational approaches that can efficiently incorporate prior information into motion prediction schemes. In this paper, we present PathRover, a general setup designed for the integration of prior information into the motion planning algorithm of rapidly exploring random trees (RRT). Each suggested motion pathway comprises a sequence of low-energy clash-free conformations that satisfy an arbitrary number of prior information constraints. These constraints can be derived from experimental data or from expert intuition about the motion. The incorporation of prior information is very straightforward and significantly narrows down the vast search in the typically high-dimensional conformational space, leading to dramatic reduction in running time. To allow the use of state-of-the-art energy functions and conformational sampling, we have integrated this framework into Rosetta, an accurate protocol for diverse types of structural modeling. The suggested framework can serve as an effective complementary tool for molecular dynamics, Normal Mode Analysis, and other prevalent techniques for predicting motion in proteins. We applied our framework to three different model systems. We show that a limited set of experimentally motivated constraints may effectively bias the simulations toward diverse predicates in an outright fashion, from distance constraints to enforcement of loop closure. In particular, our analysis sheds light on mechanisms of protein domain swapping and on the role of different residues in the motion.     

Translocation channel gating kinetics balances protein translocation efficiency with signal sequence recognition fidelity

The transition between the closed and open conformations of the Sec61 complex permits nascent protein insertion into the translocation channel. A critical event in this structural transition is the opening of the lateral translocon gate that is formed by four transmembrane (TM) spans (TM2, TM3, TM7, and TM8 in Sec61p) to expose the signal sequence-binding site. To gain mechanistic insight into lateral gate opening, mutations were introduced into a lumenal loop (L7) that connects TM7 and TM8. The sec61 L7 mutants were found to have defects in both the posttranslational and cotranslational translocation pathways due to a kinetic delay in channel gating. The translocation defect caused by L7 mutations could be suppressed by the prl class of sec61 alleles, which reduce the fidelity of signal sequence recognition. The prl mutants are proposed to act by destabilizing the closed conformation of the translocation channel. Our results indicate that the equilibrium between the open and closed conformations of the protein translocation channel maintains a balance between translocation activity and signal sequence recognition fidelity.     

The structure of the iron-binding protein, FutA1, from Synechocystis 6803

Cyanobacteria account for a significant percentage of aquatic primary productivity even in areas where the concentrations of essential micronutrients are extremely low. To better understand the mechanism of iron selectivity and transport, the structure of the solute binding domain of an ATP binding cassette iron transporter, FutA1, was determined in the presence and absence of iron. The iron ion is bound within the "C-clamp" structure via four tyrosine and one histidine residues. There are extensive interactions between these ligating residues and the rest of the protein such that the conformations of the side chains remain relatively unchanged as the iron is released by the opening of the metal binding cleft. This is in stark contrast to the zinc-binding protein, ZnuA, where the domains of the metal-binding protein remain relatively fixed, whereas the ligating residues rotate out of the binding pocket upon metal release. The rotation of the domains in FutA1 is facilitated by two flexible beta-strands running along the back of the protein that act like a hinge during domain motion. This motion may require relatively little energy since total contact area between the domains is the same whether the protein is in the open or closed conformation. Consistent with the pH dependence of iron binding, the main trigger for iron release is likely the histidine in the iron-binding site. Finally, neither FutA1 nor FutA2 binds iron as a siderophore complex or in the presence of anions, and both preferentially bind ferrous over ferric ions.     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

Cooperative transition between open and closed conformations in potassium channels

Potassium (K+) ion channels switch between open and closed conformations. The nature of this important transition was revealed by comparing the X-ray crystal structures of the MthK channel from Methanobacterium thermoautotrophicum, obtained in its open conformation, and the KcsA channel from Streptomyces lividans, obtained in its closed conformation. We analyzed the dynamic characteristics and energetics of these homotetrameric structures in order to study the role of the intersubunit cooperativity in this transition. For this, elastic models and in silico alanine-scanning mutagenesis were used, respectively. Reassuringly, the calculations manifested motion from the open (closed) towards the closed (open) conformation. The calculations also revealed a network of dynamically and energetically coupled residues. Interestingly, the network suggests coupling between the selectivity filter and the gate, which are located at the two ends of the channel pore. Coupling between these two regions was not observed in calculations that were conducted with the monomer, which emphasizes the importance of the intersubunit interactions within the tetrameric structure for the cooperative gating behavior of the channel.     

Using motion planning to map protein folding landscapes and analyze folding kinetics of known native structures.

We investigate a novel approach for studying the kinetics of protein folding. Our framework has evolved from robotics motion planning techniques called probabilistic roadmap methods (PRMs) that have been applied in many diverse fields with great success. In our previous work, we presented our PRM-based technique and obtained encouraging results studying protein folding pathways for several small proteins. In this paper, we describe how our motion planning framework can be used to study protein folding kinetics. In particular, we present a refined version of our PRM-based framework and describe how it can be used to produce potential energy landscapes, free energy landscapes, and many folding pathways all from a single roadmap which is computed in a few hours on a desktop PC. Results are presented for 14 proteins. Our ability to produce large sets of unrelated folding pathways may potentially provide crucial insight into some aspects of folding kinetics, such as proteins that exhibit both two-state and three-state kinetics that are not captured by other theoretical techniques.