The critical role of the loops of triosephosphate isomerase for its oligomerization, dynamics, and functionality

Anti-Leishmania activity and molecular docking of unusual flavonoids-rich fraction from Arrabidaea brachypoda (Bignoniaceae)

Leishmaniases comprise a group of infectious parasitic diseases caused by various species of Leishmania and are considered a significant public health problem worldwide. Only a few medications, including miltefosine, amphotericin B, and meglumine antimonate, are used in current therapy. These medications are associated with severe side effects, low efficacy, high cost, and the need for hospital support. Additionally, there have been occurrences of drug resistance. Additionally, only a limited number of drugs, such as meglumine antimonate, amphotericin B, and miltefosine, are available, all of which are associated with severe side effects. In this context, the need for new effective drugs with fewer adverse effects is evident. Therefore, this study investigated the anti-Leishmania activity of a dichloromethane fraction (DCMF) extracted from Arrabidaea brachypoda roots. This fraction inhibited the viability of L. infantum, L. braziliensis, and L. Mexicana promastigotes, with IC50 values of 10.13, 11.44, and 11.16 µg/mL, respectively, and against L. infantum amastigotes (IC50 = 4.81 µg/mL). Moreover, the DCMF exhibited moderate cytotoxicity (CC50 = 25.15) towards RAW264.7 macrophages, with a selectivity index (SI) of 5.2. Notably, the DCMF caused damage to the macrophage genome only at 40 µg/mL, which is greater than the IC50 found for all Leishmania species. The results suggest that DCMF demonstrates similar antileishmanial effectiveness to isolated brachydin B, without causing genotoxic effects on mammalian cells. This finding is crucial because the isolation of the compounds relies on several steps and is very costly while obtaining the DCMF fraction is a simple and cost-effective process. Furthermore, In addition, the potential mechanisms of action of brachydins were also investigated. The computational analysis indicates that brachydin compounds bind to the Triosephosphate isomerase (TIM) enzyme via two main mechanisms: destabilizing the interface between the homodimers and interacting with catalytic residues situated at the site of binding. Based on all the results, DCMF exhibits promise as a therapeutic agent for leishmaniasis due to its significantly reduced toxicity in comparison to the adverse effects associated with current reference treatments.

Ribosome External Electric Field Regulates Metabolic Enzyme Activity: The RAMBO Effect

Ribosomes bind to many metabolic enzymes and change their activity. A general mechanism for ribosome-mediated amplification of metabolic enzyme activity, RAMBO, was formulated and elucidated for the glycolytic enzyme triosephosphate isomerase, TPI. The RAMBO effect results from a ribosome-dependent electric field-substrate dipole interaction energy that can increase or decrease the ground state of the reactant and product to regulate catalytic rates. NMR spectroscopy was used to determine the interaction surface of TPI binding to ribosomes and to measure the corresponding kinetic rates in the absence and presence of intact ribosome particles. Chemical cross-linking and mass spectrometry revealed potential ribosomal protein binding partners of TPI. Structural results and related changes in TPI energetics and activity show that the interaction between TPI and ribosomal protein L11 mediate the RAMBO effect.

Motif-VI loop acts as a nucleotide valve in the West Nile Virus NS3 Helicase

The Orthoflavivirus NS3 helicase (NS3h) is crucial in virus replication, representing a potential drug target for pathogenesis. NS3h utilizes nucleotide triphosphate (ATP) for hydrolysis energy to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. Intermediate states along the ATP hydrolysis cycle and conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. Extensive molecular dynamics simulations of West Nile virus NS3h+ssRNA in the apo, ATP, ADP+Pi and ADP bound states were used to model the conformational ensembles along this cycle. Energetic and structural clustering analyses depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). Based on these results, MVIL mutants (D471L, D471N and D471E) were found to have a substantial reduction in ATPase activity and RNA replication compared to the wild-type. Simulations of the mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

Motif-VI Loop Acts as a Nucleotide Valve in the West Nile Virus NS3 Helicase

