Biomedical subjects

I Bahar

Publications and source records attributed to I Bahar.

At least 19 recordsLinked to original sources

PubMed2003

Allosteric changes in protein structure computed by a simple mechanical model: hemoglobin T<-->R2 transition.

Information on protein dynamics has been usually inferred from spectroscopic studies of parts of the proteins, or indirectly from the comparison of the conformations assumed in the presence of different substrates or ligands. While molecular simulations also provide information on protein dynamics, they usually suffer from incomplete sampling of conformational space, and become prohibitively expensive when exploring the collective dynamics of large macromolecular structures. Here, we explore the dynamics of a well-studied allosteric protein, hemoglobin (Hb), to show that a simple mechanical model based on Gaussian fluctuations of residues can efficiently predict the transition between the tense (T, unliganded) and relaxed (R or R2, O(2) or CO-bound) forms of Hb. The passage from T into R2 is shown to be favored by the global mode of motion, which, in turn is driven by entropic effects. The major difference between the dynamics of the T and R2 forms is the loss of the hinge-bending role of alpha(1)-beta(2) (or alpha(2)-beta(1)) interfacial residues at alpha Phe36-His45 and beta Thr87-Asn102 in the R2 form, which implies a decreased cooperativity in the higher affinity (R2) form of Hb, consistent with many experimental studies. The involvement of the proximal histidine beta His92 in this hinge region suggests that the allosteric propagation of the local structural changes (induced upon O(2) binding) into global ones occur via hinge regions. This is the first demonstration that there is an intrinsic tendency of Hb to undergo T-->R2 transition, induced by purely elastic forces of entropic origin that are uniquely defined for the particular contact topology of the T form.

Animals↗

PubMed2002

Molecular mechanisms of chaperonin GroEL-GroES function.

The dynamics of the GroEL-GroES complex is investigated with a coarse-grained model. This model is one in which single-residue points are connected to other such points, which are nearby, by identical springs, forming a network of interactions. The nature of the most important (slowest) normal modes reveals a wide variety of motions uniquely dependent upon the central cavity of the structure, including opposed torsional rotation of the two GroEL rings accompanied by the alternating compression and expansion of the GroES cap binding region, bending, shear, opposed radial breathing of the cis and trans rings, and stretching and contraction along the protein assembly's long axis. The intermediate domains of the subunits are bifunctional due to the presence of two hinges, which are alternatively activated or frozen by an ATP-dependent mechanism. ATP binding stabilizes a relatively open conformation (with respect to the central cavity) and hinders the motion of the hinge site connecting the intermediate and equatorial domains, while enhancing the flexibility of the second hinge that sets in motion the apical domains. The relative flexibilities of the hinges are reversed in the nucleotide-free form. Cooperative cross-correlations between subunits provide information about the mechanism of action of the protein. The mechanical motions driven by the different modes provide variable binding surfaces and variable sized cavities in the interior to enable accommodation of a broad range of protein substrates. These modes of motion could be used to manipulate the substrate's conformations.

Adenosine Triphosphate↗

PubMed2001

Anisotropy of fluctuation dynamics of proteins with an elastic network model.

Fluctuations about the native conformation of proteins have proven to be suitably reproduced with a simple elastic network model, which has shown excellent agreement with a number of different properties for a wide variety of proteins. This scalar model simply investigates the magnitudes of motion of individual residues in the structure. To use the elastic model approach further for developing the details of protein mechanisms, it becomes essential to expand this model to include the added details of the directions of individual residue fluctuations. In this paper a new tool is presented for this purpose and applied to the retinol-binding protein, which indicates enhanced flexibility in the region of entry to the ligand binding site and for the portion of the protein binding to its carrier protein.

Anisotropy↗

PubMed2001

Transition states and the meaning of Phi-values in protein folding kinetics.

