A new GROMOS force field for hexopyranose-based carbohydrates

The Conformation of Glycosidic Linkages According to Various Force Fields: Monte Carlo Modeling of Polysaccharides Based on Extrapolation of Short-Chain Properties

The conformational features of the glycosidic linkage are the most important variable to consider when studying di-, oligo-, and polysaccharide molecules using molecular dynamics (MD) simulations. The accuracy of the theoretical model describing this degree of freedom influences the quality of the results obtained from MD calculations based on this model. This article focuses on the following two issues related to the conformation of the glycosidic linkage. First, we describe the results of a comparative analysis of the predictions of three carbohydrate-dedicated classical force fields for MD simulations, namely, CHARMM, GLYCAM, and GROMOS, in the context of different parameters of structural and energetic nature related to the conformation of selected types of glycosidic linkages, α(1 → 4), β(1 → 3), and β(1 → 4), connecting glucopyranose units. This analysis revealed several differences, mainly concerning the energy levels of the secondary and tertiary conformers and the linkage flexibility within the dominant exo-syn conformation for α(1 → 4) and β(1 → 3) linkages. Some aspects of the comparative analysis also included the newly developed, carbohydrate-dedicated Martini 3 coarse-grained force field. Second, to overcome the time-scale problem associated with sampling slow degrees of freedom in polysaccharide chains during MD simulations, we developed a coarse-grained (CG) model based on the data from MD simulations and designed for Monte Carlo modeling. This model (CG MC) is based on information from simulations of short saccharide chains, effectively sampled in atomistic MD simulations, and is capable of extrapolating local conformational properties to the case of polysaccharides of arbitrary length. The CG MC model has the potential to estimate the conformations of very long polysaccharide chains, taking into account the influence of secondary and tertiary conformations of glycosidic linkages. With respect to the comparative analysis of force fields, the application of CG MC modeling showed that relatively small differences in the predictions of individual force fields with respect to a single glycosidic linkage accumulate when considering their effect on the structure of longer chains, leading to drastically different predictions with respect to parameters describing the polymer conformation, such as the persistence length.

Influence of degree of substitution on the hydroxypropyl-β-cyclodextrin complexation with rifampicin in water solution: a molecular simulation

CONTEXT

Hydroxypropyl-β-cyclodextrin (HPβCD) is one of the derivatized cyclodextrins most widely used as an excipient in the pharmaceutical industry, for its capacity to improve certain drugs properties. Different configurations of HPβCD are possible depending on the number and location of the 2-hydroxypropyl groups substituted on the glucose rings. Rifampicin has become the most commonly clinically used antibiotic against tuberculosis in recent years, despite its low solubility and variable bioavailability. Different techniques and materials have been proposed to enhance the properties of rifampicin: cyclodextrin complexation is one of them. The van der Waals term was the main contribution to the interaction energy, which then decisively conditioned the complex configurations. The size of rifampicin did not allow the whole molecule to fit into the host. Moreover, interaction energy was much greater when the guest was located near each rim of HPβCD, where rifampicin was partially included in the cavity and formed inclusion complexes. The piperazine tail of rifampicin was included inside the host in minimum energy structures and the guest was situated near the primary rim of HPβCD in most cases, although the complex configurations depended on the degree of substitution.

METHODS

A molecular mechanics simulation based on the GROMOS 53A6 force field was applied in this work to study the inclusion complexes formed by twelve configurations of HPβCD, with different degrees of substitution and rifampicin in water solution. We determined the penetration potential, the complex structures with minimum energies, the possibility of forming inclusion complexes other than those of minimum energies and potential energy surfaces.

The Power of Molecular Dynamics Simulations and Their Applications to Discover Cysteine Protease Inhibitors

A large family of enzymes with the function of hydrolyzing peptide bonds, called peptidases or cysteine proteases (CPs), are divided into three categories according to the peptide chain involved. CPs catalyze the hydrolysis of amide, ester, thiol ester, and thioester peptide bonds. They can be divided into several groups, such as papain-like (CA), viral chymotrypsin-like CPs (CB), papainlike endopeptidases of RNA viruses (CC), legumain-type caspases (CD), and showing active residues of His, Glu/Asp, Gln, Cys (CE). The catalytic mechanism of CPs is the essential cysteine residue present in the active site. These mechanisms are often studied through computational methods that provide new information about the catalytic mechanism and identify inhibitors. The role of computational methods during drug design and development stages is increasing. Methods in Computer-Aided Drug Design (CADD) accelerate the discovery process, increase the chances of selecting more promising molecules for experimental studies, and can identify critical mechanisms involved in the pathophysiology and molecular pathways of action. Molecular dynamics (MD) simulations are essential in any drug discovery program due to their high capacity for simulating a physiological environment capable of unveiling significant inhibition mechanisms of new compounds against target proteins, especially CPs. Here, a brief approach will be shown on MD simulations and how the studies were applied to identify inhibitors or critical information against cysteine protease from several microorganisms, such as Trypanosoma cruzi (cruzain), Trypanosoma brucei (rhodesain), Plasmodium spp. (falcipain), and SARS-CoV-2 (Mpro). We hope the readers will gain new insights and use our study as a guide for potential compound identifications using MD simulations.

Effects of sunflower pectin on thermal stability of purple sweet potato anthocyanins at different pH

