Effects of limiting extension at the alphaIIb genu on ligand binding to integrin alphaIIbbeta3.

Structural data of integrin alphaIIbbeta3 have been interpreted as supporting a model in which: 1) the receptor exists primarily in a "bent," low affinity conformation on unactivated platelets and 2) activation induces an extended, high affinity conformation prior to, or following, ligand binding. Previous studies found that "clasping" the alphaIIb head domain to the beta3 tail decreased fibrinogen binding. To study the role of alphaIIb extension about the genu, we introduced a disulfide "clamp" between the alphaIIb thigh and calf-1 domains.

Targeted molecular dynamics reveals overall common conformational changes upon hybrid domain swing-out in beta3 integrins.

The beta3 integrin family members alphaIIbeta3 and alphaVbeta3 signal bidirectionally through long-range allosteric changes, including a transition from a bent unliganded-closed low-affinity state to an extended liganded-open high-affinity state. To obtain an atomic-level description of this transition in an explicit solvent, we carried out targeted molecular dynamics simulations of the headpieces of alphaIIbeta3 and alphaVbeta3 integrins. Although minor differences were observed between these receptors, our results suggest a common transition pathway in which the hybrid domain swing-out is accompanied by conformational changes within the beta3 betaA (I-like) domain that propagate through the alpha7 helix C-terminus, and are followed by the alpha7 helix downward motion and the opening of the beta6-alpha7 loop.

Structural and therapeutic insights from the species specificity and in vivo antithrombotic activity of a novel alphaIIb-specific alphaIIbbeta3 antagonist.

We previously reported on a novel compound (Compound 1; RUC-1) identified by high-throughput screening that inhibits human alphaIIbbeta3. RUC-1 did not inhibit alphaVbeta3, suggesting that it interacts with alphaIIb, and flexible ligand/rigid protein molecular docking studies supported this speculation. We have now studied RUC-1's effects on murine and rat platelets, which are less sensitive than human to inhibition by Arg-Gly-Asp (RGD) peptides due to differences in the alphaIIb sequences contributing to the binding pocket.

Functional and computational studies of the ligand-associated metal binding site of beta3 integrins.

A combination of experimental and computational approaches was used to provide a structural context for the role of the beta3 integrin subunit ligand-associated metal binding site (LIMBS) in the binding of physiological ligands to beta3 integrins. Specifically, we have carried out (1) adhesion assays on cells expressing normal alphaIIbeta3, normal alphaVbeta3, or the corresponding beta3 D217A LIMBS mutants; and (2) equilibrium and nonequilibrium (steered) molecular dynamics (MD) simulations of eptifibatide in complex with either a fully hydrated normal alphaIIbeta3 integrin fragment (alphaIIb beta-propeller and the beta3 betaA (I-like), hybrid, and PSI domains) or the equivalent beta3 D217A mutant. Normal alphaIIbeta3 expressing cells adhered to immobilized fibrinogen and echistatin, whereas cells expressing the alphaIIbeta3 D217A LIMBS mutant failed to adhere to either ligand.

Application of high-throughput screening to identify a novel alphaIIb-specific small- molecule inhibitor of alphaIIbbeta3-mediated platelet interaction with fibrinogen.

Small-molecule alphaIIbbeta3 antagonists competitively block ligand binding by spanning between the D224 in alphaIIb and the MIDAS metal ion in beta3. They variably induce conformational changes in the receptor, which may have undesirable consequences. To identify alphaIIbbeta3 antagonists with novel structures, we tested 33 264 small molecules for their ability to inhibit the adhesion of washed platelets to immobilized fibrinogen at 16 muM.

Requirements and ontology for a G protein-coupled receptor oligomerization knowledge base.

Background: G Protein-Coupled Receptors (GPCRs) are a large and diverse family of membrane proteins whose members participate in the regulation of most cellular and physiological processes and therefore represent key pharmacological targets. Although several bioinformatics resources support research on GPCRs, most of these have been designed based on the traditional assumption that monomeric GPCRs constitute the functional receptor unit. The increase in the frequency and number of reports about GPCR dimerization/oligomerization and the implication of oligomerization in receptor function makes necessary the ability to store and access information about GPCR dimers/oligomers electronically.

Mapping early conformational changes in alphaIIb and beta3 during biogenesis reveals a potential mechanism for alphaIIbbeta3 adopting its bent conformation.

