Conformational selection and induced changes along the catalytic cycle of Escherichia coli dihydrofolate reductase

Protein function often involves changes between different conformations. Central questions are how these conformational changes are coupled to the binding or catalytic processes during which they occur, and how they affect the catalytic rates of enzymes. An important model system is the enzyme dihydrofolate reductase (DHFR) from Escherichia coli, which exhibits characteristic conformational changes of the active‐site loop during the catalytic step and during unbinding of the product. In this article, we present a general kinetic framework that can be used (1) to identify the ordering of events in the coupling of conformational changes, binding, and catalysis and (2) to determine the rates of the substeps of coupled processes from a combined analysis of nuclear magnetic resonance R₂ relaxation dispersion experiments and traditional enzyme kinetics measurements. We apply this framework to E. coli DHFR and find that the conformational change during product unbinding follows a conformational‐selection mechanism, that is, the conformational change occurs predominantly prior to unbinding. The conformational change during the catalytic step, in contrast, is an induced change, that is, the change occurs after the chemical reaction. We propose that the reason for these conformational changes, which are absent in human and other vertebrate DHFRs, is robustness of the catalytic rate against large pH variations and changes to substrate/product concentrations in E. coli. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.     

Evaluating conformational changes in protein structures binding RNA

Many protein-RNA recognition events are known to exhibit conformational changes from qualitative observations of individual complexes. However, a quantitative estimation of conformational changes is required if protein-RNA docking and template-based methods for RNA binding site prediction are to be developed. This study presents the first quantitative evaluation of conformational changes that occur when proteins bind RNA. The analysis of twelve RNA-binding proteins in the bound and unbound states using error-scaled difference distance matrices is presented. The binding site residues are mapped to each structure, and the conformational changes that affect these residues are evaluated. Of the twelve proteins four exhibit greater movements in nonbinding site residues, and a further four show the greatest movements in binding site residues. The remaining four proteins display no significant conformational change. When interface residues are found to be in conformationally variable regions of the protein they are typically seen to move less than 2 A between the bound and unbound conformations. The current data indicate that conformational changes in the binding site residues of RNA binding proteins may not be as significant as previously suggested, but a larger data set is required before wider conclusions may be drawn. The implications of the observed conformational changes for protein function prediction are discussed.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Requirements for peptidyl-prolyl isomerization activity: a comprehensive mutational analysis of the substrate-binding cavity of FK506-binding protein 12

Peptidyl-prolyl isomerase (PPIase) activity is exhibited by many proteins belonging to the PPIase family. However, the catalytic mechanism of this activity remains to be completely elucidated. Here, we selected human FK506-binding protein 12 (FKBP12) as the model PPIase and investigated the nature of amino acid residues essential for the activity. The crystal structures of several complexes of PPIase with short peptides revealed that the residues Asp37, Arg42, Phe46, Val55, Trp59, and Tyr82 in the substrate-binding cavity of FKBP12 appear to play key roles in the PPIase activity. Each of these six residues was substituted by 20 common amino acid residues. The activity of each mutant protein was measured using a peptide analog by the chymotrypsin digestion assay and then compared with wild-type FKBP12. It was found that site-specific interactions by the side chains of amino acid residues constituting the substrate-binding cavity were not essential for the PPIase activity, although the 37th, 55th, and 82nd amino acid residues significantly contributed to the activity. This suggests that the PPIase activity requires only the hydrophobic cavity that captures the Pro-containing peptide.     

A comparative study of the backbone dynamics of two closely related lipid binding proteins: Bovine heart fatty acid binding protein and porcine ileal lipid binding protein

The backbone dynamics of bovine heart fatty acid binding protein (H-FABP) and porcine ileal lipid binding protein (ILBP) were studied by 15N NMR relaxation (T1 and T2) and steady state heteronuclear 15N{1H} NOE measurements. The microdynamic parameters characterizing the backbone mobility were determined using the model-free approach. For H-FABP, the non-terminal backbone amide groups display a rather compact protein structure of low flexibility. In contrast, for ILBP an increased number of backbone amide groups display unusually high internal mobility. Furthermore, the data indicate a higher degree of conformational exchange processes in the sec-msec time range for ILBP compared to H-FABP. These results suggest significant differences in the conformational stability for these two structurally highly homologous members of the fatty acid binding protein family.     

