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.

Cell-cell communication in three dimensional microenvironments

Cellular behavior is dependent on a variety of extracellular cues required for normal tissue function, wound healing, and activation of the immune system. Removed from their in vivo microenvironment and cultured in vitro, cells lose many environmental cues and that may result in abberant behavior, making it difficult to study cellular processes. In order to mimic native tissue environments, optical tweezer and microfluidic technologies were used to place cells within defined areas of the culture environment. To provide three dimensional supports found in natural tissues, hydrogel scaffolds of poly (ethylene glycol) diacrylate and the basement membrane matrix Matrigel were used. Optical tweezer technology allowed precision placement and formation of homotypic and heterotypic arrays of human U937, HEK 293, and porcine mesenchymal stem cells. Alternatively, two microfluidic devices were designed to pattern Matrigel scaffolds. The first microfluidic device utilized laminar flow to spatially pattern multiple cell types within the device. Gradients of soluble molecules were then be formed and manipulated across the Matrigel scaffolds. Patterning Matrigel using laminar flow techniques require microfluidic expertise and do not produce consistent patterning conditions, limiting their use difficult in most cell culture laboratories. Thus, a buried Matrigel polydimethylsiloxane (PDMS) device was developed for spatial patterning of biological scaffolds. Matrigel is injected into micron sized channels of PDMS fabricated by soft lithography and allowed to thermally cure. Following curing, a second PDMS device was placed on top of the buried Matrigel channels to support media flow. In order to validate these systems, a cell-cell communication model system was developed utilizing LPS and TNFα signaling with fluorescent reporter systems to monitor communication in real time. We demonstrated the utility of microfluidic devices to support the cell-cell communication model system by co culturing three cell types within Matrigel scaffolds and monitoring signaling activity via fluorescent reporters.

Core-level photoemission investigations of the CDTE(100) and INSB(100) surface and interfacial structures

Photoemission techniques, utilizing a synchrotron light source, were used to analyze the clean (100) surfaces of the zinc-blende semiconductor materials CdTe and InSb. Several interfacial systems involving the surfaces of these materials were also studied, including the CdTe(lOO)-Ag interface, the CdTe(lOO)-Sb system, and the InSb(lOO)-Sn interface. High-energy electron diffraction was also employed to acquire information about of surface structure. A one-domain (2xl) structure was observed for the CdTe(lOO) surface. Analysis of photoemission spectra of the Cd 4d core level for this surface structure revealed two components resulting from Cd surface atoms. The total intensity of these components accounts for a full monolayer of Cd atoms on the surface. A structural model is discussed commensurate with these results. Photoemission spectra of the Cd and Te 4d core levels indicate that Ag or Sb deposited on the CdTe(l00)-(2xl) surface at room temperature do not bound strongly to the surface Cd atoms. The room temperature growth characteristics for these two elements on the CdTe(lOO)-(2xl) are discussed. The growth at elevated substrate temperatures was also studied for Sb deposition. The InSb(lOO) surface differed from the CdTe(lOO) surface. Using molecular beam epitaxy, several structures could be generated for the InSb(lOO) surface, including a c(8x2), a c(4x4), an asymmetric (lx3), a symmetric (lx3), and a (lxl). Analysis of photoemission intensities and line shapes indicates that the c(4x4) surface is terminated with 1-3/4 monolayers of Sb atoms. The c(8x2) surface is found to be terminated with 3/4 monolayer of In atoms. Structural models for both of these surfaces are proposed based upon the photoemission results and upon models of the similar GaAs(lOO) structures. The room temperature growth characteristics of grey Sn on the lnSb(lOO)-c(4x4) and InSb(l00)-c(8x2) surfaces were studied with photoemission. The discontinuity in the valence band maximum for this semiconductor heterojunction system is measured to be 0.40 eV, independent of the starting surface structure and stoichiometry. This result is reconciled with theoretical predictions for heterostructure behavior.