Interaction of peptides, peptidomimetics and other drug candidates with the DI/tripeptide transport system in intestinal epithelial cells, using the in vitro caco-2 cell culture system

An uptake system was developed using Caco-2 cell monolayers and the dipeptide, glycyl-[3H]L-proline, as a probe compound. Glycyl-[3H]L-proline uptake was via the di-/tripeptide transport system (DTS) and, exhibited concentration-, pH- and temperature-dependency. Dipeptides inhibited uptake of the probe, and the design of the system allowed competitors to be ranked against one another with respect to affinity for the transporter. The structural features required to ensure or increase interaction with the DTS were defined by studying the effect of a series of glycyl-L-proline and angiotensin-converting enzyme (ACE)-inhibitor (SQ-29852) analogues on the uptake of the probe. The SQ-29852 structure was divided into six domains (A-F) and competitors were grouped into series depending on structural variations within specific regions. Domain A was found to prefer a hydrophobic function, such as a phenyl group, and was intolerant to positive charges and H+ -acceptors and donors. SQ-29852 analogues were more tolerant of substitutions in the C domain, compared to glycyl-L-proline analogues, suggesting that interactions along the length of the SQ-29852 molecule may override the effects of substitutions in the C domain. SQ-29852 analogues showed a preference for a positive function, such as an amine group in this region, but dipeptide structures favoured an uncharged substitution. Lipophilic substituents in domain D increased affinity of SQ-29852 analogues with the DTS. A similar effect was observed for ACE-NEP inhibitor analogues. Domain E, corresponding to the carboxyl group was found to be tolerant of esterification for SQ-29852 analogues but not for dipeptides. Structural features which may increase interaction for one series of compounds, may not have the same effect for another series, indicating that the presence of multiple recognition sites on a molecule may override the deleterious effect of anyone change. Modifying current, poorly absorbed peptidomimetic structures to fit the proposed hypothetical model may improve oral bioavailability by increasing affinity for the DTS. The stereochemical preference of the transporter was explored using four series of compounds (SQ-29852, lysylproline, alanylproline and alanylalanine enantiomers). The L, L stereochemistry was the preferred conformation for all four series, agreeing with previous studies. However, D, D enantiomers were shown in some cases to be substrates for the DTS, although exhibiting a lower affinity than their L, L counterparts. All the ACE-inhibitors and β-lactam antibiotics investigated, produced a degree of inhibition of the probe, and thus show some affinity for the DTS. This contrasts with previous reports that found several ACE inhibitors to be absorbed via a passive process, thus suggesting that compounds are capable of binding to the transporter site and inhibiting the probe without being translocated into the cell. This was also shown to be the case for oligodeoxynucleotide conjugated to a lipophilic group (vitamin E), and highlights the possibility that other orally administered drug candidates may exert non-specific effects on the DTS and possibly have a nutritional impact. Molecular modelling of selected ACE-NEP inhibitors revealed that the three carbonyl functions can be oriented in a similar direction, and this conformation was found to exist in a local energy-minimised state, indicating that the carbonyls may possibly be involved in hydrogen-bond formation with the binding site of the DTS.

Studies on the mode of action and beta-lactamase inhibitory properties of novel oxapenem compounds

Four novel oxapenem compounds were evaluated for their ß-lactamase inhibitory and antibacterial properties. Two (AM-112 and AM-113) displayed intrinsic antibacterial activity with MICs of between 2 to 16µg/ml and 0.5-2µg/ml against Escherichia coli and methicillin-sensitive and -resistant Staphylococcus aureus, respectively. The isomers of these compounds, AM-115 and AM-114 did not display significant antibacterial activity. Combination of the oxapenems with ceftazidime afforded protection against ß-lactamase-producing strains, including hyperproducers of class C enzymes and extended-spectrum ß-lactamase enzymes. A fixed 4µg/ml concentration of AM-112 protected a panel of eight cephalosporins against hydrolysis by class A and class C ß-lactamase producers. In vivo studies confirmed the protective effect of AM-112 for ceftazidime against ß-lactamase producing S. aureus, Enterobacter cloacae and E. coli strains in a murine intraperitoneal infection model. Each of the oxapenems inhibited class A, class C and class D ß-lactamases isolated from whole cells and purified by isoelectric focusing. AM-114 and AM-115 were as effective as clavulanic acid against class A enzymes. AM-112 and AM-113 were less potent against these enzymes. Class C and class D enzymes proved very susceptible to inhibition by the oxapenems. Molecular modelling of the oxapenems in the active site of the class A. TEM-1 and class C P99 enzymes identified a number of potential sites of interaction. The modelling suggested that Ser-130 in TEM-1 and Tyr-150 in P99 were likely candidates for cross-linking of the inhibitor, leading to inhibition of the enzyme. Morphology studies indicated that sub-inhibitory concentrations of the oxapenems caused the formation of round-shaped cells in E. coli DC0, indicating inhibition of penicillin-binding protein 2 (PBP2). The PBP affinity profile of AM-112 was examined in isolated cell membranes of E. coli DC0, S. aureus NCTC 6571, Enterococcus faecalis SFZ and E. faecalis ATCC 29213, in competition with a radiolabelled penicillin. PBP2 was identified as the primary target for AM-112 in E. coli DC0. Studies on S. aureus NCTC 6571 failed to identify a binding target. AM-112 bound to all the PBPs of both E. faecalis strains, and a concentration of 10µg/ml inhibited all the PBPs except PBP3.

Automated synthesis and evaluation of potential new anti-microbial agents

A series of N1-benzylideneheteroarylcarboxamidrazones was prepared in an automated fashion, and tested against Mycobacterium fortuitum in a rapid screen for antimycobacterial activity. Many of the compounds from this series were also tested against Mycobacterium tuberculosis, and the usefulness as M.fortuitum as a rapid, initial screen for anti-tubercular activity evaluated. Various deletions were made to the N1-benzylideneheteroarylcarboxamidrazone structure in order to establish the minimum structural requirements for activity. The N1-benzylideneheteroarylcarbox-amidrazones were then subjected to molecular modelling studies and their activities against M.fortuitum and M.tuberculosis were analysed using quantitative structure-analysis relationship (QSAR) techniques in the computational package TSAR (Oxford Molecular Ltd.). A set of equations predictive of antimycobacterial activity was hereby obtained. The series of N1-benzylidenehetero-arylcarboxamidrazones was also tested against a multidrug-resistant strain of Staphylococcus aureus (MRSA), followed by a panel of Gram-positive and Gram-negative bacteria, if activity was observed for MRSA. A set of antimycobacterial N1-benzylideneheteroarylcarboxamidrazones was hereby discovered, the best of which had MICs against m. fortuitum in the range 4-8μgml-1 and displayed 94% inhibition of M.tuberculosis at a concentration of 6.25μgml-1. The antimycobacterial activity of these compounds appeared to be specific, since the same compounds were shown to be inactive against other classes of organisms. Compounds which were found to be sufficiently active in any screen were also tested for their toxicity against human mononuclear leucocytes. Polyethylene glycol (PEG) was used as a soluble polymeric support for the synthesis of some fatty acid derivatives, containing an isoxazoline group, which may inhibit mycolic acid synthesis in mycobacteria. Both the PEG-bound products and the cleaved, isolated products themselves were tested against M.fortuitum and some low levels of antimycobacterial activity were observed, which may serve as lead compounds for further studies.

Synthesis and biological properties of 2-phenylbenzothiazole derivatives

2-Phenylbenzothiazoles have structural similarities to the antioestrogenic 2-phenylindole, zindoxifene and to the oestrogenic isoflavone, genistein which also inhibits tyrosine kinases. Hydroxylated 2-phenylbenzothiazole derivatives were therefore produced and tested for oestrogenic and tyrosine kinase inhibitory activity. Synthesis of methoxy substituted 2-phenylbenzothiazoles was via the Jacobson method, demethylation being effected by boron tribromide at -70oC. Three amino substituted 2-phenylbenzothiazoles were also synthesised and tested for activity. Data is presented for oestrogen receptor binding activity, aromatase inhibitory activity, epidermal growth factor receptor tyrosine kinase (EGFRTK) inhibitory activity and cytotoxicity to ANN-1, 3T3, MCF-7 and WIDR cells. Oestrogen receptor binding affinity (RBA) was shown by five of the nine compounds tested. 2-(4-hydroxy)-6-hydroxybenzo-thiazole was the most active of the benzothiazoles tested (RBA 0.7). This is low but comparable to that of genistein. EGFRTK inhibitory activity was shown by four of the six benzothiazole derivatives tested; activity was comparable to that of genistein. Cytotoxicity assays have shown no selective toxicity of 2-phenylbenzothiazoles to any of the cell lines tested. Toxicity to MCF-7 cells was similar to that for other cell lines despite some compounds showing oestrogen receptor binding capacity. Amino-substituted 2-phenylbenzothiazoles showed selective toxicity towards transformed ANN-1 cells compared to normal 3T3 cells but the mechanism of this selectivity has not been established. Molecular modelling techniques, including CHEM-X, QUANTA and MOPAC were used to compare known ATP-competitive tyrosine kinase inhibitors with a model of ATP built from the crystal structure of the ATP-phosphoglycerate kinase complex. Structural features thought to be important to kinase inhibition were found and used to suggest further 2-phenylbenzothiazole analogues which may have improved activity.