The flavivirus NS3 helicase (NS3h), a highly conserved protein, plays a pivotal role in virus replication and thus represents a potential drug target for flavivirus pathogenesis. NS3h utilizes nucleotide triphosphate, such as ATP, for hydrolysis energy (ATPase) to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. The intermediate states along the ATP binding and hydrolysis cycle, as well as the conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. We use extensive molecular dynamics simulations of apo, ATP, ADP+Pi, and ADP bound to WNV NS3h+ssRNA to model the conformational ensembles along this cycle. Energetic and structural clustering analyses on these trajectories depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). These findings were experimentally corroborated using viral replicons encoding three mutations at the D471 position. Replication assays using these mutants demonstrated a substantial reduction in viral replication compared to the wild-type. Molecular simulations of the D471 mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Predicting allosteric pockets in protein biological assemblages

MOTIVATION

Allostery enables changes to the dynamic behavior of a protein at distant positions induced by binding. Here, we present APOP, a new allosteric pocket prediction method, which perturbs the pockets formed in the structure by stiffening pairwise interactions in the elastic network across the pocket, to emulate ligand binding. Ranking the pockets based on the shifts in the global mode frequencies, as well as their mean local hydrophobicities, leads to high prediction success when tested on a dataset of allosteric proteins, composed of both monomers and multimeric assemblages.

RESULTS

Out of the 104 test cases, APOP predicts known allosteric pockets for 92 within the top 3 rank out of multiple pockets available in the protein. In addition, we demonstrate that APOP can also find new alternative allosteric pockets in proteins. Particularly interesting findings are the discovery of previously overlooked large pockets located in the centers of many protein biological assemblages; binding of ligands at these sites would likely be particularly effective in changing the protein's global dynamics.

AVAILABILITY AND IMPLEMENTATION

