Biomedical subjects

Ivet Bahar

Publications and source records attributed to Ivet Bahar.

At least 19 recordsLinked to original sources

PubMed2006

Optimal design of protein docking potentials: efficiency and limitations.

Protein-protein docking is a challenging computational problem in functional genomics, particularly when one or both proteins undergo conformational change(s) upon binding. The major challenge is to define scoring function soft enough to tolerate these changes and specific enough to distinguish between near-native and "misdocked" conformations. Using a linear programming technique, we derived protein docking potentials (PDPs) that comply with this requirement. We considered a set of 63 nonredundant complexes to this aim, and generated 400,000 putative docked complexes (decoys) based on shape complementarity criterion for each complex. The PDPs were required to yield for the native (correctly docked) structure a potential energy lower than those of all the nonnative (misdocked) structures. The energy constraints applied to all complexes led to ca. 25 million inequalities, the simultaneous solution of which yielded an optimal set of PDPs that discriminated the correctly docked (up to 4.0 A root-mean-square deviation from known complex structure) structure among the 85 top-ranking (0.02%) decoys in 59/63 examined bound-bound cases. The high performance of the potentials was further verified in jackknife tests and by ranking putative docked conformation submitted to CAPRI. In addition to their utility in identifying correctly folded complexes, the PDPs reveal biologically meaningful features that distinguish docking potentials from folding potentials.

Antigens↗

PubMed2005

Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state.

Protein-protein binding usually involves structural changes that may extend beyond the rearrangements on a local scale, and cannot be explained by a classical lock-and-key mechanism. Several models have been advanced to explain the flexible binding of proteins such as the induced fit mechanism where the ligand is postulated to induce a conformational change at the interaction site upon binding, or the preexisting equilibrium hypothesis that assumes that protein samples an ensemble of conformations at equilibrium conditions and that the ligand binds selectively to an active conformation. We explored the equilibrium motions of proteins that exhibit relatively large (nonlocal) conformational changes upon protein binding using the Gaussian network model and the anisotropic network model of protein dynamics. For four complexes, LIR-1/HLA-A2, Actin/DNase I, CDK2/cyclin, and CDK6/p16(INK4a), the motions calculated for the monomer exhibiting the largest conformational change, in its unbound (free) form, correlate with the experimentally observed structural changes upon binding. This study emphasizes the preexisting equilibrium/conformational selection as a mechanism for protein-protein interaction and lends support the concept that proteins, in their native conformation, are predisposed to undergo conformational fluctuations that are relevant to, or even required for, their biological functions.

Actins↗

PubMed2005

Elastic network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies.

With advances in structure genomics, it is now recognized that knowledge of structure alone is insufficient to understand and control the mechanisms of biomolecular function. Additional information in the form of dynamics is needed. As demonstrated in a large number of studies, the machinery of proteins and their complexes can be understood to a good approximation by adopting Gaussian (or elastic) network models (GNM) for simplified normal mode analyses. While this approximation lacks chemical details, it provides us with a means for assessing the collective motions of large structures/assemblies and perform a comparative analysis of a series of proteins, thus providing insights into the mechanical aspects of biomolecular dynamics. In this paper, we discuss recent applications of GNM to a series of enzymes as well as large structures such as the HK97 bacteriophage viral capsids. Understanding the dynamics of large protein structures can be computationally challenging. To this end, we introduce a new approach for building a hierarchical, reduced rank representation of the protein topology and consequently the fluctuation dynamics.

Bacteriophages↗

PubMed2005

iGNM: a database of protein functional motions based on Gaussian Network Model.

MOTIVATION: The knowledge of protein structure is not sufficient for understanding and controlling its function. Function is a dynamic property. Although protein structural information has been rapidly accumulating in databases, little effort has been invested to date toward systematically characterizing protein dynamics. The recent success of analytical methods based on elastic network models, and in particular the Gaussian Network Model (GNM), permits us to perform a high-throughput analysis of the collective dynamics of proteins. RESULTS: We computed the GNM dynamics for 20 058 structures from the Protein Data Bank, and generated information on the equilibrium dynamics at the level of individual residues. The results are stored on a web-based system called iGNM and configured so as to permit the users to visualize or download the results through a standard web browser using a simple search engine. Static and animated images for describing the conformational mobility of proteins over a broad range of normal modes are accessible, along with an online calculation engine available for newly deposited structures. A case study of the dynamics of 20 non-homologous hydrolases is presented to illustrate the utility of the iGNM database for identifying key residues that control the cooperative motions and revealing the connection between collective dynamics and catalytic activity.