The present study aimed to examine the impact of sunflower pectin (SFP) on the thermal stability and antioxidant activity of purple sweet potato anthocyanins (PSPA) at varying pH levels. It was observed that the pH value significantly influenced the ability of pectin to protect anthocyanins from thermal degradation, which was found to be associated with the rate of binding between PSPA and SFP. The binding rate of PSPA-SFP was observed to be highest at pH 4.0, primarily due to the influence of electrostatic interaction and hydrogen bonding. Monoacylated anthocyanins exhibited a binding rate approximately 2-4 % higher than that of diacylated anthocyanins. The PSPA-SFP demonstrated its highest thermal stability at pH 4.0, with a corresponding half-life of 14.80 h at 100 °C. Molecular dynamics simulations indicated that pectin had a greater affinity for the flavylium cation and hemiketal form of anthocyanins. The antioxidant activity of anthocyanins in PSPA and PSPA-SFP increased with increasing pH, suggesting that anthocyanins at high pH had higher antioxidant activity than anthocyanins at low pH.

Probing Disaccharide Binding to Triplatin as Models for Tumor Cell Heparan Sulfate (GAG) Interactions

In this study, we have used [1H, 15N] NMR spectroscopy to investigate the interactions of the trinuclear platinum anticancer drug triplatin (1) (1,0,1/t,t,t or BBR3464) with site-specific sulfated and carboxylated disaccharides. Specifically, the disaccharides GlcNS(6S)-GlcA (I) and GlcNS(6S)-IdoA(2S) (II) are useful models of longer-chain glycosaminoglycans (GAGs) such as heparan sulfate (HS). For both the reactions of 15N-1 with I and II, equilibrium conditions were achieved more slowly (65 h) compared to the reaction with the monosaccharide GlcNS(6S) (9 h). The data suggest both carboxylate and sulfate binding of disaccharide I to the Pt with the sulfato species accounting for <1% of the total species at equilibrium. The rate constant for sulfate displacement of the aqua ligand (kL2) is 4 times higher than the analogous rate constant for carboxylate displacement (kL1). There are marked differences in the equilibrium concentrations of the chlorido, aqua, and carboxy-bound species for reactions with the two disaccharides, notably a significantly higher concentration of carboxylate-bound species for II, where sulfate-bound species were barely detectable. The trend mirrors that reported for the corresponding dinuclear platinum complex 1,1/t,t, where the rate constant for sulfate displacement of the aqua ligand was 3 times higher than that for acetate. Also similar to what we observed for the reactions of 1,1/t,t with the simple anions, aquation of the sulfato group is rapid, and the rate constant k-L2 is 3 orders of magnitude higher than that for displacement of the carboxylate (k-L1). Molecular dynamics calculations suggest that extra hydrogen-bonding interactions with the more sulfated disaccharide II may prevent or diminish sulfate binding of the triplatin moiety. The overall results suggest that Pt-O donor interactions should be considered in any full description of platinum complex cellular chemistry.

Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

Investigation of the effects of N-Acetylglucosamine on the stability of the spike protein in SARS-CoV-2 by molecular dynamics simulations

A lot of effort has been made in developing vaccine and therapeutic agents against the SARS-CoV-2, concentrating on the Spike protein that binds angiotensin-converting enzyme 2 on human cells. Nowadays, some researches study the role of the N-linked glycans as potential targets for vaccines and new agents. Due to the flexibility and diversity of the N-linked glycans, in this work, we focus on the N-Acetylglucosamine moiety, which is the precursor of nearly all eukaryotic glycans. We performed molecular dynamics simulations to study the effects of the N-Acetylglucosamine on the stability of the spike glycoprotein in SARS-CoV-2. After a 100 ns of simulation on the spike proteins without and with the N-Acetylglucosamine molecules, we found that the presence of N-Acetylglucosamine increases the local stability in their vicinity; even though their effect on the full structure is negligible. Thus; it can be inferred that the N-Acetylglucosamine moieties can potentially affect the interaction of the S protein with the ACE2 receptor. We also found that the S1 domain is more flexible than the S2 domain. We propose which of the experimentally observed glycans found on the spike may be more functional than the others. Detailed understanding of glycans is key for the development of new therapeutic strategies.

Molecular Dynamics Study of the Curvature-Driven Interactions between Carbon-Based Nanoparticles and Amino Acids

We researched the interaction between six representative carbon-based nanoparticles (CBNs) and 20 standard amino acids through all-atom molecular dynamics simulations. The six carbon-based nanoparticles are fullerene(C60), CNT55L3, CNT1010L3, CNT1515L3, CNT2020L3, and two-dimensional graphene (graphene33). Their curvatures decrease sequentially, and all of the CNTs are single-walled carbon nanotubes. We observed that as the curvature of CBNs decreases, the adsorption effect of the 20 amino acids with them has an increasing trend. In addition, we also used multi-dimensional clustering to analyze the adsorption effects of 20 amino acids on six carbon-based nanoparticles. We observed that the π-π interaction still plays an extremely important role in the adsorption of amino acids on carbon-based nanoparticles. Individual long-chain amino acids and "Benzene-like" Pro also have a strong adsorption effect on carbon-based nanoparticles.

Research and study of 2-((4,6 dimethyl pyrimidine-2-yle) thio)-N-phenyl acetamide derivatives as inhibitors of sirtuin 2 protein for the treatment of cancer using QSAR, molecular docking and molecular dynamic simulation

