Quantitative structure-property relationships of potentially bioactive fluoro phospho-silicate glasses.

In this work, the glass transition temperature and chemical durability of bioactive phospho-silicate glasses were experimentally determined and correlated to the structural descriptor Fnet derived from classical molecular dynamics simulations. The replacement of CaF2 for Na2O in the parent glass 45S5 enhances both chemical durability and density, while the replacement of CaF2 for CaO lowers chemical durability. The proposed descriptor, Fnet, provides satisfactorily correlations with glass transition temperature and chemical durability over a wide range of compositions.

Elucidation of the structural role of fluorine in potentially bioactive glasses by experimental and computational investigation.

Glasses belonging to the Na(2)O-CaO-P(2)O(5)-SiO(2) system and modified by CaF(2) substitution for CaO and Na(2)O alternatively, were synthesized and characterized experimentally and computationally. The results of molecular dynamics simulations show that fluorine is almost exclusively bonded to modifier cations (Ca and Na) with coordination number close to 4. A similar mean coordination number value is found in the crystal phases obtained by means of thermal treatment at fixed temperature. Addition of fluorine increases the polymerization of silicate tetrahedra by removing modifiers from the siliceous matrix. No appreciable amount of Si-F bonds are detected.

Properties of zinc releasing surfaces for clinical applications.

Two series of glasses of general formula (2-p) SiO2.1.1Na2O.CaO.pP2O5.xZnO (p=0.10, 0.20; x=0.0, 0.16, 0.35, and 0.78) have been analyzed for physico-chemical surface features before and after contact with simulated body fluid, morphological characteristics, and osteoblast-like cells behavior when cultured on them. The resulted good cell adhesion and growth, along with nonsignificant changes of the focal contacts, allow the authors to indicate HZ5 and HP5Z5 glasses as the ones having optimal ratio of Zn/P to maintain acceptable cell behavior, comparable to the bioactive glass (Bioglass) used as a control; results are also rationalized by means of three-dimensional models derived by molecular dynamic simulations, with decomposition and conversion rates optimized with respect to the parent Hench's Bioglass.

Theoretical quantitative structure-activity relationships of flavone ligands interacting with cytochrome P450 1A1 and 1A2 isozymes.

Theoretical descriptors obtained from quantum mechanical calculations on isolated ligands in different media and molecular dynamics simulations of ligand-enzyme complexes have been used to obtain a quantitative rationalization of the inhibition of CYP1A2 and CYP1A2 by three series of flavonoids. Predictive models obtained through one-descriptor QSAR studies and mechanistic explanations have been obtained for recognition and selectivity.

Computational approaches to structural and functional analysis of plastocyanin and other blue copper proteins.

Computational techniques are becoming increasingly important in structural and functional biology, in particular as tools to aid the interpretation of experimental results and the design of new systems. This review reports on recent studies employing a variety of computational approaches to unravel the microscopic details of the structure-function relationships in plastocyanin and other proteins belonging to the blue copper superfamily. Aspects covered include protein recognition, electron transfer and protein-solvent interaction properties of the blue copper protein family. The relevance of integrating diverse computational approaches to address the analysis of a complex protein system, such as a cupredoxin metalloprotein, is emphasized.

Theoretical descriptors for the quantitative rationalisation of plastocyanin mutant functional propertiess.

A quantitative rationalisation of the effect of specific amino acids on the recognition process and redox characteristics of plastocyanin towards cytochrome f, as determined by point mutation experiments, has been attempted in this study. To achieve this goal we derived theoretical descriptors directly from the three-dimensional structure of the plastocyanin mutants, in the same manner as it is usually done for small drug-like molecules. The protein descriptors computed can be related to: (a) the electrostatic and dipole-dipole interactions, effective at long distance; (b) polar interactions whose features are encoded by charged partial surface area descriptors; (c) the propensity of the surface residues to form hydrogen bonding interactions; and (d) dispersion and repulsive interactions. Moreover, an estimation of mutation-dependent variation of redox potential observed has been obtained by electrostatic free energy calculations. The quantitative structure-activity relationship (QSAR) models offer structural interpretation of the point mutation experiment responses and can be of help in the design of new protein engineering experiments.

Electrostatic analysis and Brownian dynamics simulation of the association of plastocyanin and cytochrome f.

