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.

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

Training civic bridge builders: Outcomes of community leadership development programs

As rural communities experience rapid economic, demographic, and political change, program interventions that focus on the development of community leadership capacity could be valuable. Community leadership development programs have been deployed in rural U.S. communities for the past 30 years by university extension units, chambers of commerce, and other nonprofit foundations. Prior research on program outcomes has largely focused on trainees’ self-reported change in individual leadership knowledge, skills, and attitudes. However, postindustrial leadership theories suggest that leadership in the community relies not on individuals but on social relationships that develop across groups akin to social bridging. The purpose of this study is to extend and strengthen prior evaluative research on community leadership development programs by examining program effects on opportunities to develop bridging social capital using more rigorous methods. Data from a quasi-experimental study of rural community leaders (n = 768) in six states are used to isolate unique program effects on individual changes in both cognitive and behavioral community leadership outcomes. Regression modeling shows that participation in community leadership development programs is associated with increased leadership development in knowledge, skills, attitudes, and behaviors that are a catalyst for social bridging. The community capitals framework is used to show that program participants are significantly more likely to broaden their span of involvement across community capital asset areas over time compared to non-participants. Data on specific program structure elements show that skills training may be important for cognitive outcomes while community development learning and group projects are important for changes in organizational behavior. Suggestions for community leadership program practitioners are presented.

Air-coupled ultrasonic tomographic imaging of concrete elements

Ultrasonic tomography is a powerful tool for identifying defects within an object or structure. This method can be applied on structures where x-ray tomography is impractical due to size, low contrast, or safety concerns. By taking many ultrasonic pulse velocity (UPV) readings through the object, an image of the internal velocity variations can be constructed. Air-coupled UPV can allow for more automated and rapid collection of data for tomography of concrete. This research aims to integrate recent developments in air-coupled ultrasonic measurements with advanced tomography technology and apply them to concrete structures. First, non-contact and semi-contact sensor systems are developed for making rapid and accurate UPV measurements through PVC and concrete test samples. A customized tomographic reconstruction program is developed to provide full control over the imaging process including full and reduced spectrum tomographs with percent error and ray density calculations. Finite element models are also used to determine optimal measurement configurations and analysis procedures for efficient data collection and processing. Non-contact UPV is then implemented to image various inclusions within 6 inch (152 mm) PVC and concrete cylinders. Although there is some difficulty in identifying high velocity inclusions, reconstruction error values were in the range of 1.1-1.7% for PVC and 3.6% for concrete. Based upon the success of those tests, further data are collected using non-contact, semi-contact, and full contact measurements to image 12 inch (305 mm) square concrete cross-sections with 1 inch (25 mm) reinforcing bars and 2 inch (51 mm) square embedded damage regions. Due to higher noise levels in collected signals, tomographs of these larger specimens show reconstruction error values in the range of 10-18%. Finally, issues related to the application of these techniques to full-scale concrete structures are discussed.