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.

The interplay between neuronal activity and actin dynamics mimic the setting of an LTD synaptic tag

Persistent forms of plasticity, such as long-term depression (LTD), are dependent on the interplay between activity-dependent synaptic tags and the capture of plasticity-related proteins. We propose that the synaptic tag represents a structural alteration that turns synapses permissive to change. We found that modulation of actin dynamics has different roles in the induction and maintenance of LTD. Inhibition of either actin depolymerisation or polymerization blocks LTD induction whereas only the inhibition of actin depolymerisation blocks LTD maintenance. Interestingly, we found that actin depolymerisation and CaMKII activation are involved in LTD synaptic-tagging and capture. Moreover, inhibition of actin polymerisation mimics the setting of a synaptic tag, in an activity-dependent manner, allowing the expression of LTD in non-stimulated synapses. Suspending synaptic activation also restricts the time window of synaptic capture, which can be restored by inhibiting actin polymerization. Our results support our hypothesis that modulation of the actin cytoskeleton provides an input-specific signal for synaptic protein capture.

STRUCTURAL ANALYSIS OF RIBOSOME BINDING ELEMENTS OF THE 3´ UTR IN TWO POSITIVE-STRAND PLANT RNA VIRUSES

Turnip crinkle virus (TCV) and Pea enation mosaic virus (PEMV) are two positive (+)-strand RNA viruses that are used to investigate the regulation of translation and replication due to their small size and simple genomes. Both viruses contain cap-independent translation elements (CITEs) within their 3´ untranslated regions (UTRs) that fold into tRNA-shaped structures (TSS) according to nuclear magnetic resonance and small angle x-ray scattering analysis (TCV) and computational prediction (PEMV). Specifically, the TCV TSS can directly associate with ribosomes and participates in RNA-dependent RNA polymerase (RdRp) binding. The PEMV kissing-loop TSS (kl-TSS) can simultaneously bind to ribosomes and associate with the 5´ UTR of the viral genome. Mutational analysis and chemical structure probing methods provide great insight into the function and secondary structure of the two 3´ CITEs. However, lack of 3-D structural information has limited our understanding of their functional dynamics. Here, I report the folding dynamics for the TCV TSS using optical tweezers (OT), a single molecule technique. My study of the unfolding/folding pathways for the TCV TSS has provided an unexpected unfolding pathway, confirmed the presence of Ψ3 and hairpin elements, and suggested an interconnection between the hairpins and pseudoknots. In addition, this study has demonstrated the importance of the adjacent upstream adenylate-rich sequence for the formation of H4a/Ψ3 along with the contribution of magnesium to the stability of the TCV TSS. In my second project, I report on the structural analysis of the PEMV kl-TSS using NMR and SAXS. This study has re-confirmed the base-pair pattern for the PEMV kl-TSS and the proposed interaction of the PEMV kl-TSS with its interacting partner, hairpin 5H2. The molecular envelope of the kl-TSS built from SAXS analysis suggests the kl-TSS has two functional conformations, one of which has a different shape from the previously predicted tRNA-shaped form. Along with applying biophysical methods to study the structural folding dynamics of RNAs, I have also developed a technique that improves the production of large quantities of recombinant RNAs in vivo for NMR study. In this project, I report using the wild-type and mutant E.coli strains to produce cost-effective, site-specific labeled, recombinant RNAs. This technique was validated with four representative RNAs of different sizes and complexity to produce milligram amounts of RNAs. The benefit of using site-specific labeled RNAs made from E.coli was demonstrated with several NMR techniques.

Data of the molecular dynamics simulations of mutations in the human connexin46 docking interface

