Elucidation of the LewisX–LewisX carbohydrate interaction with molecular dynamics simulations: A glycosynapse model

Theoretical and Experimental Elucidation of the Adsorption Process of a Bioinspired Peptide on Mineral Surfaces

Inorganic materials used for biomedical applications such as implants generally induce the adsorption of proteins on their surface. To control this phenomenon, the bioinspired peptidomimetic polymer 1 (PMP1), which aims to reproduce the adhesion of mussel foot proteins, is commonly used to graft specific proteins on various surfaces and to regulate the interfacial mechanism. To date and despite its wide application, the elucidation at the atomic scale of the PMP1 mechanism of adsorption on surfaces is still unknown. The purpose of the present work was thus to unravel this process through experimental and computational investigations of adsorption of PMP1 on gold, TiO₂, and SiO₂ surfaces. A common mechanism of adsorption is identified for the adsorption of PMP1 which emphasizes the role of electrostatics to approach the peptide onto the surface followed by a full adhesion process where the entropic desolvation step plays a key role. Besides, according to the fact that mussel naturally controls the oxidation states of its proteins, further investigations were performed for two distinct redox states of PMP1, and we conclude that even if both states are able to allow interaction of PMP1 with the surfaces, the oxidation of PMP1 leads to a stronger interaction.

Computational Studies of a DNA-Based Aptasensor: toward Theory-Driven Transduction Improvement

Aptamers are a class of bioreceptors intensively used in current analytical tools dedicated to molecular diagnostics due to their ability to perform large structural reorganization upon target binding. However, there is a lack of methodologies allowing us to rationalize their structure in order to improve the transduction efficiency in aptamer sensors. We choose here, as a model system, a three-strand DNA structure as the probe, composed of two DNA strands anchored on a gold surface and partially hybridized with an aptamer sequence sensitive to ampicillin (AMP). The DNA structure has been designed to show strong structural change upon AMP binding to its aptamer. Using a set of computational techniques including molecular dynamics simulations, we deeply investigated the structure change upon analyte binding, taking into account the grafting on the surface. Original analyses of ion distributions along the trajectories unveil a distinct pattern between both states which can be related to changes in capacitance of the interface between these states. To our knowledge, this work demonstrates the ability of computational investigations for the first time to drive, in silico, the design of aptasensors.

Glycan-to-Glycan Binding: Molecular Recognition through Polyvalent Interactions Mediates Specific Cell Adhesion

Glycan-to-glycan binding was shown by biochemical and biophysical measurements to mediate xenogeneic self-recognition and adhesion in sponges, stage-specific cell compaction in mice embryos, and in vitro tumor cell adhesion in mammals. This intermolecular recognition process is accepted as the new paradigm accompanying high-affinity and low valent protein-to-protein and protein-to-glycan binding in cellular interactions. Glycan structures in sponges have novel species-specific sequences. Their common features are the large size >100 kD, polyvalency >100 repeats of the specific self-binding oligosaccharide, the presence of fucose, and sulfated and/or pyruvylated hexoses. These structural and functional properties, different from glycosaminoglycans, inspired their classification under the glyconectin name. The molecular mechanism underlying homophilic glyconectin-to-glyconectin binding relies on highly polyvalent, strong, and structure-specific interactions of small oligosaccharide motifs, possessing ultra-weak self-binding strength and affinity. Glyconectin localization at the glycocalyx outermost cell surface layer suggests their role in the initial recognition and adhesion event during the complex and multistep process. In mammals, Le^x-to-Le^x homophilic binding is structure-specific and has ultra-weak affinity. Cell adhesion is achieved through highly polyvalent interactions, enabled by clustering of small low valent structure in plasma membranes.

An in-silico layer-by-layer adsorption study of the interaction between Rebaudioside A and the T1R2 human sweet taste receptor: modelling and biosensing perspectives

The human sweet taste receptor (T1R2) monomer-a member of the G-protein coupled receptor family that detects a wide variety of chemically and structurally diverse sweet tasting molecules, is known to pose a significant threat to human health. Protein that lack crystal structure is a challenge in structure-based protein design. This study focused on the interaction of the T1R2 monomer with rebaudioside A (Reb-A), a steviol glycoside with potential use as a natural sweetener using in-silico and biosensing methods. Herein, homology modelling, docking studies, and molecular dynamics simulations were applied to elucidate the interaction between Reb-A and the T1R2 monomer. In addition, the electrochemical sensing of the immobilised T1R2-Reb-A complex with zinc oxide nanoparticles (ZnONPs) and graphene oxide (GO) were assessed by testing the performance of multiwalled carbon nanotube (MWCNT) as an adsorbent experimentally. Results indicate a strong interaction between Reb-A and the T1R2 receptor, revealing the stabilizing interaction of the amino acids with the Reb-A by hydrogen bonds with the hydroxyl groups of the glucose moieties, along with a significant amount of hydrophobic interactions. Moreover, the presence of the MWCNT as an anchor confirms the adsorption strength of the T1R2-Reb-A complex onto the GO nanocomposite and supported with electrochemical measurements. Overall, this study could serve as a cornerstone in the development of electrochemical immunosensor for the detection of Reb-A, with applications in the food industry.

