Structural mechanisms for domain movements in proteins

Comparative dynamics of tropomyosin in vertebrates and invertebrates

Tropomyosin (Tpm) is an extended α-helical coiled-coil homodimer that regulates actinomyosin interactions in muscle. Molecular simulations of four Tpms, two from the vertebrate class Mammalia (rat and pig), and two from the invertebrate class Malacostraca (shrimp and lobster), showed that despite extensive sequence and structural homology across metazoans, dynamic behavior-particularly long-range structural fluctuations-were clearly distinct. Vertebrate Tpms were more flexible and sampled complex, multi-state conformational landscapes. Invertebrate Tpms were more rigid, sampling a highly constrained harmonic landscape. Filtering of trajectories by principle component analysis into essential subspaces showed significant overlap within but not between phyla. In vertebrate Tpms, hinge-regions decoupled long-range interhelical motions and suggested distinct domains. In contrast, crustacean Tpms did not exhibit long-range dynamic correlations-behaving more like a single rigid rod on the nanosecond time scale. These observations suggest there may be divergent mechanisms for Tpm binding to actin filaments, where conformational flexibility in mammalian Tpm allows a preorganized shape complementary to the filament surface, and where rigidity in the crustacean Tpm requires concerted bending and binding.

Structure, interactions and self-assembly of ASC-dependent inflammasomes

The inflammasome is a multi-protein platform that assembles upon the presence of cues derived from infection or tissue damage, and triggers the inflammatory response. Inflammasome components include sensor proteins that detect danger signals, procaspase 1 and the adapter ASC (apoptosis-associated speck-like protein containing a CARD) tethering these molecules together. Upon inflammasome assembly, procaspase 1 self-activates and renders functional cytokines to arbitrate in the defense mechanism. This assembly is mediated by self-association and protein interactions via Death Domains. The inflammasome plays a critical role in innate immunity and its dysregulation is the culprit of many autoimmune disorders. An in-depth understanding of the factors involved in inflammasome assembly could help fight these conditions. This review describes our current knowledge on the biophysical aspects of inflammasome formation from the perspective of ASC. The specific characteristics of the three-dimensional solution structure and interdomain dynamics of ASC are explained in relation to its function in inflammasome assembly. Additionally, the review elaborates on the identification of ASC interacting surfaces at the amino acid level using NMR techniques. Finally, the macrostructures formed by full-length ASC and its two Death Domains studied with Transmission Electron Microscopy are compared in the context of a directional model for inflammasome assembly.

Global response of diacylglycerol kinase towards substrate binding observed by 2D and 3D MAS NMR

Escherichia coli diacylglycerol kinase (DGK) is an integral membrane protein, which catalyses the ATP-dependent phosphorylation of diacylglycerol (DAG) to phosphatic acid (PA). It is a unique trimeric enzyme, which does not share sequence homology with typical kinases. It exhibits a notable complexity in structure and function despite of its small size. Here, chemical shift assignment of wild-type DGK within lipid bilayers was carried out based on 3D MAS NMR, utilizing manual and automatic analysis protocols. Upon nucleotide binding, extensive chemical shift perturbations could be observed. These data provide evidence for a symmetric DGK trimer with all of its three active sites concurrently occupied. Additionally, we could detect that the nucleotide substrate induces a substantial conformational change, most likely directing DGK into its catalytic active form. Furthermore, functionally relevant interprotomer interactions are identified by DNP-enhanced MAS NMR in combination with site-directed mutagenesis and functional assays.

Crystal structure and biochemical properties of msed_0281, the citrate synthase from Metallosphaera sedula

Metallosphaera sedula is a thermoacidophilic archaeon that has carbon fixation ability using the 3-hydroxypropionate/4-hydroxybutyrate(3-HP/4-HB) cycle, and has an incomplete TCA cycle to produce necessary biosynthetic precursors. The citrate synthase from M. sedula (MsCS) is an enzyme involved in the first step of the incomplete TCA cycle, catalyzing the conversion of oxaloacetate and acetyl-CoA into citrate and coenzyme A. To investigate the molecular mechanism of MsCS, we determined its crystal structure at 1.8 Å resolution. As other known CSs, MsCS functions as a dimer, and each monomer consists of two domains, a large domain and a small domain. We also determined the structure of the complex with acetyl-CoA and revealed the acetyl-CoA binding mode of MsCS. Structural comparison of MsCS with another CS in complex with oxaloacetate enabled us to predict the oxaloacetate binding site. Moreover, we performed inhibitory kinetic analyses of MsCS, and showed that the protein is inhibited by citrate and ATP by competitive and non-competitive inhibition modes, respectively, but not by NADH. Based on these results, we suggest that MsCS belongs to the type-I CS with structural and biochemical properties similar to those of CSs involved in the conventional TCA cycle.

Exploring Protein Conformational Diversity

The native state of proteins is composed of conformers in dynamical equilibrium. In this chapter, different issues related to conformational diversity are explored using a curated and experimentally based database called CoDNaS (Conformational Diversity in the Native State). This database is a collection of redundant structures for the same sequence. CoDNaS estimates the degree of conformational diversity using different global and local structural similarity measures. It allows the user to explore how structural differences among conformers change as a function of several structural features providing further biological information. This chapter explores the measurement of conformational diversity and its relationship with sequence divergence. Also, it discusses how proteins with high conformational diversity could affect homology modeling techniques.

