Browse
Submit Data
Databases
API
Help

Dataset Information

32 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Accurate, conformation-dependent predictions of solvent effects on protein ionization constants.

ABSTRACT: Predicting how aqueous solvent modulates the conformational transitions and influences the pKa values that regulate the biological functions of biomolecules remains an unsolved challenge. To address this problem, we developed FDPB_MF, a rotamer repacking method that exhaustively samples side chain conformational space and rigorously calculates multibody protein-solvent interactions. FDPB_MF predicts the effects on pKa values of various solvent exposures, large ionic strength variations, strong energetic couplings, structural reorganizations and sequence mutations. The method achieves high accuracy, with root mean square deviations within 0.3 pH unit of the experimental values measured for turkey ovomucoid third domain, hen lysozyme, Bacillus circulans xylanase, and human and Escherichia coli thioredoxins. FDPB_MF provides a faithful, quantitative assessment of electrostatic interactions in biological macromolecules.

SUBMITTER: Barth P

PROVIDER: S-EPMC1829236 | biostudies-literature | 2007 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Accurate, conformation-dependent predictions of solvent effects on protein ionization constants.

Barth P P Alber T T Harbury P B PB

Proceedings of the National Academy of Sciences of the United States of America 20070314 12

Predicting how aqueous solvent modulates the conformational transitions and influences the pKa values that regulate the biological functions of biomolecules remains an unsolved challenge. To address this problem, we developed FDPB_MF, a rotamer repacking method that exhaustively samples side chain conformational space and rigorously calculates multibody protein-solvent interactions. FDPB_MF predicts the effects on pKa values of various solvent exposures, large ionic strength variations, strong e ...[more]

PMID: 17360348

Similar Datasets

Ionization constants pKa of cardiolipin.

Project description:Cardiolipin is a phospholipid found in the inner mitochondrial membrane and in bacteria, and it is associated with many physiological functions. Cardiolipin has a dimeric structure consisting of two phosphatidyl residues connected by a glycerol bridge and four acyl chains, and therefore it can carry two negative charges. The pKa values of the phosphate groups have previously been reported to differ widely with pKa1 = 2.8 and pKa2 = 7.5-9.5. Still, there are several examples of experimental observations from cardiolipin-containing systems that do not fit with this dissociation behavior. Therefore, we have carried out pH-titration and titration calorimetric experiments on two synthetic cardiolipins, 1,1',2,2'-tetradecanoyl cardiolipin, CL (C14:0), and 1,1',2,2'-tetraoctadecenoyl cardiolipin, CL (C18:1). Our results show that both behave as strong dibasic acids with pKa1 about the same as the first pKa of phosphoric acid, 2.15, and pKa2 about one unit larger. The characterization of the acidic properties of cardiolipin is crucial for the understanding of the molecular organization in self-assembled systems that contain cardiolipin, and for their biological function.

| S-EPMC3772843 | biostudies-literature

Inductive effects in amino acids and peptides: Ionization constants and tryptophan fluorescence.

Project description:Although inductive effects in organic compounds are known to influence chemical properties such as ionization constants, their specific contribution to the properties/behavior of amino acids and functional groups in peptides remains largely unexplored. In this study we developed a computationally economical algorithm for ab initio calculation of the magnitude of inductive effects for non-aromatic molecules. The value obtained by the algorithm is called the Inductive Index and we observed a high correlation (R2 = 0.9427) between our calculations and the pKa values of the alpha-amino groups of amino acids with non-aromatic side-chains. Using a series of modified amino acids, we also found similarly high correlations (R2 > 0.9600) between Inductive Indexes and two wholly independent chemical properties: i) the pKa values of ionizable side-chains and, ii) the fluorescence response of the indole group of tryptophan. After assessing the applicability of the method of calculation at the amino acid level, we extended our study to tryptophan-containing peptides and established that inductive contributions of neighboring side-chains are transmitted through peptide bonds. We discuss possible contributions to the study of proteins.

| S-EPMC7498751 | biostudies-literature

Solvent-dependent conductance decay constants in single cluster junctions.

