scholarly journals Computational Study of the Interaction of a PEGylated Hyperbranched Polymer/Doxorubicin Complex with a Bilipid Membrane

Fluids ◽  
2019 ◽  
Vol 4 (1) ◽  
pp. 17 ◽  
Author(s):  
Prodromos Arsenidis ◽  
Kostas Karatasos

Fully atomistic molecular dynamics simulations are employed to study in detail the interactions between a complex comprised by a PEGylated hyperbranched polyester (HBP) and doxorubicin molecules, with a model dipalmitoylphosphatidylglycerol membrane in an aqueous environment. The effects of the presence of the lipid membrane in the drug molecules’ spatial arrangement were examined in detail and the nature of their interaction with the latter were discussed and quantified where possible. It was found that a partial migration of the drug molecules towards the membrane’s surface takes place, driven either by hydrogen-bonding (for the protonated drugs) or by hydrophobic interactions (for the neutral drug molecules). The clustering behavior of the drug molecules appeared to be enhanced in the presence of the membrane, while the development of a charge excess close to the surface of the hyperbranched polymer and of the lipid membrane was observed. The uneven charge distribution created an effective overcharging of the HBP/drug complex and the membrane/drug surface. The translational motion of the drug molecules was found to be strongly affected by the presence of the membrane. The extent of the observed changes depended on the charge of the drug molecule. The build-up of the observed charge excesses close to the surface of the polymeric host and the membrane, together with the changes in the diffusional behavior of the drug molecules are of particular interest. Both phenomena could be important at the latest stages of the liposomal disruption and the release of the drug cargo in formulations based on relevant liposomal carriers.

Cellulose ◽  
2022 ◽  
Author(s):  
Eivind Bering ◽  
Jonathan Ø. Torstensen ◽  
Anders Lervik ◽  
Astrid S. de Wijn

Abstract We investigate the dissolution mechanism of cellulose using molecular dynamics simulations in both water and a mixture solvent consisting of water with Na$$^+$$ + , OH$$^-$$ - and urea. As a first computational study of its kind, we apply periodic external forces that mimic agitation of the suspension. Without the agitation, the bundles do not dissolve, neither in water nor solvent. In the solvent mixture the bundle swells with significant amounts of urea entering the bundle, as well as more water than in the bundles subjected to pure water. We also find that the mixture solution stabilizes cellulose sheets, while in water these immediately collapse into bundles. Under agitation the bundles dissolve more easily in the solvent mixture than in water, where sheets of cellulose remain that are bound together through hydrophobic interactions. Our findings highlight the importance of urea in the solvent, as well as the hydrophobic interactions, and are consistent with experimental results. Graphical abstract


Membranes ◽  
2021 ◽  
Vol 11 (8) ◽  
pp. 580
Author(s):  
Irena Roterman ◽  
Katarzyna Stapor ◽  
Piotr Fabian ◽  
Leszek Konieczny

β-barrel membrane proteins have several important biological functions, including transporting water and solutes across the membrane. They are active in the highly hydrophobic environment of the lipid membrane, as opposed to soluble proteins, which function in a more polar, aqueous environment. Globular soluble proteins typically have a hydrophobic core and a polar surface that interacts favorably with water. In the fuzzy oil drop (FOD) model, this distribution is represented by the 3D Gauss function (3DG). In contrast, membrane proteins expose hydrophobic residues on the surface, and, in the case of ion channels, the polar residues face inwards towards a central pore. The distribution of hydrophobic residues in membrane proteins can be characterized by means of 1–3DG, a complementary 3D Gauss function. Such an analysis was carried out on the transmembrane proteins of bacteria, which, despite the considerable similarities of their super-secondary structure (β-barrel), have highly differentiated properties in terms of stabilization based on hydrophobic interactions. The biological activity and substrate specificity of these proteins are determined by the distribution of the polar and nonpolar amino acids. The present analysis allowed us to compare the ways in which the different proteins interact with antibiotics and helped us understand their relative importance in the development of the resistance mechanism. We showed that beta barrel membrane proteins with a hydrophobic core interact less strongly with the molecules they transport.


