Biomedical subjects

Ricardo Bicca de Alencastro

Publications and source records attributed to Ricardo Bicca de Alencastro.

11 recordsLinked to original sources

PubMed2006

Free-energy force-field three-dimensional quantitative structure-activity relationship analysis of a set of p38-mitogen activated protein kinase inhibitors.

The p38-mitogen-activated protein kinases (p38-MAPKs) belong to a family of serine-threonine kinases activated by pro-inflammatory or stressful stimuli that are known to be involved in several diseases. Their biological importance, related to the release of inflammatory pro-cytokines such as tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1), has generated many studies aiming at the development of selective inhibitors for the treatment of inflammatory diseases. In this work, we developed receptor-based three dimensional (3D) quantitative structure-activity relationship (QSAR) models for a series of 33 pyridinyl imidazole compounds [Liverton et al. (1999) 42:2180], using a methodology named free-energy force-field (FEFF) [Tokarski and Hopfinger (1997) 37:792], in which scaled intra- and intermolecular energy terms of the Assisted Model Building Energy Refinement (AMBER) force field combined with a hydration-shell solvation model are the independent variables used in the QSAR studies. Multiple temperature molecular-dynamics simulations (MDS) of ligand-protein complexes and genetic-function approximation (GFA) were employed using partial least squares (PLS) as the fitting functions to develop FEFF-3D-QSAR models for the binding process. The best model obtained in the FEFF-3D-QSAR receptor-dependent (RD) method shows the importance of the van der Waals energy change upon binding and the electrostatic energy in the interaction of ligands with the receptor. The QSAR equations described here show good predictability and may be regarded as representatives of the binding process of ligands to p38-MAPK. Additionally, we have compared the top FEFF-3D-QSAR model with receptor independent (RI) 4D-QSAR models developed in a recent study [Romeiro et al. (2005) 19:385].

Hydrophobic and Hydrophilic Interactions↗

PubMed2005

Construction of 4D-QSAR models for use in the design of novel p38-MAPK inhibitors.

The p38-mitogen-activated protein kinase (p38-MAPK) plays a key role in lipopolysaccharide-induced tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1) release during the inflammatory process, emerging as an attractive target for new anti-inflammatory agents. Four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis [Hopfinger et al., J. Am. Chem. Soc., 119 (1997) 10509] was applied to a series of 33 (a training set of 28 and a test set of 5) pyridinyl-imidazole and pyrimidinyl-imidazole inhibitors of p38-MAPK, with IC50 ranging from 0.11 to 2100 nM [Liverton et al., J. Med. Chem., 42 (1999) 2180]. Five thousand conformations of each analogue were sampled from a molecular dynamics simulation (MDS) during 50 ps at a constant temperature of 303 K. Each conformation was placed in a 2 angstroms grid cell lattice for each of three trial alignments. 4D-QSAR models were constructed by genetic algorithm (GA) optimization and partial least squares (PLS) fitting, and evaluated by leave-one-out cross-validation technique. In the best models, with three to six terms, the adjusted cross-validated squared correlation coefficients, Q2adj, ranged from 0.67 to 0.85. Model D (Q2adj = 0.84) was identified as the most robust model from alignment 1, and it is representative of the other best models. This model encompasses new molecular regions as containing pharmacophore sites, such as the amino-benzyl moiety of pyrimidine analogs and the N1-substituent in the imidazole ring. These regions of the ligands should be further explored to identify better anti-inflammatory inhibitors of p38-MAPK.

Anti-Inflammatory Agents, Non-Steroidal↗

PubMed2004

Theoretical evaluation of adiabatic and vertical electron affinity of some radiosensitizers in solution using FEP, ab initio and DFT methods.

The biological activity of radiosensitizers is associated to their electron affinity (EA), which can be divided in two main processes: vertical and adiabatic. In this work, we calculated the EAs of nitrofurans and nitroimidazoles (Fig. 2) using Hartree-Fock (HF) and density functional theory (DFT) methods and evaluated solvent effects (water and carbon tetrachloride) on EAs. For water, we combined the polarized continuum model (PCM) and free energy perturbation (FEP) (finite difference thermodynamic integration, FDTI) methods. For carbon tetrachloride, we used the FDTI method. The values of adiabatic EA obtained are in agreement with experimental data (deviations of 0.013 eV). The vertical EAs were calculated according to Cederbaum's outer valence Green function (OVGF) method. This methodology, which relies on theoretical aspects of free energy calculations on charged molecules in solution, was used to select potential selective radiosensitizers from recently reported compounds and could be helpful in the rational design of new and more selective bioreductive anticancer drugs.