What is the mechanism of two-state protein folding? The rate-limiting step is typically explored through a Phi-value, which is the mutation-induced change in the transition state free energy divided by the change in the equilibrium free energy of folding. Phi-values ranging from 0 to 1 have been interpreted as meaning the transition state is denatured-like (0), native-like (1) or in-between. But there is no classical interpretation for the experimental Phi-values that are negative or >1. Using a rigorous method to identity transition states via an exact lattice model, we find that nonclassical Phi-values can arise from parallel microscopic flow processes, such as those in funnel-shaped energy landscapes. Phi < 0 results when a mutation destabilizes a slow flow channel, causing a backflow into a faster flow channel. Phi > 1 implies the reverse: a backflow from a fast channel into a slow one. Using a 'landscape mapping' method, we find that Phi correlates with the acceleration/deceleration of folding induced by mutations, rather than with the degree of nativeness of the transition state.

Kinetics↗

PubMed2000

Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor.

The dynamics of alpha-amylase inhibitors has been investigated using molecular dynamics (MD) simulations and two analytical approaches, the Gaussian network model (GNM) and anisotropic network model (ANM). MD simulations use a full atomic approach with empirical force fields, while the analytical approaches are based on a coarse-grained single-site-per-residue model with a single-parameter harmonic potential between sufficiently close (r </= 7 A) residue pairs. The major difference between the GNM and the ANM is that no directional preferences can be obtained in the GNM, all residue fluctuations being theoretically isotropic, while ANM does incorporate directional preferences. The dominant modes of motions are identified by (i) the singular value decomposition (SVD) of the MD trajectory matrices, and (ii) the similarity transformation of the Kirchhoff matrices of inter-residue contacts in the GNM or ANM. The mean-square fluctuations of individual residues and the cross-correlations between domain movements retain the same characteristics, in all approaches-although the dispersion of modes and detailed amplitudes of motion obtained in the ANM conform more closely with MD results. The major weakness of the analytical approaches appears, on the other hand, to be their inadequacy to account for the anharmonic motions or multimeric transitions driven by the slowest collective mode observed in MD. Such motions usually suffer, however, from MD sampling inefficiencies, and multiple independent runs should be tested before making conclusions about their validity and detailed mechanisms. Overall this study invites attention to (i) the robustness of the average properties (mean-square fluctuations, cross-correlations) controlled by the low frequency motions, which are invariably reproduced in all approaches, and (ii) the utility and efficiency of the ANM, the computational time cost of which is of the order of "minutes" (real time), as opposed to "days" for MD simulations. Proteins 2000;40:512-524.

Computer Simulation↗

PubMed2000

Proteins with similar architecture exhibit similar large-scale dynamic behavior.

We have investigated the similarities and differences in the computed dynamic fluctuations exhibited by six members of a protein fold family with a coarse-grained Gaussian network model. Specifically, we consider the cofactor binding fragment of CysB; the lysine/arginine/ornithine-binding protein (LAO); the enzyme porphobilinogen deaminase (PBGD); the ribose-binding protein (RBP); the N-terminal lobe of ovotransferrin in apo-form (apo-OVOT); and the leucine/isoleucine/valine-binding protein (LIVBP). All have domains that resemble a Rossmann fold, but there are also some significant differences. Results indicate that similar global dynamic behavior is preserved for the members of a fold family, and that differences usually occur in regions only where specific function is localized. The present work is a computational demonstration that the scaffold of a protein fold may be utilized for diverse purposes. LAO requires a bound ligand before it conforms to the large-scale fluctuation behavior of the three other members of the family, CysB, PBGD, and RBP, all of which contain a substrate (cofactor) at the active site cleft. The dynamics of the ligand-free enzymes LIVBP and apo-OVOT, on the other hand, concur with that of unliganded LAO. The present results suggest that it is possible to construct structure alignments based on dynamic fluctuation behavior.

Apoproteins↗

PubMed2000

Characterization of anticancer agents by their growth inhibitory activity and relationships to mechanism of action and structure.

