Probing allosteric regulations with coevolution-driven molecular simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Colizzi F
PROVIDER: S-EPMC8442858 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The receptor recognition of the novel coronavirus SARS-CoV-2 relies on the "down-to-up" conformational change in the receptor-binding domain (RBD) of the spike (S) protein. Therefore, understanding the process of this change at the molecular level facilitates the design of therapeutic agents. With the help of coarse-grained molecular dynamic simulations, we provide evidence showing that the conformational dynamics of the S protein are globally cooperative. Importantly, an allosteric path was discovered that correlates the motion of the RBD with the motion of the junction between the subdomain 1 (SD1) and the subdomain 2 (SD2) of the S protein. Building on this finding, we designed non-RBD binding modulators to inhibit SARS-CoV-2 by prohibiting the conformational change of the S protein. Their inhibition effect and function stages at inhibiting SARS-CoV-2 were evaluated experimentally. In summary, our studies establish a molecular basis for future therapeutic agent design through allosteric effects.
| S-EPMC8409148 | biostudies-literature
Project description:Tubulin heterodimers are the building blocks of microtubules, a major component of the cytoskeleton, whose mechanical properties are fundamental for the life of the cell. We uncover the microscopic origins of the mechanical response in microtubules by probing features of the energy landscape of the tubulin monomers and tubulin heterodimer. To elucidate the structures of the unfolding pathways and reveal the multiple unfolding routes, we performed simulations of a self-organized polymer (SOP) model of tubulin. The SOP representation, which is a coarse-grained description of chains, allows us to perform force-induced simulations at loading rates and time scales that closely match those used in single-molecule experiments. We show that the forced unfolding of each monomer involves a bifurcation in the pathways to the stretched state. After the unfolding of the C-term domain, the unraveling continues either from the N-term domain or from the middle domain, depending on the monomer and the pathway. In contrast to the unfolding complexity of the monomers, the dimer unfolds according to only one route corresponding to the unraveling of the C-term domain and part of the middle domain of beta-tubulin. We find that this surprising behavior is due to the viscoelastic properties of the interface between the monomers. We map precise features of the complex energy landscape of tubulin by surveying the structures of the various metastable intermediates, which, in the dimer case, are characterized only by changes in the beta-tubulin monomer.
| S-EPMC2572946 | biostudies-literature
Project description:We introduce a numerical scheme to evolve functional elastic materials that can accomplish a specified mechanical task. In this scheme, the number of solutions, their spatial architectures, and the correlations among them can be computed. As an example, we consider an "allosteric" task, which requires the material to respond specifically to a stimulus at a distant active site. We find that functioning materials evolve a less-constrained trumpet-shaped region connecting the stimulus and active sites, and that the amplitude of the elastic response varies nonmonotonically along the trumpet. As previously shown for some proteins, we find that correlations appearing during evolution alone are sufficient to identify key aspects of this design. Finally, we show that the success of this architecture stems from the emergence of soft edge modes recently found to appear near the surface of marginally connected materials. Overall, our in silico evolution experiment offers a window to study the relationship between structure, function, and correlations emerging during evolution.
| S-EPMC5347607 | biostudies-literature
Project description:Using replica exchange molecular dynamics simulations and an all-atom implicit solvent model, we probed the energetics of Abeta(10-40) fibril growth. The analysis of the interactions between incoming Abeta peptides and the fibril led us to two conclusions. First, considerable variations in fibril binding propensities are observed along the Abeta sequence. The peptides in the fibril and those binding to its edge interact primarily through their N-termini. Therefore, the mutations affecting the Abeta positions 10-23 are expected to have the largest impact on fibril elongation compared with those occurring in the C-terminus and turn. Second, we performed weak perturbations of the binding free energy landscape by scanning partial deletions of side-chain interactions at various Abeta sequence positions. The results imply that strong side-chain interactions--in particular, hydrophobic contacts--impede fibril growth by favoring disordered docking of incoming peptides. Therefore, fibril elongation may be promoted by moderate reduction of Abeta hydrophobicity. The comparison with available experimental data is presented.
| S-EPMC2711452 | biostudies-literature
Project description:Investigating ligand-regulated allosteric coupling between protein domains is fundamental to understand cell-life regulation. The Hsp70 family of chaperones represents an example of proteins in which ATP binding and hydrolysis at the Nucleotide Binding Domain (NBD) modulate substrate recognition at the Substrate Binding Domain (SBD). Herein, a comparative analysis of an allosteric (Hsp70-DnaK) and a non-allosteric structural homolog (Hsp110-Sse1) of the Hsp70 family is carried out through molecular dynamics simulations, starting from different conformations and ligand-states. Analysis of ligand-dependent modulation of internal fluctuations and local deformation patterns highlights the structural and dynamical changes occurring at residue level upon ATP-ADP exchange, which are connected to the conformational transition between closed and open structures. By identifying the dynamically responsive protein regions and specific cross-domain hydrogen-bonding patterns that differentiate Hsp70 from Hsp110 as a function of the nucleotide, we propose a molecular mechanism for the allosteric signal propagation of the ATP-encoded conformational signal.
| S-EPMC3531320 | biostudies-literature
Project description:Increasing data in allostery are requiring analysis of coupling relationships among different allosteric sites on a single protein. Here, based on our previous efforts on reversed allosteric communication theory, we have developed AlloReverse, a web server for multiscale analysis of multiple allosteric regulations. AlloReverse integrates protein dynamics and machine learning to discover allosteric residues, allosteric sites and regulation pathways. Especially, AlloReverse could reveal hierarchical relationships between different pathways and couplings among allosteric sites, offering a whole map of allostery. The web server shows a good performance in re-emerging known allostery. Moreover, we applied AlloReverse to explore global allostery on CDC42 and SIRT3. AlloReverse predicted novel allosteric sites and allosteric residues in both systems, and the functionality of sites was validated experimentally. It also suggests a possible scheme for combined therapy or bivalent drugs on SIRT3. Taken together, AlloReverse is a novel workflow providing a complete regulation map and is believed to aid target identification, drug design and understanding of biological mechanisms. AlloReverse is freely available to all users at https://mdl.shsmu.edu.cn/AlloReverse/ or http://www.allostery.net/AlloReverse/.
| S-EPMC10320067 | biostudies-literature
Project description:The solution dynamics of antibodies are critical to antibody function. We explore the internal solution dynamics of antibody molecules through the combination of time-resolved fluorescence anisotropy experiments on IgG1 with more than two microseconds of all-atom molecular dynamics (MD) simulations in explicit water, an order of magnitude more than in previous simulations. We analyze the correlated motions with a mutual information entropy quantity, and examine state transition rates in a Markov-state model, to give coarse-grained descriptors of the motions. Our MD simulations show that while there are many strongly correlated motions, antibodies are highly flexible, with F(ab) and F(c) domains constantly forming and breaking contacts, both polar and non-polar. We find that salt bridges break and reform, and not always with the same partners. While the MD simulations in explicit water give the right time scales for the motions, the simulated motions are about 3-fold faster than the experiments. Overall, the picture that emerges is that antibodies do not simply fluctuate around a single state of atomic contacts. Rather, in these large molecules, different atoms come in contact during different motions.
| S-EPMC3893241 | biostudies-literature
Project description:Adenosine deaminase (ADA) catalyzes the deamination of adenosine, which is important in purine metabolism. ADA is ubiquitous to almost all human tissues, and ADA abnormalities have been reported in various diseases, including rheumatoid arthritis. ADA can be divided into two conformations based on the inhibitor that it binds to: open and closed forms. Here, we chose three ligands, namely, FR117016 (FR0), FR221647 (FR2) (open form), and HDPR (PRH, closed form), to investigate the inhibition mechanism of ADA and its effect on ADA through molecular dynamics simulations. In open forms, Egap and electrostatic potential (ESP) indicated that electron transfer might occur more easily in FR0 than in FR2. Binding free energy and hydrogen bond occupation revealed that the ADA-FR0 complex had a more stable structure than ADA-FR2. The probability of residues Pro159 to Lys171 of ADA-FR0 and ADA-FR2 to form a helix moderately increased compared with that in nonligated ADA. In comparison with FR0 and FR2 PRH could maintain ADA in a closed form to inhibit the function of ADA. The α7 helix (residues Thr57 to Ala73) of ADA in the closed form was mostly unfastened because of the effect of PRH. The number of H bonds and the relative superiority of the binding free energy indicated that the binding strength of PRH to ADA was significantly lower than that of an open inhibitor, thereby supporting the comparison of the inhibitory activities of the three ligands. Alanine scanning results showed that His17, Gly184, Asp295, and Asp296 exerted the greatest effects on protein energy, suggesting that they played crucial roles in binding to inhibitors. This study served as a theoretical basis for the development of new ADA inhibitors.
| S-EPMC6239307 | biostudies-literature
Project description:Rab5 is a small GTPase and a key regulator in early endosomal trafficking. Rab5 and its effectors are involved in a large number of infectious diseases and certain types of cancer. We performed µs atomistic molecular dynamics simulations of inactive and active full-length Rab5 anchored to a complex model bilayer with composition of the early endosome membrane. Direct interactions between the Rab5 G domain and the bilayer were observed. We found two dominant nucleotide-dependent orientations characterised by a different accessibility of the switch regions. The "buried switch" orientation was mainly associated with inactive Rab5 accompanied with a rather extended structure of the hypervariable C-terminal region. Active Rab5 preferred an orientation in which the switch regions are accessible to effector proteins. These structural differences may provide an opportunity to selectively target one Rab5 state and lead to new approaches in the development of Rab5-specific therapies.
| S-EPMC6010109 | biostudies-literature
Project description:Guanosine and related derivatives self-assemble in the presence of cations like potassium into supramolecular G-quadruplexes (SGQs), where four guanine moieties form planar tetrads (T) that coaxially stack into columnar aggregates with broad size distributions. However, SGQs made from 8-aryl-2'-deoxyguanosine derivatives (8ArGs), form mostly octamers, or two-tetrad (2T)-SGQs, while some form dodecamers (3T-SGQs), or hexadecamers (4T-SGQs), and none reported to date form higher assemblies. A theoretical model that addresses the configurational space available for the multiple pathways available for 8ArGs to self-assemble into SGQs is used to frame a series of molecular dynamics simulations (MDS) with selected SGQs. Some key insights from this work include: (a) The predicted entropic costs are not significantly higher for SGQs with more subunits due to their hierarchical assembly pathways; (b) The multiple isomeric SGQs vary in the interfacial contacts between consecutive tetrads, due to their two distinct sides (head, h; tail, t), with the MDS supporting the predicted order of stability of hh > ht > tt for octamers. (c) Such order also applies to dodecamers and hexadecamers, but with context-dependent exceptions due to strong allosteric effects. (d) The main factor disfavoring the tt interface is the repulsive dipolar interactions between the O4' from ribose moieties on adjacent tetrads. (e) SGQs with 5 or more tetrads are disfavored because the attractive interactions are not large or strong enough to overcome the many repulsive forces resulting from the addition of further tetrads. We expect these findings provide some guidelines to enable the further development of SGQs into functional materials.
| S-EPMC5649357 | biostudies-literature