Browse
Submit Data
Databases
API
Help

Dataset Information

19 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Absolute protein-protein association rate constants from flexible, coarse-grained Brownian dynamics simulations: the role of intermolecular hydrodynamic interactions in barnase-barstar association.

ABSTRACT: Theory and computation have long been used to rationalize the experimental association rate constants of protein-protein complexes, and Brownian dynamics (BD) simulations, in particular, have been successful in reproducing the relative rate constants of wild-type and mutant protein pairs. Missing from previous BD studies of association kinetics, however, has been the description of hydrodynamic interactions (HIs) between, and within, the diffusing proteins. Here we address this issue by rigorously including HIs in BD simulations of the barnase-barstar association reaction. We first show that even very simplified representations of the proteins--involving approximately one pseudoatom for every three residues in the protein--can provide excellent reproduction of the absolute association rate constants of wild-type and mutant protein pairs. We then show that simulations that include intermolecular HIs also produce excellent estimates of association rate constants, but, for a given reaction criterion, yield values that are decreased by ∼35-80% relative to those obtained in the absence of intermolecular HIs. The neglect of intermolecular HIs in previous BD simulation studies, therefore, is likely to have contributed to the somewhat overestimated absolute rate constants previously obtained. Consequently, intermolecular HIs could be an important component to include in accurate modeling of the kinetics of macromolecular association events.

SUBMITTER: Frembgen-Kesner T

PROVIDER: S-EPMC2965997 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Diffusional encounter of barnase and barstar.

Project description:We present an analysis of trajectories from Brownian dynamics simulations of diffusional protein-protein encounter for the well-studied system of barnase and barstar. This analysis reveals details about the optimal association pathways, the regions of the encounter complex, possible differences of the pathways for dissociation and association, the coupling of translational and rotation motion, and the effect of mutations on the trajectories. We found that a small free-energy barrier divides the energetically most favorable region into a region of the encounter complex above the barnase binding interface and a region around a second energy minimum near the RNA binding loop. When entering the region of the encounter complex from the region near the RNA binding loop, barstar has to change its orientation to increase the electrostatic attraction between the proteins. By concentrating the analysis on the successful binding trajectories, we found that the region of the second minimum is not essential for the binding of barstar to barnase. Nevertheless, this region may be helpful to steer barstar into the region of the encounter complex. When applying the same analysis to several barnase mutants, we found that single mutations may drastically change the free-energy landscape and may significantly alter the population of the two minima. Therefore, certain protein-protein pairs may require careful adaptation of the positions of encounter and transition states when interpreting mutation effects on kinetic rates of association and/or dissociation.

| S-EPMC1386772 | biostudies-literature

Coarse-Grained Brownian Dynamics Simulations of the 10-23 DNAzyme.

Project description:Deoxyribozymes (DNAzymes) are single-stranded DNA that catalyze nucleic acid biochemistry. Although a number of DNAzymes have been discovered by in vitro selection, the relationship between their tertiary structure and function remains unknown. We focus here on the well-studied 10-23 DNAzyme, which cleaves mRNA with a catalytic efficiency approaching that of RNase A. Using coarse-grained Brownian dynamics simulations, we find that the DNAzyme bends its substrate away from the cleavage point, exposing the reactive site and buckling the DNAzyme catalytic core. This hypothesized transition state provides microscopic insights into experimental observations concerning the size of the DNAzyme/substrate complex, the impact of the recognition arm length, and the sensitivity of the enzymatic activity to point mutations of the catalytic core. Upon cleaving the pertinent backbone bond in the substrate, we find that the catalytic core of the DNAzyme unwinds and the overall complex rapidly extends, in agreement with experiments on the related 8-17 DNAzyme. The results presented here provide a starting point for interpreting experimental data on DNAzyme kinetics, as well as developing more detailed simulation models. The results also demonstrate the limitations of using a simple physical model to understand the role of point mutations.

| S-EPMC2776252 | biostudies-literature

How optimal are the binding energetics of barnase and barstar?

