doGlycans-Tools for Preparing Carbohydrate Structures for Atomistic Simulations of Glycoproteins, Glycolipids, and Carbohydrate Polymers for GROMACS.

Carbohydrates constitute a structurally and functionally diverse group of biological molecules and macromolecules. In cells they are involved in, e.g., energy storage, signaling, and cell-cell recognition. All of these phenomena take place in atomistic scales, thus atomistic simulation would be the method of choice to explore how carbohydrates function. However, the progress in the field is limited by the lack of appropriate tools for preparing carbohydrate structures and related topology files for the simulation models. Here we present tools that fill this gap. Applications where the tools discussed in this paper are particularly useful include, among others, the preparation of structures for glycolipids, nanocellulose, and glycans linked to glycoproteins. The molecular structures and simulation files generated by the tools are compatible with GROMACS.

Glycosylation and Lipids Working in Concert Direct CD2 Ectodomain Orientation and Presentation.

Proteins embedded in the plasma membrane mediate interactions with the cell environment and play decisive roles in many signaling events. For cell-cell recognition molecules, it is highly likely that their structures and behavior have been optimized in ways that overcome the limitations of membrane tethering. In particular, the ligand binding regions of these proteins likely need to be maximally exposed. Here we show by means of atomistic simulations of membrane-bound CD2, a small cell adhesion receptor expressed by human T-cells and natural killer cells, that the presentation of its ectodomain is highly dependent on membrane lipids and receptor glycosylation acting in apparent unison. Detailed analysis shows that the underlying mechanism is based on electrostatic interactions complemented by steric interactions between glycans in the protein and the membrane surface. The findings are significant for understanding the factors that render membrane receptors accessible for binding and signaling.

Design of cholesterol arabinogalactan anchored liposomes for asialoglycoprotein receptor mediated targeting to hepatocellular carcinoma: In silico modeling, in vitro and in vivo evaluation.

We have developed active targeting liposomes to deliver anticancer agents to ASGPR which will contribute to effective treatment of hepatocellular carcinoma. Active targeting is achieved through polymeric ligands on the liposome surface. The liposomes were prepared using reverse phase evaporation method and doxorubicin hydrocholoride, a model drug, was loaded using the ammonium sulphate gradient method. Liposomes loaded with DOX were found to have a particle size of 200nm with more than 90% entrapment efficiency. Systems were observed to release the drug in a sustained manner in acidic pH in vitro. Liposomes containing targeting ligands possessed greater and selective toxicity to ASGPR positive HepG2 cell lines due to specific ligand receptor interaction. Bio-distribution studies revealed that liposomes were concentrated in the liver even after 3h of administration, thus providing conclusive evidence of targeting potential for formulated nanosystems. Tumor regression studies indicated greater tumor suppression with targeted liposomes thereby establishing superiority of the liposomal system. In this work, we used a novel methodology to guide the determination of the optimal composition of the targeting liposomes: molecular dynamics (MD) simulation that aided our understanding of the behaviour of the ligand within the bilayer. This can be seen as a demonstration of the utility of this methodology as a rational design tool for active targeting liposome formulation.

N-Glycosylation as determinant of epidermal growth factor receptor conformation in membranes.

The epidermal growth factor receptor (EGFR) regulates several critical cellular processes and is an important target for cancer therapy. In lieu of a crystallographic structure of the complete receptor, atomistic molecular dynamics (MD) simulations have recently shown that they can excel in studies of the full-length receptor. Here we present atomistic MD simulations of the monomeric N-glycosylated human EGFR in biomimetic lipid bilayers that are, in parallel, also used for the reconstitution of full-length receptors. This combination enabled us to experimentally validate our simulations, using ligand binding assays and antibodies to monitor the conformational properties of the receptor reconstituted into membranes. We find that N-glycosylation is a critical determinant of EGFR conformation, and specifically the orientation of the EGFR ectodomain relative to the membrane. In the absence of a structure for full-length, posttranslationally modified membrane receptors, our approach offers new means to structurally define and experimentally validate functional properties of cell surface receptors in biomimetic membrane environments.

Molecular dynamics, crystallography and mutagenesis studies on the substrate gating mechanism of prolyl oligopeptidase.

Altered prolyl oligopeptidase (PREP) activity is found in many common neurological and other genetic disorders, and in some cases PREP inhibition may be a promising treatment. The active site of PREP resides in an internal cavity; in addition to the direct interaction between active site and substrate or inhibitor, the pathway to reach the active site (the gating mechanism) must be understood for more rational inhibitor design and understanding PREP function. The gating mechanism of PREP has been investigated through molecular dynamics (MD) simulation combined with crystallographic and mutagenesis studies. The MD results indicate the inter-domain loop structure, comprised of 3 loops at residues, 189-209 (loop A), 577-608 (loop B), and 636-646 (loop C) (porcine PREP numbering), are important components of the gating mechanism. The results from enzyme kinetics of PREP variants also support this hypothesis: When loop A is (1) locked to loop B through a disulphide bridge, all enzyme activity is halted, (2) nicked, enzyme activity is increased, and (3) removed, enzyme activity is only reduced. Limited proteolysis study also supports the hypothesis of a loop A driven gating mechanism. The MD results show a stable network of H-bonds that hold the two protein domains together. Crystallographic study indicates that a set of known PREP inhibitors inhabit a common binding conformation, and this H-bond network is not significantly altered. Thus the domain separation, seen to occur in lower taxa, is not involved in the gating mechanism for mammalian PREP. In two of the MD simulations we observed a conformational change that involved the breaking of the H-bond network holding loops A and B together. We also found that this network was more stable when the active site was occupied, thus decreasing the likelihood of this transition.

Study of PEGylated lipid layers as a model for PEGylated liposome surfaces: molecular dynamics simulation and Langmuir monolayer studies.

We have combined Langmuir monolayer film experiments and all-atom molecular dynamics (MD) simulation of a bilayer to study the surface structure of a PEGylated liposome and its interaction with the ionic environment present under physiological conditions. Lipids that form both gel and liquid-crystalline membranes have been used in our study. By varying the salt concentration in the Langmuir film experiment and including salt at the physiological level in the simulation, we have studied the effect of salt ions present in the blood plasma on the structure of the poly(ethylene glycol) (PEG) layer. We have also studied the interaction between the PEG layer and the lipid bilayer in both the liquid-crystalline and gel states. The MD simulation shows two clear results: (a) The Na(+) ions form close interactions with the PEG oxygens, with the PEG chains forming loops around them and (b) PEG penetrates the lipid core of the membrane for the case of a liquid-crystalline membrane but is excluded from the tighter structure of the gel membrane. The Langmuir monolayer results indicate that the salt concentration affects the PEGylated lipid system, and these results can be interpreted in a fashion that is in agreement with the results of our MD simulation. We conclude that the currently accepted picture of the PEG surface layer acting as a generic neutral hydrophilic polymer entirely outside the membrane, with its effect explained through steric interactions, is not sufficient. The phenomena we have observed may affect both the interaction between the liposome and bloodstream proteins and the liquid-crystalline-gel transition and is thus relevant to nanotechnological drug delivery device design.