Project description:Single-molecule conductance measurements have focused primarily on organic molecular systems. Here, we carry out scanning tunneling microscope-based break-junction measurements on a series of metal chalcogenide Co6Se8 clusters capped with conducting ligands of varying lengths. We compare these measurements with those of individual free ligands and find that the conductance of these clusters and the free ligands have different decay constants with increasing ligand length. We also show, through measurements in two different solvents, 1-bromonaphthalene and 1,2,4-trichlorobenzene, that the conductance decay of the clusters depends on the solvent environment. We discuss several mechanisms to explain our observations.

| S-EPMC5477014 | biostudies-other

Accurate Single-Molecule Indicator of Solvent Effects.

Project description:The study of the microscopic structure of solvents is of significant importance for deciphering the essential solvation in chemical reactions and biological processes. Yet conventional technologies, such as neutron diffraction, have an inherent averaging effect as they analyze a group of molecules. In this study, we report a method to analyze the microstructure and interaction in solvents from a single-molecule perspective. A single-molecule electrical nanocircuit is used to directly analyze the dynamic microscopic structure of solvents. Through a single-molecule model reaction, the heterogeneity or homogeneity of solvents is precisely detected at the molecular level. Both the thermodynamics and the kinetics of the model reaction demonstrate the microscopic heterogeneity of alcohol-water and alcohol-n-hexane solutions and the microscopic homogeneity of alcohol-carbon tetrachloride solutions. In addition, a real-time event spectroscopy has been developed to study the dynamic characteristics of the segregated phase and the internal intermolecular interaction in microheterogeneous solvents. The development of such a unique high-resolution indicator with single-molecule and single-event accuracy provides infinite opportunities to decipher solvent effects in-depth and optimizes chemical reactions and biological processes in solution.

| S-EPMC8715489 | biostudies-literature

Accurate protein stability predictions from homology models.

Project description:Calculating changes in protein stability (ΔΔG) has been shown to be central for predicting the consequences of single amino acid substitutions in protein engineering as well as interpretation of genomic variants for disease risk. Structure-based calculations are considered most accurate, however the tools used to calculate ΔΔGs have been developed on experimentally resolved structures. Extending those calculations to homology models based on related proteins would greatly extend their applicability as large parts of e.g. the human proteome are not structurally resolved. In this study we aim to investigate the accuracy of ΔΔG values predicted on homology models compared to crystal structures. Specifically, we identified four proteins with a large number of experimentally tested ΔΔGs and templates for homology modeling across a broad range of sequence identities, and selected three methods for ΔΔG calculations to test. We find that ΔΔG-values predicted from homology models compare equally well to experimental ΔΔGs as those predicted on experimentally established crystal structures, as long as the sequence identity of the model template to the target protein is at least 40%. In particular, the Rosetta cartesian_ddg protocol is robust against the small perturbations in the structure which homology modeling introduces. In an independent assessment, we observe a similar trend when using ΔΔGs to categorize variants as low or wild-type-like abundance. Overall, our results show that stability calculations performed on homology models can substitute for those on crystal structures with acceptable accuracy as long as the model is built on a template with sequence identity of at least 40% to the target protein.

| S-EPMC9729920 | biostudies-literature

A fast and accurate computational approach to protein ionization.

Project description:We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.

| S-EPMC2578799 | biostudies-literature

Protein-protein binding pathways and calculations of rate constants using fully-continuous, explicit-solvent simulations.

Project description:A grand challenge in the field of biophysics has been the complete characterization of protein-protein binding processes at atomic resolution. This characterization requires the direct simulation of binding pathways starting from the initial, unbound state and proceeding through states that are too transient to be captured by experiment. Here, we applied the weighted ensemble path sampling strategy to orchestrate atomistic simulation of protein-protein binding pathways. Our simulation generated 203 fully-continuous and independent pathways along with rate constants for the binding process involving the barnase and barstar proteins. Results reveal multiple binding pathways along a "funnel-like" free energy landscape in which the formation of the "encounter complex" intermediate is rate-limiting followed by a relatively rapid rearrangement of the encounter complex to the bound state. Among all diffusional collisions, only ∼11% were productive. In the most probable binding pathways, the proteins rotated to a large extent (likely via electrostatic steering) in order to collide productively followed by "rolling" of the proteins along each other's binding interfaces to reach the bound state. Consistent with experiment, R59 was identified as the most kinetically important barnase residue for the binding process. Furthermore, protein desolvation occurs late in the binding process during the rearrangement of the encounter complex to the bound state. Notably, the positions of crystallographic water molecules that bridge hydrogen bonds between barnase and barstar are occupied in the bound-state ensemble. Our simulation was completed in a month using 1600 CPU cores at a time, demonstrating that it is now practical to carry out atomistic simulations of protein-protein binding.