The differences in the microenvironment of the two tryptophan residues of the glutamine-binding protein from Escherichia coli shed light on the binding properties and the structural dynamics of the protein

Glutamine-binding protein (GlnBP) from Escherichia coli is a monomer (26 kDa) that is responsible for the first step in the active transport of L-glutamine across the cytoplasmic membrane. GlnBP consists of two domains (termed large and small) linked by two antiparallel beta-strands. The large domain is similar to the small domain but it contains two additional alpha-helices and three more short antiparallel beta-strands. The deep cleft formed between the two domains contains the ligand-binding site. The binding of L-glutamine leads to cleft closing and a significant structural change with the formation of the so-called "closed form" structure. The protein contains two tryptophan residues (W32 and W220) and 10 tyrosine residues. We used phosphorescence spectroscopy measurements to characterize the role of the two tryptophan residues in the protein structure in the absence and the presence of glutamine. Our results pointed out that the phosphorescence of GlnBP is easily detected in fluid solutions where the emission of the two tryptophan residues is readily discriminated by the drastic difference in the phosphorescence lifetime allowing the assignments of the short lifetime to W220 and the long lifetime to W32. In addition, our results showed that the triplet lifetime of the superficial W220 is unusually short because of intramolecular quenching by the proximal Y163. On the contrary, the lifetime of W32 is several hundred milliseconds long, implicating a well-ordered, compact fold of the surrounding polypeptide. The spectroscopic data were analyzed and discussed together with a detailed inspection of the 3D structure of GlnBP.     

Thermodynamic analysis of halide binding to haloalkane dehalogenase suggests the occurrence of large conformational changes

Haloalkane dehalogenase (DhlA) hydrolyzes short-chain haloalkanes to produce the corresponding alcohols and halide ions. Release of the halide ion from the active-site cavity can proceed via a two-step and a three-step route, which both contain slow enzyme isomerization steps. Thermodynamic analysis of bromide binding and release showed that the slow unimolecular isomerization steps in the three-step bromide export route have considerably larger transition state enthalpies and entropies than those in the other route. This suggests that the three-step route involves different and perhaps larger conformational changes than the two-step export route. We propose that the three-step halide export route starts with conformational changes that result in a more open configuration of the active site from which the halide ion can readily escape. In addition, we suggest that the two-step route for halide release involves the transfer of the halide ion from the halide-binding site in the cavity to a binding site somewhere at the protein surface, where a so-called collision complex is formed in which the halide ion is only weakly bound. No large structural rearrangements are necessary for this latter process.     

Structure and functional analysis of the siderophore periplasmic binding protein from the fuscachelin gene cluster of Thermobifida fusca

Iron acquisition is a complex, multicomponent process critical for most organisms' survival and virulence. Small iron chelating molecules, siderophores, mediate transport as key components of common pathways for iron assimilation in many microorganisms. The chemistry and biology of the extraordinary tight and specific metal binding siderophores is of general interest in terms of host/guest chemistry and is a potential target toward the development of therapeutic treatments for microbial virulence. The siderophore pathway of the moderate thermophile, Thermobifida fusca, is an excellent model system to study the process in Gram‐positive bacteria. Here we describe the structure and characterization of the siderophore periplasmic binding protein, FscJ from the fuscachelin gene cluster of T. fusca. The structure shows a di‐domain arrangement connected with a long α‐helix hinge. Several X‐ray structures detail ligand‐free conformational changes at different pH values, illustrating complex interdomain flexibility of the siderophore receptors. We demonstrated that FscJ has a unique recognition mechanism and details the binding interaction with ferric‐fuscachelin A through ITC and docking analysis. The presented work provides a structural basis for the complex molecular mechanisms of siderophore recognition and transportation. Proteins 2016; 84:118–128. © 2015 Wiley Periodicals, Inc.     

Carnitine palmitoyltransferase activities: Effects of serum albumin,acyl-CoA binding protein and fatty acid binding protein