Current evidence supports a model in which the low-affinity state of the platelet integrin alphaIIbbeta3 results from alphaIIbbeta3 adopting a bent conformation. To assess alphaIIbbeta3 biogenesis and how alphaIIbbeta3 initially adopts the bent conformation, we mapped the conformational states occupied by alphaIIb and beta3 during biogenesis using conformation-specific monoclonal antibodies (mAbs). We found that alphaIIbbeta3 complex formation was not limited by the availability of either free pro-alphaIIb or free beta3, suggesting that other molecules, perhaps chaperones, control complex formation.

Comparative binding energy analysis considering multiple receptors: a step toward 3D-QSAR models for multiple targets.

Comparative binding energy analysis, a technique to derive receptor-based three-dimensional quantitative structure-activity relationships (3D-QSAR), is herein extended to consider both affinity and selectivity in the derivation of the QSAR model. The extension is based on allowing multiple structurally related receptors to enter the X-matrix employed in the derivation of the structure-activity model. As a result, a single model common to all of them is obtained that considers both intra- and inter-receptor affinity differences for a given congeneric series.

Modeling the structure of the StART domains of MLN64 and StAR proteins in complex with cholesterol.

Steroidogenic acute regulatory protein-related lipid transfer (StART) domains are ubiquitously involved in intracellular lipid transport and metabolism and other cell-signaling events. In this work, we use a flexible docking algorithm, comparative modeling, and molecular dynamics (MD) simulations to generate plausible three-dimensional atomic models of the StART domains of human metastatic lymph node 64 (MLN64) and steroidogenic acute regulatory protein (StAR) proteins in complex with cholesterol. Our results show that cholesterol can adopt a similar conformation in the binding cavity in both cases and that the main contribution to the protein-ligand interaction energy derives from hydrophobic contacts.

Computational approaches to model ligand selectivity in drug design.

To be effective, a designed drug must discriminate successfully the macromolecular target from alternative structures present in the organism. The last few years have witnessed the emergence of different computational tools aimed to the understanding and modeling of this process at molecular level. Although still rudimentary, these methods are shaping a coherent approach to help in the design of molecules with high affinity and specificity, both in lead discovery and in lead optimization.

Gaussian mapping of chemical fragments in ligand binding sites.

We present a new approach to automatically define a quasi-optimal minimal set of pharmacophoric points mapping the interaction properties of a user-defined ligand binding site. The method is based on a fitting algorithm where a grid of sampled interaction energies of the target protein with small chemical fragments in the binding site is approximated by a linear expansion of Gaussian functions. A heuristic approximation selects from this expansion the smallest possible set of Gaussians required to describe the interaction properties of the binding site within a prespecified accuracy.

Virtual screening with flexible docking and COMBINE-based models. Application to a series of factor Xa inhibitors.

A two-step, fully automatic virtual screening procedure consisting of flexible docking followed by activity prediction by COMparative BINding Energy (COMBINE) analysis is presented. This novel approach has been successfully applied, as an example with medicinal chemistry interest, to a recently reported series of 133 factor Xa (fXa)(1) inhibitors whose activities encompass 4 orders of magnitude. The docking algorithm is linked to the COMBINE analysis program and used to derive independent regression models of the 133 inhibitors docked within three different fXa structures (PDB entries 1fjs, 1f0r, and 1xka), so as to explore the effect of receptor conformation on the overall results.

Benzimidazole derivatives. 4. The recognition of the voluminous substituent attached to the basic amino group of 5-HT4 receptor antagonists.

We report a structure-affinity analysis of an important element in the pharmacophore model for the recognition of 5-HT4 receptor antagonists: the voluminous substituent attached to the basic nitrogen of the ligand. We have designed, synthesized and pharmacologically evaluated a series of benzimidazole derivatives 1 containing a common molecular skeleton formed by N-[(4-piperidyl)methyl]-6-chlorobenzimidazole-4-carboxamide and four different substituents (R = butyl, 2-[(methylsulfonyl)amino]ethyl, 5-[(phenylacetyl)amino]pentyl, and 5-[(benzylsulfonyl)amino]pentyl). These compounds possess binding affinities in the nM range (Ki = 0.

Benzimidazole derivatives. 3. 3D-QSAR/CoMFA model and computational simulation for the recognition of 5-HT(4) receptor antagonists.

A three-dimensional quantitative structure-affinity relationship study (3D-QSAR), using the comparative molecular field analysis (CoMFA) method, and subsequent computational simulation of ligand recognition have been successfully applied to explain the binding affinities for the 5-HT(4) receptor (5-HT(4)R) of a series of benzimidazole-4-carboxamides and carboxylates derivatives 1-24. The K(i) values of these compounds are in the range from 0.11 to 10 000 nM.