Phenyl acetamide derivatives have a wide range of biological activities, so their research and development can be useful and effective for the design production of new drugs. In this project, quantitative structure-activity relationship (QSAR) was performed. For modeling two methods of multiple linear regression (MLR) and nonlinear regression of support vector machine (SVR) were used. In the MLR stage, the best model with the values of R²_train = 0.913 and R²_test = 0.881 was selected by stepwise method. In this model, 4 descriptors of BELV2, GATS8p, GATS6e and RDF080m were included, which were used as input for the nonlinear support vector regression method. In the SVR model, the best results were obtained using the radial Gaussian kernel function (RBF) with R²_train = 0.978 and R²_test = 0.990. In the next step, using molecular docking and molecular dynamic simulation methods, the interaction between phenyl acetamide derivatives and the sirtuin 2 protein was investigated. Examining the results of molecular docking, it was observed that these derivatives formed complexes by forming hydrogen and hydrophobic bonds with the sirtuin 2 protein. Also, the results of molecular dynamic simulation show that phenyl acetamide compounds form stable complex with the sirtuin 2 protein, and it was found that the compounds with more activity have formed a number of hydrogen bonds with the protein.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

GM2/GM3 controls the organizational status of CD82/Met microdomains: further studies in GM2/GM3 complexation

At cell surface gangliosides might associate with signal transducers proteins, grown factor receptors, integrins, small G-proteins and tetraspanins establishing microdomains, which play important role in cell adhesion, cell activation, motility, and growth. Previously, we reported that GM2 and GM3 form a heterodimer that interacts with the tetraspanin CD82, controlling epithelial cell mobility by inhibiting integrin-hepatocyte growth factor-induced cMet tyrosine kinase signaling. By using molecular dynamics simulations to study the molecular basis of GM2/GM3 interaction we demonstrate, here, that intracellular levels of Ca²⁺ mediate GM2/GM3 complexation via electrostatic interaction with their carboxyl groups, while hydrogen bonds between the ceramide groups likely aid stabilizing the complex. The presence of GM2/GM3 complex alters localization of CD82 on cell surface and therefore downstream signalization. These data contribute for the knowledge of how glycosylation may control signal transduction and phenotypic changes.

Carbohydrate structure database (CSDB) oligosaccharide conformation tool

Population analysis in terms of glycosidic torsion angles is frequently used to reveal preferred conformers of glycans. However, due to high structural diversity and flexibility of carbohydrates, conformational characterization of complex glycans can be a challenging task. Herein we present a conformation module of oligosaccharide fragments occurring in natural glycan structures developed on the platform of the Carbohydrate Structure Database (CSDB). Currently, this module deposits free energy surface and conformer abundance maps plotted as a function of glycosidic torsions for 194 inter-residue bonds. Data are automatically and continuously derived from explicit-solvent molecular dynamics (MD) simulations. The module was also supplemented with high-temperature MD data of saccharides (2403 maps) provided by GlycoMapsDB (hosted by GLYCOSCIENCES.de project). Conformational data defined by up to four torsional degrees of freedom can be freely explored using a web interface of the module available at http://csdb.glycoscience.ru/database/core/search_conf.html.

Aqueous Molecular Dynamics for Understanding Glycosaminoglycan Recognition by Proteins

Glycosaminoglycans (GAGs) are biopolymers that exist in most organisms. GAGs are known to bind to hundreds of proteins and partake in multiple biological processes such as growth, morphogenesis, inflammation, infection, and others. Their intrinsic structural heterogeneity and conformational variability introduce major challenges in experimental studies. On the other hand, recent advances in force field development and computational technology have yielded phenomenal opportunity to study thousands of GAG sequences simultaneously to understand recognition of target protein(s). Here, we describe experimental setup for conventional molecular dynamics simulations of GAGs to position an experimental biologist favorably in performance, analysis and interpretation of stability, specificity, and conformational properties of GAGs, while also elucidating their interactions with amino acid residues of a protein at an atomistic level in presence of water.

Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins

The conformational properties of carbohydrates can contribute to protein structure directly through covalent conjugation in the cases of glycoproteins and proteoglycans and indirectly in the case of transmembrane proteins embedded in glycolipid-containing bilayers. However, there continue to be significant challenges associated with experimental structural biology of such carbohydrate-containing systems. All-atom explicit-solvent molecular dynamics simulations provide a direct atomic resolution view of biomolecular dynamics and thermodynamics, but the accuracy of the results depends on the quality of the force field parametrization used in the simulations. A key determinant of the conformational properties of carbohydrates is ring puckering. Here, we applied extended system adaptive biasing force (eABF) all-atom explicit-solvent molecular dynamics simulations to characterize the ring puckering thermodynamics of the ten common pyranose monosaccharides found in vertebrate biology (as represented by the CHARMM carbohydrate force field). The results, along with those for idose, demonstrate that the CHARMM force field reliably models ring puckering across this diverse set of molecules, including accurately capturing the subtle balance between 4C1 and 1C4 chair conformations in the cases of iduronate and of idose. This suggests the broad applicability of the force field for accurate modeling of carbohydrate-containing vertebrate biomolecules such as glycoproteins, proteoglycans, and glycolipids.

Structural dynamics and vibrational feature of N-Acetyl-d-glucosamine in aqueous solution