Quantifying the Intrinsic Conformation Energy Landscape Topography of Proteins with Large-Scale Open-Closed Transition

Large-scale conformational changes of proteins, including the open-closed transitions, are crucial for a variety of protein functions. These open-closed transitions are often associated with ligand binding. However, the understandings of the underlying mechanisms of the conformational changes within proteins during the open-closed transitions are still challenging at present. In this study, we quantified the intrinsic underlying conformational energy landscapes of five different proteins with large-scale open-closed transitions. This is realized by exploring the underlying density of states and the intrinsic conformational energy landscape topography measure Λ. Λ is a dimensionless ratio of conformational energy gap δE versus conformational energy roughness δE and configurational entropy S or size of the intrinsic conformational energy landscape. By quantifying the Λ of intrinsic open-closed conformational (Λ^oc) and intrinsic global folding (Λ^global) energy landscapes, we show that both intrinsic open-closed conformation energy and entropy landscapes are funneled toward the closed state. Furthermore, our results indicate the strong correlations between Λ and thermodynamics (conformational state transition temperature against trapping temperature) as well as between Λ and kinetics (open-closed kinetic time) of these proteins. This shows that the intrinsic conformational landscape topography determines both the conformational thermodynamic stability and kinetic speed of the conformational dynamics. Our investigations provide important insights for understanding the fundamental mechanisms of the protein conformational dynamics in a physical and global way.

Green function of correlated genes in a minimal mechanical model of protein evolution

The function of proteins arises from cooperative interactions and rearrangements of their amino acids, which exhibit large-scale dynamical modes. Long-range correlations have also been revealed in protein sequences, and this has motivated the search for physical links between the observed genetic and dynamic cooperativity. We outline here a simplified theory of protein, which relates sequence correlations to physical interactions and to the emergence of mechanical function. Our protein is modeled as a strongly coupled amino acid network with interactions and motions that are captured by the mechanical propagator, the Green function. The propagator describes how the gene determines the connectivity of the amino acids and thereby, the transmission of forces. Mutations introduce localized perturbations to the propagator that scatter the force field. The emergence of function is manifested by a topological transition when a band of such perturbations divides the protein into subdomains. We find that epistasis-the interaction among mutations in the gene-is related to the nonlinearity of the Green function, which can be interpreted as a sum over multiple scattering paths. We apply this mechanical framework to simulations of protein evolution and observe long-range epistasis, which facilitates collective functional modes.

Morphing and docking visualisation of biomolecular structures using Multi-Dimensional Scaling

Protein structures are often solved at atomic resolution in two states defining a functional movement but intervening conformations are usually unknown. Morphing methods generate intervening conformations between two known structures. When viewed as an animation using molecular graphics, a smooth, direct morph enables the eye to track changes in structure that might be otherwise missed. We present a morphing method that aims to linearly interpolate interatomic distances and which uses SMACOF (Scaling by MAjorisation of COmplicated Function) and multigrid techniques with a cut-off distance based weighting that optimizes the MolProbity score of intervening structures. The all-atom morphs are smooth, move directly between the two structures, and are shown, in general, to pass closer to a set of known intermediates than those generated using other methods. The techniques are also used for docking by putting the unbound structures in a "near-approach pose" and then morphing to the bound complex. The resulting GPU-accelerated tools are available on a webserver, Morphit_Pro, at http://morphit-pro.cmp.uea.ac.uk/ and more than 5000 domains movements available at the DynDom website can now be viewed as morphs http://morphit-pro.cmp.uea.ac.uk/dyndom/.

A probabilistic network model for structural transitions in biomolecules

Biological macromolecules often undergo large conformational rearrangements during a functional cycle. To simulate these structural transitions with full atomic detail typically demands extensive computational resources. Moreover, it is unclear how to incorporate, in a principled way, additional experimental information that could guide the structural transition. This article develops a probabilistic model for conformational transitions in biomolecules. The model can be viewed as a network of anharmonic springs that break, if the experimental data support the rupture of bonds. Hamiltonian Monte Carlo in internal coordinates is used to infer structural transitions from experimental data, thereby sampling large conformational transitions without distorting the structure. The model is benchmarked on a large set of conformational transitions. Moreover, we demonstrate the use of the probabilistic network model for integrative modeling of macromolecular complexes based on data from crosslinking followed by mass spectrometry.

Structural and Kinetic Studies of Asp632 Mutants and Fully Reduced NADPH-Cytochrome P450 Oxidoreductase Define the Role of Asp632 Loop Dynamics in the Control of NADPH Binding and Hydride Transfer