The structure of hCx26 derived from the X-ray analysis was used to generate a homology model for hCx46. Interacting connexin molecules were used as starting model for the molecular dynamics (MD) simulation using NAMD and allowed us to predict the dynamic behavior of hCx46wt and the cataract related mutant hCx46N188T as well as two artificial mutants hCx46N188Q and hCx46N188D. Within the 50 ns simulation time the docked complex composed of the mutants dissociate while hCx46wt remains stable. The data indicates that one hCx46 molecule forms 5-7 hydrogen bonds (HBs) with the counterpart connexin of the opposing connexon. These HBs appear essential for a stable docking of the connexons as shown by the simulation of an entire gap junction channel and were lost for all the tested mutants. The data described here are related to the research article entitled "The cataract related mutation N188T in human connexin46 (hCx46) revealed a critical role for residue N188 in the docking process of gap junction channels" (Schadzek et al., 2015) [1].

STRUCTURAL AND FUNCTIONAL STUDIES OF CYCLIC K48-LINKED DIUBIQUITIN

K48-linked di-ubiquitin exists in a dynamic equilibrium between open and closed states. The structure of K48-Ub2 in the closed conformation features a hydrophobic interface formed between the two Ub domains. The same hydrophobic residues at the interface are involved in binding to ubiquitin-associated (UBA) domains. Cyclization of K48-Ub2 should limit the range of conformations available for such interactions. Interestingly, cyclic K48-linked Ub2 (cycUb2) has been found in vivo and can be isolated in vitro to study its structure and dynamics. In this study, a crystal structure of cycUb2 was obtained, and the dynamics of cycUb2 were characterized by solution NMR. The crystal structure of cycUb2, which is in agreement with solution NMR data, is closed with the hydrophobic patches of each Ub domain buried at the interface. Despite its structural constraints, cycUb2 was still able to interact with UBA domains, albeit with lower affinity.

Structural and functional interaction between polyamine related molecules and biological membranes

Changes induced by PA on nucleic acid (NA) conformation and synthesis is proven to be a major reason for PA essentiality (1-3). However, PA interactions with other polyanions, for instance polyanionic membrane lipid bilayers and glyosaminoglycans have received less attention (3-4). The functional importance of these interactions still is an obscure but interesting area of cell and molecular biology, especially in mammalian cells for which specific PA transport systems are not fully characterized (5). In mammals, activity and turnover of the polyamine (PA) synthesis key enzyme is controlled by a set of proteins: Antizymes (OAZ1-3) and antizyme inhibitors (AZIN1 and 2). It is demonstrated that AOZ modulate polyamine uptake (6), and that PA transport to mitochondria is linked to the respiratory chain state and modulates mitochondrial permeability transition (7). Antizyme expression variants have been located in mitochondria, being proposed as a proapoptotic factor (7-8). AZIN 2 is only expressed in a reduced set of tissues that includes mast cells, where it is associated to mast cell granules membrane (9). This fact, together to the abnormalities observed in bone marrow derived mast cell granules when they are differentiated under restricted PA synthesis conditions (10 and unpublished results), point out to important roles of PA and their related proteins in structure and function of mast cell granules. We will also present novel biophysical results on tripartite interactions of PA that remark the interest of the characterization of PA interactions with lipid bilayers for biomedicine and biotechnology. Thus, the information reported in this paper integrates previously reported information with our still unpublished results, all indicating that PA and their related proteins also are important factors for structure and dynamics of biological membranes and their associated functions essential in human physiology; for instance, solute interchange with the environment (uptake and secretion), oxidative metabolism and apoptosis. The importance of these involved processes for human homeostasis claim for further research efforts. 1. Ruiz-Chica J, Medina MA, Sánchez-Jiménez F and Ramírez FJ (2001) Fourier Transform Raman study of the structural specificities on the interaction between DNA and biogenic polyamines. Biophysical J. 80:443-454. 2. Lightfoot HL, Hall J (2014) Endogenous polyamine function--the RNA perspective. Nucleic Acids Res. 42:11275-11290. 3. Igarashi K, Kashiwagi K (2010) Modulation of cellular function by polyamines. Int J Biochem Cell Biol. 42:39-51. 4. Finger S, Schwieger C, Arouri A, Kerth A, Blume A (2014) Interaction of linear polyamines with negatively charged phospholipids: the effect of polyamine charge distance. Biol Chem. 395:769-778. 5. Poulin R, Casero RA, Soulet D. (2012) Recent advances in the molecular biology of metazoan polyamine transport. Amino Acids. 42:711-723. 6. Kahana C (2009) Regulation of cellular polyamine levels and cellular proliferation by antizyme and antizyme inhibitor. Essays Biochem. 4:47-61. 7. Agostinelli E, Marques MP, Calheiros R, Gil FP, Tempera G, Viceconte N, Battaglia V, Grancara S, Toninello A (2010) Polyamines: fundamental characters in chemistry and biology. Amino Acids 38:393-403. 8. Liu GY, Liao YF, Hsu PC, Chang WH, Hsieh MC, Lin CY, Hour TC, Kao MC, Tsay GJ, Hung HC (2006) Antizyme, a natural ornithine decarboxylase inhibitor, induces apoptosis of haematopoietic cells through mitochondrial membrane depolarization and caspases' cascade. Apoptosis 11:1773-1788. 9. Kanerva K, Lappalainen J, Mäkitie LT, Virolainen S, Kovanen PT, Andersson LC (2009). Expression of antizyme inhibitor 2 in mast cells and role of polyamines as selective regulators of serotonin secretion. PLoS One 31:e6858. 10. García-Faroldi G, Rodríguez CE, Urdiales JL, Pérez-Pomares JM, Dávila JC, Pejler G, Sánchez-Jiménez F, Fajardo I (2010) Polyamines are present in mast cell secretory granules and are important for granule homeostasis. PLoS One 30:e15071.

