scholarly journals Phospholipid-Cellulose Interactions: Insight from Atomistic Computer Simulations for Understanding the Impact of Cellulose-Based Materials on Plasma Membranes

2018 ◽  
Author(s):  
Andrey A. Gurtovenko ◽  
Evgenii I. Mukhamadiarov ◽  
Andrei Yu. Kostritskii ◽  
Mikko Karttunen

AbstractCellulose is an important biocompatible and nontoxic polymer widely used in numerous biomedical applications. The impact of cellulose-based materials on cells and, more specifically, on plasma membranes that surround cells, however, remains poorly understood. To this end, here we performed atomic-scale molecular dynamics (MD) simulations of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) bilayers interacting with the surface of a cellulose crystal. Both biased umbrella sampling and unbiased simulations clearly show the existence of strong attractive interactions between phospholipids and cellulose: the free energy of the cellulose-bilayer binding was found to be −1.89 and −1.96 kJ/mol per cellulose dimer for PC and PE bilayers, respectively. Although the values are similar, there is a pronounced difference between PC and PE bilayers. The driving force in both cases is the formation of hydrogen bonds. There are two distinct types of hydrogen bonds: 1) between the lipid head groups and the hydroxyl (hydroxymethyl) groups of cellulose, and 2) lipid-water and cellulose-water bonds. The former is the dominant component for PE systems whereas the latter dominates in PC systems. This suggests that achieving controlled binding via new cellulose modifications must pay close attention to the lipid head groups involved. The observed attractive phospholipid-cellulose interactions have a significant impact on bilayer properties: a cellulose crystal induces noticeable structural perturbations on the bilayer leaflet next to the crystal. Given that such perturbations can be undesirable when it comes to the interactions of cellulose-based materials with cell membranes, our computational studies suggest that the impact of cellulose could be reduced through chemical modification of the cellulose surface which prevents cellulose-phospholipid hydrogen bonding.

2021 ◽  
Vol 12 (2) ◽  
pp. 1660-1669

In this work, structural analyses of vanillin (Vanl) and eleven of its derivatives based on the modification of the aldehyde group were investigated using density functional theory (DFT) calculations. In this regard, molecular orbital features and atomic-scale quadrupole coupling constants were evaluated for geometrically optimized structures to see the impact of structural modification on the whole structure. The results indicated that the main impact of such modification was significant only for the modification region, whereas the impact on the rest of the structure was almost negligible. However, electronic features indicated a different tendency for Vanl derivatives for involving in interaction with enzymatic targets. Because of the importance of innovating medication for COVID-19, main protease (MPro) and RNA-dependent RNA polymerase (RdRp) were chosen for the enzymatic target of Vanl ligands for the formation of ligand-target complexes through performing molecular docking (MD) simulations. The results indicated that among the complexes, Vanl 9 (–NHNH2) and 11 (–CH2Cl) could work as the best ligands for interacting with each of RdRp and MPro, respectively. Consequently, optimization of Vanl derivatives could help innovate new compounds for the possible medication of the COVID-19 pandemic.


2020 ◽  
Vol 21 (17) ◽  
pp. 6222 ◽  
Author(s):  
Veronika Zeindlhofer ◽  
Phillip Hudson ◽  
Ádám Márk Pálvölgyi ◽  
Matthias Welsch ◽  
Mazin Almarashi ◽  
...  

In this study, we investigate the influence of chiral and achiral cations on the enantiomerization of biphenylic anions in n-butylmethylether and water. In addition to the impact of the cations and solvent molecules on the free energy profile of rotation, we also explore if chirality transfer between a chiral cation and the biphenylic anion is possible, i.e., if pairing with a chiral cation can energetically favour one conformer of the anion via diastereomeric complex formation. The quantum-mechanical calculations are accompanied by polarizable MD simulations using umbrella sampling to study the impact of solvents of different polarity in more detail. We also discuss how accurate polarizable force fields for biphenylic anions can be constructed from quantum-mechanical reference data.


Author(s):  
N. D. Browning ◽  
M. M. McGibbon ◽  
M. F. Chisholm ◽  
S. J. Pennycook

The recent development of the Z-contrast imaging technique for the VG HB501 UX dedicated STEM, has added a high-resolution imaging facility to a microscope used mainly for microanalysis. This imaging technique not only provides a high-resolution reference image, but as it can be performed simultaneously with electron energy loss spectroscopy (EELS), can be used to position the electron probe at the atomic scale. The spatial resolution of both the image and the energy loss spectrum can be identical, and in principle limited only by the 2.2 Å probe size of the microscope. There now exists, therefore, the possibility to perform chemical analysis of materials on the scale of single atomic columns or planes.In order to achieve atomic resolution energy loss spectroscopy, the range over which a fast electron can cause a particular excitation event, must be less than the interatomic spacing. This range is described classically by the impact parameter, b, which ranges from ~10 Å for the low loss region of the spectrum to <1Å for the core losses.


