Theoretical evidence of a salt bridge disruption as the initiating process for the alpha1d-adrenergic receptor activation: a molecular dynamics and docking study.

This study reports the building of the three-dimensional structure of the rat alpha1d-adrenergic receptor through a topology approach based on the structure of the rhodopsin receptor from cryoelectron microscopy. The validity and reliability of the receptor model were assessed through exhaustive molecular dynamics and docking studies. Some interesting ligand-receptor interactions were identified along with significant differences between the binding mode of agonists and antagonists. The importance of the disruption of a salt bridge as a possible initial event leading to receptor activation is discussed on the basis of data from mutagenesis and molecular dynamics studies.

Neuronal nicotinic receptor agonists: a multi-approach development of the pharmacophore.

Based on the results obtained with different automated computational approaches as applied to the study of eleven high-affinity agonists of the neuronal nicotine acetylcholine receptor (nAChR), belonging to different chemical classes, new relevant features were detected which complement the existing pharmacophores. Convergent results from DISCO (Distance Comparison), QXP (Quick Explore), Catalyst/HipHop, and MIPSIM (Molecular Interaction Potential Similarity) allowed us to identify and locate, in a well defined spatial arrangement, three geometrically independent key structural features: (i) a positively charged nitrogen atom for ionic or hydrogen bond interactions, (ii) a lone pair of the pyridine nitrogen or a specific lone pair of a carbonyl oxygen, as a hydrogen bond acceptor, and (iii) a centre of a hydrophobic area generally occupied by aliphatic cycles. The pharmacophore presented herein, along with predictive 2D and 3D QSAR models recently developed in our group, could represent valuable computational tools for the design of new nAChR agonists having therapeutical potential.

3D-QSAR methods on the basis of ligand-receptor complexes. Application of COMBINE and GRID/GOLPE methodologies to a series of CYP1A2 ligands.

Many heterocyclic amines (HCA) present in cooked food exert a genotoxic activity when they are metabolised (N-oxidated) by the human cytochrome P450 1A2 (CYP1A2h). In order to rationalize the observed differences in activity of this enzyme on a series of 12 HCA, 3D-QSAR methods were applied on the basis of models of HCA-CYP1A2h complexes. The CYP1A2h enzyme model has been previously reported and was built by homology modeling based on cytochrome P450 BM3. The complexes were automatically generated applying the AUTODOCK software and refined using AMBER. A COMBINE analysis on the complexes identified the most important enzyme-ligand interactions that account for the differences in activity within the series. A GRID/GOLPE analysis was then performed on just the ligands, in the conformations and orientations found in the modeled complexes. The results from both methods were concordant and confirmed the advantages of incorporating structural information from series of ligand-receptor complexes into 3D-QSAR methodologies.

A proposed bioactive conformation of peptide T.

The conformational profiles of Peptide T, (5-8)Peptide T, [Abu5](4-8)Peptide T and (4-8)Peptide T were computed independently to assess the geometrical characteristics of the bioactive conformation of Peptide T. The conformational profiles of the peptides were computed within the molecular mechanics framework using an effective dielectric constant of 80. The conformational space was thoroughly sampled using an iterative simulated annealing protocol. The bioactive conformation was assessed by pairwise cross comparisons of each of the unique low energy conformations found for each of the different analogs studied. After a putative bioactive conformation was selected, in order to further validate our hypothesis the conformational profile of the potent compound cyclo(Thr-Thr-Asn-Tyr-Thr-Asp) was computed and the putative bioactive conformation was found. The conformation exhibits a pseudo beta-turn involving the side chain of Thr5 and the carbonyl oxygen of Tyr7 forming a C12 ring.

Conformational analysis of the highly potent bradykinin antagonist Hoe-140 by means of two different computational methods.

The AMBER 4.0 force field was used to perform the characterization of the conformational profile of the highly potent bradykinin antagonist Hoe-140 (D-Arg0-Arg1-Pro2-Hyp3-Gly4-Thi5-Ser6-D-++ +Tic7-Oic8-Arg9). The structural features of the peptide were assessed using two different computational methods, both capable to provide a good sampling of the low-energy conformations of the molecule. Specifically, the conformational space of the peptide was explored: i) computing molecular dynamics trajectories in cycles of high (900 K) and low (300 K) temperature and ii) using simulated annealing (SA) in an iterative fashion. Analysis of the structures characterized indicates that most of the low-energy conformations of the peptide exhibit a betaII'-turn motif at its C-terminus, in agreement with previous experimental and theoretical studies. On the other hand, about a 50% of the low-energy conformations characterized also exhibit different beta-turn type motifs at the N-terminus, whereas the rest of the conformations can be described as bends. Finally, in order to get new insights into the structural requirements necessary to design more potent and selective antagonists of bradykinin, present results were compared with those previously reported by this laboratory on the conformational preferences of the native nonapeptide and its DPhe7 analog.