An analysis of the growth inhibitory potency of 122 anticancer agents available from the National Cancer Institute anticancer drug screen is presented. Methods of singular value decomposition (SVD) were applied to determine the matrix of distances between all compounds. These SVD-derived dissimilarity distances were used to cluster compounds that exhibit similar tumor growth inhibitory activity patterns against 60 human cancer cell lines. Cluster analysis divides the 122 standard agents into 25 statistically distinct groups. The first eight groups include structurally diverse compounds with reactive functionalities that act as DNA-damaging agents while the remaining 17 groups include compounds that inhibit nucleic acid biosynthesis and mitosis. Examination of the average activity patterns across the 60 tumor cell lines reveals unique 'fingerprints' associated with each group. A diverse set of structural features are observed for compounds within these groups, with frequent occurrences of strong within-group structural similarities. Clustering of cell types by their response to the 122 anticancer agents divides the 60 cell types into 21 groups. The strongest within-panel groupings were found for the renal, leukemia and ovarian cell panels. These results contribute to the basis for comparisons between log(GI(50)) screening patterns of the 122 anticancer agents and additional tested compounds.

Animals↗

PubMed1999

Structure-based analysis of protein dynamics: comparison of theoretical results for hen lysozyme with X-ray diffraction and NMR relaxation data.

An analytical approach based on Gaussian network model (GNM) is proposed for predicting the rotational dynamics of proteins. The method, previously shown to successfully reproduce X-ray crystallographic temperature factors for a series of proteins is extended here to predict bond torsional mobilities and reorientation of main chain amide groups probed by 15N-H nuclear magnetic resonance (NMR) relaxation. The dynamics of hen egg-white lysozyme (HEWL) in the folded state is investigated using the proposed approach. Excellent agreement is observed between theoretical results and experimental (X-ray diffraction and NMR relaxation) data. The analysis reveals the important role of coupled rotations, or cross-correlations between dihedral angle librations, in defining the relaxation mechanism on a local scale. The crystal and solution structures exhibit some differences in their local motions, but their global motions are identical. Hinge residues mediating the cooperative movements of the alpha- and beta-domains are identified, which comprise residues in helix C, Glu35 and Ser36 on the loop succeeding helix B, Ile55 and Leu56 at the turn between strands II and III. The central part of the beta-domain long loop and the turn between strands I and II display an enhanced mobility. Finally, kinetically hot residues and key interactions are identified, which point at helix B and beta-strand III as the structural elements underlying the stability of the tertiary structure.

Animals↗

PubMed1999

Cooperative fluctuations and subunit communication in tryptophan synthase.

Tryptophan synthase (TRPS), with linearly arrayed subunits alphabetabetaalpha, catalyzes the last two reactions in the biosynthesis of L-tryptophan. The two reactions take place in the respective alpha- and beta-subunits of the enzyme, and the intermediate product, indole, is transferred from the alpha- to the beta-site through a 25 A long hydrophobic tunnel. The occurrence of a unique ligand-mediated long-range cooperativity for substrate channeling, and a quest to understand the mechanism of allosteric control and coordination in metabolic cycles, have motivated many experimental studies on the structure and catalytic activity of the TRPS alpha2beta2 complex and its mutants. The dynamics of these complexes are analyzed here using a simple but rigorous theoretical approach, the Gaussian network model. Both wild-type and mutant structures, in the unliganded and various liganded forms, are considered. The substrate binding site in the beta-subunit is found to be closely coupled to a group of hinge residues (beta77-beta89 and beta376-beta379) near the beta-beta interface. These residues simultaneously control the anticorrelated motion of the two beta-subunits, and the opening or closing of the hydrophobic tunnel. The latter process is achieved by the large amplitude fluctuations of the so-called COMM domain in the same subunit. Intersubunit communications are strengthened in the presence of external aldimines bound to the beta-site. The motions of the COMM core residues are coordinated with those of the alpha-beta hinge residues beta174-beta179 on the interfacial helix betaH6 at the entrance of the hydrophobic tunnel. And the motions of betaH6 are coupled, via helix betaH1 and alphaL6, to those of the loop alphaL2 that includes the alpha-subunit catalytically active residue Asp60. Overall, our analysis sheds light on the molecular machinery underlying subunit communication, and identifies the residues playing a key role in the cooperative transmission of conformational motions across the two reaction sites.

Allosteric Regulation↗

PubMed1999

Collective motions in HIV-1 reverse transcriptase: examination of flexibility and enzyme function.