MEASURING HEAD IMPACT CONTACT PRESSURE IN COLLEGIATE FOOTBALL GAMES TO CORRELATE HEAD KINEMATICS TO BRAIN KINETICS ELUCIDATING BRAIN INJURY DYNAMICS

Does a brain store thoughts and memories the way a computer saves its files? How can a single hit or a fall erase all those memories? Brain Mapping and traumatic brain injuries (TBIs) have become widely researched fields today. Many researchers have been studying TBIs caused to adult American football players however youth athletes have been rarely considered for these studies, contradicting to the fact that American football enrolls highest number of collegiate and high-school children than adults. This research is an attempt to contribute to the field of youth TBIs. Earlier studies have related head kinematics (linear and angular accelerations) to TBIs. However, fewer studies have dealt with brain kinetics (impact pressures and stresses) occurring during head-on collisions. The National Operating Committee on Standards for Athletic Equipment (NOCSAE) drop tests were conducted for linear impact accelerations and the Head Impact Contact Pressures (HICP) calculated from them were applied to a validated FE model. The results showed lateral region of the head as the most vulnerable region to damage from any drop height or impact distance followed by posterior region. The TBI tolerance levels in terms of Von-Mises and Maximum Principal Stresses deduced for lateral impact were 30 MPa and 18 MPa respectively. These levels were corresponding to 2.625 feet drop height. The drop heights beyond this value will result in TBI causing stress concentrations in human head without any detectable structural damage to the brain tissue. This data can be utilized for designing helmets that provide cushioning to brain along with providing a resistance to shear.

Pressure Induced Structural Changes and Gas Diffusion Pathways in Monomeric Fluorescent Proteins

