Search

Angiotensin I-converting enzyme (ACE), a key enzyme in cardiovascular pathophysiology, consists of two homologous domains (N and C), each bearing a Zn-dependent active site. We modeled the 3D-structure of the ACE N-domain using known structures of the C-domain of human ACE and the ACE homologue, ACE2, as templates. Two monoclonal antibodies (mAb), 3A5 and i2H5, developed against the human N-domain of ACE, demonstrated anticatalytic activity. N-domain modeling and mutagenesis of 21 amino acid residues allowed us to define the epitopes for these mAbs. Their epitopes partially overlap: amino acid residues K407, E403, Y521, E522, G523, P524, D529 are present in both epitopes. Mutation of 4 amino acid residues within the 3A5 epitope, N203E, R550A, D558L, and K557Q, increased the apparent binding of mAb 3A5 with the mutated N-domain 3-fold in plate precipitation assay, but abolished the inhibitory potency of this mAb. Moreover, mutation D558L dramatically decreased 3A5-induced ACE shedding from the surface of CHO cells expressing human somatic ACE. The inhibition of N-domain activity by mAbs 3A5 and i2H5 obeys similar kinetics. Both mAbs can bind to the free enzyme and enzyme-substrate complex, forming E.mAb and E.S.mAb complexes, respectively; however, only complex E.S can form a product. Kinetic analysis indicates that both mAbs bind better with the ACE N-domain in the presence of a substrate, which, in turn, implies that binding of a substrate causes conformational adjustments in the N-domain structure. Independent experiments with ELISA demonstrated better binding of mAbs 3A5 and i2H5 in the presence of the inhibitor lisinopril as well. This effect can be attributed to better binding of both mAbs with the "closed" conformation of ACE, therefore, disturbing the hinge-bending movement of the enzyme, which is necessary for catalysis.

A molecular model of the human A(2B) adenosine receptor containing seven transmembrane alpha helices connected by three intracellular and three extracellular hydrophilic loops had been constructed. A molecular docking of seven structurally diverse xanthine antagonists of the A(2B) receptor was performed, and the differences in their binding modes were investigated. The 1 ns molecular dynamics (MD) simulations of several obtained ligand-receptor complexes inserted into the phospholipid bilayer were carried out. The conformational changes of the A(2B) receptor occurring during MD simulations were explored, and the stable binding modes of the studied antagonists were determined. According to the models presented in this work, the involvement of the His251, Asn282, Ser92, Thr89, and some aromatic residues in ligand recognition was determined. The obtained binding modes of the A(2B) antagonists demonstrate good agreement with the site-directed mutagenesis data.

Modeling the structure of the C-domain of bovine angiotensin-converting enzyme revealed two putative chloride-binding sites. The kinetic parameters, K(m) and k(cat), of hydrolysis of the substrate Cbz-Phe-His-Leu catalyzed by the testicular (C-domain) enzyme were determined over a wide range of chloride concentrations. Chloride anions were found to be enzyme activators at relatively low concentrations, but they inhibit enzymatic activity at high concentrations. A general scheme for the effect of chloride anions on activity of the C-domain of bovine angiotensin-converting enzyme accounting for binding the "activating" and "inhibiting" anions is suggested.

An approach for evaluation of binding selectivity was suggested and exemplified using glycine/NMDA and AMPA receptors. For analyzing the pairwise selectivity, we propose to use the difference between biological activities (expressed as -log Ki) of ligands with respect to different receptor subtypes as a dependent variable for building comparative molecular field analysis (CoMFA) models. The resulting fields (which will be referred to as the "selectivity fields") indicate the ways of increasing selectivity of binding, inhibition, etc. As an example, CoMFA of a set of pyrazolo[1,5-c]quinazolines and triazolo[1,5-c]quinazolines was used for considering the binding selectivity with respect to glycine/NMDA and AMPA receptors. In addition, the mapping of these fields onto the molecular models of the corresponding receptors makes it possible to reveal the reasons for experimentally observed selectivity as well as to suggest additional ways of increasing selectivity.

Somatic angiotensin-converting enzyme (ACE) consists of two homologous domains, each of them containing an active site. Differences in substrate specificities and affinity to inhibitors of the active sites of the two domains of bovine ACE are described. The ACE domains demonstrate different thermostability, and the reasons for this difference are analyzed. A structural model of the ACE domains is suggested, which allows us to reveal the structural subdomain important for the protein stability and localize the hydrophobic and the carbohydrate-binding sites.