In order to study the inferences of structure for mechanism, the collective motions of the retroviral reverse transcriptase HIV-1 RT (RT) are examined using the Gaussian network model (GNM) of proteins. This model is particularly suitable for elucidating the global dynamic characteristics of large proteins such as the presently investigated heterodimeric RT comprising a total of 982 residues. Local packing density and coordination order of amino acid residues is inspected by the GNM to determine the type and range of motions, both at the residue level and on a global scale, such as the correlated movements of entire subdomains. Of the two subunits, p66 and p51, forming the RT, only p66 has a DNA-binding cleft and a functional polymerase active site. This difference in the structure of the two subunits is shown here to be reflected in their dynamic characteristics: only p66 has the potential to undergo large-scale cooperative motions in the heterodimer, while p51 is essentially rigid. Taken together, the global motion of the RT heterodimer is comprised of movements of the p66 thumb subdomain perpendicular to those of the p66 fingers, accompanied by anticorrelated fluctuations of the RNase H domain and p51 thumb, thus providing information about the details of one processivity mechanism. A few clusters of residues, generally distant in sequence but close in space, are identified in the p66 palm and connection subdomains, which form the hinge-bending regions that control the highly concerted motion of the subdomains. These regions include the catalytically active site and the non-nucleoside inhibitor binding pocket of p66 polymerase, as well as sites whose mutations have been shown to impair enzyme activity. It is easily conceivable that this hinge region, indicated by GNM analysis to play a critical role in modulating the global motion, is locked into an inactive conformation upon binding of an inhibitor. Comparative analysis of the dynamic characteristics of the unliganded and liganded dimers indicates severe repression of the mobility of the p66 thumb in RT's global mode, upon binding of non-nucleoside inhibitors.

Binding Sites↗

PubMed1998

RNA bulge entropies in the unbound state correlate with peptide binding strengths for HIV-1 and BIV TAR RNA because of improved conformational access.

For the binding of peptides to wild-type HIV-1 and BIV TAR RNA and to mutants with bulges of various sizes, changes in the DeltaDelta G values of binding were determined from experimental K d values. The corresponding entropies of these bulges are estimated by enumerating all possible RNA bulge conformations on a lattice and then applying the Boltzmann relationship. Independent calculations of entropies from fluctuations are also carried out using the Gaussian network model (GNM) recently introduced for analyzing folded structures. Strong correlations are seen between the changes in free energy determined for binding and the two different unbound entropy calculations. The fact that the calculated entropy increase with larger bulge size is correlated with the enhanced experimental binding free energy is unusual. This system exhibits a dependence on the entropy of the unbound form that is opposite to usual binding models. Instead of a large initial entropy being unfavorable since it would be reduced upon binding, here the larger entropies actually favor binding. Several interpretations are possible: (i) the higher conformational freedom implies a higher competence for binding with a minimal strain, by suitable selection amongst the set of already accessible conformations; (ii) larger bulge entropies enhance the probability of the specific favorable conformation of the bound state; (iii) the increased freedom of the larger bulges contri-butes more to the bound state than to the unbound state; (iv) indirectly the large entropy of the bound state might have an unfavorable effect on the solvent structure. Nonetheless, this unusual effect is interesting.

Animals↗

PubMed1998

Vibrational dynamics of transfer RNAs: comparison of the free and synthetase-bound forms.