Conformational changes in NADPH-cytochrome P450 oxidoreductase (CYPOR) associated with electron transfer from NADPH to electron acceptors via FAD and FMN have been investigated via structural studies of the four-electron-reduced NADP⁺-bound enzyme and kinetic and structural studies of mutants that affect the conformation of the mobile Gly631-Asn635 loop (Asp632 loop). The structure of four-electron-reduced, NADP⁺-bound wild type CYPOR shows the plane of the nicotinamide ring positioned perpendicular to the FAD isoalloxazine with its carboxamide group forming H-bonds with N1 of the flavin ring and the Thr535 hydroxyl group. In the reduced enzyme, the C8-C8 atoms of the two flavin rings are ∼1 Å closer than in the fully oxidized and one-electron-reduced structures, which suggests that flavin reduction facilitates interflavin electron transfer. Structural and kinetic studies of mutants Asp632Ala, Asp632Phe, Asp632Asn, and Asp632Glu demonstrate that the carboxyl group of Asp632 is important for stabilizing the Asp632 loop in a retracted position that is required for the binding of the NADPH ribityl-nicotinamide in a hydride-transfer-competent conformation. Structures of the mutants and reduced wild type CYPOR permit us to identify a possible pathway for NADP(H) binding to and release from CYPOR. Asp632 mutants unable to form stable H-bonds with the backbone amides of Arg634, Asn635, and Met636 exhibit decreased catalytic activity and severely impaired hydride transfer from NADPH to FAD, but leave interflavin electron transfer intact. Intriguingly, the Arg634Ala mutation slightly increases the cytochrome P450 2B4 activity. We propose that Asp632 loop movement, in addition to facilitating NADP(H) binding and release, participates in domain movements modulating interflavin electron transfer.

Large scale analysis of protein conformational transitions from aqueous to non-aqueous media

Background

Biocatalysis in organic solvents is nowadays a common practice with a large potential in Biotechnology. Several studies report that proteins which are co-crystallized or soaked in organic solvents preserve their fold integrity showing almost identical arrangements when compared to their aqueous forms. However, it is well established that the catalytic activity of proteins in organic solvents is much lower than in water. In order to explain this diminished activity and to further characterize the behaviour of proteins in non-aqueous environments, we performed a large-scale analysis (1737 proteins) of the conformational diversity of proteins crystallized in aqueous and co-crystallized or soaked in non-aqueous media.

Results

Using proteins’ experimentally determined conformational diversity taken from CoDNaS database, we found that proteins in non-aqueous media display much lower conformational diversity when compared to the corresponding conformers obtained in water. When conformational diversity is compared between conformers obtained in different non-aqueous media, their structural differences are larger and mostly independent of the presence of cognate ligands. We also found that conformers corresponding to non-aqueous media have larger but less flexible cavities, lower number of disordered regions and lower active-site residue mobility.

Conclusions

Our results show that non-aqueous media conformers have specific structural features and that they do not adopt extreme conformations found in aqueous media. This makes them clearly different from their corresponding aqueous conformers.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2044-2) contains supplementary material, which is available to authorized users.

Alternative substrate-bound conformation of bacterial solute-binding protein involved in the import of mammalian host glycosaminoglycans

Glycosaminoglycans (GAGs), constituted by repeating uronate and amino sugar units, are major components of mammalian extracellular matrices. Some indigenous and pathogenic bacteria target GAGs for colonization to and/or infection of host mammalian cells. In Gram-negative pathogenic Streptobacillus moniliformis, the solute-binding protein (Smon0123)-dependent ATP-binding cassette (ABC) transporter incorporates unsaturated GAG disaccharides into the cytoplasm after depolymerization by polysaccharide lyase. Smon0123, composed of N and C domains, adopts either a substrate-free open or a substrate-bound closed form by approaching two domains at 47° in comparison with the open form. Here we show an alternative 39°-closed conformation of Smon0123 bound to unsaturated chondroitin disaccharide sulfated at the C-4 and C-6 positions of N-acetyl-d-galactosamine residue (CΔ4S6S). In CΔ4S6S-bound Smon0123, Arg204 and Lys210 around the two sulfate groups were located at different positions from those at other substrate-bound 47°-closed conformations. Therefore, the two sulfate groups in CΔ4S6S shifted substrate-binding residue arrangements, causing dynamic conformational change. Smon0123 showed less affinity with CΔ4S6S than with non-sulfated and monosulfated substrates. ATPase activity of the Smon0123-dependent ABC transporter in the presence of CΔ4S6S was lower than that in the presence of other unsaturated chondroitin disaccharides, suggesting that CΔ4S6S-bound Smon0123 was unpreferable for docking with the ABC transporter.

Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Carbonic Anhydrase II Apo-Enzyme

Human carbonic anhydrase II (HCA II) is an enzyme that catalyzes the reversible hydration of CO2 into bicarbonate (HCO3-) and a proton (H+) as well as other reactions at an extremely high rate. This enzyme plays fundamental roles in human physiology/pathology, such as controlling the pH level in cells and so on. However, the binding mechanism between apo-HCA II and CO2 or other ligands as well as related conformational changes remains poorly understood, and atomic investigation into it could promote our understanding of related internal physiological/pathological mechanisms. In this study, long-time atomic molecular dynamics simulations as well as the clustering and free-energy analysis were performed to reveal the dynamics of apo-HCA II as well as the mechanism upon ligand binding. Our simulations indicate that the crystallographic B-factors considerably underestimate the loop dynamics: multiple conformations can be adopted by loops 1 and 2, especially for loop 1 because loop 1 is one side of the binding pocket, and its left-to-right movement can compress or extend the binding pocket, leading to one inactive (closed) state, three intermediate (semiopen) states, and one active (open) state; CO2 cannot get into the binding pocket of the inactive state but can get into those of intermediate and active states. The coexistence of multiple conformational states proposes a possible conformational selection model for the binding mechanism between apo-HCA II and CO2 or other ligands, revising our previous view of its functional mechanism of conformational change upon ligand binding and offering valuable structural insights into the workings of HCA II.

