Pathological conformations of disease mutant Ryanodine Receptors revealed by cryo-EM.

Ryanodine Receptors (RyRs) are massive channels that release Ca2+ from the endoplasmic and sarcoplasmic reticulum. Hundreds of mutations are linked to malignant hyperthermia (MH), myopathies, and arrhythmias. Here, we explore the first MH mutation identified in humans by providing cryo-EM snapshots of the pig homolog, R615C, showing that it affects an interface between three solenoid regions. We also show the impact of apo-calmodulin (apoCaM) and how it can induce opening by bending of the bridging solenoid, mediated by its N-terminal lobe. For R615C RyR1, apoCaM binding abolishes a pathological 'intermediate' conformation, distributing the population to a mixture of open and closed channels, both different from the structure without apoCaM. Comparisons show that the mutation primarily affects the closed state, inducing partial movements linked to channel activation. This shows that disease mutations can cause distinct pathological conformations of the RyR and facilitate channel opening by disrupting interactions between different solenoid regions.

A Two-Ended Data-Driven Accelerated Sampling Method for Exploring the Transition Pathways between Two Known States of Protein.

Conformational transitions of protein between different states are often associated with their biological functions. These dynamic processes, however, are usually not easy to be well characterized by experimental measurements, mainly because of inadequate temporal and spatial resolution. Meantime, sampling of configuration space with molecular dynamics (MD) simulations is still a challenge. Here we proposed a robust two-ended data-driven accelerated (teDA2) conformational sampling method, which drives the structural change in an adaptively updated feature space without introducing a bias potential. teDA2 was applied to explore adenylate kinase (ADK), a model with well characterized "open" and "closed" states. A single conformational transition event of ADK could be achieved within only a few or tens of nanoseconds sampled with teDA2. By analyzing hundreds of transition events, we reproduced different mechanisms and the associated pathways for domain motion of ADK reported in the literature. The multiroute characteristic of ADK was confirmed by the fact that some metastable states identified with teDA2 resemble available crystal structures determined at different conditions. This feature was further validated with Markov state modeling with independent MD simulations. Therefore, our work provides strong evidence for the conformational plasticity of protein, which is mainly due to the inherent degree of flexibility. As a reliable and efficient enhanced sampling protocol, teDA2 could be used to study the dynamics between functional states of various biomolecular machines.

Reconciling Simulated Ensembles of Apomyoglobin with Experimental Hydrogen/Deuterium Exchange Data Using Bayesian Inference and Multiensemble Markov State Models.

Hydrogen/deuterium exchange (HDX) is a powerful technique to investigate protein conformational dynamics at amino acid resolution. Because HDX provides a measurement of solvent exposure of backbone hydrogens, ensemble-averaged over potentially slow kinetic processes, it has been challenging to use HDX protection factors to refine structural ensembles obtained from molecular dynamics simulations. This entails dual challenges: (1) identifying structural observables that best correlate with backbone amide protection from exchange and (2) restraining these observables in molecular simulations to model ensembles consistent with experimental measurements. Here, we make significant progress on both fronts. First, we describe an improved predictor of HDX protection factors from structural observables in simulated ensembles, parametrized from ultralong molecular dynamics simulation trajectory data, with a Bayesian inference approach used to retain the full posterior distribution of model parameters. We next present a new method for obtaining simulated ensembles in agreement with experimental HDX protection factors, in which molecular simulations are performed at various temperatures and restraint biases and used to construct multiensemble Markov State Models (MSMs). Finally, the BICePs (Bayesian Inference of Conformational Populations) algorithm is then used with our HDX protection factor predictor to infer which thermodynamic ensemble agrees best with the experiment and estimate populations of each conformational state in the MSM. To illustrate the approach, we use a combination of HDX protection factor restraints and chemical shift restraints to model the conformational ensemble of apomyoglobin at pH 6. The resulting ensemble agrees well with the experiment and gives insight into the all-atom structure of disordered helices F and H in the absence of heme.

The molten globule state is unusually deformable under mechanical force.

Recently, the role of force in cellular processes has become more evident, and now with advances in force spectroscopy, the response of proteins to force can be directly studied. Such studies have found that native proteins are brittle, and thus not very deformable. Here, we examine the mechanical properties of a class of intermediates referred to as the molten globule state. Using optical trap force spectroscopy, we investigated the response to force of the native and molten globule states of apomyoglobin along different pulling axes. Unlike natively folded proteins, the molten globule state of apomyoglobin is compliant (large distance to the transition state); this large compliance means that the molten globule is more deformable and the unfolding rate is more sensitive to force (the application of force or tension will have a more dramatic effect on the unfolding rate). Our studies suggest that these are general properties of molten globules and could have important implications for mechanical processes in the cell.

A simple simulation model can reproduce the thermodynamic folding intermediate of apoflavodoxin.

Flavodoxins are single domain proteins with an alpha/beta structure, whose function and folding have been well studied. Detailed experiments have shown that several members of this protein family present a stable intermediate, which accumulates along the folding process. In this work, we use a coarse-grained model for protein folding, whose interactions are based on the topology of the native state, to analyze the thermodynamic characteristics of the folding of Anabaena apoflavodoxin. Our model shows evidence for the existence of a thermodynamic folding intermediate, which reaches a significant population along the thermal transition. According to our simulation results, the intermediate is compact, well packed, and involves distortions of the native structure similar to those experimentally found. These mainly affect the long loop in the protein surface comprising residues 120-139. Although the agreement between simulation and experiment is not perfect, something impossible for a crude model, our results show that the topology of the native state is able to dictate a folding process which includes the presence of an intermediate for this protein.

