Investigation of the Nucleation and Growth of Colloidal Indium Phosphide: from Molecular Precursors to Semiconductor Nanocrystals through In37P20(O2CR)51 as a Magic-Size Intermediate

Thesis (Ph.D.)--University of Washington, 2016-06

Erythrocyte Phospholipid Molecular Species and Fatty Acids of Down Syndrome Children Compared with Non-affected Siblings

The majority of children with Down syndrome (DS) develop Alzheimer's disease (AD) at an early age. Although long-chain n-3 fatty acids (FA) are protective of neurodegeneration, little is known about the FA status in DS. In the present study, we aimed to investigate whether children with DS presented altered plasma and erythrocyte membrane phospholipids (PL) FA composition, when compared with their non-affected siblings. Venous blood samples were analysed for plasma and erythrocyte membrane FA composition by TLC followed by GC techniques. Lipid molecular species were determined by electrospray ionisation/tandem MS (ESI-MS/MS). FA analysis measured by standard GC showed an increased concentration of MUFA and a decreased concentration of plasmalogens in major PL fractions, but there were no differences in the concentrations of arachidonic acid or DHA. However, as identified by ESI-MS/MS, children with DS had increased levels of the following erythrocyte PL molecular species: 16 : 0–16 : 0, 16 : 0–18 : 1 and 16 : 0–18 : 2n-6, with reduced levels of 16 : 0–20 : 4n-6 species. Children with DS presented significantly higher levels of MUFA in both plasma and erythrocyte membrane, as well as higher levels of saturated and monounsaturated molecular species. Of interest was the almost double proportion of 16 : 0–18 : 2n-6 and nearly half the proportion of 16 : 0–20 : 4n-6 of choline phosphoacylglycerol species in children with DS compared with their non-affected siblings. These significant differences were only revealed by ESI-MS/MS and were not observed in the GC analysis. Further investigations are needed to explore molecular mechanisms and to test the association between the pathophysiology of DS and the risk of AD.

The lipid- and lipoprotein- [LDL-Lp(a)] apheresis techniques. Updating

Therapeutic plasmapheresis allows the extracorporeal removal of plasmatic lipoproteins (Lipid-apheresis) (LA). It can be non selective (non specific), semi - selective or selective low density lipoprotein-lipoprotein(a) (specific [LDL- Lp(a)] apheresis) (Lipoprotein apheresis, LDLa). The LDL removal rate is a perfect parameter to assess the system efficiency. Plasma-Exchange (PEX) cannot be considered either specific nor, selective. In PEX the whole blood is separated into plasma and its corpuscular components usually through centrifugation or rather filtration. The corpuscular components mixed with albumin solution plus saline (NaCl 0.9%) solution at 20%-25%, are then reinfused to the patient, to substitute the plasma formerly removed. PEX eliminates atherogenic lipoproteins, but also other essential plasma proteins, such as albumin, immunoglobulins, and hemocoagulatory mediators. Cascade filtration (CF) is a method based on plasma separation and removal of plasma proteins through double filtration. During the CF two hollow–fiber filters with pores of different diameter are used to eliminate the plasma components of different weight and molecular diameter. A CF system uses a first polypropylene filter with 0.55 µm diameter pores and a second one of diacetate of cellulose with 0.02 µm pores. The first filter separates the whole blood, and the plasma is then perfused through a second filter which allows the recovery of molecules with a diameter lower than 0.02 µm, and the removal of molecules larger in diameter as apoB100–containing lipoproteins. Since both albumin and immunoglobulins are not removed, or to a negligible extent, plasma-expanders, substitution fluids, and in particular albumin, as occurs in PEX are not needed. CF however, is characterized by lower selectivity since removes also high density lipoprotein (HDL) particles which have an antiatherogenic activity. In the 80’s, a variation of Lipid-apheresis has been developed which allows the LDL-cholesterol (LDLC) (-61%) and Lp(a) (-60%) removal from plasma through processing 3 liters of filtered plasma by means of lipid-specific thermofiltration, LDL immunoadsorption, heparin-induced LDL precipitation, LDL adsorption through dextran sulphate. More recently (90’s) the DALI®, and the Liposorber D® hemoperfusion systems, effective for apoB100- containing lipoproteins removal have been developed. All the above mentioned systems are established LDL-apheresis techniques referable to the generic definition of LDLa. However, this last definition cannot describe in an appropriate manner the removal of another highly atherogenic lipoprotein particle: the Lp(a). Thus it would be better to refer the above mentioned techniques to the wider scientific and technical concept of lipoprotein apheresis. Lipid apheresis - Lipoprotein apheresis - LDL-apheresis - Severe Dyslipidemia.

Analysis of Molecular Dynamics Simulations of Protein Folding

Microsecond long Molecular Dynamics (MD) trajectories of biomolecular processes are now possible due to advances in computer technology. Soon, trajectories long enough to probe dynamics over many milliseconds will become available. Since these timescales match the physiological timescales over which many small proteins fold, all atom MD simulations of protein folding are now becoming popular. To distill features of such large folding trajectories, we must develop methods that can both compress trajectory data to enable visualization, and that can yield themselves to further analysis, such as the finding of collective coordinates and reduction of the dynamics. Conventionally, clustering has been the most popular MD trajectory analysis technique, followed by principal component analysis (PCA). Simple clustering used in MD trajectory analysis suffers from various serious drawbacks, namely, (i) it is not data driven, (ii) it is unstable to noise and change in cutoff parameters, and (iii) since it does not take into account interrelationships amongst data points, the separation of data into clusters can often be artificial. Usually, partitions generated by clustering techniques are validated visually, but such validation is not possible for MD trajectories of protein folding, as the underlying structural transitions are not well understood. Rigorous cluster validation techniques may be adapted, but it is more crucial to reduce the dimensions in which MD trajectories reside, while still preserving their salient features. PCA has often been used for dimension reduction and while it is computationally inexpensive, being a linear method, it does not achieve good data compression. In this thesis, I propose a different method, a nonmetric multidimensional scaling (nMDS) technique, which achieves superior data compression by virtue of being nonlinear, and also provides a clear insight into the structural processes underlying MD trajectories. I illustrate the capabilities of nMDS by analyzing three complete villin headpiece folding and six norleucine mutant (NLE) folding trajectories simulated by Freddolino and Schulten [1]. Using these trajectories, I make comparisons between nMDS, PCA and clustering to demonstrate the superiority of nMDS. The three villin headpiece trajectories showed great structural heterogeneity. Apart from a few trivial features like early formation of secondary structure, no commonalities between trajectories were found. There were no units of residues or atoms found moving in concert across the trajectories. A flipping transition, corresponding to the flipping of helix 1 relative to the plane formed by helices 2 and 3 was observed towards the end of the folding process in all trajectories, when nearly all native contacts had been formed. However, the transition occurred through a different series of steps in all trajectories, indicating that it may not be a common transition in villin folding. The trajectories showed competition between local structure formation/hydrophobic collapse and global structure formation in all trajectories. Our analysis on the NLE trajectories confirms the notion that a tight hydrophobic core inhibits correct 3-D rearrangement. Only one of the six NLE trajectories folded, and it showed no flipping transition. All the other trajectories get trapped in hydrophobically collapsed states. The NLE residues were found to be buried deeply into the core, compared to the corresponding lysines in the villin headpiece, thereby making the core tighter and harder to undo for 3-D rearrangement. Our results suggest that the NLE may not be a fast folder as experiments suggest. The tightness of the hydrophobic core may be a very important factor in the folding of larger proteins. It is likely that chaperones like GroEL act to undo the tight hydrophobic core of proteins, after most secondary structure elements have been formed, so that global rearrangement is easier. I conclude by presenting facts about chaperone-protein complexes and propose further directions for the study of protein folding.

A new age of fuel performance code criteria studied through advanced atomistic simulation techniques

A fundamental step in understanding the effects of irradiation on metallic uranium and uranium dioxide ceramic fuels, or any material, must start with the nature of radiation damage on the atomic level. The atomic damage displacement results in a multitude of defects that influence the fuel performance. Nuclear reactions are coupled, in that changing one variable will alter others through feedback. In the field of fuel performance modeling, these difficulties are addressed through the use of empirical models rather than models based on first principles. Empirical models can be used as a predictive code through the careful manipulation of input variables for the limited circumstances that are closely tied to the data used to create the model. While empirical models are efficient and give acceptable results, these results are only applicable within the range of the existing data. This narrow window prevents modeling changes in operating conditions that would invalidate the model as the new operating conditions would not be within the calibration data set. This work is part of a larger effort to correct for this modeling deficiency. Uranium dioxide and metallic uranium fuels are analyzed through a kinetic Monte Carlo code (kMC) as part of an overall effort to generate a stochastic and predictive fuel code. The kMC investigations include sensitivity analysis of point defect concentrations, thermal gradients implemented through a temperature variation mesh-grid, and migration energy values. In this work, fission damage is primarily represented through defects on the oxygen anion sublattice. Results were also compared between the various models. Past studies of kMC point defect migration have not adequately addressed non-standard migration events such as clustering and dissociation of vacancies. As such, the General Utility Lattice Program (GULP) code was utilized to generate new migration energies so that additional non-migration events could be included into kMC code in the future for more comprehensive studies. Defect energies were calculated to generate barrier heights for single vacancy migration, clustering and dissociation of two vacancies, and vacancy migration while under the influence of both an additional oxygen and uranium vacancy.

Characterization of nanoparticle mediated laser transfection by femtosecond laser pulses for applications in molecular medicine

Background: In molecular medicine, the manipulation of cells is prerequisite to evaluate genes as therapeutic targets or to transfect cells to develop cell therapeutic strategies. To achieve these purposes it is essential that given transfection techniques are capable of handling high cell numbers in reasonable time spans. To fulfill this demand, an alternative nanoparticle mediated laser transfection method is presented herein. The fs-laser excitation of cell-adhered gold nanoparticles evokes localized membrane permeabilization and enables an inflow of extracellular molecules into cells. Results: The parameters for an efficient and gentle cell manipulation are evaluated in detail. Efficiencies of 90% with a cell viability of 93% were achieved for siRNA transfection. The proof for a molecular medical approach is demonstrated by highly efficient knock down of the oncogene HMGA2 in a rapidly proliferating prostate carcinoma in vitro model using siRNA. Additionally, investigations concerning the initial perforation mechanism are conducted. Next to theoretical simulations, the laser induced effects are experimentally investigated by spectrometric and microscopic analysis. The results indicate that near field effects are the initial mechanism of membrane permeabilization. Conclusion: This methodical approach combined with an automated setup, allows a high throughput targeting of several 100,000 cells within seconds, providing an excellent tool for in vitro applications in molecular medicine. NIR fs lasers are characterized by specific advantages when compared to lasers employing longer (ps/ns) pulses in the visible regime. The NIR fs pulses generate low thermal impact while allowing high penetration depths into tissue. Therefore fs lasers could be used for prospective in vivo applications.

A influência da incorporação de ácido α-eleosteárico, purificado e isolado do óleo de tungue, na dinâmica molecular de lipossomos

Este trabalho descreve o isolamento e purificação do ácido α-eleosteárico (α-ESA) a partir do óleo de tungue e sua caracterização por espectroscopia de infravermelho com transformada de Fourier (FTIR), cromatografia gasosa acoplada com espectrometria de massas (GC-MS) e espectroscopia de ressonância magnética nuclear (RMN) de 1H e 13C. O α-ESA apresenta atividades biológicas (antitumorais, anti-inflamatórias e antioxidantes), tornando-se importante compreender sua interação com membranas lipídicas. Assim, este trabalho também descreve resultados referentes ao efeito da incorporação de α-ESA na dinâmica molecular de lipossomos compostos por fosfatidilcolina. O sistema lipossomal puro e contendo α-ESA foi caracterizado através do uso de espectroscopia de UV-visível, FTIR, RMN e calorimetria de varredura diferencial (DSC). Como resultados da purificação do α-ESA, obtivemos uma pureza de 95,9% utilizando acetona como solvente de recristalização em detrimento dos 92,2% em solução etanólica. Na incorporação em lipossomos, observou-se uma maior interação do α-ESA com a parte polar, de interface e os primeiros metilenos da região apolar da fosfatidilcolina. Além disso, α-ESA apresentou um efeito de redução da fluidez de lipossomos. Os resultados contribuem para a geração de conhecimento para o desenvolvimento de novos sistemas farmacológicos.

Bayesian alignment of continuous molecular shapes using random fields

Statistical methodology is proposed for comparing molecular shapes. In order to account for the continuous nature of molecules, classical shape analysis methods are combined with techniques used for predicting random fields in spatial statistics. Applying a modification of Procrustes analysis, Bayesian inference is carried out using Markov chain Monte Carlo methods for the pairwise alignment of the resulting molecular fields. Superimposing entire fields rather than the configuration matrices of nuclear positions thereby solves the problem that there is usually no clear one--to--one correspondence between the atoms of the two molecules under consideration. Using a similar concept, we also propose an adaptation of the generalised Procrustes analysis algorithm for the simultaneous alignment of multiple molecular fields. The methodology is applied to a dataset of 31 steroid molecules.

Evolution of genome sequencing techniques

The quality and the speed for genome sequencing has advanced at the same time that technology boundaries are stretched. This advancement has been divided so far in three generations. The first-generation methods enabled sequencing of clonal DNA populations. The second-generation massively increased throughput by parallelizing many reactions while the third-generation methods allow direct sequencing of single DNA molecules. The first techniques to sequence DNA were not developed until the mid-1970s, when two distinct sequencing methods were developed almost simultaneously, one by Alan Maxam and Walter Gilbert, and the other one by Frederick Sanger. The first one is a chemical method to cleave DNA at specific points and the second one uses ddNTPs, which synthesizes a copy from the DNA chain template. Nevertheless, both methods generate fragments of varying lengths that are further electrophoresed. Moreover, it is important to say that until the 1990s, the sequencing of DNA was relatively expensive and it was seen as a long process. Besides, using radiolabeled nucleotides also compounded the problem through safety concerns and prevented the automation. Some advancements within the first generation include the replacement of radioactive labels by fluorescent labeled ddNTPs and cycle sequencing with thermostable DNA polymerase, which allows automation and signal amplification, making the process cheaper, safer and faster. Another method is Pyrosequencing, which is based on the “sequencing by synthesis” principle. It differs from Sanger sequencing, in that it relies on the detection of pyrophosphate release on nucleotide incorporation. By the end of the last millennia, parallelization of this method started the Next Generation Sequencing (NGS) with 454 as the first of many methods that can process multiple samples, calling it the 2º generation sequencing. Here electrophoresis was completely eliminated. One of the methods that is sometimes used is SOLiD, based on sequencing by ligation of fluorescently dye-labeled di-base probes which competes to ligate to the sequencing primer. Specificity of the di-base probe is achieved by interrogating every 1st and 2nd base in each ligation reaction. The widely used Solexa/Illumina method uses modified dNTPs containing so called “reversible terminators” which blocks further polymerization. The terminator also contains a fluorescent label, which can be detected by a camera. Now, the previous step towards the third generation was in charge of Ion Torrent, who developed a technique that is based in a method of “sequencing-by-synthesis”. Its main feature is the detection of hydrogen ions that are released during base incorporation. Likewise, the third generation takes into account nanotechnology advancements for the processing of unique DNA molecules to a real time synthesis sequencing system like PacBio; and finally, the NANOPORE, projected since 1995, also uses Nano-sensors forming channels obtained from bacteria that conducts the sample to a sensor that allows the detection of each nucleotide residue in the DNA strand. The advancements in terms of technology that we have nowadays have been so quick, that it makes wonder: ¿How do we imagine the next generation?

Aplicación de técnicas de cocina molecular en la elaboración de propuestas de repostería utilizando mora, uvilla y babaco

Lasiguiente monografíarecapitulalosconceptosy aplicaciones dela cocinamolecularintroducidoseneláreadelarepostería cuandoes utilizadamediantela ejecución dediversastécnicasy lautilizaciónde aditivosposibilitando cambiarlatexturadeunalimentosinalterarsu sabor, yasílograr la presentacióndeun postre sin ocasionar cambios en sus características organolépticas. Estamonografíatiene comoobjetivolaaplicacióndetécnicasdecocina molecularenlaelaboracióndequincepropuestasdereposteríautilizando frutascomolamora,uvillaybabaco,al mismotiempoconocerlas propiedadesdeestasfrutasysuutilizacióndentrodelareposteríaparael desarrollodequincenuevaspropuestas,querespondan alas necesidadesdeinnovaciónen larepostería,laexperimentaciónque permitiráapreciarlautilizacióndelafrutay surespuestaalsometerlasa laejecución de las distintastécnicas decocina molecular propuestas,así comola elaboracióndeuntest de aceptación desabory agradotanto visual comogustativo dentrodeun grupodepersonas. Repostería, cocina molecular, mora, uvilla, babaco.

Designing Hydrogels that Transform their Shape in Response to Molecular Cues

There has been considerable interest in developing shape-changing soft materials for potential applications in drug delivery, microfluidics and biosensing. These shape- changing materials are inspired by the morphological changes exhibited by plants in nature, such as the Venus flytrap. One specific class of shape-change is that from a flat sheet to a folded structure (e.g., a tube). Such “self-folding” materials are usually composed of polymer hydrogels, and these typically fold in response to external stimuli such as pH and temperature. In order to develop these hydrogels for the previously described applications, it is necessary to expand the range of triggers. The focus of this dissertation is the advancement of shape-changing polymer hydrogels that are sensitive to uncommon cues such as specific biomolecules (enzymes), the substrates for such enzymes, or specific multivalent cations. First, we describe a hybrid gel that responds to the presence of low concentrations of a class of enzymes known as matrix metalloproteinases (MMPs). The hybrid gel was created by utilizing photolithographic techniques to combine two or more gels with distinct chemical composition into the same material. Certain portions of the hybrid gel are composed of a biopolymer derivative with crosslinkable groups. The hybrid gel is flat in water; however, in the presence of MMPs, the regions containing the biopolymer are degraded and the flat sheet folds to form a 3D structure. We demonstrate that hydrogels with different patterns can transform into different 3D structures such as tubes, helices and pancakes. Furthermore, this shape change can be made to occur at physiological concentrations of enzymes. Next, we report a gel with two layers that undergoes a shape change in the presence of glucose. The enzyme glucose oxidase (GOx) is immobilized in one of the layers. GOx catalyzes the conversion of glucose to gluconic acid. The production of gluconic acid decreases the local pH. The decrease in local pH causes one of the layers to swell. As a result, the flat sheet folds to form a tube. The tube unfolds to form a flat sheet when it is transferred to a solution with no glucose present. Therefore, this biomolecule- triggered shape transformation is reversible, meaning the glucose sensing gel is reusable. Furthermore, this shape change only occurs in the presence of glucose and it does not occur in the presence of other small sugars such as fructose. In our final study, we report the shape change of a gel with two layers in the presence of multivalent ions such as Ca2+ and Sr2+. The gel consists of a passive layer and an active layer. The passive layer is composed of dimethylyacrylamide (DMAA), which does not interact with multivalent ions. The active layer consists of DMAA and the biopolymer alginate. In the presence of Ca2+ ions, the alginate chains crosslink and the active layer shrinks. As a result, the gel converts from a flat sheet to a folded tube. What is particularly unusual is the direction of folding. In most cases, when flat rectangular gels fold, they do so about their short-side. However, our gels typically fold about their long-side. We hypothesize that non-homogeneous swelling determines the folding axis.

Genetic diversity assessment of the almond (Prunus dulcis (Mill.) DA Webb) traditional germplasm of Algarve, Portugal, using molecular markers

In this study, 123 almond (Prunus dulcis (Mill.) D. A. Webb) trees identified among traditional orchards in the Algarve region and 53 trees of the local field collection managed by the regional office of the Portuguese Ministry of Agriculture (DRAALG) were assessed using isozyme, inter- single sequence repeat and simple sequence repeat or microsatellite techniques for the evaluation of genetic diversity and genetic relatedness and identification of new accessions for the field collection. The isozyme analysis allowed the distribution of the 176 plants into 13 different classes of enzyme similarity, while the use of DNA markers increased the distribution of the analysed trees among 140 discriminating DNA patterns. Multiple cases of homonymy and synonymy were identified in the local germplasm. Some traditional varieties, such as Lourencinha, appeared to be relatively homogeneous, while other local denominations, e.g. Galamba, included diverse genotypes. Of the 13 commercial varieties analysed in this study, 11 assembled in one major cluster clearly differentiated from the majority of the local genotypes. These results reinforced the perception that the Algarve traditional germplasm constitutes an important repository of genetic diversity, eventually carrying alleles of high agricultural interest such as the recently identified Phomopsis resistance in the traditional variety Barrinho Grado.