Structural studies of domain movement in active-site mutants of porphobilinogen deaminase from Bacillus megaterium

The enzyme porphobilinogen deaminase (PBGD) is one of the key enzymes in tetrapyrrole biosynthesis. It catalyses the formation of a linear tetrapyrrole from four molecules of the substrate porphobilinogen (PBG). It has a dipyrromethane cofactor (DPM) in the active site which is covalently linked to a conserved cysteine residue through a thioether bridge. The substrate molecules are linked to the cofactor in a stepwise head-to-tail manner during the reaction, which is catalysed by a conserved aspartate residue: Asp82 in the B. megaterium enzyme. Three mutations have been made affecting Asp82 (D82A, D82E and D82N) and their crystal structures have been determined at resolutions of 2.7, 1.8 and 1.9 Å, respectively. These structures reveal that whilst the D82E mutant possesses the DPM cofactor, in the D82N and D82A mutants the cofactor is likely to be missing, incompletely assembled or disordered. Comparison of the mutant PBGD structures with that of the wild-type enzyme shows that there are significant domain movements and suggests that the enzyme adopts `open' and `closed' conformations, potentially in response to substrate binding.

Homology modeling in a dynamical world

A key concept in template-based modeling (TBM) is the high correlation between sequence and structural divergence, with the practical consequence that homologous proteins that are similar at the sequence level will also be similar at the structural level. However, conformational diversity of the native state will reduce the correlation between structural and sequence divergence, because structural variation can appear without sequence diversity. In this work, we explore the impact that conformational diversity has on the relationship between structural and sequence divergence. We find that the extent of conformational diversity can be as high as the maximum structural divergence among families. Also, as expected, conformational diversity impairs the well-established correlation between sequence and structural divergence, which is nosier than previously suggested. However, we found that this noise can be resolved using a priori information coming from the structure-function relationship. We show that protein families with low conformational diversity show a well-correlated relationship between sequence and structural divergence, which is severely reduced in proteins with larger conformational diversity. This lack of correlation could impair TBM results in highly dynamical proteins. Finally, we also find that the presence of order/disorder can provide useful beforehand information for better TBM performance.

A probabilistic model for detecting rigid domains in protein structures

MOTIVATION

Large-scale conformational changes in proteins are implicated in many important biological functions. These structural transitions can often be rationalized in terms of relative movements of rigid domains. There is a need for objective and automated methods that identify rigid domains in sets of protein structures showing alternative conformational states.

RESULTS

We present a probabilistic model for detecting rigid-body movements in protein structures. Our model aims to approximate alternative conformational states by a few structural parts that are rigidly transformed under the action of a rotation and a translation. By using Bayesian inference and Markov chain Monte Carlo sampling, we estimate all parameters of the model, including a segmentation of the protein into rigid domains, the structures of the domains themselves, and the rigid transformations that generate the observed structures. We find that our Gibbs sampling algorithm can also estimate the optimal number of rigid domains with high efficiency and accuracy. We assess the power of our method on several thousand entries of the DynDom database and discuss applications to various complex biomolecular systems.

AVAILABILITY AND IMPLEMENTATION

The Python source code for protein ensemble analysis is available at: https://github.com/thachnguyen/motion_detection

CONTACT

: mhabeck@gwdg.de.

A Network of Conformational Transitions in the Apo Form of NDM-1 Enzyme Revealed by MD Simulation and a Markov State Model

New Delhi metallo-β-lactamase-1 (NDM-1) is a novel β-lactamase enzyme that confers enteric bacteria with nearly complete resistance to all β-lactam antibiotics, so it raises a formidable and global threat to human health. However, the binding mechanism between apo-NDM-1 and antibiotics as well as related conformational changes remains poorly understood, which largely hinders the overcoming of its antibiotic resistance. In our study, long-time conventional molecular dynamics simulation and Markov state models were applied to reveal both the dynamical and conformational landscape of apo-NDM-1: the MD simulation demonstrates that loop L3, which is responsible for antibiotic binding, is the most flexible and undergoes dramatic conformational changes; moreover, the Markov state model built from the simulation maps four metastable states including open, semiopen, and closed conformations of loop L3 as well as frequent transitions between the states. Our findings propose a possible conformational selection model for the binding mechanism between apo-NDM-1 and antibiotics, which facilitates the design of novel inhibitors and antibiotics.

Bayesian Modeling of Biomolecular Assemblies with Cryo-EM Maps