The carnitine palmitoyltransferase activity of various subcellular preparations measured with octanoyl-CoA as substrate was markedly increased by bovine serum albumin at low M concentrations of octanoyl-CoA. However, even a large excess (500 M) of this acyl-CoA did not inhibit the activity of the mitochondrial outer carnitine palmitoyltransferase, a carnitine palmitoyltransferase isoform that is particularly sensitive to inhibition by low M concentrations of palmitoyl-CoA. This bovine serum albumin stimulation was independent of the salt activation of the carnitine palmitoyltransferase activity. The effects of acyl-CoA binding protein (ACBP) and the fatty acid binding protein were also examined with palmitoyl-CoA as substrate. The results were in line with the findings of stronger binding of acyl-CoA to ACBP but showed that fatty acid binding protein also binds acyl-CoA esters. Although the effects of these proteins on the outer mitochondrial carnitine palmitoyltransferase activity and its malonyl-CoA inhibition varied with the experimental conditions, they showed that the various carnitine palmitoyltransferase preparations are effectively able to use palmitoyl-CoA bound to ACBP in a near physiological molar ratio of 1:1 as well as that bound to the fatty acid binding protein. It is suggested that the three proteins mentioned above effect the carnitine palmitoyltransferase activities not only by binding of acyl-CoAs, preventing acyl-CoA inhibition, but also by facilitating the removal of the acylcarnitine product from carnitine palmitoyltransferase. These results support the possibility that the acyl-CoA binding ability of acyl-CoA binding protein and of fatty acid binding protein have a role in acyl-CoA metabolismin vivo.Abbreviations ACBP acyl-CoA binding protein - BSA bovine serum albumin - CPT carnitine palmitoyltransferase - CPT₀ malonyl-CoA sensitive CPT of the outer mitochondrial membrane - CPT malonyl-CoA insensitive CPT of the inner mitochondrial membrane - OG octylglucoside - OMV outer membrane vesicles - IMV inner membrane vesicles Affiliated to the Department of Experimental Medicine, University of Montreal     

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may “select” protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational‐selection and induced‐change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational‐selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced‐change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on‐ and off‐rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single‐molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational‐selection and induced‐change processes in both binding and unbinding direction.     

A method for determining binding kinetics applied to thiamine-binding protein

A rapid and simple method for assaying the binding activity of thiamine-binding protein is described. By this assay method, the binding characteristics of rice bran thiamine-binding protein have been evaluated with [14C]thiamine as ligand. Analysis of these data by Scatchard plot resulted in linear plots giving a dissociation constant (Kd) for thiamine of 0.55 microM and a maximum binding (Bmax) of 14.5 pmol of ligand bound/microgram of protein. Thiamine binding to the binding protein was time dependent and reached equilibrium at approximately 20 min. The Kob was 0.18 min-1 and the k1 was 1.25 X 10(5) min-1 M-1. Reversibility of thiamine binding at equilibrium was completed at 60 min with a k2 value of 0.052 min-1. The Kd calculated from the reverse rate constant was 0.42 microM. These results indicated that this binding assay method was substantially reliable and accurate.     

Probing kinetic drug binding mechanism in voltage-gated sodium ion channel: open state versus inactive state blockers

The kinetics and nonequilibrium thermodynamics of open state and inactive state drug binding mechanisms have been studied here using different voltage protocols in sodium ion channel. We have found that for constant voltage protocol, open state block is more efficient in blocking ionic current than inactive state block. Kinetic effect comes through peak current for mexiletine as an open state blocker and in the tail part for lidocaine as an inactive state blocker. Although the inactivation of sodium channel is a free energy driven process, however, the two different kinds of drug affect the inactivation process in a different way as seen from thermodynamic analysis. In presence of open state drug block, the process initially for a long time remains entropy driven and then becomes free energy driven. However in presence of inactive state block, the process remains entirely entropy driven until the equilibrium is attained. For oscillating voltage protocol, the inactive state blocking is more efficient in damping the oscillation of ionic current. From the pulse train analysis it is found that inactive state blocking is less effective in restoring normal repolarisation and blocks peak ionic current. Pulse train protocol also shows that all the inactive states behave differently as one inactive state responds instantly to the test pulse in an opposite manner from the other two states.     

Probing kinetic drug binding mechanism in voltage-gated sodium ion channel: open state versus inactive state blockers