The oxidation of cytochrome f by the soluble cupredoxin plastocyanin is a central reaction in the photosynthetic electron transfer chain of all oxygenic organisms. Here, two different computational approaches are used to gain new insights into the role of molecular recognition and protein-protein association processes in this redox reaction. First, a comparative analysis of the computed molecular electrostatic potentials of seven single and multiple point mutants of spinach plastocyanin (D42N, E43K, E43N, E43Q/D44N, E59K/E60Q, E59K/E60Q/E43N, Q88E) and the wt protein was carried out. The experimentally determined relative rates (k(2)) for the set of plastocyanin mutants are found to correlate well (r(2) = 0.90 - 0.97) with the computed measure of the similarity of the plastocyanin electrostatic potentials. Second, the effects on the plastocyanin/cytochrome f association rate of these mutations in the plastocyanin "eastern site" were evaluated by simulating the association of the wild type and mutant plastocyanins with cytochrome f by Brownian dynamics. Good agreement between the computed and experimental relative rates (k(2)) (r(2) = 0.89 - 0.92) was achieved for the plastocyanin mutants. The results obtained by applying both computational techniques provide support for the fundamental role of the acidic residues at the plastocyanin eastern site in the association with cytochrome f and in the overall electron-transfer process.

Control of metalloprotein reduction potential: the role of electrostatic and solvation effects probed on plastocyanin mutants.

The changes in the thermodynamics of Cu(II) reduction for spinach plastocyanin induced by point mutations altering the electrostatic potential in proximity of the copper center were determined through variable temperature direct electrochemistry experiments. In particular, the functionally important surface residues Leu12 and Gln88 were replaced with charged and polar residues, and Asn38 was substituted with Asp. The mutational variations of the reduction enthalpy and entropy were analyzed with a QSPR (quantitative structure-property relationships) approach, employing global and local theoretical descriptors defined and computed on the three-dimensional protein structure. The correlations found are informative on how electrostatic and solvation effects control the E degrees ' values in this species through the combined effects on the reduction enthalpy and entropy. The changes in reduction enthalpy can be justified with electrostatic considerations. Most notably, enthalpy-entropy compensation phenomena play a significant role: the entropic effects due to the insertion of charged residues determine E degrees ' changes that are invariably opposite to those induced by the concomitant enthalpic effects. Therefore, the resulting E degrees ' changes are small or even opposite to those expected on simple electrostatic grounds. The mutational variation in the reduction entropy appears to be linked to the hydrogen bonding donor/acceptor character of the northern part of the protein, above the metal site, and to the electrostatic potential distribution around the copper site. Both properties influence the reduction-induced reorganization of the water molecules on the protein surface in the same region.

Mapping and fitting the peripheral benzodiazepine receptor binding site by carboxamide derivatives. Comparison of different approaches to quantitative ligand-receptor interaction modeling.

The synthetic-computational approach to the study of the binding site of peripheral benzodiazepine receptor (PBR) ligands related to 1-(2-chlorophenyl)-N-methyl-N-(1-methylpropyl)-3-isoquinolinecarboxamide (PK11195, 1) within their receptor (Cappelli et al. J. Med. Chem. 1997, 40, 2910-2921) has been extended. A series of carboxamide derivatives endowed with differently substituted planar aromatic or heteroaromatic systems was designed with the aim of getting further information on the topological requisites of the carbonyl and aromatic moieties for interaction with the PBR binding site. The synthesis of most of these compounds involves Weinreb amidation of the appropriate lactone as the key step. The most potent compound, among the newly synthesized ones, shows a nanomolar PBR affinity similar to that shown by 1 and the presence of a basic N-ethyl-N-benzylaminomethyl group in 3-position of the quinoline nucleus. Thus, it may be considered the first example of a new class of water soluble derivatives of 1. Several computational methods were used to furnish descriptors of the isolated ligands (indirect approaches) able to rationalize the variation in the binding affinity of the enlarged series of compounds. Sound QSAR models are obtained by size and shape descriptors (volume approach) which codify for the short-range contributions to ligand-receptor interactions. Molecular descriptors which explicitly account for the electrostatic contribution to the interaction (CoMFA, CoMSIA, and surface approaches) perform well, but they do not improve the quantitative models. Moreover, useful hints for the identification of the antagonist binding site in the three-dimensional modeling of the receptor (direct approach) were provided by the receptor hypothesis derived by the pharmacophoric approach. The ligand-receptor complexes obtained provided a detailed description of the modalities of the interaction and interesting suggestions for further experiments.