Algorithms↗

PubMed2005

Structural dynamics of nucleosome core particle: comparison with nucleosomes containing histone variants.

The present study provides insights on the dominant mechanisms of motions of the nucleosome core particle and the changes in its functional dynamics in response to histone variants. Comparative analysis of the global dynamics of nucleosomes with native and variant H2A histones, using normal mode analysis revealed that the dynamics of the nucleosome is highly symmetric, and its interaction with the nucleosomal DNA plays a vital role in its regulation. The collective dynamics of nucleosomes are predicted to be dominated by two types of large-scale motions: (1) a global stretching-compression of nucleosome along the dyad axis by which the nucleosome undergoes a breathing motion with a massive distortion of nucleosomal DNA, modulated by histone-DNA interactions; and (2) the flipping (or bending) of both the sides of the nucleosome in an out-of-plane fashion with respect to the dyad axis, originated by the highly dynamic N-termini of H3 and (H2A.Z-H2B) dimer in agreement with the experimentally observed perturbed dynamics of the particular N-terminus under physiological conditions. In general, the nucleosomes with variant histones exhibit higher mobilities and weaker correlations between internal motions compared to the nucleosome containing ordinary histones. The differences are more pronounced at the L1 and L2 loops of the respective monomers H2B and H2A, and at the N-termini of the monomers H3 and H4, all of which closely interact with the wrapping DNA.

Amino Acid Motifs↗

PubMed2005

Coarse-grained normal mode analysis in structural biology.

The realization that experimentally observed functional motions of proteins can be predicted by coarse-grained normal mode analysis has renewed interest in applications to structural biology. Notable applications include the prediction of biologically relevant motions of proteins and supramolecular structures driven by their structure-encoded collective dynamics; the refinement of low-resolution structures, including those determined by cryo-electron microscopy; and the identification of conserved dynamic patterns and mechanically key regions within protein families. Additionally, hybrid methods that couple atomic simulations with deformations derived from coarse-grained normal mode analysis are able to sample collective motions beyond the range of conventional molecular dynamics simulations. Such applications have provided great insight into the underlying principles linking protein structures to their dynamics and their dynamics to their functions.

Binding Sites↗

PubMed2005

Maturation dynamics of bacteriophage HK97 capsid.

Maturation of the bacteriophage HK97 capsid requires a large conformational change of the virus capsid. Experimental studies have identified several intermediates along this maturation pathway. To gain insights into the molecular mechanisms of capsid maturation, we examined the fluctuation dynamics of the procapsid and mature capsid using a residue-level computational approach. The most cooperative motions of the procapsid are found to be consistent with the observed change in configuration that takes place during maturation. A few dominant modes of motion are sufficient to describe the anisotropic expansion that accompanies maturation. Based upon these modes, maturation is proposed to occur via an overall expansion and reconfiguration of the capsid initiated by puckering of the pentamers, followed by flattening and crosslinking of the hexameric subunits, and finally crosslinking of the pentameric subunits. The highly mobile E loops are stabilized by anchoring to highly stable residues belonging to neighboring subunits.

Capsid↗

PubMed2005

Coupling between catalytic site and collective dynamics: a requirement for mechanochemical activity of enzymes.

Growing evidence supports the view that enzymatic activity results from a subtle interplay between chemical kinetics and molecular motions. A systematic analysis is performed here to delineate the type and level of coupling between catalysis and conformational mechanics. The dynamics of a set of 98 enzymes representative of different EC classes are analyzed with the Gaussian network model (GNM) and compared with experimental data. In more than 70% of the examined enzymes, the global hinge centers predicted by the GNM are found to be colocalized with the catalytic sites experimentally identified. Low translational mobility (< 7%) is observed for the catalytic residues, consistent with the fine-tuned design of enzymes to achieve precise mechanochemical activities. Ligand binding sites, while closely neighboring catalytic sites, enjoy a moderate flexibility to accommodate the ligand binding. These findings could serve as additional criteria for assessing drug binding residues and could lessen the computational burden of substrate docking searches.

Binding Sites↗

PubMed2004

Escherichia coli adenylate kinase dynamics: comparison of elastic network model modes with mode-coupling (15)N-NMR relaxation data.