Fluorescent proteins (FPs) are extremely valuable biochemical markers which have found a wide range of applications in cellular and molecular biology research. The monomeric variants of red fluorescent proteins (RFPs), known as mFruits, have been especially valuable for in vivo applications in mammalian cell imaging. Fluorescent proteins consist of a chromophore caged in the beta-barrel protein scaffold. The photophysical properties of an FP is determined by its chromophore structure and its interactions with the protein barrel. Application of hydrostatic pressure on FPs results in the modification of the chromophore environment which allows a systematic study of the role of the protein-chromophore interactions on photophysical properties of FPs. Using Molecular Dynamics (MD) computer simulations, I investigated the pressure induced structural changes in the monomeric variants mCherry, mStrawberry, and Citrine. The results explain the molecular basis for experimentally observed pressure responses among FP variants. It is found that the barrel flexibility, hydrogen bonding interactions and chromophore planarity of the FPs can be correlated to their contrasting photophysical properties at vaious pressures. I also investigated the oxygen diffusion pathways in mOrange and mOrange2 which exhibit marked differences in oxygen sensitivities as well as photostability. Such computational identifications of structural changes and oxygen diffusion pathways are important in guiding mutagenesis efforts to design fluorescent proteins with improved photophysical properties.

Influence of kaolin addition on the dynamics of oxygen mass transport in polyvinyl alcohol dispersion coatings

The permeability of dispersion barriers produced from polyvinyl alcohol (PVOH) and kaolin clay blends coated onto polymeric supports has been studied by employing two different measurement methods: the oxygen transmission rate (OTR) and the ambient oxygen ingress rate (AOIR). Coatings with different thicknesses and kaolin contents were studied. Structural information of the dispersion-barrier coatings was obtained by Fourier transform infrared spectroscopy (FTIR) spectroscopy and scanning electron microscopy (SEM). These results showed that the kaolin content influences both the orientation of the kaolin and the degree of crystallinity of the PVOH coating. Increased kaolin content increased the alignment of the kaolin platelets to the basal plane of the coating. Higher kaolin content was accompanied by higher degree of crystallinity of the PVOH. The barrier thickness proved to be less important in the early stages of the mass transport process, whereas it had a significant influence on the steady-state permeability. The results from this study demonstrate the need for better understanding of how permeability is influenced by (chemical and physical) structure.

Molecular dynamics study of a dual-cation ionomer electrolyte

The poly(N1222)xLi1-x[AMPS] ionomer system with dual cations has previously shown decoupled Li ion dynamics from polymer segmental motions, characterized by the glass transition temperature, which can result in a conductive electrolyte material whilst retaining an appropriate modulus (ie. stiffness) so that it can suppress dendrite formation, thereby improving safety when used in lithium metal batteries. To understand this ion dynamics behavior, molecular dynamics techniques have been used in this work to simulate structure and dynamics in these materials. These simulations confirm that the Li ion transport is decoupled from the polymer particularly at intermediate N1222+ concentrations. At 50 mol% N1222+ concentration the polymer backbone is more rigid than for higher N1222+ concentrations, but with increasing temperature Li ion transport is more significant than polymer or quaternary ammonium cation motions. Here we suggest an ion hopping mechanism for Li+, arising from structural rearrangement of ionic clusters that could explain its decoupled behavior. Higher temperatures favor an aggregated ionic structure as well as enhancing these hopping motions. The simulations discussed here provide an atomic-level understanding of ion dynamics that could contribute to designing an improved ionomer with fast ion transport and mechanical robustness.

Exploring first-year academic achievement through structural equation modelling

The purpose of this research was to develop and test a multicausal model of the individual characteristics associated with academic success in first-year Australian university students. This model comprised the constructs of: previous academic performance, achievement motivation, self-regulatory learning strategies, and personality traits, with end-of-semester grades the dependent variable of interest. The study involved the distribution of a questionnaire, which assessed motivation, self-regulatory learning strategies and personality traits, to 1193 students at the start of their first year at university. Students' academic records were accessed at the end of their first year of study to ascertain their first and second semester grades. This study established that previous high academic performance, use of self-regulatory learning strategies, and being introverted and agreeable, were indicators of academic success in the first semester of university study. Achievement motivation and the personality trait of conscientiousness were indirectly related to first semester grades, through the influence they had on the students' use of self-regulatory learning strategies. First semester grades were predictive of second semester grades. This research provides valuable information for both educators and students about the factors intrinsic to the individual that are associated with successful performance in the first year at university.