APOP is freely available as an open-source code (https://github.com/Ambuj-UF/APOP) and as a web server at https://apop.bb.iastate.edu/.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)8-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)8-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Highly Similar Sequence and Structure Yet Different Biophysical Behavior: A Computational Study of Two Triosephosphate Isomerases

Homodimeric triosephosphate isomerases (TIMs) from Trypanosoma cruzi (TcTIM) and Trypanosoma brucei (TbTIM) have markedly similar amino-acid sequences and three-dimensional structures. However, several of their biophysical parameters, such as their susceptibility to sulfhydryl agents and their reactivation speed after being denatured, have significant differences. The causes of these differences were explored with microsecond-scale molecular dynamics (MD) simulations of three different TIM proteins: TcTIM, TbTIM, and a chimeric protein, Mut1. We examined their electrostatic interactions and explored the impact of simulation length on them. The same salt bridge between catalytic residues Lys 14 and Glu 98 was observed in all three proteins, but key differences were found in other interactions that the catalytic amino acids form. In particular, a cation-π interaction between catalytic amino acids Lys 14 and His 96 and both a salt bridge and a hydrogen bond between catalytic Glu 168 and residue Arg 100 were only observed in TcTIM. Furthermore, although TcTIM forms less hydrogen bonds than TbTIM and Mut1, its hydrogen bond network spans almost the entire protein, connecting the residues in both monomers. This work provides new insight into the mechanisms that give rise to the different behavior of these proteins. The results also show the importance of long simulations.

Catalytically active holo Homo sapiens adenosine deaminase I adopts a closed conformation

Homo sapiens adenosine deaminase 1 (HsADA1; UniProt P00813) is an immunologically relevant enzyme with roles in T-cell activation and modulation of adenosine metabolism and signaling. Patients with genetic deficiency in HsADA1 suffer from severe combined immunodeficiency, and HsADA1 is a therapeutic target in hairy cell leukemias. Historically, insights into the catalytic mechanism and the structural attributes of HsADA1 have been derived from studies of its homologs from Bos taurus (BtADA) and Mus musculus (MmADA). Here, the structure of holo HsADA1 is presented, as well as biochemical characterization that confirms its high activity and shows that it is active across a broad pH range. Structurally, holo HsADA1 adopts a closed conformation distinct from the open conformation of holo BtADA. Comparison of holo HsADA1 and MmADA reveals that MmADA also adopts a closed conformation. These findings challenge previous assumptions gleaned from BtADA regarding the conformation of HsADA1 that may be relevant to its immunological interactions, particularly its ability to bind adenosine receptors. From a broader perspective, the structural analysis of HsADA1 presents a cautionary tale for reliance on homologs to make structural inferences relevant to applications such as protein engineering or drug development.

hdANM: a new comprehensive dynamics model for protein hinges

Hinge motions are essential for many protein functions, and their dynamics are important to understand underlying biological mechanisms. The ways that these motions are represented by various computational methods differ significantly. By focusing on a specific class of motion, we have developed a new hinge-domain anisotropic network model (hdANM) that is based on the prior identification of flexible hinges and rigid domains in the protein structure and the subsequent generation of global hinge motions. This yields a set of motions in which the relative translations and rotations of the rigid domains are modulated and controlled by the deformation of the flexible hinges, leading to a more restricted, specific view of these motions. hdANM is the first model, to our knowledge, that combines information about protein hinges and domains to model the characteristic hinge motions of a protein. The motions predicted with this new elastic network model provide important conceptual advantages for understanding the underlying biological mechanisms. As a matter of fact, the generated hinge movements are found to resemble the expected mechanisms required for the biological functions of diverse proteins. Another advantage of this model is that the domain-level coarse graining makes it significantly more computationally efficient, enabling the generation of hinge motions within even the largest molecular assemblies, such as those from cryo-electron microscopy. hdANM is also comprehensive as it can perform in the same way as the well-known protein dynamics models (anisotropic network model, rotations-translations of blocks, and nonlinear rigid block normal mode analysis), depending on the definition of flexible and rigid parts in the protein structure and on whether the motions are extrapolated in a linear or nonlinear fashion. Furthermore, our results indicate that hdANM produces more realistic motions as compared to the anisotropic network model. hdANM is an open-source software, freely available, and hosted on a user-friendly website.

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

The role of ligand-gated conformational changes in enzyme catalysis

Structural and biochemical studies on diverse enzymes have highlighted the importance of ligand-gated conformational changes in enzyme catalysis, where the intrinsic binding energy of the common phosphoryl group of their substrates is used to drive energetically unfavorable conformational changes in catalytic loops, from inactive open to catalytically competent closed conformations. However, computational studies have historically been unable to capture the activating role of these conformational changes. Here, we discuss recent experimental and computational studies, which can remarkably pinpoint the role of ligand-gated conformational changes in enzyme catalysis, even when not modeling the loop dynamics explicitly. Finally, through our joint analyses of these data, we demonstrate how the synergy between theory and experiment is crucial for furthering our understanding of enzyme catalysis.

Native aggregation is a common feature among triosephosphate isomerases of different species

Triosephosphate isomerase (TIM) is an enzyme of the glycolysis pathway which exists in almost all types of cells. Its structure is the prototype of a motif called TIM-barrel or (α/β)₈ barrel, which is the most common fold of all known enzyme structures. The simplest form in which TIM is catalytically active is a homodimer, in many species of bacteria and eukaryotes, or a homotetramer in some archaea. Here we show that the purified homodimeric TIMs from nine different species of eukaryotes and one of an extremophile bacterium spontaneously form higher order aggregates that can range from 3 to 21 dimers per macromolecular complex. We analysed these aggregates with clear native electrophoresis with normal and inverse polarity, blue native polyacrylamide gel electrophoresis, liquid chromatography, dynamic light scattering, thermal shift assay and transmission electron and fluorescence microscopies, we also performed bioinformatic analysis of the sequences of all enzymes to identify and predict regions that are prone to aggregation. Additionally, the capacity of TIM from Trypanosoma brucei to form fibrillar aggregates was characterized. Our results indicate that all the TIMs we studied are capable of forming oligomers of different sizes. This is significant because aggregation of TIM may be important in some of its non-catalytic moonlighting functions, like being a potent food allergen, or in its role associated with Alzheimer’s disease.

Discovery of Functional Motifs from the Interface Region of Oligomeric Proteins Using Frequent Subgraph Mining

Modeling the interface region of a protein complex paves the way for understanding its dynamics and functionalities. Existing works model the interface region of a complex by using different approaches, such as, the residue composition at the interface region, the geometry of the interface residues, or the structural alignment of interface regions. These approaches are useful for ranking a set of docked conformation or for building scoring function for protein-protein docking, but they do not provide a generic and scalable technique for the extraction of interface patterns leading to functional motif discovery. In this work, we model the interface region of a protein complex by graphs and extract interface patterns of the given complex in the form of frequent subgraphs. To achieve this, we develop a scalable algorithm for frequent subgraph mining. We show that a systematic review of the mined subgraphs provides an effective method for the discovery of functional motifs that exist along the interface region of a given protein complex. In our experiments, we use three PDB protein structure datasets. The first two datasets are composed of PDB structures from different conformations of two dimeric protein complexes: HIV-1 protease (329 structures), and triosephosphate isomerase (TIM) (86 structures). The third dataset is a collection of different enzyme structures protein structures from the six top-level enzyme classes, namely: Oxydoreductase, Transferase, Hydrolase, Lyase, Isomerase, and Ligase. We show that for the first two datasets, our method captures the locking mechanism at the dimeric interface by taking into account the spatial positioning of the interfacial residues through graphs. Indeed, our frequent subgraph mining based approach discovers the patterns representing the dimerization lock which is formed at the base of the structure in 323 of the 329 HIV-1 protease structures. Similarly, for 86 TIM structures, our approach discovers the dimerization lock formation in 50 structures. For the enzyme structures, we show that we are able to capture the functional motifs (active sites) that are specific to each of the six top-level classes of enzymes through frequent subgraphs.

Structural Basis for the Limited Response to Oxidative and Thiol-Conjugating Agents by Triosephosphate Isomerase From the Photosynthetic Bacteria Synechocystis

In plants, the ancestral cyanobacterial triosephosphate isomerase (TPI) was replaced by a duplicated version of the cytosolic TPI. This isoform acquired a transit peptide for chloroplast localization and functions in the Calvin-Benson cycle. To gain insight into the reasons for this gene replacement in plants, we characterized the TPI from the photosynthetic bacteria Synechocystis (SyTPI). SyTPI presents typical TPI enzyme kinetics profiles and assembles as a homodimer composed of two subunits that arrange in a (β-α)₈ fold. We found that oxidizing agents diamide (DA) and H₂O₂, as well as thiol-conjugating agents such as oxidized glutathione (GSSG) and methyl methanethiosulfonate (MMTS), do not inhibit the catalytic activity of SyTPI at concentrations required to inactivate plastidic and cytosolic TPIs from the plant model Arabidopsis thaliana (AtpdTPI and AtcTPI, respectively). The crystal structure of SyTPI revealed that each monomer contains three cysteines, C47, C127, and C176; however only the thiol group of C176 is solvent exposed. While AtcTPI and AtpdTPI are redox-regulated by chemical modifications of their accessible and reactive cysteines, we found that C176 of SyTPI is not sensitive to redox modification in vitro. Our data let us postulate that SyTPI was replaced by a eukaryotic TPI, because the latter contains redox-sensitive cysteines that may be subject to post-translational modifications required for modulating TPI's enzymatic activity.

Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case

Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 "grips" the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple "open and shut" case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.

Protein Conformational Transitions from All-Atom Adaptively Biased Path Optimization

Simulation methods are valuable for elucidating protein conformational transitions between functionally diverse states given that transition pathways are difficult to capture experimentally. Nonetheless, specific computational algorithms are required because of the high free energy barriers between these different protein conformational states. Adaptively biased path optimization (ABPO) is an unrestrained, transition-path optimization method that works in a reduced-variable space to construct an adaptive biasing potential to aid convergence. ABPO was previously applied using a coarse-grained Go̅-model to study conformational activation of Lyn, a Src family tyrosine kinase. How effectively ABPO can be applied at the higher resolution of an all-atom model to explore protein conformational transitions is not yet known. Here, we report the all-atom conformational transition paths of three protein systems constructed using the ABPO methodology. Two systems, triose phosphate isomerase and dihydrofolate reductase, undergo local flipping of a short loop that promotes ligand binding. The third system, estrogen receptor α ligand binding domain, has a helix that adopts different conformations when the protein is bound to an agonist or an antagonist. For each protein, distance-based or torsion-angle reduced variables were identified from unbiased trajectories. ABPO was computed in this reduced variable space to obtain the transition path between the two states. The all-atom ABPO is shown to successfully converge an optimal transition path for each of the three systems.

New strategies for targeting glucose metabolism-mediated acidosis for colorectal cancer therapy

Colorectal cancer (CRC) is a heterogeneous group of diseases that are the result of abnormal glucose metabolism alterations with high lactate production by pyruvate to lactate conversion, which remodels acidosis and offers an evolutional advantage for tumor cells, even enhancing their aggressive phenotype. This review summarizes recent findings that involve multiple genes, molecules, and downstream signaling in the dysregulated glycolytic pathway, which can allow a tumor to initiate acid byproducts and to progress, thereby resulting in acidosis commonly found in the tumor microenvironment of CRC. Moreover, the relationship between CRC cells and the tumor acidic microenvironment, especially for regulating lactate production and lactate dehydrogenase A levels, is also discussed, as well as comprehensively defining different aspects of glycolytic pathways that affect cancer cell proliferation, invasion, and migration. Furthermore, this review concentrates on glucose metabolism-mediated transduction factors in CRC, which include acid-sensing ion channels, triosephosphate isomerase and key glycolysis-related enzymes that regulate glycolytic metabolites, coupled with the effect on tumor cell glycolysis as well as signaling pathways. In conclusion, glucose metabolism mediated by glycolytic pathways that are integral to tumor acidosis in CRC is demonstrated. Therefore, selective metabolic inhibitors or agents against these targets in glucose metabolism through glycolytic pathways may be clinically useful to regulate the tumor's acidic microenvironment for CRC treatment and to identify specific targets that regulate tumor acidosis through a cancer patient-personalized approach. Furthermore, strategies for modifying the metabolic processes that effectively inhibit cancer cell growth and tumor progression and activate potent anticancer effects may provide more effective antitumor prospects for CRC therapy.

Conformational Sampling of the Intrinsically Disordered C-Terminal Tail of DERA Is Important for Enzyme Catalysis

2-Deoxyribose-5-phosphate aldolase (DERA) catalyzes the reversible conversion of acetaldehyde and glyceraldehyde-3-phosphate into deoxyribose-5-phosphate. DERA is used as a biocatalyst for the synthesis of drugs such as statins and is a promising pharmaceutical target due to its involvement in nucleotide catabolism. Despite previous biochemical studies suggesting the catalytic importance of the C-terminal tyrosine residue found in several bacterial DERAs, the structural and functional basis of its participation in catalysis remains elusive because the electron density for the last eight to nine residues (i.e., the C-terminal tail) is absent in all available crystal structures. Using a combination of NMR spectroscopy and molecular dynamics simulations, we conclusively show that the rarely studied C-terminal tail of E. coli DERA (ecDERA) is intrinsically disordered and exists in equilibrium between open and catalytically relevant closed states, where the C-terminal tyrosine (Y259) enters the active site. Nuclear Overhauser effect distance restraints, obtained due to the presence of a substantial closed state population, were used to derive the solution-state structure of the ecDERA closed state. Real-time NMR hydrogen/deuterium exchange experiments reveal that Y259 is required for efficiency of the proton abstraction step of the catalytic reaction. Phosphate titration experiments show that, in addition to the phosphate-binding residues located near the active site, as observed in the available crystal structures, ecDERA contains previously unknown auxiliary phosphate-binding residues on the C-terminal tail which could facilitate in orienting Y259 in an optimal position for catalysis. Thus, we present significant insights into the structural and mechanistic importance of the ecDERA C-terminal tail and illustrate the role of conformational sampling in enzyme catalysis.

Differential effects on enzyme stability and kinetic parameters of mutants related to human triosephosphate isomerase deficiency

Human triosephosphate isomerase (TIM) deficiency is a very rare disease, but there are several mutations reported to be causing the illness. In this work, we produced nine recombinant human triosephosphate isomerases which have the mutations reported to produce TIM deficiency. These enzymes were characterized biophysically and biochemically to determine their kinetic and stability parameters, and also to substitute TIM activity in supporting the growth of an Escherichia coli strain lacking the tim gene. Our results allowed us to rate the deleteriousness of the human TIM mutants based on the type and severity of the alterations observed, to classify four "unknown severity mutants" with altered residues in positions 62, 72, 122 and 154 and to explain in structural terms the mutation V231M, the most affected mutant from the kinetic point of view and the only homozygous mutation reported besides E104D.

Identification of the critical residues responsible for differential reactivation of the triosephosphate isomerases of two trypanosomes

The reactivation of triosephosphate isomerase (TIM) from unfolded monomers induced by guanidine hydrochloride involves different amino acids of its sequence in different stages of protein refolding. We describe a systematic mutagenesis method to find critical residues for certain physico-chemical properties of a protein. The two similar TIMs of Trypanosoma brucei and Trypanosoma cruzi have different reactivation velocities and efficiencies. We used a small number of chimeric enzymes, additive mutants and planned site-directed mutants to produce an enzyme from T. brucei with 13 mutations in its sequence, which reactivates fast and efficiently like wild-type (WT) TIM from T. cruzi, and another enzyme from T. cruzi, with 13 slightly altered mutations, which reactivated slowly and inefficiently like the WT TIM of T. brucei. Our method is a shorter alternative to random mutagenesis, saturation mutagenesis or directed evolution to find multiple amino acids critical for certain properties of proteins.

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BACKGROUND

Accumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.

METHODS/RESULTS

We employ two computationally efficient approaches - graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.

CONCLUSIONS

Contact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.

GENERAL SIGNIFICANCE

This study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.

Altered dynamics upon oligomerization corresponds to key functional sites

It is known that over half of the proteins encoded by most organisms function as oligomeric complexes. Oligomerization confers structural stability and dynamics changes in proteins. We investigate the effects of oligomerization on protein dynamics and its functional significance for a set of 145 multimeric proteins. Using coarse-grained elastic network models, we inspect the changes in residue fluctuations upon oligomerization and then compare with residue conservation scores to identify the functional significance of these changes. Our study reveals conservation of about ½ of the fluctuations, with ¼ of the residues increasing in their mobilities and ¼ having reduced fluctuations. The residues with dampened fluctuations are evolutionarily more conserved and can serve as orthosteric binding sites, indicating their importance. We also use triosephosphate isomerase as a test case to understand why certain enzymes function only in their oligomeric forms despite the monomer including all required catalytic residues. To this end, we compare the residue communities (groups of residues which are highly correlated in their fluctuations) in the monomeric and dimeric forms of the enzyme. We observe significant changes to the dynamical community architecture of the catalytic core of this enzyme. This relates to its functional mechanism and is seen only in the oligomeric form of the protein, answering why proteins are oligomeric structures. Proteins 2017; 85:1422-1434. © 2017 Wiley Periodicals, Inc.

A Rapid Analysis of Variations in Conformational Behavior during Dihydrofolate Reductase Catalysis

Protein flexibility is central to enzyme catalysis, yet it remains challenging both to predict conformational behavior on the basis of analysis of amino acid sequence and protein structure and to provide the necessary breadth of experimental support to any such predictions. Here a generic and rapid procedure for identifying conformational changes during dihydrofolate reductase (DHFR) catalysis is described. Using DHFR from Escherichia coli (EcDHFR), selective side-chain ¹³C labeling of methionine and tryptophan residues is shown to be sufficient to detect the closed-to-occluded conformational transition that follows the chemical step in the catalytic cycle, with clear chemical shift perturbations found for both methionine methyl and tryptophan indole groups. In contrast, no such perturbations are seen for the DHFR from the psychrophile Moritella profunda, where the equivalent conformational change is absent. Like EcDHFR, Salmonella enterica DHFR shows experimental evidence of a large-scale conformational change following hydride transfer that relies on conservation of a key hydrogen bonding interaction between the M20 and GH loops, directly comparable to the closed-to-occluded conformational change observed in EcDHFR. For the hyperthermophile Thermotoga maritima, no chemical shift perturbations were observed, suggesting that no major conformational change occurs during the catalytic cycle. In spite of their conserved tertiary structures, DHFRs display variations in conformational sampling that occurs concurrently with catalysis.

Is the catalytic activity of triosephosphate isomerase fully optimized? An investigation based on maximization of entropy production

Triosephosphate isomerase (TIM) is often described as a fully evolved housekeeping enzyme with near-maximal possible reaction rate. The assumption that an enzyme is perfectly evolved has not been easy to confirm or refute. In this paper, we use maximization of entropy production within known constraints to examine this assumption by calculating steady-state cyclic flux, corresponding entropy production, and catalytic activity in a reversible four-state scheme of TIM functional states. The maximal entropy production (MaxEP) requirement for any of the first three transitions between TIM functional states leads to decreased total entropy production. Only the MaxEP requirement for the product (R-glyceraldehyde-3-phosphate) release step led to a 30% increase in enzyme activity, specificity constant kcat/KM, and overall entropy production. The product release step, due to the TIM molecular machine working in the physiological direction of glycolysis, has not been identified before as the rate-limiting step by using irreversible thermodynamics. Together with structural studies, our results open the possibility for finding amino acid substitutions leading to an increased frequency of loop six opening and product release.

Similarity in Shape Dictates Signature Intrinsic Dynamics Despite No Functional Conservation in TIM Barrel Enzymes

The conservation of the intrinsic dynamics of proteins emerges as we attempt to understand the relationship between sequence, structure and functional conservation. We characterise the conservation of such dynamics in a case where the structure is conserved but function differs greatly. The triosephosphate isomerase barrel fold (TBF), renowned for its 8 β-strand-α-helix repeats that close to form a barrel, is one of the most diverse and abundant folds found in known protein structures. Proteins with this fold have diverse enzymatic functions spanning five of six Enzyme Commission classes, and we have picked five different superfamily candidates for our analysis using elastic network models. We find that the overall shape is a large determinant in the similarity of the intrinsic dynamics, regardless of function. In particular, the β-barrel core is highly rigid, while the α-helices that flank the β-strands have greater relative mobility, allowing for the many possibilities for placement of catalytic residues. We find that these elements correlate with each other via the loops that link them, as opposed to being directly correlated. We are also able to analyse the types of motions encoded by the normal mode vectors of the α-helices. We suggest that the global conservation of the intrinsic dynamics in the TBF contributes greatly to its success as an enzymatic scaffold both through evolution and enzyme design.

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)'s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM's catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand's effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop's dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.