Molecular dynamics simulations and DFT calculations were performed for the demonstration of the structural dynamics and vibrational feature of N-Acetyl-d-glucosamine (NAG) in solution phase. The interactions between NAG and solvent molecules were evaluated through spatial distribution function and radial distribution function, and the preferred conformations of NAG in aqueous solution were revealed by cluster analysis. Results from normal mode analysis show that the solvent induced structural fluctuation of NAG could be reflected in the vibrational feature of specific chromophores, thus we can evaluate the molecular structure with the help of its vibrational signature based on the built correlation between molecular structure and vibrational frequencies of specific groups.

Fusion pores with low conductance are cation selective

Many neurotransmitters are organic ions that carry a net charge, and their release from secretory vesicles is therefore an electrodiffusion process. The selectivity of early exocytotic fusion pores is investigated by combining electrodiffusion theory, measurements of amperometric foot signals from chromaffin cells with anion substitution, and molecular dynamics simulation. The results reveal that very narrow fusion pores are cation selective, but more dilated fusion pores become anion permeable. The transition occurs around a fusion pore conductance of ~300 pS. The cation selectivity of a narrow fusion pore accelerates the release of positively charged transmitters such as dopamine, noradrenaline, adrenaline, serotonin, and acetylcholine, while glutamate release may require a more dilated fusion pore.

For transmission, a fusion pore forms when vesicle and target membranes are brought together by SNARE proteins. Delacruz et al. demonstrate that selectivity of the pore accelerates release of positively charged transmitters such as dopamine, noradrenaline, adrenaline, serotonin, and acetylcholine, while glutamate release may require a more dilated fusion pore.

Preparing and Analyzing Polarizable Molecular Dynamics Simulations with the Classical Drude Oscillator Model

Molecular dynamics (MD) simulations performed with force fields that include explicit electronic polarization are becoming more prevalent in the field. The increasing emergence of these simulations is a result of continual refinement against a range of theoretical and empirical target data, optimization of software algorithms for higher performance, and availability of graphical processing unit hardware to further accelerate the simulations. Polarizable MD simulations are likely to be most impactful in biomolecular systems in which heterogeneous environments or unique microenvironments exist that would lead to inaccuracies in simulations performed with fixed-charge, nonpolarizable force fields. The further adoption of polarizable MD simulations will benefit from tutorial material that specifically addresses preparing and analyzing their unique features. In this chapter, we introduce common protocols for preparing routine biomolecular systems containing proteins, including both a globular protein in aqueous solvent and a transmembrane model peptide in a phospholipid bilayer. Details and example input files are provided for preparation of the simulation system using CHARMM, performing the simulations with OpenMM, and analyzing interesting dipole moment properties in CHARMM.

The tug of war between Al3+ and Na+ for order-disorder transitions in lipid-A membranes

Cations play a critical role in the stability and morphology of lipid-A aggregates by neutralizing, hydrating and cross-linking these glycolipid molecules. Monophosphorylated lipid-A is the major immunostimulatory principle in commercially available adjuvants containing Al3+ such as adjuvant system 04 (AS04). The antagonist/agonist immunomodulatory properties of lipid-A are associated with chemical variations (e.g. the number of acyl chains and phosphate groups) and their aggregate arrangements (e.g. lamellar, nonlamellar or mixed). Therefore, the identification of the active form of lipid-A can provide valuable guidance in the development of vaccine adjuvants capable of boosting the immune system with decreased reactogenicity. Although the effect of mono and divalent cations on the structural polymorphism and endotoxicity of LPS has been previously investigated, much less is known about the effect of trivalent cations. We have investigated the effect of NaCl and AlCl3 salt solutions on the structural dynamics and stability of mono and diphosphorylated lipid-A membranes via atomistic MD simulations. The Al3+ ion exerts two major effects on the structural dynamics of lipid-A membranes. It acts as an efficient cross-linker of mono or diphosphorylated lipid-A molecules, thus stabilizing the lamellar arrangement of these glycolipids. It also alters the lipid-A packing and membrane fluidity, inducing disorder → order structural transitions of the membrane. This effect is promptly reversed upon the addition of NaCl solution, which promotes a nearly threefold increase in the amount of water in the carbohydrate moiety of the Al3+-containing lipid-A membranes. The exchange dynamics and residence times of cation-coordinated water molecules in these membranes provide insights into the molecular mechanism for the Na+-induced transition from a densely packed ordered phase to a disordered one. Al3+ counter-ions favor ordered lamellar aggregates, which has been previously associated with the lack of endotoxic activity and cytokine-inducing action. The resulting microscopic understanding of the structure and dynamics of lipid-A aggregates in the presence of Al3+ and Na+ salts can provide valuable guidance in the development of vaccine adjuvants capable of boosting the immune system with decreased reactogenicity.

In silico and in vitro design of cordycepin encapsulation in liposomes for colon cancer treatment

