Dataset Information

Efficient Ensemble Refinement by Reweighting.

ABSTRACT: Ensemble refinement produces structural ensembles of flexible and dynamic biomolecules by integrating experimental data and molecular simulations. Here we present two efficient numerical methods to solve the computationally challenging maximum-entropy problem arising from a Bayesian formulation of ensemble refinement. Recasting the resulting constrained weight optimization problem into an unconstrained form enables the use of gradient-based algorithms. In two complementary formulations that differ in their dimensionality, we optimize either the log-weights directly or the generalized forces appearing in the explicit analytical form of the solution. We first demonstrate the robustness, accuracy, and efficiency of the two methods using synthetic data. We then use NMR J-couplings to reweight an all-atom molecular dynamics simulation ensemble of the disordered peptide Ala-5 simulated with the AMBER99SB*-ildn-q force field. After reweighting, we find a consistent increase in the population of the polyproline-II conformations and a decrease of ?-helical-like conformations. Ensemble refinement makes it possible to infer detailed structural models for biomolecules exhibiting significant dynamics, such as intrinsically disordered proteins, by combining input from experiment and simulation in a balanced manner.

SUBMITTER: Kofinger J

PROVIDER: S-EPMC6727217 | biostudies-literature | 2019 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Efficient Ensemble Refinement by Reweighting.

Köfinger Jürgen J Stelzl Lukas S LS Reuter Klaus K Allande César C Reichel Katrin K Hummer Gerhard G

Journal of chemical theory and computation 20190417 5

Ensemble refinement produces structural ensembles of flexible and dynamic biomolecules by integrating experimental data and molecular simulations. Here we present two efficient numerical methods to solve the computationally challenging maximum-entropy problem arising from a Bayesian formulation of ensemble refinement. Recasting the resulting constrained weight optimization problem into an unconstrained form enables the use of gradient-based algorithms. In two complementary formulations that diff ...[more]

PMID: 30939006

Similar Datasets

Project description:The cytoplasmic heme binding protein from Pseudomonas aeruginosa, PhuS, plays two essential roles in regulating heme uptake and iron homeostasis. First, PhuS shuttles exogenous heme to heme oxygenase (HemO) for degradation and iron release. Second, PhuS binds DNA and modulates the transcription of the prrF/H small RNAs (sRNAs) involved in the iron-sparing response. Heme binding to PhuS regulates this dual function, as the unliganded form binds DNA, whereas the heme-bound form binds HemO. Crystallographic studies revealed nearly identical structures for apo- and holo-PhuS, and yet numerous solution-based measurements indicate that heme binding is accompanied by large conformational rearrangements. In particular, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of apo- versus holo-PhuS revealed large differences in deuterium uptake, notably in α-helices 6, 7, and 8 (α6,7,8), which contribute to the heme binding pocket. These helices were mostly labile in apo-PhuS but largely protected in holo-PhuS. In contrast, in silico-predicted deuterium uptake levels of α6,7,8 from molecular dynamics (MD) simulations of the apo- and holo-PhuS structures are highly similar, consistent only with the holo-PhuS HDX-MS data. To rationalize this discrepancy between crystal structures, simulations, and observed HDX-MS, we exploit a recently developed computational approach (HDXer) that fits the relative weights of conformational populations within an ensemble of structures to conform to a target set of HDX-MS data. Here, a combination of enhanced sampling MD, HDXer, and dimensionality reduction analysis reveals an apo-PhuS conformational landscape in which α6, 7, and 8 are significantly rearranged compared to the crystal structure, including a loss of secondary structure in α6 and the displacement of α7 toward the HemO binding interface. Circular dichroism analysis confirms the loss of secondary structure, and the extracted ensembles of apo-PhuS and of heme-transfer-impaired H212R mutant, are consistent with known heme binding and transfer properties. The proposed conformational landscape provides structural insights into the modulation by heme of the dual function of PhuS.

Project description:A critical step in injury-induced initiation of blood coagulation is the formation of the complex between the trypsin-like protease coagulation factor VIIa (FVIIa) and its cofactor tissue factor (TF), which converts FVIIa from an intrinsically poor enzyme to an active protease capable of activating zymogens of downstream coagulation proteases. Unlike its constitutively active ancestor trypsin, FVIIa is allosterically activated (by TF). Here, ensemble refinement of crystallographic structures, which uses multiple copies of the entire structure as a means of representing structural flexibility, is applied to explore the impacts of inhibitor binding to trypsin and FVIIa, as well as cofactor binding to FVIIa. To assess the conformational flexibility and its role in allosteric pathways in these proteases, main-chain hydrogen bond networks are analyzed by calculating the hydrogen-bond propensity. Mapping pairwise propensity differences between relevant structures shows that binding of the inhibitor benzamidine to trypsin has a minor influence on the protease flexibility. For FVIIa, in contrast, the protease domain is "locked" into the catalytically competent trypsin-like configuration upon benzamidine binding as indicated by the stabilization of key structural features: the nonprime binding cleft and the oxyanion hole are stabilized, and the effect propagates from the active site region to the calcium-binding site and to the vicinity of the disulphide bridge connecting with the light chain. TF binding to FVIIa furthermore results in stabilization of the 170 loop, which in turn propagates an allosteric signal from the TF-binding region to the active site. Analyses of disulphide bridge energy and flexibility reflect the striking stability difference between the unregulated enzyme and the allosterically activated form after inhibitor or cofactor binding. The ensemble refinement analyses show directly, for the first time to our knowledge, whole-domain structural footprints of TF-induced allosteric networks present in x-ray crystallographic structures of FVIIa, which previously only have been hypothesized or indirectly inferred.

Dataset Information

Efficient Ensemble Refinement by Reweighting.

Publications

Efficient Ensemble Refinement by Reweighting.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets