A new strategy to improve the predictive ability of the local lazy regression and its application to the QSAR study of melanin‐concentrating hormone receptor 1 antagonists

In the quantitative structure‐activity relationship (QSAR) study, local lazy regression (LLR) can predict the activity of a query molecule by using the information of its local neighborhood without need to produce QSAR models a priori. When a prediction is required for a query compound, a set of local models including different number of nearest neighbors are identified. The leave‐one‐out cross‐validation (LOO‐CV) procedure is usually used to assess the prediction ability of each model, and the model giving the lowest LOO‐CV error or highest LOO‐CV correlation coefficient is chosen as the best model. However, it has been proved that the good statistical value from LOO cross‐validation appears to be the necessary, but not the sufficient condition for the model to have a high predictive power. In this work, a new strategy is proposed to improve the predictive ability of LLR models and to access the accuracy of a query prediction. The bandwidth of k neighbor value for LLR is optimized by considering the predictive ability of local models using an external validation set. This approach was applied to the QSAR study of a series of thienopyrimidinone antagonists of melanin‐concentrating hormone receptor 1. The obtained results from the new strategy shows evident improvement compared with the commonly used LOO‐CV LLR methods and the traditional global linear model. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Carbonyl insertion reaction into the Pt?C bond in heterobimetallic Pt(SnCl3)(PH3)2(CO)(CH3) compound: Theoretical study

Quantum‐mechanical calculations were carried out at the MP4(SDQ)//MP2 level of theory to determine the energies and reaction mechanism for the carbonyl insertion reaction (second step in the olefin hydroformylation catalytic cycle), using a heterobimetallic Pt(SnCl₃)(PH₃)₂(CO)(CH₃) compound as a model catalytic species. The results show that this reaction proceeds through a three‐center transition state, with an activation energy of 26.4 kcal/mol, followed by an intramolecular rearrangement to the square‐planar cis‐Pt(SnCl₃)(PH₃)₂(MeCO) metal–acyl product. Analysis of the nature of the bonds shows that there is a negligible participation of the tin d‐orbitals in the formation of the Pt Sn bond. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 668–674, 2000     

Connectivity and binding‐site recognition: Applications relevant to drug design

Here, we describe a family of methods based on residue–residue connectivity for characterizing binding sites and apply variants of the method to various types of protein–ligand complexes including proteases, allosteric‐binding sites, correctly and incorrectly docked poses, and inhibitors of protein–protein interactions. Residues within ligand‐binding sites have about 25% more contact neighbors than surface residues in general; high‐connectivity residues are found in contact with the ligand in 84% of all complexes studied. In addition, a k‐means algorithm was developed that may be useful for identifying potential binding sites with no obvious geometric or connectivity features. The analysis was primarily carried out on 61 protein–ligand structures from the MEROPS protease database, 250 protein–ligand structures from the PDBSelect (25%), and 30 protein–protein complexes. Analysis of four proteases with crystal structures for multiple bound ligands has shown that residues with high connectivity tend to have less variable side‐chain conformation. The relevance to drug design is discussed in terms of identifying allosteric‐binding sites, distinguishing between alternative docked poses and designing protein interface inhibitors. Taken together, this data indicate that residue–residue connectivity is highly relevant to medicinal chemistry. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Quantitative structure–activity relationship model for prediction study of corrosion inhibition efficiency using two‐stage sparse multiple linear regression

A new quantitative structure–activity relationship (QSAR) of the inhibition of mild steel corrosion in 1 M hydrochloric acid using furan derivatives was developed by proposing two‐stage sparse multiple linear regression. The sparse multiple linear regression using ridge penalty and sparse multiple linear regression using elastic net (SMLRE) were used to develop the QSAR model. The results show that the SMLRE‐based model possesses high predictive power compared with sparse multiple linear regression using ridge penalty‐based model according to the mean‐squared errors for both training and test datasets, leave‐one‐out internal validation (Q²_int = 0.98), and external validation (Q²_ext = 0.95). In addition, the results of applicability domain assessment using the leverage approach reveal a reliable and robust SMLRE‐based model. In conclusion, the developed QSAR model using SMLRE can be efficiently used in the studies of corrosion inhibition efficiency. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessment and acceleration of binding energy calculations for protein–ligand complexes by the fragment molecular orbital method

In the field of drug discovery, it is important to accurately predict the binding affinities between target proteins and drug applicant molecules. Many of the computational methods available for evaluating binding affinities have adopted molecular mechanics‐based force fields, although they cannot fully describe protein–ligand interactions. A noteworthy computational method in development involves large‐scale electronic structure calculations. Fragment molecular orbital (FMO) method, which is one of such large‐scale calculation techniques, is applied in this study for calculating the binding energies between proteins and ligands. By testing the effects of specific FMO calculation conditions (including fragmentation size, basis sets, electron correlation, exchange‐correlation functionals, and solvation effects) on the binding energies of the FK506‐binding protein and 10 ligand complex molecule, we have found that the standard FMO calculation condition, FMO2‐MP2/6‐31G(d), is suitable for evaluating the protein–ligand interactions. The correlation coefficient between the binding energies calculated with this FMO calculation condition and experimental values is determined to be R = 0.77. Based on these results, we also propose a practical scheme for predicting binding affinities by combining the FMO method with the quantitative structure–activity relationship (QSAR) model. The results of this combined method can be directly compared with experimental binding affinities. The FMO and QSAR combined scheme shows a higher correlation with experimental data (R = 0.91). Furthermore, we propose an acceleration scheme for the binding energy calculations using a multilayer FMO method focusing on the protein–ligand interaction distance. Our acceleration scheme, which uses FMO2‐HF/STO‐3G:MP2/6‐31G(d) at R_int = 7.0 Å, reduces computational costs, while maintaining accuracy in the evaluation of binding energy. © 2015 Wiley Periodicals, Inc.     

QSID Tool: a new three-dimensional QSAR environmental tool

QSID Tool (Quantitative structure–activity relationship tool for Innovative Discovery) was developed to provide an easy-to-use, robust and high quality environmental tool for 3D QSAR. Predictive models developed with QSID Tool can accelerate the discovery of lead compounds by enabling researchers to formulate and test hypotheses for optimizing efficacy and increasing drug safety and bioavailability early in the process of drug discovery. QSID Tool was evaluated by comparison with SYBYL^® using two different datasets derived from the inhibitors of Trypsin (Böhm et al., J Med Chem 42:458, 1999) and p38-MAPK (Liverton et al., J Med Chem 42:2180, 1999; Romeiro et al., J Comput Aided Mol Des 19:385, 2005; Romeiro et al., J Mol Model 12:855, 2006). The results suggest that QSID Tool is a useful model for the prediction of new analogue activities.     

Quantum-chemical foundations of the topological substructural molecular design

The topological substructural molecular design (TOPS-MODE) approach is formulated as a tight-binding quantum-chemical method. The approach is based on certain postulates that permit to express any molecular property as a function of the spectral moments of certain types of molecular and environment-dependent energies. We use several empirical potentials to account for these intrinsic and external molecular energies. We prove that any molecular property expressed in terms of a quantitative structure-property and structure-activity relationships (QSPR/QSAR) model developed by using the TOPS-MODE method can be expressed as a bond additivity function. In addition, such a property can also be expressed as a substructural cluster expansion function. The conditions for such bond contributions being transferable are also analyzed here. Several new statistical-mechanical electronic functions are introduced as well as a bond-bond thermal Green's function or a propagator accounting for the electronic hopping between pairs of bonds. All these new concepts are applied to the development and application of a new QSAR model for describing the toxicity of polyhalogenated-dibenzo-1,4-dioxins. The QSAR model obtained displays a significant robustness and predictability. It permits an easy structural interpretation of the structure-activity relationship in terms of bond additivity functions, which display some resemblances with other theoretical parameters obtained from first principle quantum-chemical methods.     

PEPCAT—A new tool for conformational analysis of peptides

We report a new technique for the efficient analysis and visualization of peptide and protein conformations and conformational relationships, which we have implemented in a computer program called PEPCAT. PEPCAT (an abbreviation for Peptide Conformational Analysis Tool) provides a simple, graphical, and flexible framework that allows the user to define a specific structural feature or juxtaposition of amino acids and to follow the fate of the motif during a molecular dynamics simulation. Here we describe the PEPCAT analysis of the effects of environmental and chemical modifications on conformational preferences of a regulator of hemopoiesis, namely the pentapeptide pyro‐EEDCK, and of a conformational transition in the immunosuppressant drug cyclosporin A. PEPCAT, however, can be applied to the conformational analysis of peptides and proteins in general. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 446–461, 2000