Carbon Tetrachloride↗

PubMed2004

LIV-3D-QSAR models for PGI2 receptor ligands using multiple conformations.

A new 3D descriptor, the local intersection volume (LIV), was developed by our group and applied to the construction of 3D-QSAR models for ligands of the PGI(2) receptor (IP). The target compounds are a set of 42 aromatic heterocyclic derivatives [Meanwell et al., J. Med. Chem. 36 (1993), 3884], which show agonist activities in the IP receptor and are inhibitors of platelet aggregation. The LIV-3D-QSAR models were obtained through the analysis of 30% of the generated conformations for each compound, using a combined Genetic Algorithm (GA) and Partial Least Square (PLS) approach [Rogers and Hopfinger, J. Inf. Comput. Sci. 34 (1994) 854]. Statistically, Model 3 is the best as well as the most comprehensive in a mechanistic sense. Furthermore, it can be applied to design new IP ligands.

Algorithms↗

PubMed2004

Interactions of 5-deazapteridine derivatives with Mycobacterium tuberculosis and with human dihydrofolate reductases.

There are major differences between the structures of human dihydrofolate reductase (hDHFR) and Mycobacterium tuberculosis dihydrofolate reductase (mtDHFR). These differences may allow us to design more selective mtDHFR inhibitors. In this paper we study the reactions of six different compounds derived from 5-deazapteridine with human and bacterial enzymes. Results suggest that the addition of hydrophobic groups to the aminophenyl ring would increase mtDHFR-inhibitor affinity and selectivity.

Antitubercular Agents↗

PubMed2003

Synthetic and theoretical studies on the reduction of electron withdrawing group conjugated olefins using the Hantzsch 1,4-dihydropyridine ester.

The Hantzsch 1,4-dihydropyridine ester (1) has been observed to be a useful selective reducing agent for the reduction of electron-withdrawing conjugated double bonds. The rate of this reaction was observed to be dependent upon the nature of the conjugated substituents and, consequently, the electronic nature of the unsaturated double bond. Theoretical calculations confirmed the importance of the HOMO-LUMO gap for this reaction and implicated a hydride transfer, agreeing with the experimentally observed reaction rate order. The calculations also revealed the importance of a boatlike structure of the 1,4-dihydropyridine nucleus as well as a trans arrangement of the ester groups to facilitate the hydride transfer.

Journal Article↗

PubMed2003

Investigation of the induced-fit mechanism and catalytic activity of the human cytomegalovirus protease homodimer via molecular dynamics simulations.

Human cytomegalovirus (HCMV) is a highly species-specific DNA virus infecting up to 80% of the general population. The viral genome contains the open reading frame UL80, which encodes the full-length 80 kDa HCMV serine protease and its substrate. Full-length HCMV protease is composed of an N-terminal 256-amino-acid proteolytic domain, called assemblin, a linker region, and a C-terminal structural domain, the assembly protein precursor. Biochemical studies have shown that dimerization activates assemblin because of an induced stabilization of the oxyanion hole (Arg166). Thus, we performed here molecular dynamics (MD) simulations on HCMV protease models to study the induced-fit mechanism of the enzyme upon the binding of substrates and peptidyl inhibitors, and structural and energetic factors that are responsible for the catalytic activity of the enzyme dimer. Long and stable trajectories were obtained for the models of the monomeric and dimeric states, free in solution and bound to a peptidyl-activated carbonyl inhibitor, with very good agreement between theoretical and experimental results. Our results suggest that HCMV protease is indeed a novel example of serine protease that operates by an induced-fit mechanism. Also, in agreement with mutagenesis studies, our MD simulations suggest that the dimeric form is necessary to activate the enzyme because of an induced stabilization of the oxyanion hole.

Arginine↗

PubMed2003

Potential of mean force calculations on an L-type calcium channel model.