Weak carbohydrate-carbohydrate interactions in membrane adhesion are fuzzy and generic

Carbohydrates such as the trisaccharide motif LeX are key constituents of cell surfaces. Despite intense research, the interactions between carbohydrates of apposing cells or membranes are not well understood. In this article, we investigate carbohydrate-carbohydrate interactions in membrane adhesion as well as in solution with extensive atomistic molecular dynamics simulations that exceed the simulation times of previous studies by orders of magnitude. For LeX, we obtain association constants of soluble carbohydrates, adhesion energies of lipid-anchored carbohydrates, and maximally sustained forces of carbohydrate complexes in membrane adhesion that are in good agreement with experimental results in the literature. Our simulations thus appear to provide a realistic, detailed picture of LeX-LeX interactions in solution and during membrane adhesion. In this picture, the LeX-LeX interactions are fuzzy, i.e. LeX pairs interact in a large variety of short-lived, bound conformations. For the synthetic tetrasaccharide Lac 2, which is composed of two lactose units, we observe similarly fuzzy interactions and obtain association constants of both soluble and lipid-anchored variants that are comparable to the corresponding association constants of LeX. The fuzzy, weak carbohydrate-carbohydrate interactions quantified in our simulations thus appear to be a generic feature of small, neutral carbohydrates such as LeX and Lac 2.

In Vitro interaction between yeast frataxin and superoxide dismutases: Influence of mitochondrial metals

BACKGROUND

Friedreich's ataxia results from a decreased expression of the nuclear gene encoding the mitochondrial protein, frataxin. Frataxin participates in the biosynthesis of iron-sulfur clusters and heme cofactors, as well as in iron storage and protection against oxidative stress. How frataxin interacts with the antioxidant defence components is poorly understood.

METHODS

Therefore, we have investigated by kinetic, thermodynamic and modelling approaches the molecular interactions between yeast frataxin (Yfh1) and superoxide dismutases, Sod1 and Sod2, and the influence of Yfh1 on their enzymatic activities.

RESULTS

Yfh1 interacts with cytosolic Sod1 with a dissociation constant, Kd = 1.3 ± 0.3 μM, in two kinetic steps. The first step occurs in the 200 ms range and corresponds to the Yfh1-Sod1 interaction, whereas the second is slow and is assumed to be a change in the conformation of the protein-protein adduct. Furthermore, computational investigations confirm the stability of the Yfh1-Sod1 complex. Yfh1 forms two protein complexes with mitochondrial Sod2 with 1:1 and 2:1 Yfh1/Sod2 stoichiometry (Kd1 = 1.05 ± 0.05 and Kd2 = 6.6 ± 0.1 μM). Furthermore, Yfh1 increases the enzymatic activity of Sod1 while slightly affecting that of Sod2. Finally, the stabilities of the protein-protein adducts and the effect of Yfh1 on superoxide dismutase activities depend on the nature of the mitochondrial metal.

CONCLUSIONS

This work confirms the participation of Yfh1 in cellular defence against oxidative stress.

Molecular Dynamics Simulation of a RNA Aptasensor

Single-stranded RNA aptamers have emerged as novel biosensor tools. However, the immobilization procedure of the aptamer onto a surface generally induces a loss of affinity. To understand this molecular process, we conducted a complete simulation study for the Flavin mononucleotide aptamer for which experimental data are available. Several molecular dynamics simulations (MD) of the Flavin in complex with its RNA aptamer were conducted in solution, linked with six thymidines (T6) and, finally, immobilized on an hexanol-thiol-functionalized gold surface. First, we demonstrated that our MD computations were able to reproduce the experimental solution structure and to provide a meaningful estimation of the Flavin free energy of binding. We also demonstrated that the T6 linkage, by itself, does not generate a perturbation of the Flavin recognition process. From the simulation of the complete biosensor system, we observed that the aptamer stays oriented parallel to the surface at a distance around 36 Å avoiding, this way, interaction with the surface. We evidenced a structural reorganization of the Flavin aptamer binding mode related to the loss of affinity and induced by an anisotropic distribution of sodium cationic densities. This means that ionic diffusion is different between the surface and the aptamer than above this last one. We suggest that these findings might be extrapolated to other nucleic acids systems for the future design of biosensors with higher efficiency and selectivity.