Styryl dyes 4a-e containing a 15-crown-5 ether unit and a quinoline residue with a sulfonatoalkyl or sulfonatobenzyl N-substituent were synthesized. The relationship between the photochemical behavior of these dyes and their aggregates derived from complexation with Mg(2+) in MeCN was studied using (1)H NMR and absorption spectroscopy. The E-isomers of 4a-e were shown to form highly stable dimeric (2:2) complexes with Mg(2+). Upon irradiation with visible light, the dimeric complexes undergo two competing photoreactions, viz., geometric E --> Z isomerization, resulting in an anion-capped 1:1 complex of the Z-isomer with Mg(2+) and stereospecific syn-head-to-tail [2+2]-cycloaddition, affording a single isomer of bis-crown-containing cyclobutane. The N-substituent in the dye has a dramatic effect on the photochemical behavior of the dimeric complex. Molecular dynamics and semiempirical quantum-chemical calculations were carried out to interpret the observed photocycloaddition in the dimer. Conformational equilibria for the dimer of (E)-4b were analyzed using (1)H NMR spectroscopy.

Homology modeling was used to build 3D models of the N-methyl-d-aspartate (NMDA) receptor glycine binding site on the basis of an X-ray structure of the water-soluble AMPA-sensitive receptor. The docking of agonists and antagonists to these models was used to reveal binding modes of ligands and to explain known structure-activity relationships. Two types of quantitative models, 3D-QSAR/CoMFA and a regression model based on docking energies, were built for antagonists (derivatives of 4-hydroxy-2-quinolone, quinoxaline-2,3-dione, and related compounds). The CoMFA steric and electrostatic maps were superimposed on the homology-based model, and a close correspondence was marked. The derived computational models have permitted the evaluation of the structural features crucial for high glycine binding site affinity and are important for the design of new ligands.

We present new homology-based models of the glutamate binding site (in closed and open forms) of the NMDA receptor NR2B subunit derived from X-ray structures of the water soluble AMPA sensitive glutamate receptor. The models were used for revealing binding modes of agonists and competitive antagonists, as well as for rationalizing known experimental facts concerning structure-activity relationships: (i) the switching between the agonist and the antagonist modes of action upon lengthening the chain between the distal acidic group and the amino acid moiety, (ii) the preference for the methyl group attached to the alpha-amino group of ligands, (iii) the preference for the D-configuration of agonists and antagonists, and (iv) the existence of "superacidic" agonists.

An approach to the interpretation of backpropagation neural network models for quantitative structure-activity and structure-property relationships (QSAR/QSPR) studies is proposed. The method is based on analyzing the first and second moments of distribution of the values of the first and the second partial derivatives of neural network outputs with respect to inputs calculated at data points. The use of such statistics makes it possible not only to obtain actually the same characteristics as for the case of traditional "interpretable" statistical methods, such as the linear regression analysis, but also to reveal important additional information regarding the non-linear character of QSAR/QSPR relationships. The approach is illustrated by an example of interpreting a backpropagation neural network model for predicting position of the long-wave absorption band of cyane dyes.

The relationship between mutagenic activity and chemical structure was studied for 54 polycyclic compounds using two approaches: multiple linear regression analysis and artificial neural networks. Structural fragments, quantum chemical indices, and hydrophobicity (octanol-water partition coefficient) were used as descriptors (properties of the molecules introduced in the model). Both linear regression equations and nonlinear relationships obtained with the help of a neural network were shown to accurately predict mutagenic activity for the compounds structurally similar to those in the training sample. The introduction of experimentally selected descriptors is substantiated to verify the proposed mechanism of related compounds mutagenic activity.

A method to build QSAR models based on substituent constants for congeneric sets of compounds having several topologically equivalent substituent positions was proposed. The approach is based on the application of artificial neural networks (learning to construct nonlinear structure-activity relationships taking into account necessary symmetry properties of training set structures) to a training set expanded by adding the copies of compounds with the same activity values but with permuted assignment of equivalent substituent positions. The better predictive power of these constructed models, as compared with the performances of neural network models for non-expanded sets was demonstrated for the calcium channel blockers of 1,4-dihydropyridine type and for hallucinogenic phenylalkylamines.

It was shown that the full-size neurotrophic factor from pigment epithelium (PEDF) induces the cell differentiation of the human promyelocyte leukemia cell line HL-60. A structural analysis of PEDF revealed in its C-terminal region a six-membered peptide fragment PEDF-(352-357) (PEDF-6) whose sequence is highly homologous to the 41-46 fragment of the active site of the human leukocyte differentiation factor HLDF (HLDF-6). The biological effect of PEDF and synthetic peptides PEDF-6 and HLDF-6 on the HL-60 cells and the early gastrula ectoderm of Xenopus laevis embryos was studied. On the basis of the structural and functional homologies of HLDF, PEDF, and their homologous peptides and the computer models of the spatial structures of the full-size PEDF and the PEDF with the C-terminal fragment split off tby the cleavage of the Leu380-Thr381 bond in the serpin loop, a hypothesis on the functional role of the serpin loop in PEDF was put forward.

We calculated 219 topological, electronic, and steric indices for each of 59 studied embryotoxic benzylalkylamines and indolamines. Statistically significant equations correlating embryotoxicity and five calculated parameters were obtained using the method of step multiple regression. The prognostic power of the equation was 88-90%. Discriminant functions obtained by step discriminant analysis allowed prediction of the superactivity of benzylalkylamines and indolamines with 83-87% probability.