The dynamics of adenylate kinase of Escherichia coli (AKeco) and its complex with the inhibitor AP(5)A, are characterized by correlating the theoretical results obtained with the Gaussian Network Model (GNM) and the anisotropic network model (ANM) with the order parameters and correlation times obtained with Slowly Relaxing Local Structure (SRLS) analysis of (15)N-NMR relaxation data. The AMPbd and LID domains of AKeco execute in solution large amplitude motions associated with the catalytic reaction Mg(+2)*ATP + AMP --> Mg(+2)*ADP + ADP. Two sets of correlation times and order parameters were determined by NMR/SRLS for AKeco, attributed to slow (nanoseconds) motions with correlation time tau( perpendicular) and low order parameters, and fast (picoseconds) motions with correlation time tau( parallel) and high order parameters. The structural connotation of these patterns is examined herein by subjecting AKeco and AKeco*AP(5)A to GNM analysis, which yields the dynamic spectrum in terms of slow and fast modes. The low/high NMR order parameters correlate with the slow/fast modes of the backbone elucidated with GNM. Likewise, tau( parallel) and tau( perpendicular) are associated with fast and slow GNM modes, respectively. Catalysis-related domain motion of AMPbd and LID in AKeco, occurring per NMR with correlation time tau( perpendicular), is associated with the first and second collective slow (global) GNM modes. The ANM-predicted deformations of the unliganded enzyme conform to the functional reconfiguration induced by ligand-binding, indicating the structural disposition (or potential) of the enzyme to bind its substrates. It is shown that NMR/SRLS and GNM/ANM analyses can be advantageously synthesized to provide insights into the molecular mechanisms that control biological function.

Adenylate Kinase↗

PubMed2004

High-efficiency bypass of DNA damage by human DNA polymerase Q.

Endogenous DNA damage arises frequently, particularly apurinic (AP) sites. These must be dealt with by cells in order to avoid genotoxic effects. DNA polymerase theta; is a newly identified enzyme encoded by the human POLQ gene. We find that POLQ has an exceptional ability to bypass an AP site, inserting A with 22% of the efficiency of a normal template, and continuing extension as avidly as with a normally paired base. POLQ preferentially incorporates A opposite an AP site and strongly disfavors C. On nondamaged templates, POLQ makes frequent errors, incorporating G or T opposite T about 1% of the time. This very low fidelity distinguishes POLQ from other A-family polymerases. POLQ has three sequence insertions between conserved motifs in its catalytic site. One insert of approximately 22 residues into the tip of the polymerase thumb subdomain is predicted to confer considerable flexibility and additional DNA contacts to affect enzyme fidelity. POLQ is the only known enzyme that efficiently carries out both the insertion and extension steps for bypass of AP sites, commonly formed as endogenous genomic lesions.

Amino Acid Motifs↗

PubMed2004

Mining frequent patterns in protein structures: a study of protease families.

MOTIVATION: Analysis of protein sequence and structure databases usually reveal frequent patterns (FP) associated with biological function. Data mining techniques generally consider the physicochemical and structural properties of amino acids and their microenvironment in the folded structures. Dynamics is not usually considered, although proteins are not static, and their function relates to conformational mobility in many cases. RESULTS: This work describes a novel unsupervised learning approach to discover FPs in the protein families, based on biochemical, geometric and dynamic features. Without any prior knowledge of functional motifs, the method discovers the FPs for each type of amino acid and identifies the conserved residues in three protease subfamilies; chymotrypsin and subtilisin subfamilies of serine proteases and papain subfamily of cysteine proteases. The catalytic triad residues are distinguished by their strong spatial coupling (high interconnectivity) to other conserved residues. Although the spatial arrangements of the catalytic residues in the two subfamilies of serine proteases are similar, their FPs are found to be quite different. The present approach appears to be a promising tool for detecting functional patterns in rapidly growing structure databases and providing insights in to the relationship among protein structure, dynamics and function. AVAILABILITY: Available upon request from the authors.

Algorithms↗

PubMed2004

Discovery and characterization of novel small molecule inhibitors of human Cdc25B dual specificity phosphatase.