Antitumour imidazotetrazines: probes for the major groove of DNA

The antitumour imidazotetrazinones are believed to act as prodrugs for the triazene series of alkylating agents, showing a marked pteference for the alkylation of the middle guanine residue in a run of three or more contiguous guanines. However, the. exact nature of the interactions of imidazotetrazinones within the micro~environment of DNA are; as yet unknown. In order to examine such interactions a three pronged approach involving molecular modelling, synthetic chemistry and biological analysis has been undertaken during the course of this project. . Molecular modelling studies have shown that for the 8-carboxamido substituted imidazotetrazinones antitumour activity is dependent upon the. presence of a free NH group which can be involved in the formation of both intramolecular and intermolecular hydrogen bonds, and the presence of a non-bulky substituent with a small negative potential . volume. Modelling studies involving the docking of .mitozolomide into the major groove of DNA in the region of a triguanine sequence has shown that a number of hydrogen bonding interactions are feasible. A series of 8-substituted carboxamide derivatives of mitozolomide have been synthesised via the 8-acid chloride and 8-carboxylic acid derivatives including a number of peptide analogues. The peptide derivatives were based upon the key structural features of the helix-turn-helix motif of DNA-binding proteins with a view to developing agents that are capable of binding to DNA with greater selectivity. An examination of the importance of intramolecular hydrogen bonding in influencing the antitumour activity:of :the imidazotetrazinones has led to the synthesis of the novel pyrimido[4',5' :4,3]pyrazolo[5,1-d]-1,2,3,5-tetrazine ring system. In general, in vitro cytotoxicity assays showed that the new derivatives were less active against the TLX5 lymphoma cell line. than the parent compound mitozolomide despite an increased potential for hydrogen bonding interactions. Due to the high reactivity of the: tetrazinone ring system it is difficult to study the interactions between the imidazotetrazinones and DNA. Consequently a number of structural analogues that are stable under physiological conditions have been. prepared based upon the 1,2,3 triazin-4(3H)-one ring system fused with both benzene and pyrazole rings. Although the 3-methylbenzotriazinones failed to antagonise the cytotoxic activity of temozolomide encouraging results with a 3-methylpyrazolotriazinone may suggest the existence of an imidazotetrazinone receptor site within DNA. The potential of guanine rich sequences to promote the alkylating selectivity of imidazotetrazinones by acting as a catalyst for ring cleavage and thereby generation of the alkylating agent was examined. Experiments involving the monitoring: of the rate of breakdown of mitozolomide incubated in the presence of synthetic oIigonucleotides did not reveal any catalytic effect resulting from the DNA. However, it was noted that the breakdown of mitozolomide was dependent upon the type of buffer used in the incubations and this may indeed mask any catalysis by the oligonucleotides.

DNA methylation at cytosine position 5