Project description:The extracellular ribonuclease barnase and its intracellular inhibitor barstar bind fast and with high affinity. Although extensive experimental and theoretical studies have been carried out on this system, it is unclear what the relative importance of different contributions to the high affinity is and whether binding can be improved through point mutations. In this work, we first applied Poisson-Boltzmann electrostatic calculations to 65 barnase-barstar complexes with mutations in both barnase and barstar. The continuum electrostatic calculations with a van der Waals surface dielectric boundary definition result in the electrostatic interaction free energy providing the dominant contribution favoring barnase-barstar binding. The results show that the computed electrostatic binding free energy can be improved through mutations at W44/barstar and E73/barnase. Furthermore, the determinants of binding affinity were quantified by applying COMparative BINding Energy (COMBINE) analysis to derive quantitative structure-activity relationships (QSARs) for the 65 complexes. The COMBINE QSAR model highlights approximately 20 interfacial residue pairs as responsible for most of the differences in binding affinity between the mutant complexes, mainly due to electrostatic interactions. Based on the COMBINE model, together with Brownian dynamics simulations to compute diffusional association rate constants, several mutants were designed to have higher binding affinities than the wild-type proteins.

| S-EPMC1304567 | biostudies-literature

CID fragmentation, H/D exchange and supermetallization of Barnase-Barstar complex.

Project description:The barnase-barstar complex is one of the most stable protein-protein complexes and has a very wide range of possible applications. Here we report the use of top-down mass spectrometry for the investigation of the structure of this complex, its ionization via ESI, isolation and fragmentation. It was found that the asymmetry of the resulting charge state distributions of the protein monomer product ions increased as the charge state of the precursor ions increased. For the investigation of the 3D structure of the complex, the gas phase H/D exchange reaction was used. In addition, supermetallized ions of the complex with Zn were produced and investigated. It was observed that an increase in the number of metals bound to the complex results in a change in complex stability and the charge distribution between protein fragment. Analysis of the fragmentation pattern of the supermetallized complex [bn-b*?+?5Zn]10+ indicated that this ion is present in different conformations with different charges and Zn distributions. Since Zn cannot migrate, such structures must be formed during ionization.

| S-EPMC5522418 | biostudies-literature

Key intermolecular interactions in the E. coli 70S ribosome revealed by coarse-grained analysis.

Project description:The ribosome is a very large complex that consists of many RNA and protein molecules and plays a central role in protein biosynthesis in all organisms. Extensive interactions between different molecules are critical to ribosomal functional dynamics. In this work, intermolecular interactions in the Escherichia coli 70S ribosome are investigated by coarse-grained (CG) analysis. CG models are defined to preserve dynamic domains in RNAs and proteins and to capture functional motions in the ribosome, and then the CG sites are connected by harmonic springs, and spring constants are obtained by matching the computed fluctuations to those of an all-atom molecular dynamics (MD) simulation. Those spring constants indicate how strong the interactions are between the ribosomal components, and they are in good agreement with various experimental data. Nearly all the bridges between the small and large ribosomal subunits are indicated by CG interactions with large spring constants. The head of the small subunit is very mobile because it has minimal CG interactions with the rest of the subunit; however, a large number of small subunit proteins bind to maintain the internal structure of the head. The results show a clear connection between the intermolecular interactions and the structural and functional properties of the ribosome because of the reduced complexity in domain-based CG models. The present approach also provides a useful strategy to map interactions between molecules within large biomolecular complexes since it is not straightforward to investigate these by either atomistic MD simulations or residue-based elastic network models.

| S-EPMC3641354 | biostudies-literature

Gated binding of ligands to HIV-1 protease: Brownian dynamics simulations in a coarse-grained model.