To understand the mechanisms of Na(+)/Li(+) permeation at submicromolar Ca(2+) concentrations, Na(+)/Li(+) blocking at higher Ca(2+) concentrations (10(-6)-10(-4) M) and Ca(2+) permeation at millimolar Ca(2+) concentrations, we used our recently described L-type calcium channel model. For this purpose, we obtained potential of mean force (pmf) curves for the position change of one Na(+) and one Ca(2+) ion inside the channel and for the position change of a second Ca(2+) ion when the EEEE locus is coordinated to Ca(2+). The pmf curves suggest that (i) at submicromolar Ca(2+) concentrations, because of the low velocity of Ca(2+) entry in the channel, monovalent ion flux occurs; (ii) at Ca(2+) concentrations between 10(-6) and 10(-4) M, thermodynamic equilibrium between the channel and Ca(2+) is achieved; as the coordination of Ca(2+) with the locus is more favorable than the coordination of Na(+), the monovalent ion flux is blocked; and (iii) to put a second Ca(2+) ion inside the channel at an appropriate rate, the Ca(2+) concentration should reach millimolar levels. Nevertheless, the entry of a second Ca(2+) is thermodynamically unfavorable, indicating that the competition of two Ca(2+) ions for the locus leads to Ca(2+) permeation.

Calcium↗

PubMed2002

Thrombin inhibition by novel benzamidine derivatives: a free-energy perturbation study.

Thrombin is a serine protease responsible for blood coagulation. Since thrombin inhibitors appear to be effective in the treatment and prevention of thrombotic and embolic disorders, considerable attention has been focused on the structure and interactions of this enzyme. In this work, to evaluate the relative free energies of hydration and binding to thrombin for some benzamidine derivatives, we used the finite difference thermodynamic integration (FDTI) algorithm within the Discover program of MSI. By this method, two possible orders of hydration for the candidates were obtained: p-amidinophenylpyruvate > p-(2-oxo-1-propyl)benzamidine > p-methylbenzamidine > p-ethylbenzamidine > p-(1-propyl)benzamidine > benzamidine and p-amidinophenylpyruvate > p-(2-oxo-1-propyl)benzamidine > p-methylbenzamidine > p-ethylbenzamidine > benzamidine > p-(1-propyl)benzamidine. We also obtained the following order for thrombin binding: p-(2-oxo-1-propyl)benzamidine > p-ethylbenzamidine > p-(1-propyl)benzamidine > p-methylbenzamidine > benzamidine > p-amidinophenylpyruvate.

Benzamidines↗

PubMed2002

Local intersection volume: a new 3D descriptor applied to develop a 3D-QSAR pharmacophore model for benzodiazepine receptor ligands.

In this work, we have developed a new descriptor, named local intersection volume (LIV), in order to compose a 3D-QSAR pharmacophore model for benzodiazepine receptor ligands. The LIV can be classified as a 3D local shape descriptor in contraposition to the global shape descriptors. We have selected from the literature 49 non-benzodiazepine compounds as a training data set and the model was obtained and evaluated by genetic algorithms (GA) and partial least-squares (PLS) methods using LIVs as descriptors. The LIV 3D-QSAR model has a good predictive capacity according the cross-validation test by "leave-one-out" procedure (Q(2)=0.72). The developed model was compared to a comprehensive and extensive SAR pharmacophore model, recently proposed by Cook and co-workers, for benzodiazepine receptor ligands [J. Med. Chem. 43 (2000) 71]. It showed a relevant correlation with the pharmacophore groups pointed out in that work. Our LIV 3D-QSAR model was also able to predict affinity values for a series of nine compounds (test data set) that was not included into the training data set.

Algorithms↗

PubMed2002

A molecular dynamics study of an L-type calcium channel model.

In this work, we propose a molecular model of the L-type calcium channel pore from the human cardiac alpha1 subunit. Four glutamic acid residues, the EEEE locus, located at highly conserved P loops (also called SS1-SS2 segments) of the alpha1 subunit, molecularly express the calcium channel selectivity. The proposed alpha-helix structure for the SS1 segment, analyzed through molecular dynamics simulations in aqueous-phase, was validated by the plotting of Ramachandran diagrams for the averaged structures and by the analysis of i and i + 4 helical hydrogen bonding between the amino acid residues. The results of the simulation of the calcium channel model with one and two Ca2+ ions at the binding site are in accordance with mutation studies which suggest that the EEEE locus in the L-type calcium channel must form a single high-affinity binding site. These results suggest that the Ca2+ permeation through the channel would be derived from competition between two ions for the only high-affinity binding site. Furthermore, the experimentally observed blocking of the Na+ flux at micromolar Ca2+ concentrations, probably due to the occupancy of the single high-affinity binding site for one Ca2+, was also reproduced by our model.

Amino Acid Sequence↗