Empirical Bayes methods enable advanced population-level analyses of single-molecule FRET experiments.

Many single-molecule experiments aim to characterize biomolecular processes in terms of kinetic models that specify the rates of transition between conformational states of the biomolecule. Estimation of these rates often requires analysis of a population of molecules, in which the conformational trajectory of each molecule is represented by a noisy, time-dependent signal trajectory. Although hidden Markov models (HMMs) may be used to infer the conformational trajectories of individual molecules, estimating a consensus kinetic model from the population of inferred conformational trajectories remains a statistically difficult task, as inferred parameters vary widely within a population. Here, we demonstrate how a recently developed empirical Bayesian method for HMMs can be extended to enable a more automated and statistically principled approach to two widely occurring tasks in the analysis of single-molecule fluorescence resonance energy transfer (smFRET) experiments: 1), the characterization of changes in rates across a series of experiments performed under variable conditions; and 2), the detection of degenerate states that exhibit the same FRET efficiency but differ in their rates of transition. We apply this newly developed methodology to two studies of the bacterial ribosome, each exemplary of one of these two analysis tasks. We conclude with a discussion of model-selection techniques for determination of the appropriate number of conformational states. The code used to perform this analysis and a basic graphical user interface front end are available as open source software.

Adaptive and Unstructured Mesh Cleaving.

We propose a new strategy for boundary conforming meshing that decouples the problem of building tetrahedra of proper size and shape from the problem of conforming to complex, non-manifold boundaries. This approach is motivated by the observation that while several methods exist for adaptive tetrahedral meshing, they typically have difficulty at geometric boundaries. The proposed strategy avoids this conflict by extracting the boundary conforming constraint into a secondary step. We first build a background mesh having a desired set of tetrahedral properties, and then use a generalized stenciling method to divide, or "cleave", these elements to get a set of conforming tetrahedra, while limiting the impacts cleaving has on element quality. In developing this new framework, we make several technical contributions including a new method for building graded tetrahedral meshes as well as a generalization of the isosurface stuffing and lattice cleaving algorithms to unstructured background meshes.

Lattice cleaving: a multimaterial tetrahedral meshing algorithm with guarantees.

We introduce a new algorithm for generating tetrahedral meshes that conform to physical boundaries in volumetric domains consisting of multiple materials. The proposed method allows for an arbitrary number of materials, produces high-quality tetrahedral meshes with upper and lower bounds on dihedral angles, and guarantees geometric fidelity. Moreover, the method is combinatoric so its implementation enables rapid mesh construction. These meshes are structured in a way that also allows grading, to reduce element counts in regions of homogeneity. Additionally, we provide proofs showing that both element quality and geometric fidelity are bounded using this approach.

Hierarchically-coupled hidden Markov models for learning kinetic rates from single-molecule data.

We address the problem of analyzing sets of noisy time-varying signals that all report on the same process but confound straightforward analyses due to complex inter-signal heterogeneities and measurement artifacts. In particular we consider single-molecule experiments which indirectly measure the distinct steps in a biomolecular process via observations of noisy time-dependent signals such as a fluorescence intensity or bead position. Straightforward hidden Markov model (HMM) analyses attempt to characterize such processes in terms of a set of conformational states, the transitions that can occur between these states, and the associated rates at which those transitions occur; but require ad-hoc post-processing steps to combine multiple signals. Here we develop a hierarchically coupled HMM that allows experimentalists to deal with inter-signal variability in a principled and automatic way. Our approach is a generalized expectation maximization hyperparameter point estimation procedure with variational Bayes at the level of individual time series that learns an single interpretable representation of the overall data generating process.

Lattice Cleaving: Conforming Tetrahedral Meshes of Multimaterial Domains with Bounded Quality.

We introduce a new algorithm for generating tetrahedral meshes that conform to physical boundaries in volumetric domains consisting of multiple materials. The proposed method allows for an arbitrary number of materials, produces high-quality tetrahedral meshes with upper and lower bounds on dihedral angles, and guarantees geometric fidelity. Moreover, the method is combinatoric so its implementation enables rapid mesh construction. These meshes are structured in a way that also allows grading, in order to reduce element counts in regions of homogeneity.

Particle systems for adaptive, isotropic meshing of CAD models.

We present a particle-based approach for generating adaptive triangular surface and tetrahedral volume meshes from computer-aided design models. Input shapes are treated as a collection of smooth, parametric surface patches that can meet non-smoothly on boundaries. Our approach uses a hierarchical sampling scheme that places particles on features in order of increasing dimensionality. These particles reach a good distribution by minimizing an energy computed in 3D world space, with movements occurring in the parametric space of each surface patch. Rather than using a pre-computed measure of feature size, our system automatically adapts to both curvature as well as a notion of topological separation. It also enforces a measure of smoothness on these constraints to construct a sizing field that acts as a proxy to piecewise-smooth feature size. We evaluate our technique with comparisons against other popular triangular meshing techniques for this domain.

Transfer RNA-mediated regulation of ribosome dynamics during protein synthesis.

Translocation of tRNAs through the ribosome during protein synthesis involves large-scale structural rearrangement of the ribosome and ribosome-bound tRNAs that is accompanied by extensive and dynamic remodeling of tRNA-ribosome interactions. How the rearrangement of individual tRNA-ribosome interactions influences tRNA movement during translocation, however, remains largely unknown. To address this question, we used single-molecule FRET to characterize the dynamics of ribosomal pretranslocation (PRE) complex analogs carrying either wild-type or systematically mutagenized tRNAs. Our data reveal how specific tRNA-ribosome interactions regulate the rate of PRE complex rearrangement into a critical, on-pathway translocation intermediate and how these interactions control the stability of the resulting configuration. Notably, our results suggest that the conformational flexibility of the tRNA molecule has a crucial role in directing the structural dynamics of the PRE complex during translocation.

