Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.

Characterization of a Coal Fly Ash-Cement Slurry by the Absolute Foam Index

This dissertation established a standard foam index: the absolute foam index test. This test characterized a wide range of coal fly ash by the absolute volume of air-entraining admixture (AEA) necessary to produce a 15-second metastable foam in a coal fly ash-cement slurry in a specified time. The absolute foam index test was used to characterize fly ash samples having loss on ignition (LOI) values that ranged from 0.17 to 23.3 %wt. The absolute foam index characterized the fly ash samples by absolute volume of AEA, defined as the amount of undiluted AEA solution added to obtain a 15-minute endpoint signified by 15-second metastable foam. Results were compared from several foam index test time trials that used different initial test concentrations to reach termination at selected times. Based on the coefficient of variation (CV), a 15-minute endpoint, with limits of 12 to 18 minutes was chosen. Various initial test concentrations were used to accomplish consistent contact times and concentration gradients for the 15-minute test endpoint for the fly ash samples. A set of four standard concentrations for the absolute foam index test were defined by regression analyses and a procedure simplifying the test process. The set of standard concentrations for the absolute foam index test was determined by analyzing experimental results of 80 tests on coal fly ashes with loss on ignition (LOI) values ranging from 0.39 to 23.3 wt.%. A regression analysis informed selection of four concentrations (2, 6, 10, and 15 vol.% AEA) that are expected to accommodate fly ashes with 0.39 to 23.3 wt.% LOI, depending on the AEA type. Higher concentrations should be used for high-LOI fly ash when necessary. A procedure developed using these standard concentrations is expected to require only 1-3 trials to meet specified endpoint criteria for most fly ashes. The AEA solution concentration that achieved the metastable foam in the foam index test was compared to the AEA equilibrium concentration obtained from the direct adsorption isotherm test with the same fly ash. The results showed that the AEA concentration that satisfied the absolute foam index test was much less than the equilibrium concentration. This indicated that the absolute foam index test was not at or near equilibrium. Rather, it was a dynamic test where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms. Equilibrium isotherm equations obtained from direct isotherm tests were used to calculate the equilibrium concentrations and capacities of fly ash from 0.17 to 10.5% LOI. The results showed that the calculated fly ash capacity was much less than capacities obtained from isotherm tests that were conducted with higher initial concentrations. This indicated that the absolute foam index was not equilibrium. Rather, the test is dynamic where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms for fly ash of 0.17 to 10.5% LOI. Several batches of mortars were mixed for the same fly ash type increasing only the AEA concentration (dosage) in each subsequent batch. Mortar air test results for each batch showed for each increase in AEA concentration, air contents increased until a point where the next increase in AEA concentration resulted in no increase in air content. This was maximum air content that could be achieved by the particular mortar system; the system reached its air capacity at the saturation limit. This concentration of AEA was compared to the critical micelle concentration (CMC) for the AEA and the absolute foam index.

STUDY OF A DIRECT SAMPLING METHOD FOR THE INVERSE MEDIUM SCATTERING PROBLEM

Direct sampling methods are increasingly being used to solve the inverse medium scattering problem to estimate the shape of the scattering object. A simple direct method using one incident wave and multiple measurements was proposed by Ito, Jin and Zou. In this report, we performed some analytic and numerical studies of the direct sampling method. The method was found to be effective in general. However, there are a few exceptions exposed in the investigation. Analytic solutions in different situations were studied to verify the viability of the method while numerical tests were used to validate the effectiveness of the method.

SAMPLING BIAS IN EVALUATING THE PROBABILITY OF SEISMICALLY INDUCED SOIL LIQUEFACTION WITH SPT & CPT CASE HISTORIES

Several deterministic and probabilistic methods are used to evaluate the probability of seismically induced liquefaction of a soil. The probabilistic models usually possess some uncertainty in that model and uncertainties in the parameters used to develop that model. These model uncertainties vary from one statistical model to another. Most of the model uncertainties are epistemic, and can be addressed through appropriate knowledge of the statistical model. One such epistemic model uncertainty in evaluating liquefaction potential using a probabilistic model such as logistic regression is sampling bias. Sampling bias is the difference between the class distribution in the sample used for developing the statistical model and the true population distribution of liquefaction and non-liquefaction instances. Recent studies have shown that sampling bias can significantly affect the predicted probability using a statistical model. To address this epistemic uncertainty, a new approach was developed for evaluating the probability of seismically-induced soil liquefaction, in which a logistic regression model in combination with Hosmer-Lemeshow statistic was used. This approach was used to estimate the population (true) distribution of liquefaction to non-liquefaction instances of standard penetration test (SPT) and cone penetration test (CPT) based most updated case histories. Apart from this, other model uncertainties such as distribution of explanatory variables and significance of explanatory variables were also addressed using KS test and Wald statistic respectively. Moreover, based on estimated population distribution, logistic regression equations were proposed to calculate the probability of liquefaction for both SPT and CPT based case history. Additionally, the proposed probability curves were compared with existing probability curves based on SPT and CPT case histories.

Adaptive Sampling with Mobile Sensor Networks

Mobile sensor networks have unique advantages compared with wireless sensor networks. The mobility enables mobile sensors to flexibly reconfigure themselves to meet sensing requirements. In this dissertation, an adaptive sampling method for mobile sensor networks is presented. Based on the consideration of sensing resource constraints, computing abilities, and onboard energy limitations, the adaptive sampling method follows a down sampling scheme, which could reduce the total number of measurements, and lower sampling cost. Compressive sensing is a recently developed down sampling method, using a small number of randomly distributed measurements for signal reconstruction. However, original signals cannot be reconstructed using condensed measurements, as addressed by Shannon Sampling Theory. Measurements have to be processed under a sparse domain, and convex optimization methods should be applied to reconstruct original signals. Restricted isometry property would guarantee signals can be recovered with little information loss. While compressive sensing could effectively lower sampling cost, signal reconstruction is still a great research challenge. Compressive sensing always collects random measurements, whose information amount cannot be determined in prior. If each measurement is optimized as the most informative measurement, the reconstruction performance can perform much better. Based on the above consideration, this dissertation is focusing on an adaptive sampling approach, which could find the most informative measurements in unknown environments and reconstruct original signals. With mobile sensors, measurements are collect sequentially, giving the chance to uniquely optimize each of them. When mobile sensors are about to collect a new measurement from the surrounding environments, existing information is shared among networked sensors so that each sensor would have a global view of the entire environment. Shared information is analyzed under Haar Wavelet domain, under which most nature signals appear sparse, to infer a model of the environments. The most informative measurements can be determined by optimizing model parameters. As a result, all the measurements collected by the mobile sensor network are the most informative measurements given existing information, and a perfect reconstruction would be expected. To present the adaptive sampling method, a series of research issues will be addressed, including measurement evaluation and collection, mobile network establishment, data fusion, sensor motion, signal reconstruction, etc. Two dimensional scalar field will be reconstructed using the method proposed. Both single mobile sensors and mobile sensor networks will be deployed in the environment, and reconstruction performance of both will be compared.In addition, a particular mobile sensor, a quadrotor UAV is developed, so that the adaptive sampling method can be used in three dimensional scenarios.