Project description:The internal motions of proteins may serve as a "gate" in some systems, which controls ligand-protein association. This study applies Brownian dynamics simulations in a coarse-grained model to study the gated association rate constants of HIV-1 proteases and drugs. The computed gated association rate constants of three protease mutants, G48V/V82A/I84V/L90M, G48V, and L90M with three drugs, amprenavir, indinavir, and saquinavir, yield good agreements with experiments. The work shows that the flap dynamics leads to "slow gating". The simulations suggest that the flap flexibility and the opening frequency of the wild-type, the G48V and L90M mutants are similar, but the flaps of the variant G48V/V82A/I84V/L90M open less frequently, resulting in a lower gated rate constant. The developed methodology is fast and provides an efficient way to predict the gated association rate constants for various protease mutants and ligands.

| S-EPMC1459512 | biostudies-literature

Calculation of therapeutic antibody viscosity with coarse-grained models, hydrodynamic calculations and machine learning-based parameters.

Project description:High viscosity presents a challenge for manufacturing and drug delivery of therapeutic antibodies. The viscosity is determined by protein-protein interactions among many antibodies. Molecular simulation is a promising method to study protein-protein interactions; however, all-atom models do not allow the simulation of multiple molecules, which is necessary to compute viscosity directly. Coarse-grained models, on the other hand can do this. In this work, a 12-bead coarse-grained model based on Swan and coworkers (J. Phys. Chem. B 2018, 122, 2867-2880) was applied to study antibody interactions. Two adjustable parameters related to the short-range interactions on the variable and constant regions were determined by fitting experimental data of 20 IgG1 monoclonal antibodies at 150 mg/mL. The root-mean-square deviation improved from 1 to 0.68, and the correlation coefficient improved from 0.63 to 0.87 compared to that of a previous model that assumed the short-range interactions were the same for all the beads. Our model is also able to calculate the viscosity over a wide range of concentrations without additional parameters. A tabulated viscosity based on our model is provided to facilitate antibody screening in early-stage design.

| S-EPMC8043186 | biostudies-literature

The effect of different force applications on the protein-protein complex Barnase-Barstar.

Project description:Steered molecular dynamics simulations are a tool to examine the energy landscape of protein-protein complexes by applying external forces. Here, we analyze the influence of the velocity and geometry of the probing forces on a protein complex using this tool. With steered molecular dynamics, we probe the stability of the protein-protein complex Barnase-Barstar. The individual proteins are mechanically labile. The Barnase-Barstar binding site is more stable than the folds of the individual proteins. By using different force protocols, we observe a variety of responses of the system to the applied tension.

| S-EPMC2741585 | biostudies-literature

Sliding of proteins non-specifically bound to DNA: Brownian dynamics studies with coarse-grained protein and DNA models.

Project description:DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.

| S-EPMC4263358 | biostudies-literature

A novel strategy for promoting homoplasmic plastid transformant production using the barnase-barstar system.

Project description:Plastid transformants form biofactories that are able to produce extra proteins in plastids when they are in a homoplasmic state. To date, plastid transformation has been reported in about twenty plant species; however, the production of homoplasmic plastid transformants is not always successful or easy. Heteroplasmic plants that contain wild-type plastids produce fewer target proteins and do not always successfully transfer transgenes to progeny. In order to promote the generation of homoplasmic plants, we developed a novel system using barnase-barster to eliminate wild-type plastids from heteroplasmic cells systematically. In this system, a chemically inducible cytotoxic barnase under a plastid transit signal was introduced into nuclear DNA and barster, which inhibits barnase, was integrated into plastid DNA with the primary selection markers aminoglycoside 3'-adenylyltransferase (aadA) and green fluorescence protein (GFP) gene. As expected, the expression of the plastid barnase was lethal to cells as seen in leaf segments, but barster expression in plastids rescued them. We then investigated the regeneration frequency of homoplasmic shoots from heteroplasmic leaf segments with or without barnase expression. The regeneration frequency of homoplasmic-like shoots expressing barnase-barster system was higher than that of shoots not expressing this. We expect that the application of this novel strategy for transformation of plastids will be supportive to generate homoplasmic plastid transformants in other plant species.

| S-EPMC7434676 | biostudies-literature

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