Blue copper proteins: a comparative analysis of their molecular interaction properties.

Blue copper proteins are type-I copper-containing redox proteins whose role is to shuttle electrons from an electron donor to an electron acceptor in bacteria and plants. A large amount of experimental data is available on blue copper proteins; however, their functional characterization is hindered by the complexity of redox processes in biological systems. We describe here the application of a semiquantitative method based on a comparative analysis of molecular interaction fields to gain insights into the recognition properties of blue copper proteins. Molecular electrostatic and hydrophobic potentials were computed and compared for a set of 33 experimentally-determined structures of proteins from seven blue copper subfamilies, and the results were quantified by means of similarity indices. The analysis provides a classification of the blue copper proteins and shows that (I) comparison of the molecular electrostatic potentials provides useful information complementary to that highlighted by sequence analysis; (2) similarities in recognition properties can be detected for proteins belonging to different subfamilies, such as amicyanins and pseudoazurins, that may be isofunctional proteins; (3) dissimilarities in interaction properties, consistent with experimentally different binding specificities, may be observed between proteins belonging to the same subfamily, such as cyanobacterial and eukaryotic plastocyanins; (4) proteins with low sequence identity, such as azurins and pseudoazurins, can have sufficient similarity to bind to similar electron donors and acceptors while having different binding specificity profiles.

Theoretical investigation of substrate specificity for cytochromes P450 IA2, P450 IID6 and P450 IIIA4.

Three-dimensional models of the cytochromes P450 IA2, P450 IID6 and P450 IIIA4 were built by means of comparative modeling using the X-ray crystallographic structures of P450 CAM, P450 BM-3, P450 TERP and P450 ERYF as templates. The three cytochromes were analyzed both in their intrinsic structural features and in their interaction properties with fifty specific and non-specific substrates. Substrate/enzyme complexes were obtained by means of both automated rigid and flexible body docking. The comparative analysis of the three cytochromes and the selected substrates, in their free and bound forms, allowed for the building of semi-quantitative models of substrate specificity based on both molecular and intermolecular interaction descriptors. The results of this study provide new insights into the molecular determinants of substrate specificity for the three different eukaryotic P450 isozymes and constitute a useful tool for predicting the specificity of new compounds.

Novel potent and selective central 5-HT3 receptor ligands provided with different intrinsic efficacy. 2. Molecular basis of the intrinsic efficacy of arylpiperazine derivatives at the central 5-HT3 receptors.

Novel 5-HT3 receptor ligands were designed and synthesized with the aim of obtaining deeper insight into the molecular basis of the intrinsic efficacy of arylpiperazines interacting with the central 5-HT3 receptor. The newly synthesized compounds and some previously published compounds belonging to the same class of heteroarylpiperazines were tested for their potential ability to displace [3H]granisetron from rat cortical membranes. These 5-HT3 receptor binding studies revealed subnanomolar affinity in several of the compounds under study. The most active ligands were quipazine derivatives bearing a phenyl group in the 4-position and various oxygenated alkyl side chains in the 3-position of the quinoline nucleus. Qualitative and theoretical quantitative structure-affinity relationship studies were carried out, and the interaction model for the 5-HT3 ligands related to quipazine with their receptor, proposed in part 1 of the present work, was updated to incorporate the latest data. The potential 5-HT3 agonist/antagonist activity of 12 selected compounds was assessed in vitro on the 5-HT3 receptor-dependent [14C]guanidinium uptake in NG 108-15 cells. Their intrinsic efficacy ranged from the 5-HT3 full agonist properties of compounds 7a and 8h, i to those of partial agonists 10a,d and antagonists 8b,d,e, and 9c, d,h,i. The comparison between these functional data and those relative to the previously described compounds suggested that in this class of 5-HT3 ligands the intrinsic efficacy is modulated in a rather subtle manner by the steric features of the heteroaryl moiety.

Synthesis, pharmacological evaluation, and structure-activity relationship and quantitative structure-activity relationship studies on novel derivatives of 2,4-diamino-6,7-dimethoxyquinazoline alpha1-adrenoceptor antagonists.