Three-dimensional modelling of human cytochrome P450 1A2 and its interaction with caffeine and MeIQ.

The three-dimensional modelling of proteins is a useful tool to fill the gap between the number of sequenced proteins and the number of experimentally known 3D structures. However, when the degree of homology between the protein and the available 3D templates is low, model building becomes a difficult task and the reliability of the results depends critically on the correctness of the sequence alignment. For this reason, we have undertaken the modelling of human cytochrome P450 1A2 starting by a careful analysis of several sequence alignment strategies (multiple sequence alignments and the TOPITS threading technique). The best results were obtained using TOPITS followed by a manual refinement to avoid unlikely gaps. Because TOPITS uses secondary structure predictions, several methods that are available for this purpose (Levin, Gibrat, DPM, NnPredict, PHD, SOPM and NNSP) have also been evaluated on cytochromes P450 with known 3D structures. More reliable predictions on alpha-helices have been obtained with PHD, which is the method implemented in TOPITS. Thus, a 3D model for human cytochrome P450 1A2 has been built using the known crystal coordinates of P450 BM3 as the template. The model was refined using molecular mechanics computations. The model obtained shows a consistent location of the substrate recognition segments previously postulated for the CYP2 family members. The interaction of caffeine and a carcinogenic aromatic amine (MeIQ), which are characteristic P450 1A2 substrates, has been investigated. The substrates were solvated taking into account their molecular electrostatic potential distributions. The docking of the solvated substrates in the active site of the model was explored with the AUTODOCK programme, followed by molecular mechanics optimisation of the most interesting complexes. Stable complexes were obtained that could explain the oxidation of the considered substrates by cytochrome P450 1A2 and could offer an insight into the role played by water molecules.

Computational study of the conformational domains of peptide T.

The conformational preferences of peptide T (ASTTTNYT) were analysed by means of computational methods. A thorough exploration of the conformational space was carried out within the framework of the molecular mechanics approach, using simulated annealing as a searching strategy. Specifically, in order to obtain a subset of low-energy conformations with energies close to the global minimum as complete as possible, a simulated annealing protocol was repeated several times in a recursive fashion. The results of the search indicate that the peptide exhibits a alpha-helical character although most of the conformations characterized, including the global minimum, can be described as bent conformations. Conformations exhibiting beta-turn motives previously proposed from NMR studies were also characterized, although they are not very predominant in the set of low-energy conformations.

On the bioactive conformation of a small peptide and its set of thermodynamically accessible conformations.

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics frame-work using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.

Docking of a human rhinovirus neutralizing antibody onto the viral capsid.

The structure of the complex between the Fab fragment of a human rhinovirus serotype 2 (HRV2) neutralizing antibody (8F5) and a cross-reactive synthetic peptide derived from the viral capsid protein VP2 has been recently determined by crystallographic methods. The conformation adopted by the peptide was very similar to and could be superimposed onto the corresponding region of the viral protein VP2 of human rhinovirus 1A (HRV1A) whose three-dimensional structure is known. The structure of the Fab fragment determined in the complex was docked onto the viral capsid using the superimposition transformation found for the peptide. In the resulting model the Fab protrudes almost radially to about 60 A from the surface of the virion without any major steric problem. The Fab fragment was then placed on each one of the 60 equivalent epitopes using the T = 1 icosahedral symmetry of the virus. The closest pairs of Fab fragments are related by viral 2-fold axes and run almost parallel to each other without clashing. These axes of symmetry from the viral particle could thus be coincident with the dyad axes of the antibodies. Furthermore, comparison of the three-dimensional structure of the Fab/peptide complex with the structure of the Fab fragment alone indicates that the flexibility of the antibody's elbow would facilitate bivalent attachment to the same viral particle. In accordance with the docking results, experimental determination of the stoichiometry of binding yielded a ratio of 30 IgG molecules per virion also suggesting bivalent attachment of antibody 8F5 onto the viral particle. The neutralization of viral infectivity, being neither aggregation (this paper) nor inhibition of receptor binding, might be mainly achieved by reducing viral spread from cell to cell and/or inhibition of uncoating.

Characterization of the conformational domains of bradykinin by computational methods.

The AMBER 4.0 force field was used to perform a characterization of the conformational profile of the nonapeptide bradykinin. A thorough conformational search was carried out using molecular dynamics as sampling technique, by computing cycles of high (900 K) and low (300 K) temperature trajectories. A total of 2400 minima were generated and subsequently clustered using the root-mean-square of the backbone dihedral angles as criterium. After the use of a tolerance value of 20 degrees, the conformations were clustered in 233 unique conformations with energies up to 40 kcal/mol above the lowest minimum. The analysis of the low-energy conformations indicate that the peptide exhibits a high tendency to adopt a beta-turn at the C-terminus and a propensity to adopt a bent structure at the N-terminus. These results are in agreement with the experimental evidence reported in the literature and provide detailed information necessary to understand the conformational preferences of the peptide.