Status of GPCR modeling and docking as reflected by community-wide GPCR Dock 2010 assessment

The community-wide GPCR Dock assessment is conducted to evaluate the status of molecular modeling and ligand docking for human G protein-coupled receptors. The present round of the assessment was based on the recent structures of dopamine D3 and CXCR4 chemokine receptors bound to small molecule antagonists and CXCR4 with a synthetic cyclopeptide. Thirty-five groups submitted their receptor-ligand complex structure predictions prior to the release of the crystallographic coordinates. With closely related homology modeling templates, as for dopamine D3 receptor, and with incorporation of biochemical and QSAR data, modern computational techniques predicted complex details with accuracy approaching experimental. In contrast, CXCR4 complexes that had less-characterized interactions and only distant homology to the known GPCR structures still remained very challenging. The assessment results provide guidance for modeling and crystallographic communities in method development and target selection for further expansion of the structural coverage of the GPCR universe.

Kinetostatic-model-based Stiffness Analysis of Exechon PKM

As a comparative newly-invented PKM with over-constraints in kinematic chains, the Exechon has attracted extensive attention from the research society. Different from the well-recognized kinematics analysis, the research on the stiffness characteristics of the Exechon still remains as a challenge due to the structural complexity. In order to achieve a thorough understanding of the stiffness characteristics of the Exechon PKM, this paper proposed an analytical kinetostatic model by using the substructure synthesis technique. The whole PKM system is decomposed into a moving platform subsystem, three limb subsystems and a fixed base subsystem, which are connected to each other sequentially through corresponding joints. Each limb body is modeled as a spatial beam with a uniform cross-section constrained by two sets of lumped springs. The equilibrium equation of each individual limb assemblage is derived through finite element formulation and combined with that of the moving platform derived with Newtonian method to construct the governing kinetostatic equations of the system after introducing the deformation compatibility conditions between the moving platform and the limbs. By extracting the 6 x 6 block matrix from the inversion of the governing compliance matrix, the stiffness of the moving platform is formulated. The computation for the stiffness of the Exechon PKM at a typical configuration as well as throughout the workspace is carried out in a quick manner with a piece-by-piece partition algorithm. The numerical simulations reveal a strong position-dependency of the PKM's stiffness in that it is symmetric relative to a work plane due to structural features. At the last stage, the effects of some design variables such as structural, dimensional and stiffness parameters on system rigidity are investigated with the purpose of providing useful information for the structural optimization and performance enhancement of the Exechon PKM. It is worthy mentioning that the proposed methodology of stiffness modeling in this paper can also be applied to other overconstrained PKMs and can evaluate the global rigidity over workplace efficiently with minor revisions.

A Lagrangian approach for modeling road collisions using second-order models of traffic flow

In [M. Herty, A. Klein, S. Moutari, V. Schleper, and G. Steinaur, IMA J. Appl. Math., 78(5), 1087–1108, 2013] and [M. Herty and V. Schleper, ZAMM J. Appl. Math. Mech., 91, 763–776, 2011], a macroscopic approach, derived from fluid-dynamics models, has been introduced to infer traffic conditions prone to road traffic collisions along highways’ sections. In these studies, the governing equations are coupled within an Eulerian framework, which assumes fixed interfaces between the models. A coupling in Lagrangian coordinates would enable us to get rid of this (not very realistic) assumption. In this paper, we investigate the well-posedness and the suitability of the coupling of the governing equations within the Lagrangian framework. Further, we illustrate some features of the proposed approach through some numerical simulations.

Modeling study of woody biomass:Interactions of cellulose, hemicellulose, and lignin

Lignocellulosic biomass pretreatment and the subsequent thermal conversion processes to produce solid, liquid, and gas biofuels are attractive solutions for today's energy challenges. The structural study of the main components in biomass and their macromolecular complexes is an active and ongoing research topic worldwide. The interactions among the three main components, cellulose, hemicellulose, and lignin, are studied in this paper using electronic structure methods, and the study includes examining the hydrogen bond network of cellulose-hemicellulose systems and the covalent bond linkages of hemicellulose-lignin systems. Several methods (semiempirical, Hartree-Fock, and density functional theory) using different basis sets were evaluated. It was shown that theoretical calculations can be used to simulate small model structures representing wood components. By comparing calculation results with experimental data, it was concluded that B3LYP/6-31G is the most suitable basis set to describe the hydrogen bond system and B3LYP/6-31G(d,p) is the most suitable basis set to describe the covalent system of woody biomass. The choice of unit model has a much larger effect on hydrogen bonding within cellulose-hemicellulose system, whereas the model choice has a minimal effect on the covalent linkage in the hemicellulose-lignin system. © 2011 American Chemical Society.