| S-EPMC6385678 | biostudies-literature

A chemogenomic analysis of ionization constants--implications for drug discovery.

Project description:Chemogenomics methods seek to characterize the interaction between drugs and biological systems and are an important guide for the selection of screening compounds. The acid/base character of drugs has a profound influence on their affinity for the receptor, on their absorption, distribution, metabolism, excretion and toxicity (ADMET) profile and the way the drug can be formulated. In particular, the charge state of a molecule greatly influences its lipophilicity and biopharmaceutical characteristics. This study investigates the acid/base profile of human small-molecule drugs, chemogenomics datasets and screening compounds including a natural products set. We estimate the acid-ionization constant (pK(a)) values of these compounds and determine the identity of the ionizable functional groups in each set. We find substantial differences in acid/base profiles of the chemogenomic classes. In many cases, these differences can be linked to the nature of the target binding site and the corresponding functional groups needed for recognition of the ligand. Clear differences are also observed between the acid/base characteristics of drugs and screening compounds. For example, the proportion of drugs containing a carboxylic acid was 20 %, in stark contrast to a value of 2.4 % for the screening set sample. The proportion of aliphatic amines was 27 % for drugs and only 3.4 % for screening compounds. This suggests that there is a mismatch between commercially available screening compounds and the compounds that are likely to interact with a given chemogenomic target family. Our analysis provides a guide for the selection of screening compounds to better target specific chemogenomic families with regard to the overall balance of acids, bases and pK(a) distributions.

| S-EPMC3777741 | biostudies-literature

Solvent interactions determine carbohydrate conformation.

Project description:The relationship between the three-dimensional structures of oligosaccharides and polysaccharides and their biological properties has been the focus of many recent studies. The overall conformation of an oligosaccharide depends primarily on the orientation of the torsion angles (phi, psi, and omega) between glycosyl residues. Numerous experimental studies have shown that in glucopyranosides the omega-torsion angle (O(6)-C(6)-C(5)-O(5)) displays a preference for gauche orientations, in disagreement with predictions based on gas-phase quantum mechanics calculations. In contrast, the omega-angle in galactopyranosides displays a high proportion of the anti-orientation. For oligosaccharides containing glycosidic linkages at the 6-position (1-->6 linked), variations in rotamer population have a direct effect on the oligosaccharides' structure and function, and yet the physical origin of these conformational preferences remains unclear. Although it is generally recognized that the gauche effect in carbohydrates is a solvent-dependent phenomenon, the mechanism through which solvent induces the gauche preference is not understood. In the present work, quantum mechanics and solvated molecular dynamics calculations were performed on two representative carbohydrates, methyl alpha-D-glucopyranoside and methyl alpha-D-galactopyranoside. We show that correct reproduction of the experimental rotamer distributions about the omega-angles is obtained only after explicit water is included in the molecular dynamics simulations. The primary role of the water appears to be to disrupt the hydrogen bonding within the carbohydrate, thereby allowing the rotamer populations to be determined by internal electronic and steric repulsions between the oxygen atoms. The results reported here provide a quantitative explanation of the conformational behavior of (1-->6)-linked carbohydrates.

| S-EPMC58501 | biostudies-literature

Determination of ionization constants by paper electrophoresis.

Project description:Dimensionless apparent ionization constants of charged low-molecular-weight species may be obtained from paper-electrophoretic data at 20-25 degrees C with buffers (I0.1-0.5) of measured pH (1.5-12.5) containing oxalate ions. Relative mobilities rather than absolute mobilities were measured by using glycerol and m-nitrobenzenesulphonate respectively as standards of zero and unit mobility. Application of the procedure to ionizations of adenine, adenosine, 2'-deoxyadenosine, 3'-deoxyadenosine, 3':5'-cyclic AMP, ADP, ADP-glucose-agrocin 84 and ATP is described.

| S-EPMC1162905 | biostudies-other

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data