A method to construct the dynamic landscape of a bio-membrane with experiment and simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Smith AA
PROVIDER: S-EPMC8748619 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Conformational changes in the envelope proteins of flaviviruses help to expose the highly conserved fusion peptide (FP), a region which is critical to membrane fusion and host cell infection, and which represents a significant target for antiviral drugs and antibodies. In principle, extended timescale atomic-resolution simulations may be used to characterize the dynamics of such peptides. However, the resultant accuracy is critically dependent upon both the underlying force field and sufficient conformational sampling. In the present study, we report a comprehensive comparison of three simulation methods and four force fields comprising a total of more than 40??s of sampling. Additionally, we describe the conformational landscape of the FP fold across all flavivirus family members. All investigated methods sampled conformations close to available X-ray structures, but exhibited differently populated ensembles. The best force field / sampling combination was sufficiently accurate to predict that the solvated peptide fold is less ordered than in the crystallographic state, which was subsequently confirmed via circular dichroism and spectrofluorometric measurements. Finally, the conformational landscape of a mutant incapable of membrane fusion was significantly shallower than wild-type variants, suggesting that dynamics should be considered when therapeutically targeting FP epitopes.
| S-EPMC4726195 | biostudies-literature
Project description:Experimental data from protein engineering studies and NMR spectroscopy have been used by theoreticians to develop algorithms for helix propensity and to benchmark computer simulations of folding pathways and energy landscapes. Molecular dynamic simulations of the unfolding of chymotrypsin inhibitor 2 (CI2) have provided detailed structural models of the transition state ensemble for unfolding/folding of the protein. We now have used the simulated transition state structures to design faster folding mutants of CI2. The models pinpoint a number of unfavorable local interactions at the carboxyl terminus of the single alpha-helix and in the protease-binding loop region of CI2. By removing these interactions or replacing them with stabilizing ones, we have increased the rate of folding of the protein up to 40-fold (tau = 0.4 ms). This correspondence, and other examples of agreement between experiment and theory in general, Phi-values and molecular dynamics simulations, in particular, suggest that significant progress has been made toward describing complete folding pathways at atomic resolution by combining experiment and simulation.
| S-EPMC21100 | biostudies-literature
Project description:The most widely-used method for detecting genome-wide protein-DNA interactions is chromatin immunoprecipitation on tiling microarrays, commonly known as ChIP-chip. Here, we conducted the first objective analysis of tiling array platforms and analysis algorithms in a simulated ChIP-chip experiment. Mixtures of human genomic DNA and "spike-ins" comprised of nearly 100 human sequences at various concentrations were hybridized to four tiling array platforms by eight independent groups. Blind to the number of spike-ins, their locations, and the range of concentrations, each group made predictions of the spike-in locations. All commercial tiling array platforms performed well, although each platform and analysis algorithm had distinct performance and cost characteristics. Simple sequence repeats and genome redundancy tend to result in false positives on oligonucleotide platforms. The spike-in DNA samples and the resulting array data presented here provide a stable benchmark against which future ChIP platforms, protocol improvements, and analysis methods can be evaluated. Keywords: chip-ChIP simulation For data usage terms and conditions, please refer to http://www.genome.gov/27528022 and http://www.genome.gov/Pages/Research/ENCODE/ENCODEDataReleasePolicyFinal2008.pdf
2007-12-01 | GSE9732 | GEO
Project description:Membrane-integral pyrophosphatases (mPPases) couple the hydrolysis of pyrophosphate (PPi) to the pumping of Na+, H+, or both these ions across a membrane. Recently solved structures of the Na+-pumping Thermotoga maritima mPPase (TmPPase) and H+-pumping Vigna radiata mPPase revealed the basis of ion selectivity between these enzymes and provided evidence for the mechanisms of substrate hydrolysis and ion-pumping. Our atomistic molecular dynamics (MD) simulations of TmPPase demonstrate that loop 5-6 is mobile in the absence of the substrate or substrate-analogue bound to the active site, explaining the lack of electron density for this loop in resting state structures. Furthermore, creating an apo model of TmPPase by removing ligands from the TmPPase:IDP:Na structure in MD simulations resulted in increased dynamics in loop 5-6, which results in this loop moving to uncover the active site, suggesting that interactions between loop 5-6 and the imidodiphosphate and its associated Mg2+ are important for holding a loop-closed conformation. We also provide further evidence for the transport-before-hydrolysis mechanism by showing that the non-hydrolyzable substrate analogue, methylene diphosphonate, induces low levels of proton pumping by VrPPase.
| S-EPMC5336470 | biostudies-literature
Project description:BackgroundMany computer studies have employed either dynamic simulation or metabolic flux analysis (MFA) to predict the behaviour of biochemical pathways. Dynamic simulation determines the time evolution of pathway properties in response to environmental changes, whereas MFA provides only a snapshot of pathway properties within a particular set of environmental conditions. However, owing to the large amount of kinetic data required for dynamic simulation, MFA, which requires less information, has been used to manipulate large-scale pathways to determine metabolic outcomes.ResultsHere we describe a simulation method based on cooperation between kinetics-based dynamic models and MFA-based static models. This hybrid method enables quasi-dynamic simulations of large-scale metabolic pathways, while drastically reducing the number of kinetics assays needed for dynamic simulations. The dynamic behaviour of metabolic pathways predicted by our method is almost identical to that determined by dynamic kinetic simulation.ConclusionThe discrepancies between the dynamic and the hybrid models were sufficiently small to prove that an MFA-based static module is capable of performing dynamic simulations as accurately as kinetic models. Our hybrid method reduces the number of biochemical experiments required for dynamic models of large-scale metabolic pathways by replacing suitable enzyme reactions with a static module.
| S-EPMC1262783 | biostudies-literature
Project description:In this work, the isomer mixture of 4,4'-diphenylmethane diisocyanate (MDI) and 2,4'-MDI was separated and purified by dynamic falling film melt crystallization, and 99.3% purity and 50.8% yield of 4,4'-MDI could be obtained under optimized conditions. The separation mechanism was simulated by density functional theory (DFT) and molecular dynamics (MD) simulation. Results showed that compared with 2,4'-MDI, 4,4'-MDI molecules could form a more stable and symmetrical crystal structure due to their stronger charge density symmetry and electrostatic potential energy. Furthermore, the separation phenomenon and the formation of the crystal structure were observed according to the radial distribution function (RDF) and orientation correlation function obtained from MD simulation. Finally, the attachment energy (AE) model was used to observe and compare different crystal surfaces; it was proposed that the aggregation of 4,4'-MDI was attributed to the polar attraction between isocyanate groups according to the results of the orientation correlation function. It was also observed that compared with 2,4'-MDI, 4,4'-MDI molecules on the (110) crystal surface were easier to form crystal structures.
| S-EPMC9244906 | biostudies-literature
Project description:Cellular transport machinery, such as channels and pumps, is working against the background of unassisted material transport through membranes. The permeation of a blocked tryptophan through a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) membrane is investigated to probe unassisted or physical transport. The transport rate is measured experimentally and modeled computationally. The time scale measured by parallel artificial membrane permeation assay (PAMPA) experiments is ~8 h. Simulations with the milestoning algorithm suggest mean first passage time (MFPT) of ~4 h and the presence of a large barrier at the center of the bilayer. A similar calculation with the solubility-diffusion model yields a MFPT of ~15 min. This permeation rate is 9 orders of magnitude slower than the permeation rate of only a tryptophan side chain (computed by us and others). This difference suggests critical dependence of transport time on permeant size and hydrophilicity. Analysis of the simulation results suggests that the permeant partially preserves hydrogen bonding of the peptide backbone to water and lipid molecules even when it is moving closer to the bilayer center. As a consequence, defects of the membrane structure are developed to assist permeation.
| S-EPMC3302722 | biostudies-literature
Project description:In this study, surfactants were used to enhance heavy oil-solid separation, and a detailed mechanism was explored by SARA (saturates, aromatics, resins, asphaltenes) analysis, element analysis, AFM measurement, and molecular dynamic simulation. Surfactants could effectively decrease oil/solid interaction force and then oil-solid separation would be enhanced. The oil-solid interactive force was in relation to surfactants concentration, pH value, asphaltene content, and salinity. The molecular dynamics simulation results show that the dissociation of saturated hydrocarbon, aromatic hydrocarbon, resin, and asphaltene (SARA) on carbonate minerals is gradually weakened for all surfactants. In the process of molecular dynamics simulation of surfactant stripping SARA, firstly, the surfactant molecules adsorb on the surface of SARA molecules. After that, the surfactant peels SARA molecules off the surface of calcite under the influence of molecular thermal motion. In this process, surfactant molecules will not be directly adsorbed on the surface of trace minerals. The results of energy/temperature balance indicated that saturates, aromatics and resins could remain stable when the molecular dynamics simulation time reached 2000 ps with the phenomenon that saturates, aromatics could liberate from minerals totally within 2000 ps. The molecular dynamics simulation of asphaltenes will not liberate from calcite surface within 6000 ps, meanwhile, they could not reach the energy balance/energy balance within 6000 ps. The functional groups of surfactant molecules would have interactions with the SARA functional group, resulting in different dissociation effects of SARA. The results of molecular dynamics simulation are consistent with the experiment results. The separation effect of saturated hydrocarbon, aromatic hydrocarbon, resin, and asphaltene in five kinds of surfactants were different. The molecular dynamic simulation results were in accordance with the SARA analysis.
| S-EPMC8308391 | biostudies-literature
Project description:The lack of a universal simulation method for triboelectric nanogenerator (TENG) makes the device design and optimization difficult before experiment, which protracts the research and development process and hinders the landing of practical TENG applications. The existing electrostatic induction models for TENGs have limitations in simulating TENGs with complex geometries and their dynamic behaviors under practical movements due to the topology change issues. Here, a dynamic finite element method (FEM) model is proposed. The introduction of air buffer layers and the moving mesh method eliminates the topology change issues during practical movement and allows simulation of dynamic and time-varying behaviors of TENGs with complex 2D/3D geometries. Systematic investigations are carried out to optimize the air buffer thickness and mesh densities, and the optimized results show excellent consistency with the experimental data and results based on other existing methods. It also shows that a 3D disk-type rotating TENG can be simulated using the model, clearly demonstrating the capability and superiority of the dynamic FEM model. Moreover, the dynamic FEM model is used to optimize the shape of the tribo-material, which is used as a preliminary example to demonstrate the possibility of designing a TENG-based sensor.
| S-EPMC7506595 | biostudies-literature
Project description:The most widely-used method for detecting genome-wide protein-DNA interactions is chromatin immunoprecipitation on tiling microarrays, commonly known as ChIP-chip. Here, we conducted the first objective analysis of tiling array platforms and analysis algorithms in a simulated ChIP-chip experiment. Mixtures of human genomic DNA and "spike-ins" comprised of nearly 100 human sequences at various concentrations were hybridized to four tiling array platforms by eight independent groups. Blind to the number of spike-ins, their locations, and the range of concentrations, each group made predictions of the spike-in locations. All commercial tiling array platforms performed well, although each platform and analysis algorithm had distinct performance and cost characteristics. Simple sequence repeats and genome redundancy tend to result in false positives on oligonucleotide platforms. The spike-in DNA samples and the resulting array data presented here provide a stable benchmark against which future ChIP platforms, protocol improvements, and analysis methods can be evaluated. Keywords: chip-ChIP simulation For data usage terms and conditions, please refer to http://www.genome.gov/27528022 and http://www.genome.gov/Pages/Research/ENCODE/ENCODEDataReleasePolicyFinal2008.pdf For each replicate array, we labeled 1µg of spike-in sample and 1µg of control sample. The Bioprime Plus Array CGH labeling Module (Invitrogen Catolog number 18095-014) was used. The spike-in was labeled with the red channel dye (Alexa Fluor-647) and the control was labeled with the green channel dye (Alexa Fluor-555) in two of the replicates. A dye swap was performed for the third replicate. These samples were then competitively hybridized to three replicate 244k arrays from Agilent, (AMADID 25150451). We hybridized the samples to the arrays using the Agilent Array CGH protocol, with some modifications. Briefly, labeled DNA was combined with Cot-1 DNA, blocking reagent, and hybridization buffer. The hybridizations were carried out at 60°C. The arrays were scanned using an Agilent dual laser scanner, and the images were processed using Agilent Feature Extraction Software.
2010-05-22 | E-GEOD-9732 | biostudies-arrayexpress