DNA methylation appears to be involved in the regulation of gene expression. Transcriptionally inactive (silenced) genes normally contain a high proportion of 5-methyl-2'-deoxycytosine residues whereas transcriptionally active genes show much reduced levels. There appears good reason to believe that chemical agents capable of methylating 2'-deoxycytosine might affect gene expression and as a result of hypermethylating promoter regions of cytosine-guanine rich oncogenic sequences, cancer related genes may be silenced. This thesis describes the synthesis of a number of `electrophilic' S-methylsulphonium compounds and assesses their ability to act as molecules capable of methylating cytosine at position 5 and also considers their potential as cytotoxic agents. DNA is methylated in vivo by DNA methyltransferase utilising S-adenoxylmethionine as the methyl donor. This thesis addresses the theory that S-adenoxylmethionine may be replaced as the methyl donor for DNA methytransferase by other sulphonium compounds. S-[3H-methyl]methionine sulphonium iodide was synthesised and experiments to assess the ability of this compounds to transfer methyl groups to cytosine in the presence of DNA methyltransferase were unsuccessful. A proline residue adjacent to a cysteine residue has been identified to a highly conserved feature of the active site region of a large number of prokaryotic DNA methyltransferases. The thesis examines the possibility that short peptides containing the Pro-Cys fragment may be able to facilitate the alkylation of cytosine position 5 by sulphonium compounds. Peptides were synthesised up to 9 amino acids in length but none were shown to exhibit significant activity. Molecular modelling techniques, including Chem-X, Quanta, BIPED and protein structure prediction programs were used to assess any structural similarities that may exist between short peptides containing a Pro-Cys fragment and similar sequences present in proteins. A number of similar structural features were observed.

The use of site-directed mutagenesis to study GPCRs

G protein coupled receptors (GPCRs) are highly flexible and dynamic proteins, which are able to interact with diverse ligands, effectors, and regulatory proteins. Site-directed mutagenesis (SDM) is a powerful tool for providing insight into how these proteins actually work, both in its own right and when used in conjunction with information provided by other techniques such as crystallography or molecular modelling. Mutagenesis has been used to identify and characterise a myriad of functionally important residues, motifs and domains within the GPCR architecture, and to identify aspects of similarity and differences between the major families of GPCRs. This chapter presents the necessary information for undertaking informative SDM of these proteins. Whilst this is relevant to protein structure/function studies in -general, specific pitfalls and protocols suited to investigating GPCRs in particular will be highlighted.

Toward the atomistic simulation of T cell epitopes automated construction of MHC:peptide structures for free energy calculations

Epitopes mediated by T cells lie at the heart of the adaptive immune response and form the essential nucleus of anti-tumour peptide or epitope-based vaccines. Antigenic T cell epitopes are mediated by major histocompatibility complex (MHC) molecules, which present them to T cell receptors. Calculating the affinity between a given MHC molecule and an antigenic peptide using experimental approaches is both difficult and time consuming, thus various computational methods have been developed for this purpose. A server has been developed to allow a structural approach to the problem by generating specific MHC:peptide complex structures and providing configuration files to run molecular modelling simulations upon them. A system has been produced which allows the automated construction of MHC:peptide structure files and the corresponding configuration files required to execute a molecular dynamics simulation using NAMD. The system has been made available through a web-based front end and stand-alone scripts. Previous attempts at structural prediction of MHC:peptide affinity have been limited due to the paucity of structures and the computational expense in running large scale molecular dynamics simulations. The MHCsim server (http://igrid-ext.cryst.bbk.ac.uk/MHCsim) allows the user to rapidly generate any desired MHC:peptide complex and will facilitate molecular modelling simulation of MHC complexes on an unprecedented scale.

Pluronic F127-modified liposome-containing tacrolimus-cyclodextrin inclusion complexes:improved solubility, cellular uptake and intestinal penetration

Objective The aim of this study was to investigate Pluronic F127-modified liposome-containing cyclodextrin (CD) inclusion complex (FLIC) for improving the solubility, cellular uptake and intestinal penetration of tacrolimus (FK 506) in the gastrointestinal (GI) tract. Methods Molecular modelling was performed to screen the optimal CD for the solubilization of FK 506. FLIC was prepared by thin-lipid film hydration with the inclusion complex solutions followed by membrane extrusion. Dilution tests were conducted in simulated gastric fluids and phosphate-buffered solution of sodium taurocholate to investigate the solubility improvement of FK506. The cellular uptake of nanocarriers was studied in Caco-2 cells, and intestinal mucous membrane penetration in the GI tract was evaluated in Sprague-Dawley rats. Key findings The results showed that β-CD had the strongest binding energy with the guest molecule among the CDs. The prepared FLIC has an average diameter of 180.8 ± 8.1 nm with a spherical shape. The solubility and cellular uptake of FK 506 was greatly improved by the incorporation of CD complexes in the Pluronic F127-modified liposomes. Intestinal mucous membrane penetration was also significantly improved by the preparation of FLIC. Conclusion With improved drug solubility and intestinal mucous membrane penetration, FLIC shows a promising oral delivery system for FK 506. © 2013 Royal Pharmaceutical Society.

The structure and dynamics of 2-dimensional fluids in swelling clays

The interlayer pores of swelling 2:1 clays provide an ideal 2-dimensional environment in which to study confined fluids. In this paper we discuss our understanding of the structure and dynamics of interlayer fluid species in expanded clays, based primarily on the outcome of recent molecular modelling and neutron scattering studies. Counterion solvation is compared with that measured in bulk solutions, and at a local level the cation-oxygen coordination is found to be remarkably similar in these two environments. However, for the monovalent ions the contribution to the first coordination shell from the clay surfaces increases with counterion radius. This gives rise to inner-sphere (surface) complexes in the case of potassium and caesium. In this context, the location of the negative clay surface charge (i.e. arising from octahedral or tetrahedral substitution) is also found to be of major importance. Divalent cations, such as calcium, eagerly solvate to form outer-sphere complexes. These complexes are able to pin adjacent clay layers together, and thereby prevent colloidal swelling. Confined water molecules form hydrogen bonds to each other and to the clays' surfaces. In this way their local environment relaxes to close to the bulk water structure within two molecular layers of the clay surface. Finally, we discuss the way in which the simple organic molecules methane, methanol and ethylene glycol behave in the interlayer region of hydrated clays. Quasi-elastic neutron scattering of isotopically labelled interlayer CH 3OD and (CH2OD)2 in deuterated clay allows us to measure the diffusion of the CH3- and CH2-groups in both clay and liquid environments. We find that in both the one-layer methanol solvates and the two-layer glycol solvates the diffusion of the most mobile organic molecules is close to that in the bulk solution.

Towards the in silico identification of class II restricted T-cell epitopes:a partial least squares iterative self-consistent algorithm for affinity prediction

Motivation: The immunogenicity of peptides depends on their ability to bind to MHC molecules. MHC binding affinity prediction methods can save significant amounts of experimental work. The class II MHC binding site is open at both ends, making epitope prediction difficult because of the multiple binding ability of long peptides. Results: An iterative self-consistent partial least squares (PLS)-based additive method was applied to a set of 66 pep- tides no longer than 16 amino acids, binding to DRB1*0401. A regression equation containing the quantitative contributions of the amino acids at each of the nine positions was generated. Its predictability was tested using two external test sets which gave r pred =0.593 and r pred=0.655, respectively. Furthermore, it was benchmarked using 25 known T-cell epitopes restricted by DRB1*0401 and we compared our results with four other online predictive methods. The additive method showed the best result finding 24 of the 25 T-cell epitopes. Availability: Peptides used in the study are available from http://www.jenner.ac.uk/JenPep. The PLS method is available commercially in the SYBYL molecular modelling software package. The final model for affinity prediction of peptides binding to DRB1*0401 molecule is available at http://www.jenner.ac.uk/MHCPred. Models developed for DRB1*0101 and DRB1*0701 also are available in MHC- Pred

An allosteric role for receptor activity-modifying proteins in defining GPCR pharmacology

G protein-coupled receptors are allosteric proteins that control transmission of external signals to regulate cellular response. Although agonist binding promotes canonical G protein signalling transmitted through conformational changes, G protein-coupled receptors also interact with other proteins. These include other G protein-coupled receptors, other receptors and channels, regulatory proteins and receptor-modifying proteins, notably receptor activity-modifying proteins (RAMPs). RAMPs have at least 11 G protein-coupled receptor partners, including many class B G protein-coupled receptors. Prototypic is the calcitonin receptor, with altered ligand specificity when co-expressed with RAMPs. To gain molecular insight into the consequences of this protein–protein interaction, we combined molecular modelling with mutagenesis of the calcitonin receptor extracellular domain, assessed in ligand binding and functional assays. Although some calcitonin receptor residues are universally important for peptide interactions (calcitonin, amylin and calcitonin gene-related peptide) in calcitonin receptor alone or with receptor activity-modifying protein, others have RAMP-dependent effects, whereby mutations decreased amylin/calcitonin gene-related peptide potency substantially only when RAMP was present. Remarkably, the key residues were completely conserved between calcitonin receptor and AMY receptors, and between subtypes of AMY receptor that have different ligand preferences. Mutations at the interface between calcitonin receptor and RAMP affected ligand pharmacology in a RAMP-dependent manner, suggesting that RAMP may allosterically influence the calcitonin receptor conformation. Supporting this, molecular dynamics simulations suggested that the calcitonin receptor extracellular N-terminal domain is more flexible in the presence of receptor activity-modifying protein 1. Thus, RAMPs may act in an allosteric manner to generate a spectrum of unique calcitonin receptor conformational states, explaining the pharmacological preferences of calcitonin receptor-RAMP complexes. This provides novel insight into our understanding of G protein-coupled receptor-protein interaction that is likely broadly applicable for this receptor class.

Design et synthèse d’inhibiteurs d’une ectonucléotide pyrophosphatase/phosphodiestérase de type 1 (ENPP1) et leur activité anticancéreuse

La calcification de la valve aortique (CVA) est une maladie cardiovasculaire de plus en plus répandue, particulièrement en Amérique du Nord. Elle cause le rétrécissement de la valve aortique et le seul traitement actuellement disponible est le remplacement chirurgical. Des études menées par le Dr Patrick Mathieu (Institut de Cardiologie et de Pneumologie de Québec) ont montré qu’une surexpression d’une ectonucléotide pyrophosphatase/phosphodiestérase de type 1 (ENPP1) est à l’origine de cette sténose. Une solution à cette maladie serait donc de trouver un inhibiteur d’ENPP1. Inspirées des travaux du groupe de Pfizer visant ENPP1 pour le traitement de la chondrocalcinose articulaire et l’ostéoarthrite, quelques familles d’inhibiteurs de type quinazoline-4-pipéridine sulfamides (QPS) ont été synthétisés et testées in vitro. Une étude en modélisation moléculaire sur le site potentiel de liaison des inhibiteurs sur ENPP1 est en cours, en collaboration avec le Pr Patrick Lagüe (Université Laval, Département de biochimie, microbiologie et bio-informatique) et son équipe pour optimiser le design de la structure des composés. Les composés d’une des familles, les QPS-pyrimidine, ont été testés in vitro sur quelques lignées cellulaires cancéreuses (HT-1080, HT-29, M21 et MCF-7) pour mesurer leur activité antiproliférative. Ces composés ont une inhibition de croissance médiane (IC50) de l’ordre du micromolaire et représentent donc un point de départ intéressant pour la mise au point de nouveaux traitements anticancéreux.

Besnoitia besnoiti protein disulfide isomerase (BbPDI): molecular characterization, expression and in silico modelling

Besnoitia besnoiti is an apicomplexan parasite responsible for bovine besnoitiosis, a disease with a high prevalence in tropical and subtropical regions and re-emerging in Europe. Despite the great economical losses associated with besnoitiosis, this disease has been underestimated and poorly studied, and neither an effective therapy nor an efficacious vaccine is available. Protein disulfide isomerase (PDI) is an essential enzyme for the acquisition of the correct three-dimensional structure of proteins. Current evidence suggests that in Neosporacaninum and Toxoplasmagondii, which are closely related to B. besnoiti, PDI play an important role in host cell invasion, is a relevant target for the host immune response, and represents a promising drug target and/or vaccine candidate. In this work, we present the nucleotide sequence of the B. besnoiti PDI gene. BbPDI belongs to the thioredoxin-like superfamily (cluster 00388) and is included in the PDI_a family (cluster defined cd02961) and the PDI_a_PDI_a'_c subfamily (cd02995). A 3D theoretical model was built by comparative homology using Swiss-Model server, using as a template the crystallographic deduced model of Tapasin-ERp57 (PDB code 3F8U chain C). Analysis of the phylogenetic tree for PDI within the phylum apicomplexa reinforces the close relationship among B. besnoiti, N. caninum and T. gondii. When subjected to a PDI-assay based on the polymerisation of reduced insulin, recombinant BbPDI expressed in E. coli exhibited enzymatic activity, which was inhibited by bacitracin. Antiserum directed against recombinant BbPDI reacted with PDI in Western blots and by immunofluorescence with B. besnoiti tachyzoites and bradyzoites.

Structural modelling of liquid ZrO2 and glassy B2O3 by using AB INITIO Molecular Dynamics method

Atomic structure of ZrO2 and B2O3 was investigated in this work. New data under extreme conditions (T = 3100 K) was obtained for the liquid ZrO2 structure. A fractional number of boron was investigated for glassy structure of B2O3. It was shown that it is possible to obtain an agreement for the fractional number between NMR and DFT techniques using a suitable initial configuration.