CryptoSite: Expanding the Druggable Proteome by Characterization and Prediction of Cryptic Binding Sites.

Many proteins have small-molecule binding pockets that are not easily detectable in the ligand-free structures. These cryptic sites require a conformational change to become apparent; a cryptic site can therefore be defined as a site that forms a pocket in a holo structure, but not in the apo structure. Because many proteins appear to lack druggable pockets, understanding and accurately identifying cryptic sites could expand the set of drug targets.

ModBase, a database of annotated comparative protein structure models and associated resources.

ModBase (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence-structure alignment, model building and model assessment (http://salilab.

Mapping polymerization and allostery of hemoglobin S using point mutations.

Hemoglobin is a complex system that undergoes conformational changes in response to oxygen, allosteric effectors, mutations, and environmental changes. Here, we study allostery and polymerization of hemoglobin and its variants by application of two previously described methods: (i) AllosMod for simulating allostery dynamics given two allosterically related input structures and (ii) a machine-learning method for dynamics- and structure-based prediction of the mutation impact on allostery (Weinkam et al. J.

All-atom ensemble modeling to analyze small-angle x-ray scattering of glycosylated proteins.

The flexible and heterogeneous nature of carbohydrate chains often renders glycoproteins refractory to traditional structure determination methods. Small-angle X-ray scattering (SAXS) can be a useful tool for obtaining structural information of these systems. All-atom modeling of glycoproteins with flexible glycan chains was applied to interpret the solution SAXS data for a set of glycoproteins.

Impact of mutations on the allosteric conformational equilibrium.

Allostery in a protein involves effector binding at an allosteric site that changes the structure and/or dynamics at a distant, functional site. In addition to the chemical equilibrium of ligand binding, allostery involves a conformational equilibrium between one protein substate that binds the effector and a second substate that less strongly binds the effector. We run molecular dynamics simulations using simple, smooth energy landscapes to sample specific ligand-induced conformational transitions, as defined by the effector-bound and effector-unbound protein structures.

Structure-based model of allostery predicts coupling between distant sites.

Allostery is a phenomenon that couples effector ligand binding at an allosteric site to a structural and/or dynamic change at a distant regulated site. To study an allosteric transition, we vary the size of the allosteric site and its interactions to construct a series of energy landscapes with pronounced minima corresponding to both the effector bound and unbound crystal structures. We use molecular dynamics to sample these landscapes.

Prediction of native-state hydrogen exchange from perfectly funneled energy landscapes.

Simulations based on perfectly funneled energy landscapes often capture many of the kinetic features of protein folding. We examined whether simulations based on funneled energy functions can also describe fluctuations in native-state protein ensembles. We quantitatively compared the site-specific local stability determined from structure-based folding simulations, with hydrogen exchange protection factors measured experimentally for ubiquitin, chymotrypsin inhibitor 2, and staphylococcal nuclease.

The folding energy landscape and free energy excitations of cytochrome c.

The covalently bound heme cofactor plays a dominant role in the folding of cytochrome c. Because of the complicated inorganic chemistry of the heme, some might consider the folding of cytochrome c to be a special case, following principles different from those used to describe the folding of proteins without cofactors. Recent investigations, however, demonstrate that common models describing folding for many proteins work well for cytochrome c when heme is explicitly introduced, generally providing results that agree with experimental observations.

The spectrum of biomolecular states and motions.

The universe of conformational substates of a protein molecule is huge. The complete energy landscape of proteins is, therefore, complex when studied at low temperature. Many experiments under physiological conditions commonly reveal a simpler spectrum of states.

Chemical frustration in the protein folding landscape: grand canonical ensemble simulations of cytochrome c.

A grand canonical formalism is developed to combine discrete simulations for chemically distinct species in equilibrium. Each simulation is based on a perturbed funneled landscape. The formalism is illustrated using the alkaline-induced transitions of cytochrome c as observed by FTIR spectroscopy and with various other experimental approaches.

Electrostatic effects on funneled landscapes and structural diversity in denatured protein ensembles.

The denatured state of proteins is heterogeneous and susceptible to general hydrophobic and electrostatic forces, but to what extent does the funneled nature of protein energy landscapes play a role in the unfolded ensemble? We simulate the denatured ensemble of cytochrome c using a series of models. The models pinpoint the efficacy of incorporating energetic funnels toward the native state in contrast with models having no native structure-seeking tendency. These models also contain varying strengths of electrostatic effects and hydrophobic collapse.

Characterization of alkaline transitions in ferricytochrome c using carbon-deuterium infrared probes.

The alkaline-induced structural transitions of ferricytochrome c have been studied intensively as a model for how changes in metal ligation contribute to protein function and folding. Previous studies have demonstrated that multiple non-native species accumulate with increasing pH. Here, we used a combination of experiments and simulations to provide a high-resolution view of the changes associated with increasing alkaline conditions.

Origins of barriers and barrierless folding in BBL.

In contrast to classical chemical phenomenology, theory suggests that proteins may undergo downhill folding without an activation barrier under certain thermodynamic conditions. Recently, the BBL protein was proposed to fold by such a downhill scenario, but discrepancies between experimental results found in different groups argue against this. After briefly reviewing the major experimental studies of the BBL folding mechanism, we show that simulations of both coarse-grained and atomistic models can reconcile the seemingly conflicting observations.

A funneled energy landscape for cytochrome c directly predicts the sequential folding route inferred from hydrogen exchange experiments.

Proteins fold through a variety of mechanisms. For a given protein, folding routes largely depend on the protein's stability and its native-state geometry, because the landscape is funneled. These ideas are corroborated for cytochrome c by using a coarse-grained topology-based model with a perfect funnel landscape that includes explicit modeling of the heme.

Building programmable jigsaw puzzles with RNA.

One challenge in supramolecular chemistry is the design of versatile, self-assembling building blocks to attain total control of arrangement of matter at a molecular level. We have achieved reliable prediction and design of the three-dimensional structure of artificial RNA building blocks to generate molecular jigsaw puzzle units called tectosquares. They can be programmed with control over their geometry, topology, directionality, and addressability to algorithmically self-assemble into a variety of complex nanoscopic fabrics with predefined periodic and aperiodic patterns and finite dimensions.