Cordycepin or 3'-deoxyadenosine is an interesting anti-cancer drug candidate that is found in abundance in the fungus Cordyceps militaris. It inhibits cellular growth of many cancers including lung carcinoma, melanoma, bladder cancer, and colon cancer by inducing apoptosis, anti-proliferation, anti-metastasis and by arresting the cell cycle. Cordycepin has, however, poor stability and low solubility in water, resulting in loss of its bioactivity. Liposomes can be used to overcome these obstacles. Our aim is to improve cordycepin's anti-colon cancer activity by liposome encapsulation. Cordycepin-encapsulated liposomes were designed and fabricated based on a combination of theoretical and experimental studies. Molecular dynamics (MD) simulations and free energy calculations suggest that phosphatidylcholine (PC) lipid environment is favorable for cordycepin adsorption. Cordycepin passively permeates into PC lipid bilayers without membrane damage and strongly binds to the lipids' polar groups by flipping its deoxyribose sugar toward the bilayer center. Our fabricated liposomes containing 10 : 1 molar ratio of egg yolk PC : cholesterol showed encapsulation efficiency (%EE) of 99% using microfluidic hydrodynamic focusing (MHF) methods. In our in vitro study using the HT-29 colon cancer cell line, cordycepin was able to inhibit growth by induction of apoptosis. Cell viability was significantly decreased below 50% at 125 μg mL-1 dosage after 48 h treatment with non-encapsulated and encapsulated cordycepin. Importantly, encapsulation provided (1) a 2-fold improvement in the inhibition of cancer cell growth at 125 μg mL-1 dosage and (2) 4-fold increase in release time. These in silico and in vitro studies indicate that cordycepin-encapsulated liposomes could be a potent drug candidate for colon cancer therapy.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

Out of Sight, Out of Mind: The Effect of the Equilibration Protocol on the Structural Ensembles of Charged Glycolipid Bilayers

Molecular dynamics (MD) simulations represent an essential tool in the toolbox of modern chemistry, enabling the prediction of experimental observables for a variety of chemical systems and processes and majorly impacting the study of biological membranes. However, the chemical diversity of complex lipids beyond phospholipids brings new challenges to well-established protocols used in MD simulations of soft matter and requires continuous assessment to ensure simulation reproducibility and minimize unphysical behavior. Lipopolysaccharides (LPS) are highly charged glycolipids whose aggregation in a lamellar arrangement requires the binding of numerous cations to oppositely charged groups deep inside the membrane. The delicate balance between the fully hydrated carbohydrate region and the smaller hydrophobic core makes LPS membranes very sensitive to the choice of equilibration protocol. In this work, we show that the protocol successfully used to equilibrate phospholipid bilayers when applied to complex lipopolysaccharide membranes occasionally leads to a small expansion of the simulation box very early in the equilibration phase. Although the use of a barostat algorithm controls the system dimension and particle distances according to the target pressure, fluctuation in the fleeting pressure occasionally enables a few water molecules to trickle into the hydrophobic region of the membrane, with spurious solvent buildup. We show that this effect stems from the initial steps of NPT equilibration, where initial pressure can be fairly high. This can be solved with the use of a stepwise-thermalization NVT/NPT protocol, as demonstrated for atomistic MD simulations of LPS/DPPE and lipid-A membranes in the presence of different salts using an extension of the GROMOS forcefield within the GROMACS software. This equilibration protocol should be standard procedure for the generation of consistent structural ensembles of charged glycolipids starting from atomic coordinates not previously pre-equilibrated. Although different ways to deal with this issue can be envisioned, we investigated one alternative that could be readily available in major MD engines with general users in mind.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Computational tools for drawing, building and displaying carbohydrates: a visual guide

Drawing and visualisation of molecular structures are some of the most common tasks carried out in structural glycobiology, typically using various software. In this perspective article, we outline developments in the computational tools for the sketching, visualisation and modelling of glycans. The article also provides details on the standard representation of glycans, and glycoconjugates, which helps the communication of structure details within the scientific community. We highlight the comparative analysis of the available tools which could help researchers to perform various tasks related to structure representation and model building of glycans. These tools can be useful for glycobiologists or any researcher looking for a ready to use, simple program for the sketching or building of glycans.

ADD Force Field for Sugars and Polyols: Predicting the Additivity of Protein-Osmolyte Interaction

The protein–osmolyte interaction has been shown experimentally to follow an additive construct, where the individual osmolyte–backbone and osmolyte–side-chain interactions contribute to the overall conformational stability of proteins. Here, we computationally reconstruct this additive relation using molecular dynamics simulations, focusing on sugars and polyols, including sucrose and sorbitol, as model osmolytes. A new set of parameters (ADD) is developed for this purpose, using the individual Kirkwood–Buff integrals for sugar–backbone and sugar–side-chain interactions as target experimental data. We show that the ADD parameters can reproduce the additivity of protein–sugar interactions and correctly predict sucrose and sorbitol self-association and their interaction with water. The accurate description of the separate osmolyte–backbone and osmolyte–side-chain contributions also automatically translates into a good prediction of preferential exclusion from the surface of ribonuclease A and α-chymotrypsinogen A. The description of sugar polarity is improved compared to previous force fields, resulting in closer agreement with the experimental data and better compatibility with charged groups, such as the guanidinium moiety. The ADD parameters are developed in combination with the CHARMM36m force field for proteins, but good compatibility is also observed with the AMBER 99SB-ILDN and the OPLS-AA force fields. Overall, exploiting the additivity of protein–osmolyte interactions is a promising approach for the development of new force fields.

Molecular dynamic simulation: Conformational properties of single-stranded curdlan in aqueous solution