A new series of novel piperazine and non-piperazine derivatives of 2, 4-diamino-6,7-dimethoxyquinazoline was synthesized and evaluated for binding affinity toward alpha1-adrenergic and other G-protein-coupled aminergic receptors. The alpha1-adrenoceptor (AR) subtype selectivity was also investigated for the most interesting compounds. Only compound 16 showed moderate selectivity toward the alpha1b-AR subtype. Selected compounds were tested in vivo in a dog model indicating activity on blood pressure and on the lower urinary tract. Compound 10 showed in vivo potency close to that of prazosin. Powerful interpretative and predictive theoretical QSAR models have been obtained. The theoretical descriptors employed in the rationalization of the alpha1-adrenergic binding affinity depict the key features for receptor binding which can be summarized in an electrostatic interaction between the protonated amine function and a primary nucleophilic site of the receptor, complemented by short-range attractive (polar and dispersive) and repulsive (steric) intermolecular interactions. Moreover, on predictive grounds, the ad hoc derived size and shape QSAR model developed in a previous paper (Rastelli, G.; et al. J. Mol. Struct. 1991, 251, 307-318) proved to be successful in predicting nanomolar alpha1-adrenergic binding affinity for compound 28.

Relevance of theoretical molecular descriptors in quantitative structure-activity relationship analysis of alpha1-adrenergic receptor antagonists.

A quantitative structure-activity relationship (QSAR) study of a wide series of structurally diverse alpha1-adrenergic receptor antagonists was performed using the CODESSA (Comprehensive Descriptors for Structural and Statistical Analysis) technique. Theoretical descriptors derived on a single structure and ad hoc defined size and shape descriptors were considered in the attempt of describing information relevant to receptor interaction. The relative effectiveness of these two classes of parameters in developing QSAR models for native (alpha1A and alpha1B) and cloned (alpha1a, alpha1b, and alpha1d) adrenergic receptor binding affinity, functional activity of vascular and lower urinary tract tissues, and in vitro and in vivo selectivity was evaluated.

Theoretical descriptors in quantitative structure-affinity and selectivity relationship study of potent N4-substituted arylpiperazine 5-HT1A receptor antagonists.

The ability of ad hoc defined size and shape descriptors and theoretical descriptors derived on a single structure to give powerful interpretative and predictive QSAR models has been compared and evaluated with respect to the quality of the pharmacological data available for a series of structurally diverse 5-HT1A receptor antagonists, displaying selectivity towards the alpha 1-adrenergic receptor.

Computer modeling of size and shape descriptors of alpha 1-adrenergic receptor antagonists and quantitative structure-affinity/selectivity relationships.

Computational chemistry and molecular modeling procedures allow us to define and compute ad hoc size and shape descriptors on the different prototropic forms assumed by drugs in biotest solutions. Together with experimental data measured on a well-identified target receptor, these descriptors are essential elements for obtaining simple, consistent, comparable, and easily interpretable theoretical quantitative structure-activity relationship (QSAR) models based on the ligand similarity-target receptor complementarity paradigm. In this context, quantitative size and shape affinity/subtype selectivity relationships have been modeled for a large set of very heterogeneous alpha 1a-, alpha 1b-, and alpha 1c-adrenergic receptor antagonists. The linear QSAR models generated have been validated by predicting both binding affinity and selectivity of a test set of noncongeneric antagonists. The satisfactory results obtained highlight both the simplicity and the versatility of the approach presented.

Novel potent and selective central 5-HT3 receptor ligands provided with different intrinsic efficacy. 1. Mapping the central 5-HT3 receptor binding site by arylpiperazine derivatives.