The vibrational dynamics of transfer RNAs, both free, and complexed with the cognate synthetase, are analyzed using a model (Gaussian network model) which recently proved to satisfactorily describe the collective motions of folded proteins. The approach is similar to a normal mode analysis, with the major simplification that no residue specificity is taken into consideration, which permits us (i) to cast the problem into an analytical form applicable to biomolecular systems including about 10(3 )residues, and (ii) to acquire information on the essential dynamics of such large systems within computational times at least two orders of magnitude shorter than conventional simulations. On a local scale, the fluctuations calculated for yeast tRNAPhe and tRNAAsp in the free state, and for tRNAGln complexed with glutaminyl-tRNA synthetase (GlnRS) are in good agreement with the corresponding crystallographic B factors. On a global scale, a hinge-bending region comprising nucleotides U8 to C12 in the D arm, G20 to G22 in the D loop, and m7G46 to C48 in the variable loop (for tRNAPhe), is identified in the free tRNA, conforming with previous observations. The two regions subject to the largest amplitude anticorrelated fluctuations in the free form, i.e. the anticodon region and the acceptor arm are, at the same time, the regions that experience the most severe suppression in their flexibilities upon binding to synthetase, suggesting that their sampling of the conformational space facilitates their recognition by the synthetase. Likewise, examination of the global mode of motion of GlnRS in the complex indicates that residues 40 to 45, 260 to 270, 306 to 314, 320 to 327 and 478 to 485, all of which cluster near the ATP binding site, form a hinge-bending region controlling the cooperative motion, and thereby the catalytic function, of the enzyme. The distal beta-barrel and the tRNA acceptor binding domain, on the other hand, are distinguished by their high mobilities in the global modes of motion, a feature typical of recognition sites, also observed for other proteins. Most of the conserved bases and residues of tRNA and GlnRS are severely constrained in the global motions of the molecules, suggesting their having a role in stabilizing and modulating the global motion.

Amino Acyl-tRNA Synthetases↗

PubMed1998

Recognition of native structure from complete enumeration of low-resolution models with constraints.

Complete sets of low-resolution conformations are generated for eight small proteins by rotating the Calpha-Calpha virtual bonds at selected flexible regions, while the remaining structural elements are assumed to move in rigid blocks. Several filtering criteria are used to reduce the ensemble size and to ensure the sampling of well-constructed conformations. These filters, based on structure and energy constraints deduced from knowledge-based studies, include the excluded volume requirement, the radius of gyration constraint, and the occurrence of sufficiently strong attractive inter-residue potentials to stabilize compact forms. About 8,000 well-constructed decoys or "probable folds" (PFs) are constructed for each protein. A correlation between root-mean-square (rms) deviations from X-ray structure and total energies is observed, revealing a decrease in energy as the rms deviation decreases. The conformation with the lowest energy exhibits an rms deviation smaller than 3.0 A, in most of the proteins considered. The results are highly sensitive to the choice of flexible regions. A strong tendency to assume native state rotational angles is revealed for some flexible bonds from the analysis of the distributions of dihedral angles in the PFs, suggesting the formation of foldons near these locally stable regions at early folding pathway.

Computer Simulation↗

PubMed1998

Coarse-grained simulations of conformational dynamics of proteins: application to apomyoglobin.

A coarse-grained dynamic Monte Carlo method is proposed for investigating the conformational dynamics of proteins. Each residue is represented by two interaction sites, one at the alpha-carbon, and the other on the amino acid sidechain. Geometry and energy parameters extracted from databank structures are used. The method is applied to the crystal structure of apomyoglobin (apo-Mb). Equilibrium and dynamic properties of apo-Mb are characterized within computation times one order of magnitude shorter than conventional molecular dynamics (MD) simulations. The calculated rms fluctuations in alpha-carbons are in good agreement with crystallographic temperature factors. Regions exhibiting enhanced conformational mobilities are identified. Among the loops connecting the eight helices A to H, the loop CD undergoes the fastest motions, leading to partial unwinding of helix D. Helix G is the most stable helix on the basis of the kinetic stability of dihedral angles, followed by the respective helices A, E, H, and B. These results, in agreement with H/D exchange and two-dimensional NMR experiments, as well as with MD simulations, lend support to the use of the proposed approach as an efficient, yet physically plausible, means of characterizing protein conformational dynamics.

Animals↗

PubMed1998

Correlation between native-state hydrogen exchange and cooperative residue fluctuations from a simple model.