Dynamic Modeling and Interaction of Hybrid Natural Gas and Electricity Supply System in Microgrid

Natural gas (NG) network and electric network are becoming tightly integrated by microturbines in the microgrid. Interactions between these two networks are not well captured by the traditional microturbine (MT) models. To address this issue, two improved models for single-shaft MT and split-shaft MT are proposed in this paper. In addition, dynamic models of the hybrid natural gas and electricity system (HGES) are developed for the analysis of their interactions. Dynamic behaviors of natural gas in pipes are described by partial differential equations (PDEs), while the electric network is described by differential algebraic equations (DAEs). So the overall network is a typical two-time scale dynamic system. Numerical studies indicate that the two-time scale algorithm is faster and can capture the interactions between the two networks. The results also show the HGES with a single-shaft MT is a weakly coupled system in which disturbances in the two networks mainly influence the dc link voltage of the MT, while the split-shaft MT is a strongly coupled system where the impact of an event will affect both networks.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

Concrete damage plasticity model for modeling FRP-to-concrete bond behavior

The technique of externally bonding fiber-reinforced polymer (FRP) composites has become very popular worldwide for retrofitting existing reinforced concrete (RC) structures. Debonding of FRP from the concrete substrate is a typical failure mode in such strengthened structures. The bond behavior between FRP and concrete thus plays a crucial role in these structures. The FRP-to-concrete bond behavior has been extensively investigated experimentally, commonly using a single or double shear test of the FRP-to-concrete bonded joint. Comparatively, much less research has been concerned with numerical simulation, chiefly due to difficulties in the accurate modeling of the complex behavior of concrete. This paper presents a simple but robust finite-element (FE) model for simulating the bond behavior in the entire debonding process for the single shear test. A concrete damage plasticity model is proposed to capture the concrete-to-FRP bond behavior. Numerical results are in close agreement with test data, validating the model. In addition to accuracy, the model has two further advantages: it only requires the basic material parameters (i.e., no arbitrary user-defined parameter such as the shear retention factor is required) and it can be directly implemented in the FE software ABAQUS.

Solutions of fourth-order parabolic equation modeling thin film growth

In this paper we study the well-posedness for a fourth-order parabolic equation modeling epitaxial thin film growth. Using Kato's Method [1], [2] and [3] we establish existence, uniqueness and regularity of the solution to the model, in suitable spaces, namelyC⁰([0,T];L^p(Ω)) where  with 1<α<2, n∈N and n≥2. We also show the global existence solution to the nonlinear parabolic equations for small initial data. Our main tools are L_p–L_q-estimates, regularization property of the linear part of e^−tΔ2 and successive approximations. Furthermore, we illustrate the qualitative behavior of the approximate solution through some numerical simulations. The approximate solutions exhibit some favorable absorption properties of the model, which highlight the stabilizing effect of our specific formulation of the source term associated with the upward hopping of atoms. Consequently, the solutions describe well some experimentally observed phenomena, which characterize the growth of thin film such as grain coarsening, island formation and thickness growth.

Assessment Of various sensors for structural health monitoring for bridge weigh-in-motion (B-WIM).

In recent years, Structural Health Monitoring (SHM) systems have been developed to monitor bridge deterioration, assess real load levels and hence extend bridge life and safety. A road bridge is only safe if the stresses caused by the passing vehicles are less than the capacity of the bridge to resist them. Conventional SHM systems can be used to improve knowledge of the bridges capacity to resist stresses but generally give no information on the causes of any increase in stresses (based on measuring strain). The concept of in Bridge Weigh-in-Motion (B-WIM) is to establish axle loads, without interruption to traffic flow, by using strain sensors at a bridge soffit and subsequently converting the data to real time axle loads or stresses. Recent studies have shown it would be most beneficial to develop a portable system which can be easily attached to existing and new bridge structures for a specified monitoring period. The sensors could then be left in place while the data acquisition can be moved for various other sites. Therefore it is necessary to find accurate sensors capable of capturing peak strains under dynamic load and suitable methods for attaching these strain sensors to existing and new bridge structures. Additionally, it is important to ensure accurate strain transfer between concrete and steel, the adhesives layer and the strain sensor. This paper describes research investigating the suitably of using various sensors for the monitoring of concrete structures under dynamic vehicle load. Electrical resistance strain (ERS) gauges, vibrating wire (VW) gauges and fibre optic sensors (FOS) are commonly used for SHM. A comparative study will be carried out to select a suitable sensor for a bridge Weigh in Motion System. This study will look at fixing methods, durability, scanning rate and accuracy range. Finite element modeling is used to predict the strains which are then validated in laboratory trials.

Switched Diversity Techniques for Indoor Off-body Communications Channels: An Experimental Analysis and Modeling

This paper investigates the potential improvement in signal reliability for indoor off-body communications channels operating at 5.8 GHz using switched diversity techniques. In particular we investigate the performance of switch-and-stay combining (SSC), switch-and-examine combining (SEC) and switch-and-examine combining with post-examining selection (SECps) schemes which utilize multiple spatially separated antennas at the base station. During the measurements a test subject, wearing an antenna on his chest, performed a number of walking movements towards and then away from a uniform linear array. It was found that all of the considered diversity schemes provided a worthwhile signal improvement. However, the performance of the diversity systems varied according to the switching threshold that was adopted. To model the fading envelope observed at the output of each of the combiners, we have applied diversity specific equations developed under the assumption of Nakagami-$m$ fading. As a measure of the goodness-of-fit, the Kullback-Leibler divergence between the empirical and theoretical probability density functions (PDFs) was calculated and found to be close to 0. To assist with the interpretation of the goodness-of-fit achieved in this study, the standard deviation, $\sigma$, of a zero-mean, $\sigma^2$ variance Gaussian PDF used to approximate a zero-mean, unit variance Gaussian PDF is also presented. These were generally quite close to 1 indicating that the theoretical models provided an adequate fit to the measured data.

Phase equilibria modeling for biofuels production

Em consequência de uma série de problemas ambientais, económicos e políticos relacionados com o uso de combustíveis convencionais, vários países estão agora a focar as suas atenções em combustíveis alternativos. O biodiesel está na linha da frente das alternativas ao petróleo no sector dos transportes, sendo considerado uma opção a curto prazo visto que o seu preço é competitivo e não são necessárias mudanças nos motores para implementar o seu uso. De entres os possíveis processos de produzir biodiesel, a reacção de transesterificação com catálise básica é o método preferido. Depois da reacção são sempre necessários processos de purificação de modo ao biodiesel produzido cumprir os standards definidos para os combustíveis alternativos, reduzindo problemas de motor e consequentemente aumentando a sua aceitação por parte dos consumidores. De entre as especificações encontram-se o conteúdo em água, em álcool e em glicerol. Ser-se capaz de descrever correctamente o equilíbrio de fases de sistemas que são de interesse para os processos de purificação de biodieseis numa gama alargada de condições termodinâmicas é uma condição necessária para uma correcta simulação do processo industrial, de modo a se atingir uma elevada produtividade a baixos custos de operação. O uso de moléculas oxigenadas como combustíveis representa uma alteração significativa em termos da termodinâmica de soluções. Para combustíveis baseados em petróleo as equações de estado cúbicas e os modelos clássicos de coeficientes de actividade mostraram ser apropriados, no entanto para combustíveis novos como o biodiesel, sendo mais complexos do ponto de vista das interacções intermoleculares com formação de dipolos e pontes de hidrogénio, são necessários modelos termodinâmicos mais complexos para descrever essas interacções. Neste trabalho a CPA EoS (Cubic-Plus-Association Equation of State) será desenvolvida de modo a permitir uma descrição adequada dos equilíbrios líquido-vapor e líquido-líquido para uma serie de sistemas binários e multicomponentes contendo água, ácidos gordos, ésteres de ácidos gordos, glicerol e álcoois.

Measurement and modeling of hyperfine properties in ferroic materials

This thesis presents the results of perturbed angular correlation (PAC) experiments , an experimental technique which measures the hyperfine interaction at probes (radioactive ions implanted in the materials to study), from which one infers local information on an atomic scale. Furthermore, abinitio calculations using density functional theory electronic obtain results that directly complement the experiments, and are also used for theoretical research. These methods were applied in two families of materials. The manganites, with the possible existence of magnetic, charge, orbital and ferroelectric orders, are of fundamental and technological interest. The experimental results are obtained in the alkaline-earth manganites (Ca, Ba, Sr), with special interest due to the structural variety of possible polymorphs. With probes of Cd and In the stability of the probe and its location in a wide temperature range is established and a comparison with calculations allows the physical interpretation of the results. Calculations of hyperfine properties in rare-earth manganites are also presented. The second type of materials in which hyperfine properties were studied are the Manganese pnictides: MnAs, MnSb, and MnBi, compounds in which magnetism is fundamental. The experimental results obtained mainly consider the MnAs compound, whose magneto-structural transition is of great interest. The transition is analyzed in detail with the local resolution characteristic of the technique, obtaining information of the character of the transition also with complementary, more conventional techniques. The last work in this thesis uses only the first principles calculations, continuing the theme of the hyperfine interactions, but this time with respect to ferroelectrics. Several transition metal oxides with perovskite or distorted structures are considered. The electric field gradient which exists due to the quadrupole interaction in nuclei is related to the spontaneous electric polarization, the main quantity measured in ferroelectrics. This study provides a fundamental theoretical basis for previous empirical studies, suggesting new directions for research in ferroelectrics and multiferroics using techniques which measure the electric field gradient.

Machine learning methods for quantitative structure-property relationship modeling

Tese de doutoramento, Informática (Bioinformática), Universidade de Lisboa, Faculdade de Ciências, 2014

Ab initio modeling and molecular dynamics simulation of the alpha 1b-adrenergic receptor activation.

This work describes the ab initio procedure employed to build an activation model for the alpha 1b-adrenergic receptor (alpha 1b-AR). The first version of the model was progressively modified and complicated by means of a many-step iterative procedure characterized by the employment of experimental validations of the model in each upgrading step. A combined simulated (molecular dynamics) and experimental mutagenesis approach was used to determine the structural and dynamic features characterizing the inactive and active states of alpha 1b-AR. The latest version of the model has been successfully challenged with respect to its ability to interpret and predict the functional properties of a large number of mutants. The iterative approach employed to describe alpha 1b-AR activation in terms of molecular structure and dynamics allows further complications of the model to allow prediction and interpretation of an ever-increasing number of experimental data.

Modeling reaction kinetics and mass transfer in ozonation in water solutions

This dissertation is based on four articles dealing with modeling of ozonation. The literature part of this considers some models for hydrodynamics in bubble column simulation. A literature review of methods for obtaining mass transfer coefficients is presented. The methods presented to obtain mass transfer are general models and can be applied to any gas-liquid system. Ozonation reaction models and methods for obtaining stoichiometric coefficients and reaction rate coefficients for ozonation reactions are discussed in the final section of the literature part. In the first article, ozone gas-liquid mass transfer into water in a bubble column was investigated for different pH values. A more general method for estimation of mass transfer and Henry’s coefficient was developed from the Beltrán method. The ozone volumetric mass transfer coefficient and the Henry’s coefficient were determined simultaneously by parameter estimation using a nonlinear optimization method. A minor dependence of the Henry’s law constant on pH was detected at the pH range 4 - 9. In the second article, a new method using the axial dispersion model for estimation of ozone self-decomposition kinetics in a semi-batch bubble column reactor was developed. The reaction rate coefficients for literature equations of ozone decomposition and the gas phase dispersion coefficient were estimated and compared with the literature data. The reaction order in the pH range 7-10 with respect to ozone 1.12 and 0.51 the hydroxyl ion were obtained, which is in good agreement with literature. The model parameters were determined by parameter estimation using a nonlinear optimization method. Sensitivity analysis was conducted using object function method to obtain information about the reliability and identifiability of the estimated parameters. In the third article, the reaction rate coefficients and the stoichiometric coefficients in the reaction of ozone with the model component p-nitrophenol were estimated at low pH of water using nonlinear optimization. A novel method for estimation of multireaction model parameters in ozonation was developed. In this method the concentration of unknown intermediate compounds is presented as a residual COD (chemical oxygen demand) calculated from the measured COD and the theoretical COD for the known species. The decomposition rate of p-nitrophenol on the pathway producing hydroquinone was found to be about two times faster than the p-nitrophenol decomposition rate on the pathway producing 4- nitrocatechol. In the fourth article, the reaction kinetics of p-nitrophenol ozonation was studied in a bubble column at pH 2. Using the new reaction kinetic model presented in the previous article, the reaction kinetic parameters, rate coefficients, and stoichiometric coefficients as well as the mass transfer coefficient were estimated with nonlinear estimation. The decomposition rate of pnitrophenol was found to be equal both on the pathway producing hydroquinone and on the path way producing 4-nitrocathecol. Comparison of the rate coefficients with the case at initial pH 5 indicates that the p-nitrophenol degradation producing 4- nitrocathecol is more selective towards molecular ozone than the reaction producing hydroquinone. The identifiability and reliability of the estimated parameters were analyzed with the Marcov chain Monte Carlo (MCMC) method. @All rights reserved. No part of the publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior permission of the author.