Genetic diversity assessment of the almond (Prunus dulcis (Mill.) DA Webb) traditional germplasm of Algarve, Portugal, using molecular markers

In this study, 123 almond (Prunus dulcis (Mill.) D. A. Webb) trees identified among traditional orchards in the Algarve region and 53 trees of the local field collection managed by the regional office of the Portuguese Ministry of Agriculture (DRAALG) were assessed using isozyme, inter- single sequence repeat and simple sequence repeat or microsatellite techniques for the evaluation of genetic diversity and genetic relatedness and identification of new accessions for the field collection. The isozyme analysis allowed the distribution of the 176 plants into 13 different classes of enzyme similarity, while the use of DNA markers increased the distribution of the analysed trees among 140 discriminating DNA patterns. Multiple cases of homonymy and synonymy were identified in the local germplasm. Some traditional varieties, such as Lourencinha, appeared to be relatively homogeneous, while other local denominations, e.g. Galamba, included diverse genotypes. Of the 13 commercial varieties analysed in this study, 11 assembled in one major cluster clearly differentiated from the majority of the local genotypes. These results reinforced the perception that the Algarve traditional germplasm constitutes an important repository of genetic diversity, eventually carrying alleles of high agricultural interest such as the recently identified Phomopsis resistance in the traditional variety Barrinho Grado.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

MOLECULAR MODELING OF EPON 862-DETDA / CARBON COMPOSITES

The thermoset epoxy resin EPON 862, coupled with the DETDA hardening agent, are utilized as the polymer matrix component in many graphite (carbon fiber) composites. Because it is difficult to experimentally characterize the interfacial region, computational molecular modeling is a necessary tool for understanding the influence of the interfacial molecular structure on bulk-level material properties. The purpose of this research is to investigate the many possible variables that may influence the interfacial structure and the effect they will have on the mechanical behavior of the bulk level composite. Molecular models are established for EPON 862-DETDA polymer in the presence of a graphite surface. Material characteristics such as polymer mass-density, residual stresses, and molecular potential energy are investigated near the polymer/fiber interface. Because the exact degree of crosslinking in these thermoset systems is not known, many different crosslink densities (degrees of curing) are investigated. It is determined that a region exists near the carbon fiber surface in which the polymer mass density is different than that of the bulk mass density. These surface effects extend ~10 Å into the polymer from the center of the outermost graphite layer. Early simulations predict polymer residual stress levels to be higher near the graphite surface. It is also seen that the molecular potential energy in polymer atoms decreases with increasing crosslink density. New models are then established in order to investigate the interface between EPON 862-DETDA polymer and graphene nanoplatelets (GNPs) of various atomic thicknesses. Mechanical properties are extracted from the models using Molecular Dynamics techniques. These properties are then implemented into micromechanics software that utilizes the generalized method of cells to create representations of macro-scale composites. Micromechanics models are created representing GNP doped epoxy with varying number of graphene layers and interfacial polymer crosslink densities. The initial micromechanics results for the GNP doped epoxy are then taken to represent the matrix component and are re-run through the micromechanics software with the addition of a carbon fiber to simulate a GNP doped epoxy/carbon fiber composite. Micromechanics results agree well with experimental data, and indicate GNPs of 1 to 2 atomic layers to be highly favorable. The effect of oxygen bonded to the surface of the GNPs is lastly investigated. Molecular Models are created for systems with varying graphene atomic thickness, along with different amounts of oxygen species attached to them. Models are created for graphene containing hydroxyl groups only, epoxide groups only, and a combination of epoxide and hydroxyl groups. Results show models of oxidized graphene to decrease in both tensile and shear modulus. Attaching only epoxide groups gives the best results for mechanical properties, though pristine graphene is still favored.