Recently, molecular dynamic simulation technique has been proved to be a powerful tool providing structural insights for better understanding the functionality of carbohydrates. Here, by using molecular dynamic simulation method we investigated the detailed conformational properties of the single-stranded curdlan with 12 glucose units. The results showed that the right-handed 6/1 helix structure was thermodynamically the most stable conformation in solution. The formation of the helix conformation was determined by many factors such as the glycosidic linkage, explicit water solvation and hydrogen bonds. When temperature was increased, the representative helix conformation was found becoming unstable giving rise to metastable conformations because when water mobility was accelerated with temperature, the hydrogen bonding strength between the curdlan chain and water went down, breaking the continuity of the hydrogen bonding network of water and hydroxyl groups. When the number of repeating glucose units varied from 6 to 24, the major helix conformation remained, but the conformational properties of longer chains were more apparently aff ;ected by chain flexibility.

Extending the Martini Coarse-Grained Force Field to N-Glycans

Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.

Rigorous analysis of free solution glycosaminoglycan dynamics using simple, new tools

Heparin/heparan sulfates (H/HS) are ubiquitous biopolymers that interact with many proteins to induce a range of biological functions. Unfortunately, how these biopolymers recognize their preferred protein targets remain poorly understood. It is suggested that computational simulations offer attractive avenues but a number of challenges, e.g., difficulty of selecting a comprehensive force field, few simple tools to interpret data, among others, remain. This work addresses several such challenges so as to help ease the implementation and analysis of computational experiments. First, this work presents a rigorous comparison of two different recent force fields, CHARMM36 and GLYCAM06, for H/HS studies. Second, it introduces two new straightforward parameters, i.e., end-to-end distance and minimum volume enclosing ellipsoid, to understand the myriad conformational forms of oligosaccharides that evolve over time in water. Third, it presents an application to elucidate the number and nature of inter and intramolecular, nondirect bridging water molecules, which help stabilize unique forms of H/HS. The results show that nonspecialists can use either CHARMM36 or GLYCAM06 force fields because both gave comparable results, albeit with small differences. The comparative study shows that the HS hexasaccharide samples a range of conformations with nearly equivalent energies, which could be the reason for its recognition by different proteins. Finally, analysis of the nondirect water bridges across the dynamics trajectory shows their importance in stabilization of certain conformational forms, which may become important for protein recognition. Overall, the work aids nonspecialists employ computational studies for understanding the solution behavior of H/HS.

Interactions of a Bacterial RND Transporter with a Transmembrane Small Protein in a Lipid Environment

The small protein AcrZ in Escherichia coli interacts with the transmembrane portion of the multidrug efflux pump AcrB and increases resistance of the bacterium to a subset of the antibiotic substrates of that transporter. It is not clear how the physical association of the two proteins selectively changes activity of the pump for defined substrates. Here, we report cryo-EM structures of AcrB and the AcrBZ complex in lipid environments, and comparisons suggest that conformational changes occur in the drug-binding pocket as a result of AcrZ binding. Simulations indicate that cardiolipin preferentially interacts with the AcrBZ complex, due to increased contact surface, and we observe that chloramphenicol sensitivity of bacteria lacking AcrZ is exacerbated when combined with cardiolipin deficiency. Taken together, the data suggest that AcrZ and lipid cooperate to allosterically modulate AcrB activity. This mode of regulation by a small protein and lipid may occur for other membrane proteins.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Molecular dynamics simulations of two double-helical hexamer fragments of iota-carrageenan in aqueous solution

The gelation of anionic carrageenans is known to occur through a coil-to-helix transition followed by further aggregation or association on which positive counterions play a significant role. In the present work, molecular dynamics (MD) simulations were performed on two double-helical iota-carrageenan hexamer fragments along with their sodium counterions using the Carbohydrate Solution Force Field (CSFF) in an aqueous (TIP4P) solution with the GROMACS molecular dynamics package. Results showed a counterion condensation between the two double helices and that the subsequent forces of interaction between them were predominantly attractive. By varying the distance separating the two double helices, the effect of distance on the counterion distribution and the forces of attraction was also investigated. In the presence of counterions, an increase in the forces of attraction was observed as the distance between the two double helices decreases which can be attributed to the greater counterion density between the two like-charged oligosaccharides.

Pairwise-additive and polarizable atomistic force fields for molecular dynamics simulations of proteins

Protein force fields have been undergoing continual development since the first complete parameter sets were introduced nearly four decades ago. The functional forms that underlie these models have many common elements for the treatment of bonded and nonbonded forces, which are reviewed here. The most widely used force fields to date use a fixed-charge convention in which electronic polarization effects are treated via a mean-field approximation during partial charge assignment. Despite success in modeling folded proteins over many years, the fixed-charge assumption has limitations that cannot necessarily be overcome within their potential energy equations. To overcome these limitations, several force fields have recently been derived that explicitly treat electronic polarization effects with straightforward extensions of the potential energy functions used by nonpolarizable force fields. Here, we review the history of the most popular nonpolarizable force fields (AMBER, CHARMM, OPLS, and GROMOS) as well as studies that have validated them and applied them to studies of protein folding and misfolding. Building upon these force fields are more recent polarizable interaction potentials, including fluctuating charge models, POSSIM, AMOEBA, and the classical Drude oscillator. These force fields differ in their implementations but all attempt to model electronic polarization in a computationally tractable manner. Despite their recent emergence in the field of protein folding, several studies have already applied these polarizable models to challenging problems in this domain, including the role of polarization in folding free energies and sequence-specific effects on the stability of α-helical structures.