WEBnm@ v2.0: Web server and services for comparing protein flexibility

BACKGROUND

Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics-function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins.

RESULTS

We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma .

CONCLUSION

WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.

Comparing the intrinsic dynamics of multiple protein structures using elastic network models

BACKGROUND

Elastic network models (ENMs) are based on the simple idea that a protein can be described as a set of particles connected by springs, which can then be used to describe its intrinsic flexibility using, for example, normal mode analysis. Since the introduction of the first ENM by Monique Tirion in 1996, several variants using coarser protein models have been proposed and their reliability for the description of protein intrinsic dynamics has been widely demonstrated. Lately an increasing number of studies have focused on the meaning of slow dynamics for protein function and its potential conservation through evolution. This leads naturally to comparisons of the intrinsic dynamics of multiple protein structures with varying levels of similarity.

SCOPE OF REVIEW

We describe computational strategies for calculating and comparing intrinsic dynamics of multiple proteins using elastic network models, as well as a selection of examples from the recent literature.

MAJOR CONCLUSIONS

The increasing interest for comparing dynamics across protein structures with various levels of similarity, has led to the establishment and validation of reliable computational strategies using ENMs. Comparing dynamics has been shown to be a viable way for gaining greater understanding for the mechanisms employed by proteins for their function. Choices of ENM parameters, structure alignment or similarity measures will likely influence the interpretation of the comparative analysis of protein motion.

GENERAL SIGNIFICANCE

Understanding the relation between protein function and dynamics is relevant to the fundamental understanding of protein structure-dynamics-function relationship. This article is part of a Special Issue entitled Recent developments of molecular dynamics.