A growing array of experimental techniques allows us to characterize the three-dimensional structure of large biological assemblies at increasingly higher resolution. In addition to X-ray crystallography and nuclear magnetic resonance in solution, new structure determination methods such cryo-electron microscopy (cryo-EM), crosslinking/mass spectrometry and solid-state NMR have emerged. Often it is not sufficient to use a single experimental method, but complementary data need to be collected by using multiple techniques. The integration of all datasets can only be achieved by computational means. This article describes Inferential structure determination, a Bayesian approach to integrative modeling of biomolecular complexes with hybrid structural data. I will introduce probabilistic models for cryo-EM maps and outline Markov chain Monte Carlo algorithms for sampling model structures from the posterior distribution. I will focus on rigid and flexible modeling with cryo-EM data and discuss some of the computational challenges of Bayesian inference in the context of biomolecular modeling.

Architecture and coevolution of allosteric materials

We introduce a numerical scheme to evolve functional elastic materials that can accomplish a specified mechanical task. In this scheme, the number of solutions, their spatial architectures, and the correlations among them can be computed. As an example, we consider an "allosteric" task, which requires the material to respond specifically to a stimulus at a distant active site. We find that functioning materials evolve a less-constrained trumpet-shaped region connecting the stimulus and active sites, and that the amplitude of the elastic response varies nonmonotonically along the trumpet. As previously shown for some proteins, we find that correlations appearing during evolution alone are sufficient to identify key aspects of this design. Finally, we show that the success of this architecture stems from the emergence of soft edge modes recently found to appear near the surface of marginally connected materials. Overall, our in silico evolution experiment offers a window to study the relationship between structure, function, and correlations emerging during evolution.

Conformational diversity analysis reveals three functional mechanisms in proteins

Protein motions are a key feature to understand biological function. Recently, a large-scale analysis of protein conformational diversity showed a positively skewed distribution with a peak at 0.5 Å C-alpha root-mean-square-deviation (RMSD). To understand this distribution in terms of structure-function relationships, we studied a well curated and large dataset of ~5,000 proteins with experimentally determined conformational diversity. We searched for global behaviour patterns studying how structure-based features change among the available conformer population for each protein. This procedure allowed us to describe the RMSD distribution in terms of three main protein classes sharing given properties. The largest of these protein subsets (~60%), which we call "rigid" (average RMSD = 0.83 Å), has no disordered regions, shows low conformational diversity, the largest tunnels and smaller and buried cavities. The two additional subsets contain disordered regions, but with differential sequence composition and behaviour. Partially disordered proteins have on average 67% of their conformers with disordered regions, average RMSD = 1.1 Å, the highest number of hinges and the longest disordered regions. In contrast, malleable proteins have on average only 25% of disordered conformers and average RMSD = 1.3 Å, flexible cavities affected in size by the presence of disordered regions and show the highest diversity of cognate ligands. Proteins in each set are mostly non-homologous to each other, share no given fold class, nor functional similarity but do share features derived from their conformer population. These shared features could represent conformational mechanisms related with biological functions.

Fingerprints of Conformational States of Human Hsp70 at Sub-THz Frequencies

Large multidomain proteins occur in different conformational states to function. Detection and monitoring of these different structural states are of crucial interest for understanding the mechanics of proteins. Using computational methods, we show that different protein conformational states of the two-domain 70 kDa human Heat-shock protein (hHsp70), with similar vibrational density of states, lead to remarkably different far-IR spectra at acoustical frequencies (ν < 300 GHz). We found that the slow damped motions of the positively charged residues of hHsp70 contribute the most to collective IR active modes at low frequencies (ν < 300 GHz). We predicted that different structural states and functional modes of large proteins, such as hHsp70, might be detected in the sub-THz frequency range by single-molecule spectroscopy similar to the recent extraordinary acoustic Raman spectroscopy (Wheaton S.; Nat. Photonics2015, 9, 68-72).

An efficient perturbation method to predict the functionally key sites of glutamine binding protein

Glutamine-Binding Protein (GlnBP) of Escherichia coli, an important member of the periplasmic binding protein family, is responsible for the first step in the active transport of glutamine across the cytoplasmic membrane. In this work, the functionally key regulation sites of GlnBP were identified by utilizing a perturbation method proposed by our group, in which the residues whose perturbations markedly change the binding free energy between GlnBP and glutamine are considered to be functionally key residues. The results show that besides the substrate binding sites, some other residues distant from the binding pocket, including the ones in the hinge regions between the two domains, the front- and back- door channels and the exposed region, are important for the function of glutamine binding and transport. The predicted results are well consistent with the theoretical and experimental data, which indicates that our method is an effective approach to identify the key residues important for both ligand binding and long-range allosteric signal transmission. This work can provide some insights into the function performance of GlnBP and the physical mechanism of its allosteric regulation.

Theoretical Insights into Sub-Terahertz Acoustic Vibrations of Proteins Measured in Single-Molecule Experiments

Proteins are an important class of nanobioparticles with acoustical modes in the sub-THz frequency range. There is considerable interest to measure and establish the role of these acoustical vibrations for biological function. So far, the technique providing the most detailed information about the acoustical modes of proteins is the very recent Extraordinary Acoustic Raman (EAR) spectroscopy. In this technique, proteins are trapped in nanoholes and excited by two optical lasers of slightly different wavelengths producing an electric field at low frequency (<100 GHz). We demonstrate that the acoustical modes of proteins studied by EAR spectroscopy are both infrared- and Raman-active modes, and we provided interpretation of the spectroscopic fingerprints measured at the single-molecule level. A combination of the present calculations with techniques based on the excitation of a single nanobioparticle by an electric field, such as EAR spectroscopy, should provide a wealth of information on the role of molecular dynamics for biological function.

A Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Binary Complex of Escherichia coli 6-Hydroxymethyl-7,8-dihydropterin Pyrophosphokinase with MgATP

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that, in its binary form, even bound by analogues of MgATP, loops 2 and 3 remain rather flexible; this raises questions about the putative large-scale induced-fit conformational change of the HPPK-MgATP binary complex. In this work, long-time all-atomic molecular dynamics simulations were conducted to investigate the loop dynamics in this complex. Our simulations show that, with loop 3 closed, multiple conformations of loop 2, including the open, semiopen, and closed forms, are all accessible to the binary complex. These results provide valuable structural insights into the details of conformational changes upon 6-hydroxymethyl-7,8-dihydropterin (HP) binding and biological activities of HPPK. Conformational network analysis and principal component analysis related to the loops are also discussed.

Single molecule optical measurements of orientation and rotations of biological macromolecules

Subdomains of macromolecules often undergo large orientation changes during their catalytic cycles that are essential for their activity. Tracking these rearrangements in real time opens a powerful window into the link between protein structure and functional output. Site-specific labeling of individual molecules with polarized optical probes and measurement of their spatial orientation can give insight into the crucial conformational changes, dynamics, and fluctuations of macromolecules. Here we describe the range of single molecule optical technologies that can extract orientation information from these probes, review the relevant types of probes and labeling techniques, and highlight the advantages and disadvantages of these technologies for addressing specific inquiries.

CoDNaS 2.0: a comprehensive database of protein conformational diversity in the native state

CoDNaS (conformational diversity of the native state) is a protein conformational diversity database. Conformational diversity describes structural differences between conformers that define the native state of proteins. It is a key concept to understand protein function and biological processes related to protein functions. CoDNaS offers a well curated database that is experimentally driven, thoroughly linked, and annotated. CoDNaS facilitates the extraction of key information on small structural differences based on protein movements. CoDNaS enables users to easily relate the degree of conformational diversity with physical, chemical and biological properties derived from experiments on protein structure and biological characteristics. The new version of CoDNaS includes ∼70% of all available protein structures, and new tools have been added that run sequence searches, display structural flexibility profiles and allow users to browse the database for different structural classes. These tools facilitate the exploration of protein conformational diversity and its role in protein function.

Database URL:http://ufq.unq.edu.ar/codnas

Evolutionary Conserved Positions Define Protein Conformational Diversity

Conformational diversity of the native state plays a central role in modulating protein function. The selection paradigm sustains that different ligands shift the conformational equilibrium through their binding to highest-affinity conformers. Intramolecular vibrational dynamics associated to each conformation should guarantee conformational transitions, which due to its importance, could possibly be associated with evolutionary conserved traits. Normal mode analysis, based on a coarse-grained model of the protein, can provide the required information to explore these features. Herein, we present a novel procedure to identify key positions sustaining the conformational diversity associated to ligand binding. The method is applied to an adequate refined dataset of 188 paired protein structures in their bound and unbound forms. Firstly, normal modes most involved in the conformational change are selected according to their corresponding overlap with structural distortions introduced by ligand binding. The subspace defined by these modes is used to analyze the effect of simulated point mutations on preserving the conformational diversity of the protein. We find a negative correlation between the effects of mutations on these normal mode subspaces associated to ligand-binding and position-specific evolutionary conservations obtained from multiple sequence-structure alignments. Positions whose mutations are found to alter the most these subspaces are defined as key positions, that is, dynamically important residues that mediate the ligand-binding conformational change. These positions are shown to be evolutionary conserved, mostly buried aliphatic residues localized in regular structural regions of the protein like β-sheets and α-helix.

E. coli metabolic protein aldehyde-alcohol dehydrogenase-E binds to the ribosome: a unique moonlighting action revealed

It is becoming increasingly evident that a high degree of regulation is involved in the protein synthesis machinery entailing more interacting regulatory factors. A multitude of proteins have been identified recently which show regulatory function upon binding to the ribosome. Here, we identify tight association of a metabolic protein aldehyde-alcohol dehydrogenase E (AdhE) with the E. coli 70S ribosome isolated from cell extract under low salt wash conditions. Cryo-EM reconstruction of the ribosome sample allows us to localize its position on the head of the small subunit, near the mRNA entrance. Our study demonstrates substantial RNA unwinding activity of AdhE which can account for the ability of ribosome to translate through downstream of at least certain mRNA helices. Thus far, in E. coli, no ribosome-associated factor has been identified that shows downstream mRNA helicase activity. Additionally, the cryo-EM map reveals interaction of another extracellular protein, outer membrane protein C (OmpC), with the ribosome at the peripheral solvent side of the 50S subunit. Our result also provides important insight into plausible functional role of OmpC upon ribosome binding. Visualization of the ribosome purified directly from the cell lysate unveils for the first time interactions of additional regulatory proteins with the ribosome.