The impact of active site mutations of South African HIV PR on drug resistance: Insight from molecular dynamics simulations, binding free energy and per-residue footprints

Molecular dynamics simulations and binding free energy calculations were used to provide an understanding of the impact of active site drug-resistant mutations of the South African HIV protease subtype C (C-SA HIV PR), V82A and V82F/I84V on drug resistance. Unique per-residue interaction energy 'footprints' were developed to map the overall drug-binding profiles for the wild type and mutants. Results confirmed that these mutations altered the overall binding landscape of the amino acid residues not only in the active site region but also in the flaps as well. Four FDA-approved drugs were investigated in this study; these include ritonavir (RTV), saquinavir (SQV), indinavir (IDV), and nelfinavir (NFV). Computational results compared against experimental findings were found to be complementary. Against the V82F/I84V variant, saquinavir, indinavir, and nelfinavir lose remarkable entropic contributions relative to both wild-type and V82A C-SA HIV PRs. The per-residue energy 'footprints' and the analysis of ligand-receptor interactions for the drug complexes with the wild type and mutants have also highlighted the nature of drug interactions. The data presented in this study will prove useful in the design of more potent inhibitors effective against drug-resistant HIV strains.

Studies on the binding modes of Lassa nucleoprotein complexed with m7GpppG and dTTP by molecular dynamic simulations and free energy calculations

Lassa virus can cause dreadful human hemorrhagic disease, for which there is no effective therapy. A recent study points out that the amino (N)-terminal domain of Lassa virus nucleoprotein (NP) plays an important role in viral RNA synthesis and firstly solved the X-ray crystal structures of NP complexed with the capped Deoxythymidine triphosphate (dTTP) analog, but the binding mode of m7GpppG to the N domain of NP, which is required for viral RNA transcription, has not been studied. In this study, molecular dynamics (MD) simulations have been carried out to investigate the characters of dTTP binding to two forms of NP, i.e. the NP without the C domain and the full-length NP model, using two different force fields, ff03 and ff99SB, respectively. Our calculated results show that the truncated model is reasonable and can replace the full protein model in the following MD simulations, and that ff99SB combined with the general AMBER force field is more suitable for sampling the structure of small molecule NP complex. From the comparisons of stability of hydrogen bonds between small molecule and protein in the dTTP and Uridine 5'-Triphosphate complexes, one finds that the stable hydrogen bonds between the second phosphate group of small molecules and two residues, Thr178 and Arg323, are critical for cap analogs binding to the N domain of NP. Additionally, docking method combined with MD simulations have been applied to predict the binding mode of m7GpppG to NP; and the hydrogen bond analysis and the binding free energy decomposition method (MM/GBSA) are conducted to study the interactions in the putative binding mode. The calculated results are expected to provide guidance for drug development.

Conformational Search on the Lewis X Structure by Molecular Dynamic: Study of Tri- and Pentasaccharide

Carbohydrates play vital roles in many biological processes, such as recognition, adhesion, and signalling between cells. The Lewis X determinant is a trisaccharide fragment implicated as a specific differentiation antigen, tumor antigen, and key component of the ligand for the endothelial leukocyte adhesion molecule, so it is necessary or essential to determine and to know their conformational and structural properties. In this work, conformational analysis was performed using molecular dynamics (MD) simulation with the AMBER10 program package in order to study the dynamic behavior of of the Lewis X trisaccharide (β-D-Gal-(1,4)-[α-L-Fuc-(1,3)]-β-D-GlcNAc-OMe) and the Lewis X pentasaccharide (β-D-Gal-(1,4)-[α-L-Fuc-(1,3)]-β-D-GlcNAc-(1,3)-β-D-Gal-(1,4)-β-D-Glu-OMe) in explicit water model at 300 K for 10 ns using the GLYCAM 06 force field.

Is inhibition process better described with MD(QM/MM) simulations? The case of urokinase type plasminogen activator inhibitors

Urokinase plasminogen activator (uPA) is an enzyme involved in cancer growth and metastasis. Therefore, the design of inhibitors of uPA is of high therapeutic value, and several chemical families have been explored, even if none has still emerged, emphasizing the need of a rationalized approach. This work represents a complete computational study of uPA complexed with five inhibitors, which present weak similarities. Molecular dynamics simulations in explicit solvent were conducted, and structural analyses, along with molecular mechanics (MM)/Poisson-Boltzmann surface area free energies estimations, yield precious structure-activity relationships of these inhibitors. Besides, we realized supplemental QM/MM computations that improved drastically the quality of our models providing original information on the hydrogen bonds and charge transfer effects, which are, most often, neglected in other studies. We suggest that these simulations and analyses could be reproduced for other systems involving protein/ligand molecular recognitions.