2021 ◽  
Author(s):  
Eivind Bering ◽  
Jonathan Torstensen ◽  
Anders Lervik ◽  
Astrid S. de Wijn

Abstract We investigate the dissolution mechanism of cellulose using molecular dynamics simulations in both water and a mixture solvent consisting of water with Na + , OH - and urea. As a first computational study of its kind, we apply periodic external forces that mimic agitation of the suspension. Without the agitation, the bundles do not dissolve, neither in water nor solvent. In the solvent mixture the bundle swells up with significant amounts of urea entering the bundle, as well as more water than in the bundles subjected to pure water. We also find that the mixture solution stabilizes cellulose sheets, while in water these immediately collapse into bundles. Under agitation the bundles dissolve more easily in the solvent mixture than in water, where sheets of cellulose remain that are bound together through hydrophobic interactions. Our findings highlight the importance of urea in the solvent, as well as the hydrophobic interactions, and are consistent with experimental results.


2018 ◽  
Vol 150 (7) ◽  
pp. 949-968 ◽  
Author(s):  
Priyanka Samanta ◽  
Yitang Wang ◽  
Shadi Fuladi ◽  
Jinjing Zou ◽  
Ye Li ◽  
...  

Tight junctions are macromolecular structures that traverse the space between adjacent cells in epithelia and endothelia. Members of the claudin family are known to determine tight junction permeability in a charge- and size-selective manner. Here, we use molecular dynamics simulations to build and refine an atomic model of claudin-15 channels and study its transport properties. Our simulations indicate that claudin-15 forms well-defined channels for ions and molecules and otherwise “seals” the paracellular space through hydrophobic interactions. Ionic currents, calculated from simulation trajectories of wild-type as well as mutant channels, reflect in vitro measurements. The simulations suggest that the selectivity filter is formed by a cage of four aspartic acid residues (D55), contributed by four claudin-15 molecules, which creates a negative electrostatic potential to favor cation flux over anion flux. Charge reversal or charge ablation mutations of D55 significantly reduce cation permeability in silico and in vitro, whereas mutations of other negatively charged pore amino acid residues have a significantly smaller impact on channel permeability and selectivity. The simulations also indicate that water and small ions can pass through the channel, but larger cations, such as tetramethylammonium, do not traverse the pore. Thus, our model provides an atomic view of claudin channels, their transport function, and a potential three-dimensional organization of its selectivity filter.


2020 ◽  
Author(s):  
Robert Stepic ◽  
Lara Jurković ◽  
Ksenia Klementyeva ◽  
Marko Ukrainczyk ◽  
Matija Gredičak ◽  
...  

In many living organisms, biomolecules interact favorably with various surfaces of calcium carbonate. In this work, we have considered the interactions of aspartate (Asp) derivatives, as models of complex biomolecules, with calcite. Using kinetic growth experiments, we have investigated the inhibition of calcite growth by Asp, Asp2 and Asp3.This entailed the determination of a step-pinning growth regime as well as the evaluation of the adsorption constants and binding free energies for the three species to calcite crystals. These latter values are compared to free energy profiles obtained from fully atomistic molecular dynamics simulations. When using a flat (104) calcite surface in the models, the measured trend of binding energies is poorly reproduced. However, a more realistic model comprised of a surface with an island containing edges and corners, yields binding energies that compare very well with experiments. Surprisingly, we find that most binding modes involve the positively charged, ammonium group. Moreover, while attachment of the negatively charged carboxylate groups is also frequently observed, it is always balanced by the aqueous solvation of an equal or greater number of carboxylates. These effects are observed on all calcite features including edges and corners, the latter being associated with dominant affinities to Asp derivatives. As these features are also precisely the active sites for crystal growth, the experimental and theoretical results point strongly to a growth inhibition mechanism whereby these sites become blocked, preventing further attachment of dissolved ions and halting further growth.


2020 ◽  
Author(s):  
Sean A. Newmister ◽  
Kinshuk Raj Srivastava ◽  
Rosa V. Espinoza ◽  
Kersti Caddell Haatveit ◽  
Yogan Khatri ◽  
...  

Biocatalysis offers an expanding and powerful strategy to construct and diversify complex molecules by C-H bond functionalization. Due to their high selectivity, enzymes have become an essential tool for C-H bond functionalization and offer complementary reactivity to small-molecule catalysts. Hemoproteins, particularly cytochromes P450, have proven effective for selective oxidation of unactivated C-H bonds. Previously, we reported the in vitro characterization of an oxidative tailoring cascade in which TamI, a multifunctional P450 functions co-dependently with the TamL flavoprotein to catalyze regio- and stereoselective hydroxylations and epoxidation to yield tirandamycin A and tirandamycin B. TamI follows a defined order including 1) C10 hydroxylation, 2) C11/C12 epoxidation, and 3) C18 hydroxylation. Here we present a structural, biochemical, and computational investigation of TamI to understand the molecular basis of its substrate binding, diverse reactivity, and specific reaction sequence. The crystal structure of TamI in complex with tirandamycin C together with molecular dynamics simulations and targeted mutagenesis suggest that hydrophobic interactions with the polyene chain of its natural substrate are critical for molecular recognition. QM/MM calculations and molecular dynamics simulations of TamI with variant substrates provided detailed information on the molecular basis of sequential reactivity, and pattern of regio- and stereo-selectivity in catalyzing the three-step oxidative cascade.<br>


2018 ◽  
Author(s):  
Anthony Nash ◽  
Nora H de Leeuw ◽  
Helen L Birch

<div> <div> <div> <p>The computational study of advanced glycation end-product cross- links remains largely unexplored given the limited availability of bonded force constants and equilibrium values for molecular dynamics force fields. In this article, we present the bonded force constants, atomic partial charges and equilibrium values of the arginine-lysine cross-links DOGDIC, GODIC and MODIC. The Hessian was derived from a series of <i>ab initio</i> quantum mechanical electronic structure calculations and from which a complete set of force constant and equilibrium values were generated using our publicly available software, ForceGen. Short <i>in vacuo</i> molecular dynamics simulations were performed to validate their implementation against quantum mechanical frequency calculations. </p> </div> </div> </div>


2021 ◽  
Vol 23 (4) ◽  
pp. 2964-2971
Author(s):  
Bernadeta Jasiok ◽  
Mirosław Chorążewski ◽  
Eugene B. Postnikov ◽  
Claude Millot

Thermophysical properties of liquid dibromomethane are investigated by molecular dynamics simulations between 268 and 328 K at pressures up to 3000 bar. Notably, the isotherms of the isobaric thermal expansivity cross around 800 bar.


Glycobiology ◽  
2021 ◽  
Author(s):  
Margrethe Gaardløs ◽  
Sergey A Samsonov ◽  
Marit Sletmoen ◽  
Maya Hjørnevik ◽  
Gerd Inger Sætrom ◽  
...  

Abstract Mannuronan C-5 epimerases catalyse the epimerization of monomer residues in the polysaccharide alginate, changing the physical properties of the biopolymer. The enzymes are utilized to tailor alginate to numerous biological functions by alginate-producing organisms. The underlying molecular mechanisms that control the processive movement of the epimerase along the substrate chain is still elusive. To study this, we have used an interdisciplinary approach combining molecular dynamics simulations with experimental methods from mutant studies of AlgE4, where initial epimerase activity and product formation were addressed with NMR spectroscopy, and characteristics of enzyme-substrate interactions were obtained with isothermal titration calorimetry and optical tweezers. Positive charges lining the substrate-binding groove of AlgE4 appear to control the initial binding of poly-mannuronate, and binding also seems to be mediated by both electrostatic and hydrophobic interactions. After the catalytic reaction, negatively charged enzyme residues might facilitate dissociation of alginate from the positive residues, working like electrostatic switches, allowing the substrate to translocate in the binding groove. Molecular simulations show translocation increments of two monosaccharide units before the next productive binding event resulting in MG-block formation, with the epimerase moving with its N-terminus towards the reducing end of the alginate chain. Our results indicate that the charge pair R343-D345 might be directly involved in conformational changes of a loop that can be important for binding and dissociation. The computational and experimental approaches used in this study complement each other, allowing for a better understanding of individual residues’ roles in binding and movement along the alginate chains.


Sign in / Sign up

Export Citation Format

Share Document