Clusters of Structurally Similar MHC I HLA-A2 Molecules, Found with a New Method, Suggest Mechanisms of T-Cell Receptor Avidity

Only α1 and α2 domains of the α-chain of the major histocompatibility complex (MHC) directly bind peptide antigens (Ag-s) and the T-cell receptor (TCR). Significant plasticity was found in the TCR but only minor in (α1 + α2). The α3-domain position variation was noted only in connection to its binding the coreceptor CD8. We apply our methods for identifying functional conformational changes in proteins to a systematic study of similarities between 43 X-ray structures of the entire α chains of MHC-I HLA-A2. Out of 903 different αHLA-A2 pairs 203 show similarities within the earlier determined uncertainty threshold and unexpectedly form three similarity clusters (SCs) with all/most structures in a cluster similar within the uncertainty threshold. Pairs from different SCs always differ above the threshold, mainly due to variations in the α3 position/structure. All structures in SC3 cannot bind the CD8 coreceptor. Strong hydrogen bonds between (α1 + α2) and α3 differ between SC1 and SC2 but are nearly invariant within each SC. Small conformational changes in the (α1 + α2), caused by Ag-s differences, act as an α3 "allosteric switch" between SC2 and SC1. Binding of CD8 to SC2-HLA-A2 (Tax-type Ag-s) changes it to SC1-HLA-A2 (HuD-type Ag-s). HuD binding to HLA-A2 is much less stable than Tax binding. CD8-liganded HLA-A2 preference for binding HuD suggests that CD8-HLA-A2 may present a weakly binding peptide for TCR recognition, supporting the hypothesis that CD8 increases TCR avidity to weak Ag-s. Other HLA-A2 functions may involve α3. TCR-A6-liganded-Tax-type-HLA-A2s form two small clusters, similar to either A6-liganded-HuD or A6-liganded-native-Tax HLA-A2s.

Structural and mechanistic insight into substrate binding from the conformational dynamics in apo and substrate-bound DapE enzyme

Conformational dynamics induced by substrate binding in DapE enzyme.

Intrinsic Localized Modes in Proteins

Protein dynamics is essential for proteins to function. Here we predicted the existence of rare, large nonlinear excitations, termed intrinsic localized modes (ILMs), of the main chain of proteins based on all-atom molecular dynamics simulations of two fast-folder proteins and of a rigid α/β protein at 300 K and at 380 K in solution. These nonlinear excitations arise from the anharmonicity of the protein dynamics. The ILMs were detected by computing the Shannon entropy of the protein main-chain fluctuations. In the non-native state (significantly explored at 380 K), the probability of their excitation was increased by a factor between 9 and 28 for the fast-folder proteins and by a factor 2 for the rigid protein. This enhancement in the non-native state was due to glycine, as demonstrated by simulations in which glycine was mutated to alanine. These ILMs might play a functional role in the flexible regions of proteins and in proteins in a non-native state (i.e. misfolded or unfolded states).

Domain Motion Enhanced (DoME) Model for Efficient Conformational Sampling of Multidomain Proteins

Large conformational changes of multidomain proteins are difficult to simulate using all-atom molecular dynamics (MD) due to the slow time scale. We show that a simple modification of the structure-based coarse-grained (CG) model enables a stable and efficient MD simulation of those proteins. "Motion Tree", a tree diagram that describes conformational changes between two structures in a protein, provides information on rigid structural units (domains) and the magnitudes of domain motions. In our new CG model, which we call the DoME (domain motion enhanced) model, interdomain interactions are defined as being inversely proportional to the magnitude of the domain motions in the diagram, whereas intradomain interactions are kept constant. We applied the DoME model in combination with the Go model to simulations of adenylate kinase (AdK). The results of the DoME-Go simulation are consistent with an all-atom MD simulation for 10 μs as well as known experimental data. Unlike the conventional Go model, the DoME-Go model yields stable simulation trajectories against temperature changes and conformational transitions are easily sampled despite domain rigidity. Evidently, identification of domains and their interfaces is useful approach for CG modeling of multidomain proteins.

Molecular dynamics simulations of the Escherichia coli HPPK apo-enzyme reveal a network of conformational transitions

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that even when confined by crystal contacts, loops 2 and 3 remain rather flexible when the enzyme is in its apo form, raising questions about the putative large-scale induced-fit conformational change of HPPK. To investigate the loop dynamics in a crystal-free environment, we performed conventional molecular dynamics simulations of the apo-enzyme at two different temperatures (300 and 350 K). Our simulations show that the crystallographic B-factors considerably underestimate the loop dynamics; multiple conformations of loops 2 and 3, including the open, semi-open, and closed conformations that an enzyme must adopt throughout its catalytic cycle, are all accessible to the apo-enzyme. These results revise our previous view of the functional mechanism of conformational change upon MgATP binding and offer valuable structural insights into the workings of HPPK. In this paper, conformational network analysis and principal component analysis related to the loops are discussed to support the presented conclusions.

Structural comparison between the open and closed forms of citrate synthase from Thermus thermophilus HB8