Synthesis and pharmacological evaluation of a series of condensed quinoline and pyridine derivatives bearing a N-methylpiperazine moiety attached to the 2-position of the quinoline or pyridine nucleus are described. 5-HT receptor binding studies revealed subnanomolar affinity for the 5-HT3 receptor subtype in some of the compounds under study. The most active compound (5b) displayed a Ki value about 1 order of magnitude higher than that of quipazine along with a higher selectivity. The potential 5-HT3 agonist/antagonist activity of four selected compounds was assessed in vitro on 5-HT3 receptor-dependent [14C]guanidinium uptake in NG 108-15 cells. Compound 5j acted as a 5-HT3 agonist in this assay with an EC50 value close to that reported for quipazine, while 5b was a partial agonist with an EC50 value of about 0.25 nM, and compound 5c possessed antagonist properties with an IC50 value (approximately 8 nM) in the same range as those of previously characterized 5-HT3 receptor antagonists. Qualitative and quantitative structure-affinity relationship studies carried out by making use of theoretical molecular descriptors allowed to elucidate the role of the main pharmacophoric components and to develop a model for the interaction of the 5-HT3 ligands related to quipazine with their receptor.

Theoretical investigation of IL-6 multiprotein receptor assembly.

Interleukin-6 (IL-6) is a multifunctional cytokine that regulates cell growth, differentiation, and cellular functions in many cell lineages. Recently, evidences for the formation of an active hexameric complex with an IL-6:IL-6R alpha:gp130 stoichiometry of 2:2:2 have been obtained by different experimental approaches. Analysis of the electrostatic potential complementarity between IL-6 and its receptors has been used, in this study, to guide the assembly of homology-based 3D models of the components. The results strongly support a mechanism whereby the active cytokine (IL-6: IL-6R alpha) associates with the signal transducing gp130 protein, and the trimeric complex formed further dimerizes to form the hexameric species. Furthermore, computational simulations of the multiprotein complexes provide a rationalization of data from mutation experiments and highlight some key protein-protein interactions which have not yet been the subject of mutagenesis studies.

Mapping the peripheral benzodiazepine receptor binding site by conformationally restrained derivatives of 1-(2-chlorophenyl)-N-methyl-N-(1-methylpropyl)-3- isoquinolinecarboxamide (PK11195).

A synthetic-computational approach to the study of the binding site of peripheral benzodiazepine receptor (PBR) ligands related to 1-(2-chlorophenyl)-N-methyl-N-(1-methylpropyl)-3-isoquinolinecarboxam ide (PK11195, 1) within their receptor has been developed. A wide series of conformationally restrained derivatives of 1 has been designed with the aim of probing the PBR binding site systematically. The synthesis of these compounds involves palladium-catalyzed coupling and amidation as the key steps. Twenty-nine rigid and semirigid derivatives of 1 were tested in binding studies using [3H]-1, and most of these showed PBR affinities in the nanomolar range. The essential role of the carbonyl moiety as a primary pharmacophoric element in the recognition by and the binding to PBR has been confirmed, and the restricted range of the carbonyl orientations, which characterizes the most potent ligands, points to a specific hydrogen-bonding interaction, mainly directed by the geometrical factors, when the electronic ones are fulfilled. Moreover, the fundamental importance of the short-range dispersive interactions in the modulation of the binding affinity and, hence, in the stabilization of the ligand-receptor complex, emerged from the QSAR models reported.

Alpha 1-adrenoceptor subtype selectivity: molecular modelling and theoretical quantitative structure-affinity relationships.

This study constitutes a preliminary rationalization, at the molecular level, of antagonist selectivity towards the three cloned alpha 1-adrenergic receptor (alpha 1-AR) subtypes. Molecular dynamics simulations allowed a structural/dynamics analysis of the seven alpha-helix-bundle models of the bovine alpha 1a-, hamster alpha 1b-, and rat alpha 1d-AR subtypes. The results showed that the transmembrane domains of these subtypes have different dynamic behaviours and different topographies of the binding sites, which are mainly constituted by conserved residues. In particular, the alpha 1a-AR binding site is more flexible and topographically different with respect to the other two subtypes. The results of the theoretical structural/dynamics analysis of the isolated receptors are consistent with the binding affinities of the 16 antagonists tested towards the three cloned alpha 1-AR subtypes. Moreover, the theoretical quantitative structure-affinity relationships obtained from the antagonist-receptor interaction models further corroborates the hypothesis that selectivity towards one preferential subtype is mainly modulated by receptor and/or ligand distortion energies. In other words, subtype selectivity seems to be mainly guided by the dynamic complementarity (induced fit) between ligand and receptor. On the basis of the quantitative models presented it is possible to predict both affinities and selectivities of putative alpha 1-AR ligands as well as to estimate the theoretical alpha 1-AR subtype affinities and selectivities of existing antagonists.