Increasing the precision of comparative models with YASARA NOVA--a self-parameterizing force field

One of the conclusions drawn at the CASP4 meeting in Asilomar was that applying various force fields during refinement of template-based models tends to move predictions in the wrong direction, away from the experimentally determined coordinates. We have derived an all-atom force field aimed at protein and nucleotide optimization in vacuo (NOVA), which has been specifically designed to avoid this problem. NOVA resembles common molecular dynamics force fields but has been automatically parameterized with two major goals: (i) not to make high resolution X-ray structures worse and (ii) to improve homology models built by WHAT IF. Force-field parameters were not required to be physically correct; instead, they were optimized with random Monte Carlo moves in force-field parameter space, each one evaluated by simulated annealing runs of a 50-protein optimization set. Errors inherent to the approximate force-field equation could thus be canceled by errors in force-field parameters. Compared with the optimization set, the force field did equally well on an independent validation set and is shown to move in silico models closer to reality. It can be applied to modeling applications as well as X-ray and NMR structure refinement. A new method to assign force-field parameters based on molecular trees is also presented. A NOVA server is freely accessible at http://www.yasara.com/servers     

基于启发式广度优先搜索的SVR参数优化方法研究

支持向量回归机(Support vector regressio,SVR)模型的拟合精度和泛化能力取决于其相关参数的选择,其参数选择实质上是一个优化搜索过程。根据启发式广度优先搜索(Heuristic Breadth first Search,HBFS)算法在求解优化问题上高效的特点,提出了一种以k-fold交叉验证的最小化误差为目标,HBFS为寻优策略的SVR参数选择方法,通过3个基准数据集对该模型进行了仿真实验,结果表明该方法在保证预测精度前提下,大幅度的缩短了训练建模时间,为大样本的SVR参数选择提供了一种新的有效解决方案。     

Use of decoys to optimize an all-atom force field including hydration

A novel method of parameter optimization is proposed. It makes use of large sets of decoys generated for six nonhomologous proteins with different architecture. Parameter optimization is achieved by creating a free energy gap between sets of nativelike and nonnative conformations. The method is applied to optimize the parameters of a physics-based scoring function consisting of the all-atom ECEPP05 force field coupled with an implicit solvent model (a solvent-accessible surface area model). The optimized force field is able to discriminate near-native from nonnative conformations of the six training proteins when used either for local energy minimization or for short Monte Carlo simulated annealing runs after local energy minimization. The resulting force field is validated with an independent set of six nonhomologous proteins, and appears to be transferable to proteins not included in the optimization; i.e., for five out of the six test proteins, decoys with 1.7- to 4.0-Å all-heavy-atom root mean-square deviations emerge as those with the lowest energy. In addition, we examined the set of misfolded structures created by Park and Levitt using a four-state reduced model. The results from these additional calculations confirm the good discriminative ability of the optimized force field obtained with our decoy sets.     

Determination method of the balance of the secondary-structure-forming tendencies of force fields

We propose a new method of optimisation of backbone torsion-energy parameters in the force field for molecular simulations of protein systems. This method is based on the idea of balancing the secondary-structure-forming tendencies, namely, those of α-helix and β-sheet structures. We perform a minimisation of the backbone dihedral angle-based root-mean-square deviation of the helix and β structure regions in many protein structures. As an example, we optimised the backbone torsion-energy parameters of AMBER parm96 force field using 100 protein molecules from the Protein Data Bank. We then performed folding simulations of α-helical and β-hairpin peptides, using the optimised force field. The results imply that the new force-field parameters give structures more consistent with the experimental implications than the original AMBER parm96 force field.     

Continuous optimization algorithms for tuning real and integer parameters of swarm intelligence algorithms

The performance of optimization algorithms, including those based on swarm intelligence, depends on the values assigned to their parameters. To obtain high performance, these parameters must be fine-tuned. Since many parameters can take real values or integer values from a large domain, it is often possible to treat the tuning problem as a continuous optimization problem. In this article, we study the performance of a number of prominent continuous optimization algorithms for parameter tuning using various case studies from the swarm intelligence literature. The continuous optimization algorithms that we study are enhanced to handle the stochastic nature of the tuning problem. In particular, we introduce a new post-selection mechanism that uses F-Race in the final phase of the tuning process to select the best among elite parameter configurations. We also examine the parameter space of the swarm intelligence algorithms that we consider in our study, and we show that by fine-tuning their parameters one can obtain substantial improvements over default configurations.     

Protein structure alignment by deterministic annealing

MOTIVATION: Protein structure alignment is one of the most important computational problems in molecular biology and plays a key role in protein structure prediction, fold family classification, motif finding, phylogenetic tree reconstruction and so on. From the viewpoint of computational complexity, a pairwise structure alignment is also a NP-hard problem, in contrast to the polynomial time algorithm for a pairwise sequence alignment. RESULTS: We propose a method for solving the structure alignment problem in an accurate manner at the amino acid level, based on a mean field annealing technique. We define the structure alignment as a mixed integer-programming (MIP) problem. By avoiding complicated combinatorial computation and exploiting the special structure of the continuous partial problem, we transform the MIP into a reduced non-linear continuous optimization problem (NCOP) with a much simpler form. To optimize the reduced NCOP, a mean field annealing procedure is adopted with a modified Potts model, whose solution is generally identical to that of the MIP. There is no 'soft constraint' in our mean field model and all constraints are automatically satisfied throughout the annealing process, thereby not only making the optimization more efficient but also eliminating many unnecessary parameters that depend on problems and usually require careful tuning. A number of benchmark examples are tested by the proposed method with comparisons to several existing approaches.     

Determination of best-fit potential parameters for a reactive force field using a genetic algorithm

The ReaxFF interatomic potential, used for organic materials, involves more than 600 adjustable parameters, the best-fit values of which must be determined for different materials. A new method of determining the set of best-fit parameters for specific molecules containing carbon, hydrogen, nitrogen and oxygen is presented, based on a parameter reduction technique followed by genetic algorithm (GA) minimization. This work has two novel features. The first is the use of a parameter reduction technique to determine which subset of parameters plays a significant role for the species of interest; this is necessary to reduce the optimization space to manageable levels. The second is the application of the GA technique to a complex potential (ReaxFF) with a very large number of adjustable parameters, which implies a large parameter space for optimization. In this work, GA has been used to optimize the parameter set to determine best-fit parameters that can reproduce molecular properties to within a given accuracy. As a test problem, the use of the algorithm has been demonstrated for nitromethane and its decomposition products.     

Simplified method for the computation of parameters of power-law rate equations from time-series

Modeling biological processes from time-series data is a resourceful procedure which has received much attention in the literature. For models established in the context of non-linear differential equations, parameter-dependent phenomenological tentative response functions are tested by comparing would-be solutions of those models to the experimental time-series. Those values of the parameters for which a tested solution is a best fit are then retained. It is done with the help of some appropriate optimization algorithm which simplifies the searching procedure within the range of variability of the parameters that are to be estimated. The procedure works well in problems with a small number of adjustable parameters or/and with narrow searching ranges. However, it may start to be problematic for models with a large number of problem parameters inasmuch as convergence to the best fit is not necessarily ensured. In this case, a reduction in size of the parameter estimation problem must be undertaken. We presently address this issue by proposing a systematic procedure that does so in problems in which the system's response to a sufficiently small pulse perturbation of steady-state can be obtained. The response is then assumed to be a solution of the linearized equations, the Jacobian of which can be retrieved by a simple multilinear regression. The calculated n(2) Jacobian entries provide as many relationships among problem parameters, thus cutting substantially the size of the starting problem. After this preliminary treatment is applied, only (kappa-n(2)) of the initial kappa adjustable parameters are left for evaluation by means of a non-linear optimization procedure. The benefits of the present variant are both in economy of computation and in accuracy in determining the parameter values. The performance of the method is established under different circumstances. It is illustrated in the context of power-law rates, although this does not preclude its applicability to more general functional responses.     

Model selection and efficiency testing for normalization of cDNA microarray data

In this study we present two novel normalization schemes for cDNA microarrays. They are based on iterative local regression and optimization of model parameters by generalized cross-validation. Permutation tests assessing the efficiency of normalization demonstrated that the proposed schemes have an improved ability to remove systematic errors and to reduce variability in microarray data. The analysis also reveals that without parameter optimization local regression is frequently insufficient to remove systematic errors in microarray data.     

A geometry force field which converts low-resolution X-ray models to structures with properties found at ultra high resolution

A geometry optimization force field was developed using ultra high-resolution structures and tested using high- and low-resolution X-ray structures. Protein and small molecule X-ray data was used. When applied to ultra high-resolution structures the force field conserves the internal geometry and local strain energy. When applied to low-resolution structures there is a small change in geometry accompanied by a large drop in local strain energy. Although optimization causes only small structural changes in low-resolution X-ray models, it dramatically modifies profiles for hydrogen bonding, Van der Waals contact, bonded geometry, and local strain energy, making them almost indistinguishable from those found at high resolution. Further insight into the effect of the force field was obtained by comparing geometries of homologous proteins before and after geometry optimization. Optimization causes homologous regions of structures to become similar in internal geometry and energies. Once again, the changes only require small atomic movements. These findings provide insights into the structure of molecular complexes. The new force field contains only short-range interatomic potential functions. Its effectiveness shows that local geometries are determined by short-range interactions which are well modeled by the force field. Potential applications of this study include detection of possible structural errors, correction of errors with minimal change in geometry, improved understanding and prediction of the effects of modifying ligands or proteins, and computational addition of structural water.     

RankAggreg, an R package for weighted rank aggregation

Background

Determining the parameters of a mathematical model from quantitative measurements is the main bottleneck of modelling biological systems. Parameter values can be estimated from steady-state data or from dynamic data. The nature of suitable data for these two types of estimation is rather different. For instance, estimations of parameter values in pathway models, such as kinetic orders, rate constants, flux control coefficients or elasticities, from steady-state data are generally based on experiments that measure how a biochemical system responds to small perturbations around the steady state. In contrast, parameter estimation from dynamic data requires time series measurements for all dependent variables. Almost no literature has so far discussed the combined use of both steady-state and transient data for estimating parameter values of biochemical systems.

Results

In this study we introduce a constrained optimization method for estimating parameter values of biochemical pathway models using steady-state information and transient measurements. The constraints are derived from the flux connectivity relationships of the system at the steady state. Two case studies demonstrate the estimation results with and without flux connectivity constraints. The unconstrained optimal estimates from dynamic data may fit the experiments well, but they do not necessarily maintain the connectivity relationships. As a consequence, individual fluxes may be misrepresented, which may cause problems in later extrapolations. By contrast, the constrained estimation accounting for flux connectivity information reduces this misrepresentation and thereby yields improved model parameters.

Conclusion

The method combines transient metabolic profiles and steady-state information and leads to the formulation of an inverse parameter estimation task as a constrained optimization problem. Parameter estimation and model selection are simultaneously carried out on the constrained optimization problem and yield realistic model parameters that are more likely to hold up in extrapolations with the model.     

Efficient measurement error correction with spatially misaligned data

Association studies in environmental statistics often involve exposure and outcome data that are misaligned in space. A common strategy is to employ a spatial model such as universal kriging to predict exposures at locations with outcome data and then estimate a regression parameter of interest using the predicted exposures. This results in measurement error because the predicted exposures do not correspond exactly to the true values. We characterize the measurement error by decomposing it into Berkson-like and classical-like components. One correction approach is the parametric bootstrap, which is effective but computationally intensive since it requires solving a nonlinear optimization problem for the exposure model parameters in each bootstrap sample. We propose a less computationally intensive alternative termed the "parameter bootstrap" that only requires solving one nonlinear optimization problem, and we also compare bootstrap methods to other recently proposed methods. We illustrate our methodology in simulations and with publicly available data from the Environmental Protection Agency.     

Mathematical and kinetic modeling of biofilm reactor based on ant colony optimization

Inverse estimation of model parameters via mathematical modeling route, known as inverse modeling (IM), is an attractive alternative approach to the experimental methods. This approach makes use of efficient optimization techniques in the course of solution of an inverse problem with the aid of measured data. In this study, a novel optimization method based on ant colony optimization (ACO), denoted by ACO-IM, is presented for inverse estimation of kinetic and film thickness parameters of biofilm models that describe an experimental fixed bed anaerobic reactor. The proposed optimization method for parameter estimation emulates the fact that ants are capable of finding the shortest path from a food source to their nest by depositing a trial of pheromone during their walk. The efficacy of the ACO-IM for numerical estimation of bio-kinetic parameters is demonstrated through its application for the anaerobic treatment of industry wastewater in a fixed bed biofilm process. The results explain the rigorousness of mathematical models, the form of kinetic and film thickness models and the type of packing to be used with the biofilm process for accurate determination of kinetic and film thickness parameters so as to ensure reliable predictive performance of the biofilm reactor models.     

Parameterization and optimization of the menthol force field for molecular dynamics simulations

Menthol’s various biological properties render it a useful component for medical and cosmetological applications, while its three centers of asymmetry mean that it can be used in a range of organic reactions. Menthol-substituted ionic liquids (ILs) have been found to exhibit promising antimicrobial and antielectrostatic properties, as well as being useful in organic catalysis and biochemical studies. However, so far, a force field designed and validated specifically for the menthol molecule has not been constructed. In the present work, the validation and optimization of force field parameters with regard to the ability to reproduce the macroscopic properties of menthol is presented. The set of optimized potentials for liquid simulations all atom (OPLS-AA) compatible parameters was tested and carefully tuned. The refinement of parameters included fitting of partial atomic charges, optimization of Lennard-Jones parameters, and recalculation of the dihedral angle parameters needed to reproduce quantum energy profiles. To validate the force field, a variety of physicochemical properties were calculated for liquid menthol. Both thermodynamic and kinetic properties were taken into account, including density, surface tension, enthalpy of vaporization, and shear viscosity. The obtained force field was proven to accurately reproduce the properties of the investigated compound while being fully compatible with the OPLS-AA force field.

     

基于混合优化算法的经颅磁刺激线圈阵列设计

为了研究对经颅磁刺激激励线圈聚焦性能的优化,利用混合优化算法与CST软件的外部通信接口,建立优化的激励线圈模型。依据多信道线圈阵列方法,利用磁场叠加原理,对影响磁场分布的线圈可调参数进行分析,结合混合优化算法对可调参数进行优化。结果对比显示,经优化的线圈阵列有良好的磁聚焦性,其刺激强度与聚焦程度都有了不同程度提高。可用于改善TMS系统聚焦性能,实验有助于进一步探索全面优化激励线圈的空间结构。     

Protein fold recognition without Boltzmann statistics or explicit physical basis.

We present a fast method for finding optimal parameters for a low-resolution (threading) force field intended to distinguish correct from incorrect folds for a given protein sequence. In contrast to other methods, the parameterization uses information from >10(7) misfolded structures as well as a set of native sequence-structure pairs. In addition to testing the resulting force field's performance on the protein sequence threading problem, results are shown that characterize the number of parameters necessary for effective structure recognition.     

Optimal Experimental Design and Accuracy of Parameter Estimation for Nonlinear Regression Models Used in Long-term Selection

An experimental design problem is considered for the analysis of long-term selection experiments with nonlinear regression models. For a 3-parametric exponential regression function whose parameters have also a reasonable biological interpretation approximate formulas for the determination of the necessary number of observations at each generation are constructed in such a way that the half expected length of an (1 — α)-confidence interval for a chosen parameter is not greater than a given value. In this sense the accuracy of the parameter estimators can be described.     

Cellular traction force recovery: An optimal filtering approach in two-dimensional Fourier space

Quantitative estimation of cellular traction has significant physiological and clinical implications. As an inverse problem, traction force recovery is essentially susceptible to noise in the measured displacement data. For traditional procedure of Fourier transform traction cytometry (FTTC), noise amplification is accompanied in the force reconstruction and small tractions cannot be recovered from the displacement field with low signal-noise ratio (SNR). To improve the FTTC process, we develop an optimal filtering scheme to suppress the noise in the force reconstruction procedure. In the framework of the Wiener filtering theory, four filtering parameters are introduced in two-dimensional Fourier space and their analytical expressions are derived in terms of the minimum-mean-squared-error (MMSE) optimization criterion. The optimal filtering approach is validated with simulations and experimental data associated with the adhesion of single cardiac myocyte to elastic substrate. The results indicate that the proposed method can highly enhance SNR of the recovered forces to reveal tiny tractions in cell-substrate interaction.     

Optimisation of OPLS–UA force-field parameters for protein systems using protein data bank

We have previously proposed a method for refining force-field parameters of protein systems, which consists of minimising the summation of the square of the force acting on each atom in the proteins with the structures from the protein data bank (PDB). The results showed that the modified force-field parameters for all-atom model gave structures more consistent with the experimental implications than the original force fields. In this work, we applied this method and a new method to the OPLS–UA force field. In the new method, we perform a minimisation of the average of the root-mean-square deviation of various protein structures from the native structure. We selected some torsion-energy parameters for this optimisation, and 100 molecules from the PDB were used. The results imply that the new force-field parameters gave structures of two peptides more consistent with the experimental implications for the secondary structure-forming tendencies than the original OPLS–UA force field.