Computational design and selections for an engineered, thermostable terpene synthase.

Terpenoids include structurally diverse antibiotics, flavorings, and fragrances. Engineering terpene synthases for control over the synthesis of such compounds represents a long sought goal. We report computational design, selections, and assays of a thermostable mutant of tobacco 5-epi-aristolochene synthase (TEAS) for the catalysis of carbocation cyclization reactions at elevated temperatures. Selection for thermostability included proteolytic digestion followed by capture of intact proteins. Unlike the wild-type enzyme, the mutant TEAS retains enzymatic activity at 65°C. The thermostable terpene synthase variant denatures above 80°C, approximately twice the temperature of the wild-type enzyme.

Graphical models for inferring single molecule dynamics.

BACKGROUND: The recent explosion of experimental techniques in single molecule biophysics has generated a variety of novel time series data requiring equally novel computational tools for analysis and inference. This article describes in general terms how graphical modeling may be used to learn from biophysical time series data using the variational Bayesian expectation maximization algorithm (VBEM). The discussion is illustrated by the example of single-molecule fluorescence resonance energy transfer (smFRET) versus time data, where the smFRET time series is modeled as a hidden Markov model (HMM) with Gaussian observables. A detailed description of smFRET is provided as well. RESULTS: The VBEM algorithm returns the model's evidence and an approximating posterior parameter distribution given the data. The former provides a metric for model selection via maximum evidence (ME), and the latter a description of the model's parameters learned from the data. ME/VBEM provide several advantages over the more commonly used approach of maximum likelihood (ML) optimized by the expectation maximization (EM) algorithm, the most important being a natural form of model selection and a well-posed (non-divergent) optimization problem. CONCLUSIONS: The results demonstrate the utility of graphical modeling for inference of dynamic processes in single molecule biophysics.

Learning rates and states from biophysical time series: a Bayesian approach to model selection and single-molecule FRET data.

Time series data provided by single-molecule Förster resonance energy transfer (smFRET) experiments offer the opportunity to infer not only model parameters describing molecular complexes, e.g., rate constants, but also information about the model itself, e.g., the number of conformational states. Resolving whether such states exist or how many of them exist requires a careful approach to the problem of model selection, here meaning discrimination among models with differing numbers of states. The most straightforward approach to model selection generalizes the common idea of maximum likelihood--selecting the most likely parameter values--to maximum evidence: selecting the most likely model. In either case, such an inference presents a tremendous computational challenge, which we here address by exploiting an approximation technique termed variational Bayesian expectation maximization. We demonstrate how this technique can be applied to temporal data such as smFRET time series; show superior statistical consistency relative to the maximum likelihood approach; compare its performance on smFRET data generated from experiments on the ribosome; and illustrate how model selection in such probabilistic or generative modeling can facilitate analysis of closely related temporal data currently prevalent in biophysics. Source code used in this analysis, including a graphical user interface, is available open source via http://vbFRET.sourceforge.net.

Allosteric collaboration between elongation factor G and the ribosomal L1 stalk directs tRNA movements during translation.

Determining the mechanism by which tRNAs rapidly and precisely transit through the ribosomal A, P, and E sites during translation remains a major goal in the study of protein synthesis. Here, we report the real-time dynamics of the L1 stalk, a structural element of the large ribosomal subunit that is implicated in directing tRNA movements during translation. Within pretranslocation ribosomal complexes, the L1 stalk exists in a dynamic equilibrium between open and closed conformations. Binding of elongation factor G (EF-G) shifts this equilibrium toward the closed conformation through one of at least two distinct kinetic mechanisms, where the identity of the P-site tRNA dictates the kinetic route that is taken. Within posttranslocation complexes, L1 stalk dynamics are dependent on the presence and identity of the E-site tRNA. Collectively, our data demonstrate that EF-G and the L1 stalk allosterically collaborate to direct tRNA translocation from the P to the E sites, and suggest a model for the release of E-site tRNA.

Transcription factor logic using chemical complementation.

Chemical complementation was used to make a transcription factor circuit capable of performing complex Boolean logic.

Molecular Dynamics Simulation of WSK-3, a Computationally Designed, Water-Soluble Variant of the Integral Membrane Protein KcsA.

Poor solubility and low expression levels often make membrane proteins difficult to study. An alternative to the use of detergents to solubilize these aggregation-prone proteins is the partial redesign of the sequence so as to confer water solubility. Recently, computationally assisted membrane protein solubilization (CAMPS) has been reported, where exposed hydrophobic residues on a protein's surface are computationally redesigned. Herein, the structure and fluctuations of a designed, water-soluble variant of KcsA (WSK-3) were studied using molecular dynamics simulations. The root mean square deviation of the protein from its starting structure, where the backbone coordinates are those of KcsA, was 1.8 angstroms. The structure of salt bridges involved in structural specificity and solubility were examined. The preferred configuration of ions and water in the selectivity filter of WSK-3 was consistent with the reported preferences for KcsA. The structure of the selectivity filter was maintained, which is consistent with WSK-3 having an affinity for agitoxin2 comparable to that of wild-type KcsA. In contrast to KcsA, the central cavity's side chains were observed to reorient, allowing water diffusion through the side of the cavity wall. These simulations provide an atomistic analysis of the CAMPS strategy and its implications for further investigations of membrane proteins.