Biomedical subjects

Ramy Farid

Publications and source records attributed to Ramy Farid.

6 recordsLinked to original sources

PubMed2006

Novel procedure for modeling ligand/receptor induced fit effects.

We present a novel protein-ligand docking method that accurately accounts for both ligand and receptor flexibility by iteratively combining rigid receptor docking (Glide) with protein structure prediction (Prime) techniques. While traditional rigid-receptor docking methods are useful when the receptor structure does not change substantially upon ligand binding, success is limited when the protein must be "induced" into the correct binding conformation for a given ligand. We provide an in-depth description of our novel methodology and present results for 21 pharmaceutically relevant examples. Traditional rigid-receptor docking for these 21 cases yields an average RMSD of 5.5 A. The average ligand RMSD for docking to a flexible receptor for the 21 pairs is 1.4 A; the RMSD is < or =1.8 A for 18 of the cases. For the three cases with RMSDs greater than 1.8 A, the core of the ligand is properly docked and all key protein/ligand interactions are captured.

Binding Sites↗

PubMed2006

New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies.

We created a homology model of the homo-tetrameric pore domain of HERG using the crystal structure of the bacterial potassium channel, KvAP, as a template. We docked a set of known blockers with well-characterized effects on channel function into the lumen of the pore between the selectivity filter and extracellular entrance using a novel docking and refinement procedure incorporating Glide and Prime. Key aromatic groups of the blockers are predicted to form multiple simultaneous ring stacking and hydrophobic interactions among the eight aromatic residues lining the pore. Furthermore, each blocker can achieve these interactions via multiple docking configurations. To further interpret the docking results, we mapped hydrophobic and hydrophilic potentials within the lumen of each refined docked complex. Hydrophilic iso-potential contours define a 'propeller-shaped' volume at the selectivity filter entrance. Hydrophobic contours define a hollow 'crown-shaped' volume located above the 'propeller', whose hydrophobic 'rim' extends along the pore axis between Tyr652 and Phe656. Blockers adopt conformations/binding orientations that closely mimic the shapes and properties of these contours. Blocker basic groups are localized in the hydrophilic 'propeller', forming electrostatic interactions with Ser624 rather than a generally accepted pi-cation interaction with Tyr652. Terfenadine, cisapride, sertindole, ibutilide, and clofilium adopt similar docked poses, in which their N-substituents bridge radially across the hollow interior of the 'crown' (analogous to the hub and spokes of a wheel), and project aromatic/hydrophobic portions into the hydrophobic 'rim'. MK-499 docks with its longitudinal axis parallel to the axis of the pore and 'crown', and its hydrophobic groups buried within the hydrophobic 'rim'.

Conserved Sequence↗

PubMed2006

Use of an induced fit receptor structure in virtual screening.

Structured-based drug design has traditionally relied on a single receptor structure as a target for docking and screening studies. However, it has become increasingly clear that in many cases where protein flexibility is an issue, it is critical to accurately model ligand-induced receptor movement in order to obtain high enrichment factors. We present a novel protein-ligand docking method that accounts for both ligand and receptor flexibility and accurately predicts the conformation of protein-ligand binding complexes. This method can generate viable receptor ensembles that can be used in virtual database screens.

Databases, Protein↗

PubMed2005

Crystal structures of two bacterial 3-hydroxy-3-methylglutaryl-CoA lyases suggest a common catalytic mechanism among a family of TIM barrel metalloenzymes cleaving carbon-carbon bonds.

The enzyme 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) lyase catalyzes the terminal steps in ketone body generation and leucine degradation. Mutations in this enzyme cause a human autosomal recessive disorder called primary metabolic aciduria, which typically kills victims because of an inability to tolerate hypoglycemia. Here we present crystal structures of the HMG-CoA lyases from Bacillus subtilis and Brucella melitensis at 2.7 and 2.3 A resolution, respectively. These enzymes share greater than 45% sequence identity with the human orthologue. Although the enzyme has the anticipated triose-phosphate isomerase (TIM) barrel fold, the catalytic center contains a divalent cation-binding site formed by a cluster of invariant residues that cap the core of the barrel, contrary to the predictions of homology models. Surprisingly, the residues forming this cation-binding site and most of their interaction partners are shared with three other TIM barrel enzymes that catalyze diverse carbon-carbon bond cleavage reactions believed to proceed through enolate intermediates (4-hydroxy-2-ketovalerate aldolase, 2-isopropylmalate synthase, and transcarboxylase 5S). We propose the name "DRE-TIM metallolyases" for this newly identified enzyme family likely to employ a common catalytic reaction mechanism involving an invariant Asp-Arg-Glu (DRE) triplet. The Asp ligates the divalent cation, while the Arg probably stabilizes charge accumulation in the enolate intermediate, and the Glu maintains the precise structural alignment of the Asp and Arg. We propose a detailed model for the catalytic reaction mechanism of HMG-CoA lyase based on the examination of previously reported product complexes of other DRE-TIM metallolyases and induced fit substrate docking studies conducted using the crystal structure of human HMG-CoA lyase (reported in the accompanying paper by Fu, et al. (2006) J. Biol. Chem. 281, 7526-7532). Our model is consistent with extensive mutagenesis results and can guide subsequent studies directed at definitive experimental elucidation of this enzyme's reaction mechanism.

2-Isopropylmalate Synthase↗

PubMed2005

Thermal equilibrium of high- and low-spin forms of cytochrome P450 BM-3: repositioning of the substrate?

We demonstrate that cytochrome P450 BM-3 in complex with N-palmitoylglycine undergoes a spin state change between room temperature, where optimal activity is seen, and low temperatures, where X-ray diffraction characterization has been carried out. On the basis of NMR measurements of the full-length protein, this spin state change is likely to be accompanied by a general structural rearrangement in the enzyme pocket. The substrate remains bound at all temperatures. We propose that the substrate may "slide" from a position directly atop the heme (thus displacing the ligating water) to the more distant position (thus restoring the ligating water) as the temperature is lowered. This proposal is evaluated on the basis of computational modeling of the protein-ligand complex, using a novel induced fit methodology. We thereby generate a structure with the ligand in close contact with the heme, similar in energy to the experimental structure. With this combination of theory and experiment we provide a specific proposal of how ligands may be positioned for chemistry for this enzyme.

Bacterial Proteins↗

PubMed2005

Integrated Modeling Program, Applied Chemical Theory (IMPACT).

We provide an overview of the IMPACT molecular mechanics program with an emphasis on recent developments and a description of its current functionality. With respect to core molecular mechanics technologies we include a status report for the fixed charge and polarizable force fields that can be used with the program and illustrate how the force fields, when used together with new atom typing and parameter assignment modules, have greatly expanded the coverage of organic compounds and medicinally relevant ligands. As we discuss in this review, explicit solvent simulations have been used to guide our design of implicit solvent models based on the generalized Born framework and a novel nonpolar estimator that have recently been incorporated into the program. With IMPACT it is possible to use several different advanced conformational sampling algorithms based on combining features of molecular dynamics and Monte Carlo simulations. The program includes two specialized molecular mechanics modules: Glide, a high-throughput docking program, and QSite, a mixed quantum mechanics/molecular mechanics module. These modules employ the IMPACT infrastructure as a starting point for the construction of the protein model and assignment of molecular mechanics parameters, but have then been developed to meet specialized objectives with respect to sampling and the energy function.

Computer Simulation↗