Entropy deepens loading chemical potentials of small alcohols by narrow carbon nanotubes

Small alcohol confinement within narrow carbon nanotubes has been extensively and systematically studied via rigorous free-energy calculations.

Molecular Dynamics with Orientational Tensorial Constraints: A New Approach to Probe the Torsional Angle Distributions of Small Rotationally Flexible Molecules

The potential of residual dipolar couplings (RDCs) in conformational studies of small molecules is now widely recognized, but current theoretical approaches for their interpretation have several limitations and there is still the need for a general method to probe the torsional angle distributions applicable to any rotationally flexible molecule. Molecular dynamics simulations with RDC-based orientational tensorial constraints (MDOC), implemented in the software COSMOS, are presented here as a conceptually new strategy. For the cases of the fluorinated anti-inflammatory drug diflunisal and the disaccharide cellobiose, we demonstrate that MDOC simulations with one-bond RDCs as tensorial constraints unveil torsion distributions and allow the determination of relative configuration in the presence of rotational flexibility. The independence of the initial structure or any a priori assumption as well as the possibility to combine different experimental constraints represent features, which make the COSMOS software a promising tool for the investigation of torsional angle distributions of flexible molecules, regardless of their size and degree of freedom.

Wood–Moisture Relationships Studied with Molecular Simulations: Methodological Guidelines

This paper aims at providing a methodological framework for investigating wood polymers using atomistic modeling, namely, molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. Atomistic simulations are used to mimic water adsorption and desorption in amorphous polymers, make observations on swelling, mechanical softening, and on hysteresis. This hygromechanical behavior, as observed in particular from the breaking and reforming of hydrogen bonds, is related to the behavior of more complex polymeric composites. Wood is a hierarchical material, where the origin of wood-moisture relationships lies at the nanoporous material scale. As water molecules are adsorbed into the hydrophilic matrix in the cell walls, the induced fluid–solid interaction forces result in swelling of these cell walls. The interaction of the composite polymeric material, that is the layer S2 of the wood cell wall, with water is known to rearrange its internal material structure, which makes it moisture sensitive, influencing its physical properties. In-depth studies of the coupled effects of water sorption on hygric and mechanical properties of different polymeric components can be performed with atomistic modeling. The paper covers the main components of knowledge and good practice for such simulations.

Structural Aspects of the O-glycosylation Linkage in Glycopeptides via MD Simulations and Comparison with NMR Experiments

A powerful conformational searching and enhanced sampling simulation method, and unbiased molecular dynamics simulations have been used along with NMR spectroscopic observables to provide a detailed structural view of O-glycosylation. For four model systems, the force-field parameters can accurately predict experimental NMR observables (J couplings and NOE's). This enables us to derive conclusions based on the generated ensembles, in which O-glycosylation affects the peptide backbone conformation by forcing it towards to an extended conformation. An exception is described for β-GalNAc-Thr where the α content is increased and stabilized via hydrogen bonding between the sugar and the peptide backbone, which was not observed in the rest of the studied systems. These observations might offer an explanation for the evolutionary preference of α-linked GalNAc glycosylation instead of a β link.

Perspective on computational simulations of glycosaminoglycans

Glycosaminoglycans (GAGs) represent a formidable frontier for chemists, biochemists, biologists, medicinal chemists and drug delivery specialists because of massive structural complexity. GAGs are arguably the most complex, natural linear biopolymers with theoretical diversity orders of magnitude higher than proteins and nucleic acids. Yet, this diversity remains generally untapped. Computational approaches offer major routes to understand GAG structure and dynamics so as to enable novel applications of these biopolymers. In fact, computational algorithms, softwares, online tools and techniques have reached a level of sophistication that help understand atomistic details of conformational variation and protein recognition of individual GAG sequences. This review describes current approaches and challenges in computational study of GAGs. It presents a history of major findings since the earliest mention of GAGs (the 1960s), the development of parameters and force fields specific for GAGs, and the application of these tools in understanding GAG structure-function relationship. This review also presents a section on how to perform simulation of GAGs, which is directed toward researchers interested in entering this promising field with potential to impact therapy.

Glucose Can Protect Membranes against Dehydration Damage by Inducing a Glassy Membrane State at Low Hydrations

The physical effects of small sugars on membranes have been studied for decades, primarily because of their membrane stabilization in cold or dehydrated environments. We studied the effects of up to 20 mol% glucose in bilayers made of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) at low hydration by combining X-ray diffraction and Molecular Dynamics (MD) simulations. In agreement with previous studies, we observe membrane thinning at low and membrane thickening at high sugar concentrations. Glucose was found to preferentially localize to the outer head region of phospholipid bilayers at all concentrations, and partitioning of sugar in the membranes was found to monotonically increase with increasing sugar concentration. While the number of gauche defects in the lipid acyl tails and the lipid packing in the presence of sugar resembled values of a fluid lipid bilayer, tail dynamics, as assessed by autocorrelation of the carbon atoms in the phospholipid tails, were slowed down significantly with increasing glucose content. Thus, our findings suggest that sugar leads to a a disordered, glassy state of the hydrophobic membrane core. The non-monotonic effect of glucose on membrane thickness was found to be an effect of fluidification at low concentrations and decreased interdigitation in the higher sugar concentration regime.