Interacting Lewis-X carbohydrates in condensed phase: a first-principles molecular dynamics study

We performed first-principles molecular dynamics calculations at finite temperature, to study the interacting conformations of Lewis-X (LeX) trisaccharides in the crystalline phase. The calculated cell parameters and detailed atomic structure of the LeX molecule compare well to the experimental data obtained by X-ray diffraction. We identify and characterize the hydrogen-bond network, responsible for the mutual interaction of the LeX pairs, whereas we find the intramolecular conformation and stability to be mainly assured by dispersion forces. The relative contributions to the crystallization energy of the hydrogen bonds and of the dispersion forces are defined and quantified. From this study, candidate configurations for the fully hydrated, in vivo structures of homotypic LeX-LeX interactions at cell surfaces can be proposed. We discuss how these configurations could also be relevant for the adhesion and self-assembly of nanostructures.

Membrane adhesion via homophilic saccharide-saccharide interactions investigated by neutron scattering

Solid-supported membrane multilayers doped with membrane-anchored oligosaccharides bearing the LewisX motif (Le(X) lipid) were utilized as a model system of membrane adhesion mediated via homophilic carbohydrate-carbohydrate interactions. Specular and off-specular neutron scattering in bulk aqueous electrolytes allowed us to study multilayer structure and membrane mechanics at full hydration at various Ca(2+) concentrations, indicating that membrane-anchored Le(X) cross-links the adjacent membranes. To estimate forces and energies required for cross-linking, we theoretically modeled the interactions between phospholipid membranes and compared this model with our experimental results on membranes doped with Le(X) lipids. We demonstrated that the bending rigidity, extracted from the off-specular scattering signals, is not significantly influenced by the molar fraction of Le(X) lipids, while the vertical compression modulus (and thus the intermembrane confinement) increases with the molar fraction of Le(X) lipids.

Prediction of zanamivir efficiency over the possible 2009 influenza A (H1N1) mutants by multiple molecular dynamics simulations and free energy calculations

As one of the most important antiviral drugs against 2009 influenza A (H1N1), will zanamivir be effective for the possible drug resistant mutants? To answer this question, we combined multiple molecular dynamics simulations and molecular mechanics generalized Born surface area (MM-GBSA) calculations to study the efficiency of zanamivir over the most frequent drug-resistant strains of neuraminidase including R293K, R152K, E119A/D and H275Y mutants. The calculated results indicate that the modeled mutants of the 2009-H1N1 strains except H275Y will be significantly resistant to zanamivir. The resistance to zanamivir is mainly caused by the loss of polar interactions. The identified potential resistance sites in this study will be useful for the development of new effective anti-influenza drugs and to avoid the occurrence of the state without effective drugs to new mutant influenza strains.

In silico identification of the potential drug resistance sites over 2009 influenza A (H1N1) virus neuraminidase

The outbreak and high speed global spread of the new strain of influenza A (H1N1) virus in 2009 poses a serious threat to the general population and governments. At present, the most effective drugs for the treatment of 2009 influenza A (H1N1) virus are neuraminidase inhibitors: mainly oseltamivir and zanamivir. The use of these two inhibitors will undoubtedly increase, and therefore it is more likely that drug-resistant influenza strains will arise. The identification of the potential resistance sites for these drugs in advance and the understanding of corresponding molecular basis to cause drug resistance are no doubt very important to fight against the new resistant influenza strains. In this study, first, the complexes of neuraminidase with the substrate sialic acid and two inhibitors oseltamivir and zanamivir were obtained by fitting them to the 3D structure of 2009 influenza A (H1N1) neuraminidase obtained by homology modeling. By using these complexes as the initial structures, molecular dynamics simulation and molecular mechanics generalized Born surface area (MM-GBSA) calculations were performed to identify the residues with significant contribution to the binding of substrate and inhibitors. By analyzing the difference of interaction profiles of substrate and inhibitors, the potential drug resistance sites for two inhibitors were identified. Parts of the identified sites have been verified to confer resistance to oseltamivir and zanamivir for influenza virus of the past flu epidemic. The identified potential resistance sites in this study will be useful for the development of new effective drugs against the drug resistance and avoid the situation of having no effective drugs to treat new mutant influenza strains.