Cdc25A and Cdc25B dual-specificity phosphatases are key regulators of cell cycle transition and proliferation. They have oncogenic properties and are overexpressed in many human tumors. Because selective Cdc25 phosphatase inhibitors would be valuable biological tools and possible therapeutic agents, we have assayed a small molecule library for in vitro inhibition of Cdc25. We now report the identification of two new structurally distinct classes of Cdc25 inhibitors with cellular activity. The cyclopentaquinoline 3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-4,8-dicarboxylic acid (5661118) and the naphthofurandione 3-benzoyl-naphtho[1,2-b]furan-4,5-dione (5169131) had in vitro IC50 values of 2.5 to 11 microM against recombinant Cdc25 and were less potent inhibitors of other phosphatases. Unlike 5661118, 5169131 caused reversible inhibition of Cdc25B and displayed competitive inhibitor kinetics. No growth inhibitory activity was seen with 5661118, whereas 10 to 30 microM 5169131 caused G1/S and G2/M arrest. We also found that 5169131 inhibited human PC-3 prostate and MDA-MB-435 breast cancer cell proliferation. Concentration-dependent Tyr15 hyperphosphorylation was seen on cyclin-dependent kinase with a 1-h 5169131 treatment, consistent with Cdc25 inhibition. Cells resistant to DNA toposiomerase II inhibitors were as sensitive to 5169131 as parental cells, indicating that this quinone compound does not inhibit topoisomerase II in vivo. Molecular modeling was used to predict a potential interaction site between the inhibitor and Cdc25B and to provide insights as to the molecular origins of the experimental observations. Based on its kinetic profile and cellular activity, we suggest that 5169131 could be an excellent tool for further studies on the cellular roles of Cdc25.

Cell Cycle↗

PubMed2004

Identification of core amino acids stabilizing rhodopsin.

Rhodopsin is the only G protein-coupled receptor (GPCR) whose 3D structure is known; therefore, it serves as a prototype for studies of the GPCR family of proteins. Rhodopsin dysfunction has been linked to misfolding, caused by chemical modifications that affect the naturally occurring disulfide bond between C110 and C187. Here, we identify the structural elements that stabilize rhodopsin by computational analysis of the rhodopsin structure and comparison with data from previous in vitro mutational studies. We simulate the thermal unfolding of rhodopsin by breaking the native-state hydrogen bonds sequentially in the order of their relative strength, using the recently developed Floppy Inclusion and Rigid Substructure Topography (FIRST) method [Jacobs, D. J., Rader, A. J., Kuhn, L. A. & Thorpe, M. F. (2001) Proteins 44, 150-165]. Residues most stable under thermal denaturation are part of a core, which is assumed to be important for the formation and stability of folded rhodopsin. This core includes the C110-C187 disulfide bond at the center of residues forming the interface between the transmembrane and the extracellular domains near the retinal binding pocket. Fast mode analysis of rhodopsin using the Gaussian network model also identifies the disulfide bond and the retinal ligand binding pocket to be the most rigid region in rhodopsin. Experiments confirm that 90% of the amino acids predicted by the FIRST method to be part of the core cause misfolding upon mutation. The observed high degree of conservation (78.9%) of this disulfide bond across all GPCR classes suggests that it is critical for the stability and function of GPCRs.

Amino Acid Sequence↗

PubMed2004

Global ribosome motions revealed with elastic network model.

The motions of large systems such as the ribosome are not fully accessible with conventional molecular simulations. A coarse-grained, less-than-atomic-detail model such as the anisotropic network model (ANM) is a convenient informative tool to study the cooperative motions of the ribosome. The motions of the small 30S subunit, the larger 50S subunit, and the entire 70S assembly of the two subunits have been analyzed using ANM. The lowest frequency collective modes predicted by ANM show that the 50S subunit and 30S subunit are strongly anti-correlated in the motion of the 70S assembly. A ratchet-like motion is observed that corresponds well to the experimentally reported ratchet motion. Other slow modes are also examined because of their potential links to the translocation steps in the ribosome. We identify several modes that may facilitate the E-tRNA exiting from the assembly. The A-site t-RNA and P-site t-RNA are found to be strongly coupled and positively correlated in these slow modes, suggesting that the translocations of these two t-RNAs occur simultaneously, while the motions of the E-site t-RNA are less correlated, and thus less likely to occur simultaneously. Overall the t-RNAs exhibit relatively large deformations. Animations of these slow modes of motion can be viewed at.

Anisotropy↗

PubMed2004

A sequence alignment-independent method for protein classification.