2019 ◽  
Vol 25 (7) ◽  
pp. 750-773 ◽  
Author(s):  
Pabitra Narayan Samanta ◽  
Supratik Kar ◽  
Jerzy Leszczynski

The rapid advancement of computer architectures and development of mathematical algorithms offer a unique opportunity to leverage the simulation of macromolecular systems at physiologically relevant timescales. Herein, we discuss the impact of diverse structure-based and ligand-based molecular modeling techniques in designing potent and selective antagonists against each adenosine receptor (AR) subtype that constitutes multitude of drug targets. The efficiency and robustness of high-throughput empirical scoring function-based approaches for hit discovery and lead optimization in the AR family are assessed with the help of illustrative examples that have led to nanomolar to sub-micromolar inhibition activities. Recent progress in computer-aided drug discovery through homology modeling, quantitative structure-activity relation, pharmacophore models, and molecular docking coupled with more accurate free energy calculation methods are reported and critically analyzed within the framework of structure-based virtual screening of AR antagonists. Later, the potency and applicability of integrated molecular dynamics (MD) methods are addressed in the context of diligent inspection of intricated AR-antagonist binding processes. MD simulations are exposed to be competent for studying the role of the membrane as well as the receptor flexibility toward the precise evaluation of the biological activities of antagonistbound AR complexes such as ligand binding modes, inhibition affinity, and associated thermodynamic and kinetic parameters.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yogeeshwar Ajjugal ◽  
Kripi Tomar ◽  
D. Krishna Rao ◽  
Thenmalarchelvi Rathinavelan

AbstractBase pair mismatches in DNA can erroneously be incorporated during replication, recombination, etc. Here, the influence of A…A mismatch in the context of 5′CAA·5′TAG sequence is explored using molecular dynamics (MD) simulation, umbrella sampling MD, circular dichroism (CD), microscale thermophoresis (MST) and NMR techniques. MD simulations reveal that the A…A mismatch experiences several transient events such as base flipping, base extrusion, etc. facilitating B–Z junction formation. A…A mismatch may assume such conformational transitions to circumvent the effect of nonisostericity with the flanking canonical base pairs so as to get accommodated in the DNA. CD and 1D proton NMR experiments further reveal that the extent of B–Z junction increases when the number of A…A mismatch in d(CAA)·d(T(A/T)G) increases (1–5). CD titration studies of d(CAA)·d(TAG)n=5 with the hZαADAR1 show the passive binding between the two, wherein, the binding of protein commences with B–Z junction recognition. Umbrella sampling simulation indicates that the mismatch samples anti…+ syn/+ syn…anti, anti…anti & + syn…+ syn glycosyl conformations. The concomitant spontaneous transitions are: a variety of hydrogen bonding patterns, stacking and minor or major groove extrahelical movements (with and without the engagement of hydrogen bonds) involving the mismatch adenines. These transitions frequently happen in anti…anti conformational region compared with the other three regions as revealed from the lifetime of these states. Further, 2D-NOESY experiments indicate that the number of cross-peaks diminishes with the increasing number of A…A mismatches implicating its dynamic nature. The spontaneous extrahelical movement seen in A…A mismatch may be a key pre-trapping event in the mismatch repair due to the accessibility of the base(s) to the sophisticated mismatch repair machinery.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Rajendra P. Koirala ◽  
Rudramani Pokhrel ◽  
Prabin Baral ◽  
Purushottam B. Tiwari ◽  
Prem P. Chapagain ◽  
...  

Abstract Methylation induced DNA base-pairing damage is one of the major causes of cancer. O6-alkylguanine-DNA alkyltransferase (AGT) is considered a demethylation agent of the methylated DNA. Structural investigations with thermodynamic properties of the AGT-DNA complex are still lacking. In this report, we modeled two catalytic states of AGT-DNA interactions and an AGT-DNA covalent complex and explored structural features using molecular dynamics (MD) simulations. We utilized the umbrella sampling method to investigate the changes in the free energy of the interactions in two different AGT-DNA catalytic states, one with methylated GUA in DNA and the other with methylated CYS145 in AGT. These non-covalent complexes represent the pre- and post-repair complexes. Therefore, our study encompasses the process of recognition, complex formation, and separation of the AGT and the damaged (methylated) DNA base. We believe that the use of parameters for the amino acid and nucleotide modifications and for the protein-DNA covalent bond will allow investigations of the DNA repair mechanism as well as the exploration of cancer therapeutics targeting the AGT-DNA complexes at various functional states as well as explorations via stabilization of the complex.


Molecules ◽  
2021 ◽  
Vol 26 (5) ◽  
pp. 1484
Author(s):  
Christopher Reinknecht ◽  
Anthony Riga ◽  
Jasmin Rivera ◽  
David A. Snyder

