Toolbox for analyzing finite two-state trajectories

In many experiments, the aim is to deduce an underlying multi-substate on-off kinetic scheme (KS) from the statistical properties of a two-state trajectory. However, the mapping of a KS into a two-state trajectory leads to the loss of information about the KS, and so, in many cases, more than one KS can be associated with the data. We recently showed that the optimal way to solve this problem is to use canonical forms of reduced dimensions (RD). RD forms are on-off networks with connections only between substates of different states, where the connections can have non-exponential waiting time probability density functions (WT-PDFs). In theory, only a single RD form can be associated with the data. To utilize RD forms in the analysis of the data, a RD form should be associated with the data. Here, we give a toolbox for building a RD form from a finite two-state trajectory. The methods in the toolbox are based on known statistical methods in data analysis, combined with statistical methods and numerical algorithms designed specifically for the current problem. Our toolbox is self-contained - it builds a mechanism based only on the information it extracts from the data, and its implementation on the data is fast (analyzing a 10^6 cycle trajectory from a thirty-parameter mechanism takes a couple of hours on a PC with a 2.66 GHz processor). The toolbox is automated and is freely available for academic research upon electronic request

Regulatory mechanisms triggered by enzyme interactions with lipid membrane surfaces

Recruitment of enzymes to intracellular membranes often modulates their catalytic activity, which can be important in cell signaling and membrane trafficking. Thus, re-localization is not only important for these enzymes to gain access to their substrates, but membrane interactions often allosterically regulate enzyme function by inducing conformational changes across different time and amplitude scales. Recent structural, biophysical and computational studies have revealed how key enzymes interact with lipid membrane surfaces, and how this membrane binding regulates protein structure and function. This review summarizes the recent progress in understanding regulatory mechanisms involved in enzyme-membrane interactions

Millisecond Timescale Motions Connect Amino Acid Interaction Networks in Alpha Tryptophan Synthase

Tryptophan synthase is a model system for understanding allosteric regulation within enzyme complexes. Amino acid interaction networks were previously delineated in the isolated alpha subunit (αTS) in the absence of the beta subunit (βTS). The amino acid interaction networks were different between the ligand-free enzyme and the enzyme actively catalyzing turnover. Previous X-ray crystallography studies indicated only minor localized changes when ligands bind αTS, and so, structural changes alone could not explain the changes to the amino acid interaction networks. We hypothesized that the network changes could instead be related to changes in conformational dynamics. As such, we conducted nuclear magnetic resonance relaxation studies on different substrate- and products-bound complexes of αTS. Specifically, we collected 15N R2 relaxation dispersion data that reports on microsecond-to-millisecond timescale motion of backbone amide groups. These experiments indicated that there are conformational exchange events throughout αTS. Substrate and product binding change specific motional pathways throughout the enzyme, and these pathways connect the previously identified network residues. These pathways reach the αTS/βTS binding interface, suggesting that the identified dynamic networks may also be important for communication with the βTS subunit

In Silico Elucidation of the Recognition Dynamics of Ubiquitin

Elucidation of the mechanism of biomacromolecular recognition events has been a topic of intense interest over the past century. The inherent dynamic nature of both protein and ligand molecules along with the continuous reshaping of the energy landscape during the binding process renders it difficult to characterize this process at atomic detail. Here, we investigate the recognition dynamics of ubiquitin via microsecond all-atom molecular dynamics simulation providing both thermodynamic and kinetic information. The high-level of consistency found with respect to experimental NMR data lends support to the accuracy of the in silico representation of the conformational substates and their interconversions of free ubiquitin. Using an energy-based reweighting approach, the statistical distribution of conformational states of ubiquitin is monitored as a function of the distance between ubiquitin and its binding partner Hrs-UIM. It is found that extensive and dense sampling of conformational space afforded by the µs MD trajectory is essential for the elucidation of the binding mechanism as is Boltzmann sampling, overcoming inherent limitations of sparsely sampled empirical ensembles. The results reveal a population redistribution mechanism that takes effect when the ligand is at intermediate range of 1–2 nm from ubiquitin. This mechanism, which may be depicted as a superposition of the conformational selection and induced fit mechanisms, also applies to other binding partners of ubiquitin, such as the GGA3 GAT domain

Vaccine-derived mutation in motif D of poliovirus RNA-dependent RNA polymerase lowers nucleotide incorporation fidelity

All viral RNA-dependent RNA polymerases (RdRps) have a conserved structural element termed motif D. Studies of the RdRp from poliovirus (PV) have shown that a conformational change of motifDleads to efficient and faithful nucleotide addition by bringing Lys-359 into the active site where it serves as a general acid. The RdRp of the Sabin I vaccine strain has Thr-362 changed to Ile. Such a drastic change so close to Lys-359 might alter RdRp function and contribute in some way to the attenuated phenotype of Sabin type I. Here we present our characterization of the T362I RdRp.Wefind that the T362I RdRp exhibits a mutator phenotype in biochemical experiments in vitro. Using NMR, we show that this change in nucleotide incorporation fidelity correlates with a change in the structural dynamics of motif D. A recombinant PV expressing the T362I RdRp exhibits normal growth properties in cell culture but expresses a mutator phenotype in cells. For example, the T362I-containing PV is more sensitive to the mutagenic activity of ribavirin than wildtype PV. Interestingly, the T362I change was sufficient to cause a statistically significant reduction in viral virulence. Collectively, these studies suggest that residues of motif D can be targeted when changes in nucleotide incorporation fidelity are desired. Given the observation that fidelity mutants can serve as vaccine candidates, it may be possible to use engineering of motif D for this purpose

Manipulating Protein Conformations By Single-molecule Afm-fret Nanoscopy

Combining atomic force microscopy and fluorescence resonance energy transfer spectroscopy (AFM-FRET), we have developed a single-molecule AFM-FRET nanoscopy approach capable of effectively pinpointing and mechanically manipulating a targeted dye-labeled single protein in a large sampling area and simultaneously monitoring the conformational changes of the targeted protein by recording single-molecule FRET time trajectories. We have further demonstrated an application of using this nanoscopy on manipulation of single-molecule protein conformation and simultaneous single-molecule FRET measurement of a Cy3-Cy5-labeled kinase enzyme, HPPK (6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase). By analyzing time-resolved FRET trajectories and correlated AFM force pulling curves of the targeted single-molecule enzyme, we are able to observe the protein conformational changes of a specific coordination by AFM mechanic force pulling

Determination of Conformational Equilibria in Proteins Using Residual Dipolar Couplings

In order to carry out their functions, proteins often undergo significant conformational fluctuations that enable them to interact with their partners. The accurate characterization of these motions is key in order to understand the mechanisms by which macromolecular recognition events take place. Nuclear magnetic resonance spectroscopy offers a variety of powerful methods to achieve this result. We discuss a method of using residual dipolar couplings as replica-averaged restraints in molecular dynamics simulations to determine large amplitude motions of proteins, including those involved in the conformational equilibria that are established through interconversions between different states. By applying this method to ribonuclease A, we show that it enables one to characterize the ample fluctuations in interdomain orientations expected to play an important functional role

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods. © 2013 American Chemical Society