Annotation of the rapidly accumulating body of sequence data relies heavily on the detection of remote homologues and functional motifs in protein families. The most popular methods rely on sequence alignment. These include programs that use a scoring matrix to compare the probability of a potential alignment with random chance and programs that use curated multiple alignments to train profile hidden Markov models (HMMs). Related approaches depend on bootstrapping multiple alignments from a single sequence. However, alignment-based programs have limitations. They make the assumption that contiguity is conserved between homologous segments, which may not be true in genetic recombination or horizontal transfer. Alignments also become ambiguous when sequence similarity drops below 40%. This has kindled interest in classification methods that do not rely on alignment. An approach to classification without alignment based on the distribution of contiguous sequences of four amino acids (4-grams) was developed. Interest in 4-grams stemmed from the observation that almost all theoretically possible 4-grams (20(4)) occur in natural sequences and the majority of 4-grams are uniformly distributed. This implies that the probability of finding identical 4-grams by random chance in unrelated sequences is low. A Bayesian probabilistic model was developed to test this hypothesis. For each protein family in Pfam-A and PIR-PSD, a feature vector called a probe was constructed from the set of 4-grams that best characterised the family. In rigorous jackknife tests, unknown sequences from Pfam-A and PIR-PSD were compared with the probes for each family. A classification result was deemed a true positive if the probe match with the highest probability was in first place in a rank-ordered list. This was achieved in 70% of cases. Analysis of false positives suggested that the precision might approach 85% if selected families were clustered into subsets. Case studies indicated that the 4-grams in common between an unknown and the best matching probe correlated with functional motifs from PRINTS. The results showed that remote homologues and functional motifs could be identified from an analysis of 4-gram patterns.

Algorithms↗

PubMed2004

Inflammatory modulation of hepatocyte apoptosis by nitric oxide: in vivo, in vitro, and in silico studies.

Nitric oxide (NO*) and its reaction products are key players in the physiology and pathophysiology of inflammatory settings such as sepsis and shock. The consequences of the expression of inducible NO* synthase (iNOS, NOS-2) can be either protective or damaging to the liver. We have delineated two distinct hepatoprotective actions of NO*: the stimulation of cyclic guanosine monophosphate and the inhibition of caspases by S-nitrosation. In contrast, iNOS/NO* promotes hepatocyte death under conditions of severe redox stress, such as hemorrhagic shock or ischemia/reperfusion. Redox stress activates an unknown molecular switch that transforms NO*, which is hepatoprotective under resting conditions, into an agent that induces hepatocyte death. We hypothesize that the magnitude of the redox stress is a major determinant for the effects of NO* on cell survival by controlling the chemical fate of NO*. To address this hypothesis, we have carried out studies in relevant in vivo and in vitro settings. Moreover, we have constructed an initial mathematical model of caspase activation and coupled it to a model describing some of the reactions of NO* in hepatocytes. Our studies suggest that modulation of iron, oxygen, and superoxide may dictate whether NO* is hepatoprotective or hepatotoxic.

Animals↗

PubMed2004

Conformational changes in HIV-1 reverse transcriptase induced by nonnucleoside reverse transcriptase inhibitor binding.

Nonnucleoside reverse transcriptase inhibitors (NNRTI) are a group of small hydrophobic compounds with diverse structures that specifically inhibit HIV-1 reverse transcriptase (RT). NNRTIs interact with HIV-1 RT by binding to a single site on the p66 subunit of the p66/p51 heterodimeric enzyme, termed the NNRTI-binding pocket (NNRTI-BP). This binding interaction results in both short-range and long-range distortions of RT structure. In this article, we review the structural, computational and experimental evidence of the NNRTI-induced conformational changes in HIV-1 RT and relate them to the mechanism by which these compounds inhibit HIV-1 reverse transcription.

Anti-HIV Agents↗

PubMed2003

The origin and extent of coarse-grained regularities in protein internal packing.

Despite the suitability of various lattice geometries for coarse-grained modeling of proteins, the actual packing geometry of residues in folded structures has remained largely unexplored. A strong tendency to assume a regular packing geometry is shown here by optimally reorienting and superimposing clusters of neighboring residues from databank structures examined on a coarse-grained (single-site-per-residue) scale. The orientation function (or order parameter) of the examined coordination clusters with respect to fcc lattice directions is found to be 0.82. The observed geometry, which may be termed an incomplete distorted face-centered cubic (fcc) packing, is apparently favored by the drive to maximize packing density, in a fashion analogous to the way identical spheres pack densely and follow fcc geometry. About 2/3 of all residues obey this packing geometry, while the remainder occupy other context-dependent positions. The preferred coordination directions show relatively small variations over the various amino acid types, consistent with uniform residue viewpoint. Both the extremes of solvent-exposed and completely buried residue neighborhoods approximate the same generic packing, the only difference being in the numbers (and not the orientations) of coordination sites that are occupied (or left void for solvent occupancy). We observe the prevalence of a rather uniform (tight) residue packing density throughout the structure, including even the residues packed near solvent-exposed regions. The observed orientation distribution reveals an underlying, intrinsic orientation lattice for proteins.

Amino Acid Sequence↗