Development of an aequorin luminescence calcium assay for high-throughput screening using a plate reader, the LumiLux.

A luminescence assay using a new plate reader, the LumiLux (PerkinElmer, Waltham, MA), has been validated for high-throughput screening (HTS). In this study, we compared the aequorin luminescence-based calcium mobilization assay to the fluorescence-based calcium assay. A cell line stably co-expressing apo-aequorin, a chimeric G-protein, and a G-protein-coupled dopamine receptor was used to screen a collection of 8,106 compounds using the Hamamatsu Photonics (Bridgewater, NJ) FDSS6000 and LumiLux as the plate readers. The assay parameters evaluated included hit rate correlation, signal-to-noise ratio, and overall assay performance calculated by Z and standard deviation. The average Z values and hit rates were comparable between assay platforms;however, the standard deviation for the agonist aequorin assay was significantly smaller. There was also a significant decrease in the number of false-positives with the aequorin assay. These results suggest that the aequorin assay in combination with the new plate reader, LumiLux, provides a simple, cost-effective, robust, and sensitive assay for HTS

Unfolding of globular proteins: monte carlo dynamics of a realistic reduced model.

Reduced lattice models of proteins and Monte Carlo dynamics were used to simulate the initial stages of the unfolding of several proteins of various structural types, and the results were compared to experiment. The models semiquantitatively reproduce the approximate order of events of unfolding as well as subtle mutation effects and effects resulting from differences in sequences of similar folds. The short-time mobility of particular residues, observed in simulations, correlates with the crystallographic temperature factor. The main factor controlling unfolding is the native state topology, with sequence playing a less important role. The correlation with various experiments, especially for sequence-specific effects, strongly suggests that properly designed reduced models of proteins can be used for qualitative studies (or prediction) of protein unfolding pathways.

Indirect assessment of small hydrophobic ligand binding to a model protein using a combination of ESI MS and HDX/ESI MS.

Direct mass spectrometric characterization of interactions between proteins and small hydrophobic ligands often poses a serious problem due to the complex instability in the gas phase. We have developed a method that probes the efficacy of ligand-protein interactions indirectly by monitoring changes in protein flexibility. The latter is assessed quantitatively using a combination of charge state distribution analysis and amide hydrogen exchange under both native and mildly denaturing conditions. The method was used to evaluate binding of a model protein cellular retinoic acid binding protein I to its natural ligand all-trans retinoic acid (RA), isomers 13-cis- and 9-cis-RA, and retinol, yielding the following order of ligand affinities: All-trans RA > 9-cis RA > 13-cis RA, with no detectable binding of retinol. This order is in agreement with the results of earlier fluorimetric titration studies. Furthermore, binding energy of the protein to each of retinoic acid isomers was determined based on the measured hydrogen exchange kinetics data acquired under native conditions.

Coarse-grained simulations of conformational dynamics of proteins: application to apomyoglobin.

A coarse-grained dynamic Monte Carlo method is proposed for investigating the conformational dynamics of proteins. Each residue is represented by two interaction sites, one at the alpha-carbon, and the other on the amino acid sidechain. Geometry and energy parameters extracted from databank structures are used. The method is applied to the crystal structure of apomyoglobin (apo-Mb). Equilibrium and dynamic properties of apo-Mb are characterized within computation times one order of magnitude shorter than conventional molecular dynamics (MD) simulations. The calculated rms fluctuations in alpha-carbons are in good agreement with crystallographic temperature factors. Regions exhibiting enhanced conformational mobilities are identified. Among the loops connecting the eight helices A to H, the loop CD undergoes the fastest motions, leading to partial unwinding of helix D. Helix G is the most stable helix on the basis of the kinetic stability of dihedral angles, followed by the respective helices A, E, H, and B. These results, in agreement with H/D exchange and two-dimensional NMR experiments, as well as with MD simulations, lend support to the use of the proposed approach as an efficient, yet physically plausible, means of characterizing protein conformational dynamics.

Conformation studies on and assessment by spectral analysis of the protein-chromophore interaction of the macromolecular antitumor antibiotic C-1027.

Characterization of the secondary structure of the antitumor antibiotic C-1027 has been made from a comparison of C-1027 and its apoprotein by various analytical means. The results indicated the antibiotic to be abundant in beta-structure by measurements of Fourier-transform infrared (FT-IR) spectroscopy and the circular dichroism (CD) spectrum, and by a prediction of the secondary structure based on the amino acid sequence of the peptide. In comparison of the IR spectra of their proteins in D2O, the apoprotein exhibited a faster H-D exchange than C-1027, indicating an increase in the "non-motile parts" of the beta-sheets formed through the protein-chromophore interaction in holo-C-1027. The prediction of hydropathic index indicated the hydrophobic residues of the apoprotein to be predominantly located in the beta-sheet structures, suggesting hydrophobic interaction in the binding between chromophore and apoprotein. Further, the interaction between chromophore and apoprotein was detected by a fluorescence method. The result showed the dissociation constant (Kd) to be 6.88 x 10(-5) M, indicating that the chromophore is tightly bound to the protein moiety.