Switch Region for Pathogenic Structural Change in Conformational Disease and Its Prediction

Many diseases are believed to be related to abnormal protein folding. In the first step of such pathogenic structural changes, misfolding occurs in regions important for the stability of the native structure. This destabilizes the normal protein conformation, while exposing the previously hidden aggregation-prone regions, leading to subsequent errors in the folding pathway. Sites involved in this first stage can be deemed switch regions of the protein, and can represent perfect binding targets for drugs to block the abnormal folding pathway and prevent pathogenic conformational changes. In this study, a prediction algorithm for the switch regions responsible for the start of pathogenic structural changes is introduced. With an accuracy of 94%, this algorithm can successfully find short segments covering sites significant in triggering conformational diseases (CDs) and is the first that can predict switch regions for various CDs. To illustrate its effectiveness in dealing with urgent public health problems, the reason of the increased pathogenicity of H5N1 influenza virus is analyzed; the mechanisms of the pandemic swine-origin 2009 A(H1N1) influenza virus in overcoming species barriers and in infecting large number of potential patients are also suggested. It is shown that the algorithm is a potential tool useful in the study of the pathology of CDs because: (1) it can identify the origin of pathogenic structural conversion with high sensitivity and specificity, and (2) it provides an ideal target for clinical treatment.

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.

Coastal community hazard mitigation and community rating system of NFIP

Storm force flooding continues to be a major concern in the hurricane season and causes considerable loss to the coastal communities. National Flood Insurance Program (NFIP) provides recovery resources for the flood disaster and dissuades uneconomic uses from locating in flood hazard area. In order to motivate flood insurance purchase and promote increased flood hazard mitigation, the Community Rating System (CRS) that is a part of NFIP, credits 18 community floodplain management activities. However, CRS has been marked by a lack of active participation since its inception limiting its potential effectiveness. As of January 2008, 1080 communities, representing only 5% of all the NFIP communities have enrolled in CRS. Little empirical evidence exists to shed light on what factors influence the establishment of local hazard mitigation projects. To fill this gap, we propose to analyze flood hazard mitigation projects in 37 North Carolina coastal counties between 2002 and 2008. Specifically, we will examine the influence of physical, risk, and socioeconomic factors on coastal community hazard mitigation decisions as reflected in the CRS score. Ultimately, our project will forge a better understanding of community decision making, as related to natural hazards. (PDF contains 4 pages)

Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.