The kinetics and nonequilibrium thermodynamics of open state and inactive state drug binding mechanisms have been studied here using different voltage protocols in sodium ion channel. We have found that for constant voltage protocol, open state block is more efficient in blocking ionic current than inactive state block. Kinetic effect comes through peak current for mexiletine as an open state blocker and in the tail part for lidocaine as an inactive state blocker. Although the inactivation of sodium channel is a free energy driven process, however, the two different kinds of drug affect the inactivation process in a different way as seen from thermodynamic analysis. In presence of open state drug block, the process initially for a long time remains entropy driven and then becomes free energy driven. However in presence of inactive state block, the process remains entirely entropy driven until the equilibrium is attained. For oscillating voltage protocol, the inactive state blocking is more efficient in damping the oscillation of ionic current. From the pulse train analysis it is found that inactive state blocking is less effective in restoring normal repolarisation and blocks peak ionic current. Pulse train protocol also shows that all the inactive states behave differently as one inactive state responds instantly to the test pulse in an opposite manner from the other two states.     

Probing flexibility and "induced-fit" phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations

Aldose reductase is a promising target for the treatment of diabetic complications, and as such, has become the focus of various drug design projects. As revealed by a survey of available crystal structures, the protein shows pronounced induced-fit effects upon ligand binding. Although helping to explain the enzyme's substrate promiscuity, phenomena of this kind are still responsible for significant complications in structure-based design efforts directed to aldose reductase. Accordingly, a deeper understanding of the principles governing conformational alterations in this enzyme would be of utmost practical importance. As a first step in addressing this issue, molecular dynamics (MD) simulations have been carried out. The ultrahigh resolution crystal structure of aldose reductase complexed with inhibitor IDD594 served as ideal starting point for a set of different simulations of nanosecond time scale: the native complexed state with bound inhibitor, the uncomplexed state (after removal of the inhibitor) at standard temperature, and the uncomplexed state at elevated temperature. The reference simulation of the complex exhibits extraordinary stability of the overall fold, whereas two distinct conformational substates are found for the binding-site region. In contrast, already at standard temperature pronounced changes are observed in the binding region during the simulation of the uncomplexed state. Leu300, for example, closes the access to the pocket opened by IDD594. On the other hand, conformations around the catalytic site are highly conserved, with the His110-Tyr48-NADP+ orientation being stabilized by a water molecule. Detailed analysis of the trajectories allows to reveal a set of distinct conformational substates that may prove useful as alternative structural templates in virtual screening for new aldose reductase inhibitors.     

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.     

Purification and partial peptide sequence analysis of the boar proacrosin binding protein

Boar proacrosin binding protein has been purified and the partial peptide sequence of the CNBr‐digested proacrosin binding protein has been determined. Proacrosin binding protein was purified as a proacrosin and proacrosin binding protein complex from the acid extracts of boar spermatozoa through gel filtration. After the proacrosin binding protein was dissociated from proacrosin by freeze‐thaw method, the proacrosin binding protein was purified through gel filtration. Fractions containing the proacrosin binding protein were pooled and were concentrated by lyophilization and then subjected to CNBr digestion. Four major CNBr‐digested peptides were subjected to N‐terminal peptide sequencing. All four showed the same N‐terminus sequence. Mol. Reprod. Dev. 54:76–80, 1999. © 1999 Wiley‐Liss, Inc.     

Calculating pKa values in the cAMP-dependent protein kinase: the effect of conformational change and ligand binding

The conformational change observed upon ligand binding and phosphorylation for the cAMP-dependent protein kinase (protein kinase A-PKA) is of high importance for the regulation of its activity. We calculate pKa values and net charges for 18 3D structures of PKA in various conformations and liganded states to examine the role of electrostatics in ligand binding and activation. We find that the conformational change of PKA takes place without any significant net proton uptake/release at all pH values, thus indicating that PKA has evolved to reduce any pH-dependent barriers to the conformational motion. We furthermore find that the binding of ligands induces large changes in the net charge of PKA at most pH values, but significantly, we find that the net charge difference at physiological pH is close to zero, thus indicating that the active-site pKa values have been preorganized for substrate binding. We are unable to unequivocally resolve the identity of the groups responsible for determining the pH-activity profile of PKA but speculate that the titration of Lys 168 or the titration of ATP itself could be responsible for the loss of activity at high pH values. Finally, we examine the effect of point mutations on the pKa values of the PKA catalytic residues and find these to be relatively insensitive to both noncharge-altering and charge-altering mutations.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma