Prediction of antibiotic activity and synthesis of new pentadecapeptides based on lactoferricins.

The antibacterial activity against Escherichia coli and Staphylococcus aureus has been studied for a number of modified pentadecapeptides based on lactoferricins of different origin. The peptides were classified by multivariate methods and quantitative structure-activity relationships (QSAR) were developed using theoretically derived variables for the amino acids. For the modified peptides based on bovine lactoferricin (LFB) a model was calculated and used for prediction of new peptides that were then tested for antibacterial activity in order to improve peptide activity and to check the validity of the model. Models were also calculated including lactoferricins of different origin. Theories of the mechanism of action of the peptides are briefly discussed.     

Variations in antigen-antibody association kinetics as a function of pH and salt concentration: a QSAR and molecular modeling study

The relationship between three environmental factors (ionic strength, pH, and temperature) and antigen-antibody binding kinetics was investigated using QSAR (quantitative structure-activity relationship) and molecular modeling approaches. The interaction used for this analysis is that between the camel antibody fragment cAbLys3 and lysozyme. Binding kinetics were measured using a Biacore 2000 instrument, at NaCl concentrations between 50 and 500 mM, at pH's between 5 and 10, and at temperatures between 15 and 30 degrees C, according to multivariate experimental designs. Variations in kinetic on- and off-rate parameters were up to 400- and 16-fold, respectively. Mathematical models that relate log k(on) to experimental conditions were developed. They indicated an influence of all three factors, with a clear dependency between pH and NaCl concentration for their effect on k(on). These models were able to predict on-rate parameters under new experimental conditions. Titration calculations using continuum electrostatics were performed on the crystallographic structures of the isolated and bound proteins to gain structural insight for the on-rate enhancement observed at pH <6.5 and low salt concentrations. These calculations rule out electrostatic steering linked to global and/or local charge variations in the molecules as the factor responsible for the on-rate enhancement at low pH. His 111 of cAbLys3, located at the binding interface, can adopt two side chain orientations with different intramolecular contacts. The results of the calculations suggest an alternative mechanism whereby the conformation of the interfacial His 111 depends on the charge, and these differences in conformation may influence the solvation energy and the subsequent binding kinetics. Our results stress the complex relationship between environmental conditions and molecular binding properties.     

Comparative quantitative structure-activity study of radical scavengers

Classic and three-dimensional (3-D) QSAR analyses of 13 radical scavengers (1-13) were performed to derive two classic, two Apex-3-D and one comparative field analysis (CoMFA) models. Two classical models with predictive cross-validated r2 (Q2) over 0.96 indicated that the activity was attributed to the electronic COH and ELUMO, steric molar refractivity (MR) and lipophilic log P. Three-dimensional quantitative structure-activity relationship (3-D-QSAR) studies were performed by 3-D pharmacophore generation (Apex-3-D) and CoMFA techniques. For Apex-3-D studies, two best models with high Q2 (0.94 and 0.97) were yielded. Structural properties contributing to the activity were not only lipophilic but also the optimum steric property and geometry of side-chain composition. For CoMFA studies, the sp3 C(+1) probe provided the best Q2 of 0.79 with steric and electrostatic contributions of 42.3 and 57.7%, respectively. The activity of four new compounds (14-17) not included in the derivation were predicted with these models. Although the derived models were from limited data, the statistic relation was predictive. The linear correlations between the experimental IC50 values and the predicted values from classical and Apex-3-D models were found to be high and significant. The predicted activity of 17 from CoMFA was much lower than the experimental value; this deviation occurred according to the missing of hydrophobic field in standard CoMFA study. In vitro and ex vivo antilipid peroxidation in mouse brain and ESR studies of 14-17 were investigated for the radical-scavenging ability. The difference between the in vitro results, antilipid peroxidation and electron spin resonance (ESR) and ex vivo results in coumarin series was found. Thus, other properties for good bioavailability besides log P should also be taken into consideration.     

Molecular docking and 3D QSAR studies on 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes based on the structural modeling of human CCR5 receptor

In the present study, we have used an approach combining protein structure modeling, molecular dynamics (MD) simulation, automated docking, and 3D QSAR analyses to investigate the detailed interactions of CCR5 with their antagonists. Homology modeling and MD simulation were used to build the 3D model of CCR5 receptor based on the high-resolution X-ray structure of bovine rhodopsin. A series of 64 CCR5 antagonists, 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes, were docked into the putative binding site of the 3D model of CCR5 using the docking method, and the probable interaction model between CCR5 and the antagonists were obtained. The predicted binding affinities of the antagonists to CCR5 correlate well with the antagonist activities, and the interaction model could be used to explain many mutagenesis results. All these indicate that the 3D model of antagonist-CCR5 interaction is reliable. Based on the binding conformations and their alignment inside the binding pocket of CCR5, three-dimensional structure-activity relationship (3D QSAR) analyses were performed on these antagonists using comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) methods. Both CoMFA and CoMSIA provide statistically valid models with good correlation and predictive power. The q(2)(r(cross)(2)) values are 0.568 and 0.587 for CoMFA and CoMSIA, respectively. The predictive ability of these models was validated by six compounds that were not included in the training set. Mapping these models back to the topology of the active site of CCR5 leads to a better understanding of antagonist-CCR5 interaction. These results suggest that the 3D model of CCR5 can be used in structure-based drug design and the 3D QSAR models provide clear guidelines and accurate activity predictions for novel antagonist design.     

Chromatographic evaluation of the toxicity in fish of pesticides

Ecotoxicity assessment is essential before placing new chemical substances on the market. An investigation of the use of the chromatographic retention (log k) in biopartitioning micellar chromatography (BMC) as an in vitro approach to evaluate the toxicity in fish of pesticides (acute toxicity levels as pLC(50)) is proposed. A heterogeneous data set of 85 pesticides from six chemical families with available experimental fish toxicity data (ECOTOX database from U.S. Environmental Protection Agency (EPA)) was used. For pesticides exhibiting non-polar narcosis mechanism in fish (non-specific toxicity), more reliable models and precise pLC(50) estimations are obtained from log k (quantitative retention-activity relationships, QRAR) than from log P (quantitative structure-activity relationships, QSAR) or ECOSAR (ECOSAR program from U.S. EPA).     

Comparative structural and energetic analysis of WW domain-peptide interactions

WW domains are small globular protein interaction modules found in a wide spectrum of proteins. They recognize their target proteins by binding specifically to short linear peptide motifs that are often proline-rich. To infer the determinants of the ligand binding propensities of WW domains, we analyzed 42 WW domains. We built models of the 3D structures of the WW domains and their peptide complexes by comparative modeling supplemented with experimental data from peptide library screens. The models provide new insights into the orientation and position of the peptide in structures of WW domain-peptide complexes that have not yet been determined experimentally. From a protein interaction property similarity analysis (PIPSA) of the WW domain structures, we show that electrostatic potential is a distinguishing feature of WW domains and we propose a structure-based classification of WW domains that expands the existent ligand-based classification scheme. Application of the comparative molecular field analysis (CoMFA), GRID/GOLPE and comparative binding energy (COMBINE) analysis methods permitted the derivation of quantitative structure-activity relationships (QSARs) that aid in identifying the specificity-determining residues within WW domains and their ligand-recognition motifs. Using these QSARs, a new group-specific sequence feature of WW domains that target arginine-containing peptides was identified. Finally, the QSAR models were applied to the design of a peptide to bind with greater affinity than the known binding peptide sequences of the yRSP5-1 WW domain. The prediction was verified experimentally, providing validation of the QSAR models and demonstrating the possibility of rationally improving peptide affinity for WW domains. The QSAR models may also be applied to the prediction of the specificity of WW domains with uncharacterized ligand-binding properties.     

QSAR and docking studies of anthraquinone derivatives by similarity cluster prediction

Forty anthraquinone derivatives have been downloaded from PubChem database and investigated in a quantitative structure-activity relationships (QSAR) study. The models describing log P and LD50 of this set were built up on the hypermolecule scheme that mimics the investigated receptor space; the models were validated by the leave-one-out procedure, in the external test set and in a new version of prediction by using similarity clusters. Molecular docking approach using Lamarckian Genetic Algorithm was made on this class of anthraquinones with respect to 3Q3B receptor. The best scored molecules in the docking assay were used as leaders in the similarity clustering procedure. It is demonstrated that the LD50 data of this set of anthraquinones are related to the binding energies of anthraquinone ligands to the 3Q3B receptor.     

Chemometric modelling of antimalarial activity of aryltriazolylhydroxamates

We have performed quantitative structure–activity relationship (QSAR) and quantitative activity–activity relationship (QAAR) studies for aryltriazolylhydroxamates having antimalarial activity data against both chloroquine-sensitive (D6 clone) and chloroquine-resistant (W2 clone) strains of Plasmodium falciparum to understand the relationships between the biological activity and molecular properties for the design of new compounds. The QSAR studies were performed using 35 compounds among which 26 molecules were taken using k-means clustering technique in the training set for the derivation of the QSAR models and nine molecules were kept as the test-set compounds to evaluate the predictive ability of the derived models. The chemometric tool used for the analysis was the genetic function approximation. The developed models were analysed in terms of their predictive ability, and comparable results were obtained for cross-validated predictive variance (Q ²) and externally predicted variance (R ² _pred) values (0.761 and 0.829, respectively, for the D6 model, 0.708 and 0.748, respectively, for the W2 model and 0.984 and 0.982, respectively for the QAAR model). The QSAR models suggest that the number of methylene groups (between the triazolyl and hydroxamate moieties) and partially negatively charged surface areas of the molecules are important parameters for the antimalarial activity.