Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Steuer J
PROVIDER: S-EPMC8373139 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The guanidine-II riboswitch, also known as mini-ykkC, is a conserved mRNA element with more than 800 examples in bacteria. It consists of two stem-loops capped by identical, conserved tetraloops that are separated by a linker region of variable length and sequence. Like the guanidine-I riboswitch, it controls the expression of guanidine carboxylases and SugE-like genes. The guanidine-II riboswitch specifically binds free guanidinium cations and functions as a translationally controlled on-switch. Here we report the structure of a P2 stem-loop from the Pseudomonas aeruginosa guanidine-II riboswitch aptamer bound to guanidine at 1.57 Å resolution. The hairpins dimerize via the conserved tetraloop, which also contains the binding pocket. Two guanidinium molecules bind near the dimerization interface, one in each tetraloop. The guanidinium cation is engaged in extensive hydrogen bonding to the RNA. Contacts include the Hoogsteen face of a guanine base and three nonbridging phosphate oxygens. Cation-π interactions and ionic interactions also stabilize ligand binding. The guanidine-II riboswitch utilizes the same recognition strategies as the guanidine-I riboswitch while adopting an entirely different and much smaller RNA fold.
| S-EPMC5558903 | biostudies-literature
Project description:Sampling multiple binding modes of a ligand in a single molecular dynamics simulation is difficult. A given ligand may have many internal degrees of freedom, along with many different ways it might orient itself in a binding site or across several binding sites, all of which might be separated by large energy barriers. We have developed a novel Monte Carlo move called molecular darting (MolDarting) to reversibly sample between predefined binding modes of a ligand. Here, we couple this with nonequilibrium candidate Monte Carlo (NCMC) to improve acceptance of moves. We apply this technique to a simple dipeptide system, a ligand binding to T4 lysozyme L99A, and ligand binding to HIV integrase to test this new method. We observe significant increases in acceptance compared to uniformly sampling the internal and rotational/translational degrees of freedom in these systems.
| S-EPMC8121195 | biostudies-literature
Project description:The guanidine-I riboswitch is a conserved RNA element with approximately 2,000 known examples across four phyla of bacteria. It exists upstream of nitrogen metabolism and multidrug resistance transporter genes and alters expression through the specific recognition of a free guanidinium cation. Here we report the structure of a guanidine riboswitch aptamer from Sulfobacillus acidophilus at 2.7 Å resolution. Helices P1, P1a, P1b, and P2 form a coaxial stack that acts as a scaffold for ligand binding. A previously unidentified P3 helix docks into P1a to form the guanidinium binding pocket, which is completely enclosed. Every functional group of the ligand is recognized through hydrogen bonding to guanine bases and phosphate oxygens. Guanidinium binding is further stabilized through cation-? interactions with guanine bases. This allows the riboswitch to recognize guanidinium while excluding other bacterial metabolites with a guanidino group, including the amino acid arginine.
| S-EPMC5317007 | biostudies-literature
Project description:Elucidating the mechanisms of specific small-molecule (ligand) recognition by proteins is a long-standing conundrum. While the structures of these molecules, proteins and ligands, have been extensively studied, protein-ligand interactions, or binding modes, have not been comprehensively analyzed. Although methods for assessing similarities of binding site structures have been extensively developed, the methods for the computational treatment of binding modes have not been well established. Here, we developed a computational method for encoding the information about binding modes as graphs, and assessing their similarities. An all-against-all comparison of 20,040 protein-ligand complexes provided the landscape of the protein-ligand binding modes and its relationships with protein- and chemical spaces. While similar proteins in the same SCOP Family tend to bind relatively similar ligands with similar binding modes, the correlation between ligand and binding similarities was not very high (R(2) = 0.443). We found many pairs with novel relationships, in which two evolutionally distant proteins recognize dissimilar ligands by similar binding modes (757,474 pairs out of 200,790,780 pairs were categorized into this relationship, in our dataset). In addition, there were an abundance of pairs of homologous proteins binding to similar ligands with different binding modes (68,217 pairs). Our results showed that many interesting relationships between protein-ligand complexes are still hidden in the structure database, and our new method for assessing binding mode similarities is effective to find them.
| S-EPMC5338237 | biostudies-literature
Project description:A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.
| S-EPMC1302794 | biostudies-literature
Project description:Accurate ligand-protein binding affinity prediction, for a set of similar binders, is a major challenge in the lead optimization stage in drug development. In general, docking and scoring functions perform unsatisfactorily in this application. Docking calculations, followed by molecular dynamics simulations and free energy calculations can be applied to improve the predictions. However, for targets with large, flexible binding sites, with no experimentally determined binding modes for a set of ligands, insufficient sampling can decrease the accuracy of the free energy calculations. Cytochrome P450s, a protein family of major importance for drug metabolism, is an example of a challenging target for binding affinity predictions. As a result, the choice of starting structure from the docking solutions becomes crucial. In this study, an iterative scheme is introduced that includes multiple independent molecular dynamics simulations to obtain weighted ensemble averages to be used in the linear interaction energy method. The proposed scheme makes the initial pose selection less crucial for further simulation, as it automatically calculates the relative weights of the various poses. It also properly takes into account the possibility that multiple binding modes contribute similarly to the overall affinity, or of similar compounds occupying very different poses. The method was applied to a set of 12 compounds binding to cytochrome P450 2C9 and it displayed a root mean-square error of 2.9 kJ/mol.
| S-EPMC2877349 | biostudies-literature
Project description:Multisite ligand binding by proteins is commonly utilized in the regulation of biological systems and exploited in a range of biochemical technologies. Aptamers, although widely utilized in many rationally designed biochemical systems, are rarely capable of multisite ligand binding. The cocaine-binding aptamer is often used for studying and developing sensor and aptamer-based technologies. Here, we use isothermal titration calorimetry (ITC) and NMR spectroscopy to demonstrate that the cocaine-binding aptamer switches from one-site to two-site ligand binding, dependent on NaCl concentration. The high-affinity site functions at all buffer conditions studied, the low-affinity site only at low NaCl concentrations. ITC experiments show the two ligand-binding sites operate independently of one another with different affinities and enthalpies. NMR spectroscopy shows the second binding site is located in stem 2 near the three-way junction. This ability to control ligand binding at the second site by adjusting the concentration of NaCl is rare among aptamers and may prove a useful in biotechnology applications. This work also demonstrates that in vitro selected biomolecules can have functions as complex as those found in nature.
| S-EPMC5388413 | biostudies-literature
Project description:Cytochrome P450 (CYP) monoxygenses utilize heme cofactors to catalyze oxidation reactions. They play a critical role in metabolism of many classes of drugs, are an attractive target for drug development, and mediate several prominent drug interactions. Many substrates and inhibitors alter the spin state of the ferric heme by displacing the heme's axial water ligand in the resting enzyme to yield a five-coordinate iron complex, or they replace the axial water to yield a nitrogen-ligated six-coordinate iron complex, which are traditionally assigned by UV-vis spectroscopy. However, crystal structures and recent pulsed electron paramagnetic resonance (EPR) studies find a few cases where molecules hydrogen bond to the axial water. The water-bridged drug-H2O-heme has UV-vis spectra similar to nitrogen-ligated, six-coordinate complexes, but are closer to "reverse type I" complexes described in older liteature. Here, pulsed and continuous wave (CW) EPR demonstrate that water-bridged complexes are remarkably common among a range of nitrogenous drugs or drug fragments that bind to CYP3A4 or CYP2C9. Principal component analysis reveals a distinct clustering of CW EPR spectral parameters for water-bridged complexes. CW EPR reveals heterogeneous mixtures of ligated states, including multiple directly-coordinated complexes and water-bridged complexes. These results suggest that water-bridged complexes are under-represented in CYP structural databases and can have energies similar to other ligation modes. The data indicates that water-bridged binding modes can be identified and distinguished from directly-coordinated binding by CW EPR.
| S-EPMC5921899 | biostudies-literature
Project description:Cytidine triphosphate synthase (CTPS), which comprises an ammonia ligase domain and a glutamine amidotransferase domain, catalyzes the final step of de novo CTP biosynthesis. The activity of CTPS is regulated by the binding of four nucleotides and glutamine. While glutamine serves as an ammonia donor for the ATP-dependent conversion of UTP to CTP, the fourth nucleotide GTP acts as an allosteric activator. Models have been proposed to explain the mechanisms of action at the active site of the ammonia ligase domain and the conformational changes derived by GTP binding. However, actual GTP/ATP/UTP binding modes and relevant conformational changes have not been revealed fully. Here, we report the discovery of binding modes of four nucleotides and a glutamine analog 6-diazo-5-oxo-L-norleucine in Drosophila CTPS by cryo-electron microscopy with near-atomic resolution. Interactions between GTP and surrounding residues indicate that GTP acts to coordinate reactions at both domains by directly blocking ammonia leakage and stabilizing the ammonia tunnel. Additionally, we observe the ATP-dependent UTP phosphorylation intermediate and determine interacting residues at the ammonia ligase. A noncanonical CTP binding at the ATP binding site suggests another layer of feedback inhibition. Our findings not only delineate the structure of CTPS in the presence of all substrates but also complete our understanding of the underlying mechanisms of the allosteric regulation and CTP synthesis.
| S-EPMC8325340 | biostudies-literature
Project description:The adenine riboswitch aptamer, the A box, positively regulates gene expression upon adenine binding. To provide insight into structure-function relationships, important for the adenine riboswitch aptamer, we have created alignments for six aptamer sequences that reveal the core requirements. In addition, 2-aminopurine (2AP) binding studies have been used to test the consensus sequence derived from the alignment. Overall, the consensus secondary structure is consistent with 2AP binding studies. However, a position in the core, previously identified as variable, shows restriction in nucleotide sequence. Furthermore, this restriction is found to be related with the ligand specificity of the riboswitch. The implications of this relationship for the riboswitch gene regulation mechanism are discussed.
| S-EPMC1800517 | biostudies-literature