The crystal structures of citrate synthase from the thermophilic eubacteria Thermus thermophilus HB8 (TtCS) were determined for an open form at 1.5 Å resolution and for closed form at 2.3 Å resolution, respectively. In the absence of ligands TtCS in the open form was crystalized into a tetragonal form with a single subunit in the asymmetric unit. TtCS was also co-crystallized with citrate and coenzyme-A to form an orthorhombic crystal with two homodimers in the asymmetric unit. Citrate and CoA are found in the active site situated between the large domain and the small domain in all subunit whereas the complex shows two distinct closed conformations, the fully closed form and partially closed form. Structural comparisons are performed to describe conformational changes associated with binding of products of TtCS. Upon binding of citrate, basic residues in the active site move toward citrate and make a hydrogen bond network in the active site, inducing a large-scale rotation of the small domain relative to the large domain. CoA is sandwiched between the small and large domains and then the cysteamine tail is inserted into the active site with a cooperative rotation around mainchain dihedrals in the hinge region connecting helices M and N. According to this rotation these helices are extended to close the active site completely. The considerable flexibility and structural rearrangements in the hinge region are crucial for an ordered bibi reaction in catalysis for microbial CSs.

Protein dynamics and the all-ferrous [Fe4 S4 ] cluster in the nitrogenase iron protein

In nitrogen fixation by Azotobacter vinelandii nitrogenase, the iron protein (FeP) binds to and subsequently transfers electrons to the molybdenum–FeP, which contains the nitrogen fixation site, along with hydrolysis of two ATPs. However, the nature of the reduced state cluster is not completely clear. While reduced FeP is generally thought to contain an [Fe₄S₄]¹⁺ cluster, evidence also exists for an all‐ferrous [Fe₄S₄]⁰ cluster. Since the former indicates a single electron is transferred per two ATPs hydrolyzed while the latter indicates two electrons could be transferred per two ATPs hydrolyzed, an all‐ferrous [Fe₄S₄]⁰ cluster in FeP is potenially two times more efficient. However, the 1+/0 reduction potential has been measured in the protein at both 460 and 790 mV, causing the biological significance to be questioned. Here, “density functional theory plus Poisson Boltzmann” calculations show that cluster movement relative to the protein surface observed in the crystal structures could account for both measured values. In addition, elastic network mode analysis indicates that such movement occurs in low frequency vibrations of the protein, implying protein dynamics might lead to variations in reduction potential. Furthermore, the different reductants used in the conflicting measurements of the reduction potential could be differentially affecting the protein dynamics. Moreover, even if the all‐ferrous cluster is not the biologically relevant cluster, mutagenesis to stabilize the conformation with the more exposed cluster may be useful for bioengineering more efficient enzymes.

Structure and Function in Homodimeric Enzymes: Simulations of Cooperative and Independent Functional Motions

Large-scale conformational change is a common feature in the catalytic cycles of enzymes. Many enzymes function as homodimers with active sites that contain elements from both chains. Symmetric and anti-symmetric cooperative motions in homodimers can potentially lead to correlated active site opening and/or closure, likely to be important for ligand binding and release. Here, we examine such motions in two different domain-swapped homodimeric enzymes: the DcpS scavenger decapping enzyme and citrate synthase. We use and compare two types of all-atom simulations: conventional molecular dynamics simulations to identify physically meaningful conformational ensembles, and rapid geometric simulations of flexible motion, biased along normal mode directions, to identify relevant motions encoded in the protein structure. The results indicate that the opening/closure motions are intrinsic features of both unliganded enzymes. In DcpS, conformational change is dominated by an anti-symmetric cooperative motion, causing one active site to close as the other opens; however a symmetric motion is also significant. In CS, we identify that both symmetric (suggested by crystallography) and asymmetric motions are features of the protein structure, and as a result the behaviour in solution is largely non-cooperative. The agreement between two modelling approaches using very different levels of theory indicates that the behaviours are indeed intrinsic to the protein structures. Geometric simulations correctly identify and explore large amplitudes of motion, while molecular dynamics simulations indicate the ranges of motion that are energetically feasible. Together, the simulation approaches are able to reveal unexpected functionally relevant motions, and highlight differences between enzymes.

Parsimony in Protein Conformational Change

Protein conformational change is analyzed by finding the minimalist backbone torsion angle rotations that superpose crystal structures within experimental error. Of several approaches for enforcing parsimony during flexible least-squares superposition, an ℓ(1)-norm restraint provided greatest consistency with independent indications of flexibility from nuclear magnetic resonance relaxation dispersion and chemical shift perturbation in arginine kinase and four previously studied systems. Crystallographic cross-validation shows that the dihedral parameterization describes conformational change more accurately than rigid-group approaches. The rotations that superpose the principal elements of structure constitute a small fraction of the raw (φ, ψ) differences that also reflect local conformation and experimental error. Substantial long-range displacements can be mediated by modest dihedral rotations, accommodated even within α helices and β sheets without disruption of hydrogen bonding at the hinges. Consistency between ligand-associated and intrinsic motions (in the unliganded state) implies that induced changes tend to follow low-barrier paths between conformational sub-states that are in intrinsic dynamic equilibrium.