Recently, we developed a simple analytical model based on local residue packing densities and the distribution of tertiary contacts for describing the conformational fluctuations of proteins in their folded state. This so-called Gaussian network model (GNM) is applied here to the interpretation of experimental hydrogen exchange (HX) behavior of proteins in their native state or under weakly denaturing conditions. Calculations are performed for five proteins: bovine pancreatic trypsin inhibitor, cytochrome c, plastocyanin, staphylococcal nuclease, and ribonuclease H. The results are significant in two respects. First, a good agreement is reached between calculated fluctuations and experimental measurements of HX despite the simplicity of the model and within computational times 2 or 3 orders of magnitude faster than earlier, more complex simulations. Second, the success of a theory, based on the coupled conformational fluctuations of residues near the native state, to satisfactorily describe the native-state HX behavior indicates the significant contribution of local, but cooperative, fluctuations to protein conformational dynamics. The correlation between the HX data and the unfolding kinetics of individual residues further suggests that local conformational susceptibilities as revealed by the GNM approach may have implications relevant to the global dynamics of proteins.

Aprotinin↗

PubMed1998

Identification of kinetically hot residues in proteins.

A number of recent studies called attention to the presence of kinetically important residues underlying the formation and stabilization of folding nuclei in proteins, and to the possible existence of a correlation between conserved residues and those participating in the folding nuclei. Here, we use the Gaussian network model (GNM), which recently proved useful in describing the dynamic characteristics of proteins for identifying the kinetically hot residues in folded structures. These are the residues involved in the highest frequency fluctuations near the native state coordinates. Their high frequency is a manifestation of the steepness of the energy landscape near their native state positions. The theory is applied to a series of proteins whose kinetically important residues have been extensively explored: chymotrypsin inhibitor 2, cytochrome c, and related C2 proteins. Most of the residues previously pointed out to underlie the folding process of these proteins, and to be critically important for the stabilization of the tertiary fold, are correctly identified, indicating a correlation between the kinetic hot spots and the early forming structural elements in proteins. Additionally, a strong correlation between kinetically hot residues and loci of conserved residues is observed. Finally, residues that may be important for the stability of the tertiary structure of CheY are proposed.

Amino Acid Sequence↗

PubMed1998

Empirical solvent-mediated potentials hold for both intra-molecular and inter-molecular inter-residue interactions.

Whether knowledge-based intra-molecular inter-residue potentials are valid to represent inter-molecular interactions taking place at protein-protein interfaces has been questioned in several studies. Differences in the chain connectivity effect and in residue packing geometry between interfaces and single chain monomers have been pointed out as possible sources of distinct energetics for the two cases. In the present study, the interfacial regions of protein-protein complexes are examined to extract inter-molecular inter-residue potentials, using the same statistical methods as those previously adopted for intra-molecular residue pairs. Two sets of energy parameters are derived, corresponding to solvent-mediation and "average residue" mediation. The former set is shown to be highly correlated (correlation coefficient 0.89) with that previously obtained for inter-residue interactions within single chain monomers, while the latter exhibits a weaker correlation (0.69) with its intra-molecular counterpart. In addition to the close similarity of intra- and inter-molecular solvent-mediated potentials, they are shown to be significantly more residue-specific and thereby discriminative compared to the residue-mediated ones, indicating that solvent-mediation plays a major role in controlling the effective inter-residue interactions, either at interfaces, or within single monomers. Based on this observation, a reduced set of energy parameters comprising 20 one-body and 3 two-body terms is proposed (as opposed to the 20 x 20 tables of inter-residue potentials), which reproduces the conventional 20 x 20 tables with a correlation coefficient of 0.99.

Amino Acid Sequence↗

PubMed1998

Packing of sidechains in low-resolution models for proteins.

BACKGROUND: Atomic level rotamer libraries for sidechains in proteins have been proposed by several groups. Conformations of side groups in coarse-grained models, on the other hand, have not yet been analyzed, although low resolution approaches are the only efficient way to explore global structural features. RESULTS: A residue-specific backbone-dependent library for sidechain isomers, compatible with a coarse-grained model, is proposed. The isomeric states are utilized in packing sidechains of known backbone structures. Sidechain positions are predicted with a root-mean-square deviation (r.m.s.d.) of 2.40 A with respect to crystal structure for 50 test proteins. The rmsd for core residues is 1.60 A and decreases to 1.35 A when conformational correlations and directional effects in inter-residue couplings are considered. CONCLUSIONS: An automated method for assigning sidechain positions in coarse-grained model proteins is proposed and made available on the internet; the method accounts satisfactorily for sidechain packing, particularly in the core.

Animals↗