Design, synthesis, and biomedical applications of synthetic sulphated polysaccharides

This review summarizes the synthetic methods to sulphated polysaccharides, describes their compositional and structural diversity in regards to activity, and showcases their biomedical applications.

Molecular Dynamics Simulation for Mechanism Elucidation of Food Processing and Safety: State of the Art

Molecular dynamics (MD) simulation is a useful technique to study the interaction between molecules and how they are affected by various processes and processing conditions. This review summarizes the application of MD simulations in food processing and safety, with an emphasis on the effects that emerging nonthermal technologies (for example, high hydrostatic pressure, pulsed electric field) have on the molecular and structural characteristics of foods and biomaterials. The advances and potential projection of MD simulations in the science and engineering aspects of food materials are discussed and focused on research work conducted to study the effects of emerging technologies on food components. It is expected by showing key case studies that it will stir novel developments as a valuable tool to study the effects of emerging food technologies on biomaterials. This review is useful to food researchers and the food industry, as well as researchers and practitioners working on flavor and nutraceutical encapsulations, dietary carbohydrate product developments, modified starches, protein engineering, and other novel food applications.

Optimizing Nonbonded Interactions of the OPLS Force Field for Aqueous Solutions of Carbohydrates: How to Capture Both Thermodynamics and Dynamics

Knowledge on thermodynamic and transport properties of aqueous solutions of carbohydrates is of great interest for process and product design in the food, pharmaceutical, and biotechnological industries. Molecular simulation is a powerful tool to calculate these properties, but current classical force fields cannot provide accurate estimates for all properties of interest. The poor performance of the force fields is mainly observed for concentrated solutions, where solute–solute interactions are overestimated. In this study, we propose a method to refine force fields, such that solute–solute interactions are more accurately described. The OPLS force field combined with the SPC/Fw water model is used as a basis. We scale the nonbonded interaction parameters of sucrose, a disaccharide. The scaling factors are chosen in such a way that experimental thermodynamic and transport properties of aqueous solutions of sucrose are accurately reproduced. Using a scaling factor of 0.8 for Lennard-Jones energy parameters (ϵ) and a scaling factor of 0.95 for partial atomic charges (q), we find excellent agreement between experiments and computed liquid densities, thermodynamic factors, shear viscosities, self-diffusion coefficients, and Fick (mutual) diffusion coefficients. The transferability of these optimum scaling factors to other carbohydrates is verified by computing thermodynamic and transport properties of aqueous solutions of d-glucose, a monosaccharide. The good agreement between computed properties and experiments suggests that the scaled interaction parameters are transferable to other carbohydrates, especially for concentrated solutions.

Biopolymer nanofibrils: structure, modeling, preparation, and applications

Biopolymer nanofibrils exhibit exceptional mechanical properties with a unique combination of strength and toughness, while also presenting biological functions that interact with the surrounding environment. These features of biopolymer nanofibrils profit from their hierarchical structures that spun angstrom to hundreds of nanometer scales. To maintain these unique structural features and to directly utilize these natural supramolecular assemblies, a variety of new methods have been developed to produce biopolymer nanofibrils. In particular, cellulose nanofibrils (CNFs), chitin nanofibrils (ChNFs), silk nanofibrils (SNFs) and collagen nanofibrils (CoNFs), as the four most abundant biopolymer nanofibrils on earth, have been the focus of research in recent years due to their renewable features, wide availability, low-cost, biocompatibility, and biodegradability. A series of top-down and bottom-up strategies have been accessed to exfoliate and regenerate these nanofibrils for versatile advanced applications. In this review, we first summarize the structures of biopolymer nanofibrils in nature and outline their related computational models with the aim of disclosing fundamental structure-property relationships in biological materials. Then, we discuss the underlying methods used for the preparation of CNFs, ChNFs, SNF and CoNFs, and discuss emerging applications for these biopolymer nanofibrils.

Statistical thermodynamics for the unexpectedly large difference between disaccharide stereoisomers in terms of solubility in water

We unravel the physical origins of the large difference between cellobiose and maltose, which consist of two β-1,4 and α-1,4 linked d-glucose units, respectively, in terms of the solubility in water. We construct a thermodynamic theory where the chemical-potential difference between disaccharides in water and in vacuum is identified as the key free-energy function. Its energetic and entropic components are calculated for cellobiose and maltose by statistical-mechanical theories for solute hydration. The disaccharide structures are taken into account at the atomic level and a molecular model is adopted for water. Molecular dynamics simulations are used to account for the conformational fluctuation of a disaccharide molecule, which also enables us to estimate the conformational entropy. We show that the cellobiose/maltose solubility ratio calculated is in good agreement with the experimental value. The solubility becomes much lower for cellobiose due to conformational-entropy and water-entropy effects. The former effect is relevant to higher stability of the intramolecular hydrogen bond between oxygen atoms in the six-membered ring and in the neighboring hydroxyl group: the hydration alters the fluctuation of a molecular conformation to a larger or less regular one, but the degree of this alteration is smaller. The latter effect is attributed to more separation of two hydroxymethyl groups in a molecule, causing lower probability of the overlap of excluded volumes generated by the groups for water molecules. We suggest that physicochemical properties of disaccharides in water become variable depending on the stereoisomerism through hydration effects and the origins of the variety are entropic.