Efficient sampling of protein structures by model hopping.
Ontology highlight
ABSTRACT:
SUBMITTER: Kwak W
PROVIDER: S-EPMC1356668 | biostudies-literature | 2005 Sep
REPOSITORIES: biostudies-literature
Publications
Efficient sampling of protein structures by model hopping.
Kwak Wooseop W Hansmann Ulrich H E UH
Physical review letters 20050922 13
We introduce a novel simulation method, model hopping, that enhances sampling of low-energy configurations in complex systems. The approach is illustrated for a protein-folding problem. Thermodynamic quantities of proteins with up to 46 residues are evaluated from all-atom simulations with this method. ...[more]
Similar Datasets
Project description:Clearly, a protein cannot sample all of its conformations (e.g., approximately 3(100) approximately 10(48) for a 100 residue protein) on an in vivo folding timescale (<1 s). To investigate how the conformational dynamics of a protein can accommodate subsecond folding time scales, we introduce the concept of the native topomer, which is the set of all structures similar to the native structure (obtainable from the native structure through local backbone coordinate transformations that do not disrupt the covalent bonding of the peptide backbone). We have developed a computational procedure for estimating the number of distinct topomers required to span all conformations (compact and semicompact) for a polypeptide of a given length. For 100 residues, we find approximately 3 x 10(7) distinct topomers. Based on the distance calculated between different topomers, we estimate that a 100-residue polypeptide diffusively samples one topomer every approximately 3 ns. Hence, a 100-residue protein can find its native topomer by random sampling in just approximately 100 ms. These results suggest that subsecond folding of modest-sized, single-domain proteins can be accomplished by a two-stage process of (i) topomer diffusion: random, diffusive sampling of the 3 x 10(7) distinct topomers to find the native topomer ( approximately 0.1 s), followed by (ii) intratopomer ordering: nonrandom, local conformational rearrangements within the native topomer to settle into the precise native state.
| S-EPMC15813 | biostudies-literature
Project description:We have investigated photoinduced hole hopping in a Pseudomonas aeruginosa azurin mutant Re126WWCuI, where two adjacent tryptophan residues (W124 and W122) are inserted between the CuI center and a Re photosensitizer coordinated to a H126 imidazole (Re = ReI(H126)(CO)3(dmp)+, dmp = 4,7-dimethyl-1,10-phenanthroline). Optical excitation of this mutant in aqueous media (≤40 μM) triggers 70 ns electron transport over 23 Å, yielding a long-lived (120 μs) ReI(H126)(CO)3(dmp•-)WWCuII product. The Re126FWCuI mutant (F124, W122) is not redox-active under these conditions. Upon increasing the concentration to 0.2-2 mM, {Re126WWCuI}2 and {Re126FWCuI}2 are formed with the dmp ligand of the Re photooxidant of one molecule in close contact (3.8 Å) with the W122' indole on the neighboring chain. In addition, {Re126WWCuI}2 contains an interfacial tryptophan quadruplex of four indoles (3.3-3.7 Å apart). In both mutants, dimerization opens an intermolecular W122' → //*Re ET channel (// denotes the protein interface, *Re is the optically excited sensitizer). Excited-state relaxation and ET occur together in two steps (time constants of ∼600 ps and ∼8 ns) that lead to a charge-separated state containing a Re(H126)(CO)3(dmp•-)//(W122•+)' unit; then (CuI)' is oxidized intramolecularly (60-90 ns) by (W122•+)', forming ReI(H126)(CO)3(dmp•-)WWCuI//(CuII)'. The photocycle is closed by ∼1.6 μs ReI(H126)(CO)3(dmp•-) → //(CuII)' back ET that occurs over 12 Å, in contrast to the 23 Å, 120 μs step in Re126WWCuI. Importantly, dimerization makes Re126FWCuI photoreactive and, as in the case of {Re126WWCuI}2, channels the photoproduced "hole" to the molecule that was not initially photoexcited, thereby shortening the lifetime of ReI(H126)(CO)3(dmp•-)//CuII. Although two adjacent W124 and W122 indoles dramatically enhance CuI → *Re intramolecular multistep ET, the tryptophan quadruplex in {Re126WWCuI}2 does not accelerate intermolecular electron transport; instead, it acts as a hole storage and crossover unit between inter- and intramolecular ET pathways. Irradiation of {Re126WWCuII}2 or {Re126FWCuII}2 also triggers intermolecular W122' → //*Re ET, and the Re(H126)(CO)3(dmp•-)//(W122•+)' charge-separated state decays to the ground state by ∼50 ns ReI(H126)(CO)3(dmp•-)+ → //(W122•+)' intermolecular charge recombination. Our findings shed light on the factors that control interfacial hole/electron hopping in protein complexes and on the role of aromatic amino acids in accelerating long-range electron transport.
| S-EPMC6384139 | biostudies-literature
Project description:Associative memory Hamiltonian structure prediction potentials are not overly rugged, thereby suggesting their landscapes are like those of actual proteins. In the present contribution we show how basin-hopping global optimization can identify low-lying minima for the corresponding mildly frustrated energy landscapes. For small systems the basin-hopping algorithm succeeds in locating both lower minima and conformations closer to the experimental structure than does molecular dynamics with simulated annealing. For large systems the efficiency of basin-hopping decreases for our initial implementation, where the steps consist of random perturbations to the Cartesian coordinates. We implemented umbrella sampling using basin-hopping to further confirm when the global minima are reached. We have also improved the energy surface by employing bioinformatic techniques for reducing the roughness or variance of the energy surface. Finally, the basin-hopping calculations have guided improvements in the excluded volume of the Hamiltonian, producing better structures. These results suggest a novel and transferable optimization scheme for future energy function development.
| S-EPMC2674628 | biostudies-literature
Project description:A new and efficient Monte Carlo algorithm for sampling protein configurations in the continuous space is presented; the efficiency of this algorithm, named Local Moves for Proteins (LMProt), was compared to other alternative algorithms. For this purpose, we used an intrachain interaction energy function that is proportional to the root mean square deviation (rmsd) with respect to alpha-carbons from native structures of real proteins. For phantom chains, the LMProt method is approximately 10(4) and 20 times faster than the algorithms Thrashing (no local moves) and Sevenfold Way (local moves), respectively. Additionally, the LMProt was tested for real chains (excluded-volume all-atoms model); proteins 5NLL (138 residues) and 1BFF (129 residues) were used to determine the folding success xi as a function of the number eta of residues involved in the chain movements, and as a function of the maximum amplitude of atomic displacement delta r(max). Our results indicate that multiple local moves associated with relative chain flexibility, controlled by appropriate adjustments for eta and delta r(max), are essential for configurational search efficiency.
| S-EPMC1304563 | biostudies-literature
Project description:RNA-protein complexes underlie numerous cellular processes including translation, splicing, and posttranscriptional regulation of gene expression. The structures of these complexes are crucial to their functions but often elude high-resolution structure determination. Computational methods are needed that can integrate low-resolution data for RNA-protein complexes while modeling de novo the large conformational changes of RNA components upon complex formation. To address this challenge, we describe RNP-denovo, a Rosetta method to simultaneously fold-and-dock RNA to a protein surface. On a benchmark set of diverse RNA-protein complexes not solvable with prior strategies, RNP-denovo consistently sampled native-like structures with better than nucleotide resolution. We revisited three past blind modeling challenges involving the spliceosome, telomerase, and a methyltransferase-ribosomal RNA complex in which previous methods gave poor results. When coupled with the same sparse FRET, crosslinking, and functional data used previously, RNP-denovo gave models with significantly improved accuracy. These results open a route to modeling global folds of RNA-protein complexes from low-resolution data.
| S-EPMC6318048 | biostudies-literature
Project description:Human decisions are based on finite information, which makes them inherently imprecise. But what determines the degree of such imprecision? Here, we develop an efficient coding framework for higher-level cognitive processes in which information is represented by a finite number of discrete samples. We characterize the sampling process that maximizes perceptual accuracy or fitness under the often-adopted assumption that full adaptation to an environmental distribution is possible, and show how the optimal process differs when detailed information about the current contextual distribution is costly. We tested this theory on a numerosity discrimination task, and found that humans efficiently adapt to contextual distributions, but in the way predicted by the model in which people must economize on environmental information. Thus, understanding decision behavior requires that we account for biological restrictions on information coding, challenging the often-adopted assumption of precise prior knowledge in higher-level decision systems.
| S-EPMC7492090 | biostudies-literature
Project description:Scientific advance is often driven by identifying conceptually simple models underlying complex phenomena. This process commonly ignores imperfections which, however, might give rise to non-trivial collective behavior. For example, already a small amount of disorder can dramatically change the transport properties of a system compared to the underlying simple model. While systems with disordered potentials were already studied in detail, experimental investigations on systems with disordered hopping are still in its infancy. To this end, we experimentally study a dipole-dipole-interacting three-dimensional Rydberg system and map it onto a simple XY model with random couplings by spectroscopic evidence. We discuss the localization-delocalization crossover emerging in the model and present experimental signatures of it. Our results demonstrate that Rydberg systems are a useful platform to study random hopping models with the ability to access the microscopic degrees of freedom. This will allow to study transport processes and localization phenomena in random hopping models with a high level of control.
| S-EPMC8632899 | biostudies-literature
Project description:High-resolution structure determination of small proteins in solution is one of the big assets of NMR spectroscopy in structural biology. Improvements in the efficiency of NMR structure determination by advances in NMR experiments and automation of data handling therefore attracts continued interest. Here, non-uniform sampling (NUS) of 3D heteronuclear-resolved [(1)H,(1)H]-NOESY data yielded two- to three-fold savings of instrument time for structure determinations of soluble proteins. With the 152-residue protein NP_372339.1 from Staphylococcus aureus and the 71-residue protein NP_346341.1 from Streptococcus pneumonia we show that high-quality structures can be obtained with NUS NMR data, which are equally well amenable to robust automated analysis as the corresponding uniformly sampled data.
| S-EPMC4576834 | biostudies-literature
Project description:Multi-step electron transfer reactions are important to the function of many cellular systems. The ways in which such systems have evolved to direct electrons along specific pathways are largely understood, but less so are the ways in which the reduction-oxidation potentials of individual redox sites are controlled. We prepared a series of three new artificial variants of Pseudomonas aeruginosa azurin where a tyrosine (Tyr109) is situated between the native Cu ion and a Ru(II) photosensitizer tethered to a histidine (His107). Arginine, glutamine, or methionine were introduced as position 122, which is near to Tyr109. We investigated the rate of CuI oxidation by a flash-quench generated Ru(III) oxidant over pH values from 5 to 9. While the identity of the residue at position 122 affects some of the physical properties of Tyr109, the rates of CuI oxidation are only weakly dependent on the identity of the residue at 122. The results highlight that more work is still needed to understand how non-covalent interactions of redox active groups are affected in redox proteins.
| S-EPMC10818705 | biostudies-literature
Project description:Both experimental and computational methods are available to gather information about a protein's conformational space and interpret changes in protein structure. However, experimentally observing and computationally modeling large proteins remain critical challenges for structural biology. Our work aims at addressing these challenges by combining computational and experimental techniques relying on each other to overcome their respective limitations. Indeed, despite its advantages, an experimental technique such as hydrogen-exchange monitoring cannot produce structural models because of its low resolution. Additionally, the computational methods that can generate such models suffer from the curse of dimensionality when applied to large proteins. Adopting a common solution to this issue, we have recently proposed a framework in which our computational method for protein conformational sampling is biased by experimental hydrogen-exchange data. In this paper, we present our latest application of this computational framework: generating an atomic-resolution structural model for an unknown protein state. For that, starting from an available protein structure, we explore the conformational space of this protein, using hydrogen-exchange data on this unknown state as a guide. We have successfully used our computational framework to generate models for three proteins of increasing size, the biggest one undergoing large-scale conformational changes.
| S-EPMC6280153 | biostudies-literature