Proteins are molecular machines requiring flexibility to function. Crystallographic B-factors and Molecular Dynamics (MD) simulations both provide insights into protein flexibility on an atomic scale. Nuclear Magnetic Resonance (NMR) lacks a universally accepted analog of the B-factor. However, a lack of convergence in atomic coordinates in an NMR-based structure calculation also suggests atomic mobility. This paper describes a pattern in the coordinate uncertainties of backbone heavy atoms in NMR-derived structural “ensembles” first noted in the development of FindCore2 (previously called Expanded FindCore: DA Snyder, J Grullon, YJ Huang, R Tejero, GT Montelione, Proteins: Structure, Function, and Bioinformatics 82 (S2), 219–230) and demonstrates that this pattern exists in coordinate variances across MD trajectories but not in crystallographic B-factors. This either suggests that MD trajectories and NMR “ensembles” capture motional behavior of peptide bond units not captured by B-factors or indicates a deficiency common to force fields used in both NMR and MD calculations.


Symmetry ◽  
2021 ◽  
Vol 13 (5) ◽  
pp. 766
Author(s):  
Barbara Bankiewicz ◽  
Marcin Palusiak

The DFT calculations have been performed on a series of two-element complexes formed by substituted 2,6-diaminopyridine (R−PDA) and pyridine (R−Pyr) with X−CY3 molecules (where X = Cl, Br and Y = H, F). The primary aim of this study was to examine the intermolecular hydrogen and halogen bonds in the condition of their mutual coexistence. Symmetry/antisymmetry of the interrelation between three individual interactions is addressed. It appears that halogen bonds play the main role in the stabilization of the structures of the selected systems. However, the occurrence of one or two hydrogen bonds was associated with the favourable geometry of the complexes. Moreover, the impact of different substituent groups attached in the para position to the aromatic ring of the 2,6-diaminopyridine and pyridine on the character of the intermolecular hydrogen and halogen bonds was examined. The results indicate that the presence of electron-donating substituents strengthens the bonds. In turn, the presence of electron-withdrawing substituents reduces the strength of halogen bonds. Additionally, when hydrogen and halogen bonds lose their leading role in the complex formation, the nonspecific electrostatic interactions between dipole moments take their place. Analysis was based on geometric, energetic, and topological parameters of the studied systems.


Crystals ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 127
Author(s):  
YongChao Wang ◽  
YinBo Zhu ◽  
HengAn Wu

The porous characteristics of disordered carbons are critical factors to their performance on hydrogen storage and electrochemical capacitors. Even though the porous information can be estimated indirectly by gas adsorption experiments, it is still hard to directly characterize the porous morphology considering the complex 3D connectivity. To this end, we construct full-atom disordered graphene networks (DGNs) by mimicking the chlorination process of carbide-derived carbons using annealing-MD simulations, which could model the structure of disordered carbons at the atomic scale. The porous characteristics, including pore volume, pore size distribution (PSD), and specific surface area (SSA), were then computed from the coordinates of carbon atoms. From the evolution of structural features, pores grow dramatically during the formation of polyaromatic fragments and sequent disordered framework. Then structure is further graphitized while the PSD shows little change. For the obtained DGNs, the porosity, pore size, and SSA increase with decreasing density. Furthermore, SSA tends to saturate in the low-density range. The DGNs annealed at low temperatures exhibit larger SSA than high-temperature DGNs because of the abundant free edges.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Arunabh Choudhury ◽  
Taj Mohammad ◽  
Nikhil Samarth ◽  
Afzal Hussain ◽  
Md. Tabish Rehman ◽  
...  

AbstractConserved telomere maintenance component 1 (CTC1) is an important component of the CST (CTC1-STN1-TEN1) complex, involved in maintaining the stability of telomeric DNA. Several non-synonymous single-nucleotide polymorphisms (nsSNPs) in CTC1 have been reported to cause Coats plus syndrome and Dyskeratosis congenital diseases. Here, we have performed sequence and structure analyses of nsSNPs of CTC1 using state-of-the-art computational methods. The structure-based study focuses on the C-terminal OB-fold region of CTC1. There are 11 pathogenic mutations identified, and detailed structural analyses were performed. These mutations cause a significant disruption of noncovalent interactions, which may be a possible reason for CTC1 instability and consequent diseases. To see the impact of such mutations on the protein conformation, all-atom molecular dynamics (MD) simulations of CTC1-wild-type (WT) and two of the selected mutations, R806C and R806L for 200 ns, were carried out. A significant conformational change in the structure of the R806C mutant was observed. This study provides a valuable direction to understand the molecular basis of CTC1 dysfunction in disease progression, including Coats plus syndrome.


Sign in / Sign up

Export Citation Format

Share Document