SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25Linked to original sources

PubMed2004

Complementarity of structure ensembles in protein-protein binding.

Protein-protein association is often accompanied by changes in receptor and ligand structure. This interplay between protein flexibility and protein-protein recognition is currently the largest obstacle both to our understanding of and to the reliable prediction of protein complexes. We performed two sets of molecular dynamics simulations for the unbound receptor and ligand structures of 17 protein complexes and applied shape-driven rigid body docking to all combinations of representative snapshots. The crossdocking of structure ensembles increased the likelihood of finding near-native solutions. The free ensembles appeared to contain multiple complementary conformations. These were in general not related to the bound structure. We suggest that protein-protein binding follows a three-step mechanism of diffusion, free conformer selection, and refolding. This model combines previously conflicting ideas and is in better agreement with the current data on interaction forces, time scales, and kinetics.

Models, Molecular↗

PubMed1995

Genetic algorithms in molecular recognition and design.

Genetic algorithms provide a novel tool for the investigation of combinatorial optimization problems. A genetic algorithm takes an initial set of possible starting solutions, and iteratively improves them by means of crossover and mutation operators that are related to those involved in Darwinian evolution. This approach is illustrated by reference to applications in molecular modelling, the docking of flexible ligands into protein active sites and de novo ligand design.

Algorithms↗

PubMed2005

Probing the effect of point mutations at protein-protein interfaces with free energy calculations.

We have studied the effect of point mutations of the primary binding residue (P1) at the protein-protein interface in complexes of chymotrypsin and elastase with the third domain of the turkey ovomucoid inhibitor and in trypsin with the bovine pancreatic trypsin inhibitor, using molecular dynamics simulations combined with the linear interaction energy (LIE) approach. A total of 56 mutants have been constructed and docked into their host proteins. The free energy of binding could be reliably calculated for 52 of these mutants that could unambiguously be fitted into the binding sites. We find that the predicted binding free energies are in very good agreement with experimental data with mean unsigned errors between 0.50 and 1.03 kcal/mol. It is also evident that the standard LIE model used to study small drug-like ligand binding to proteins is not suitable for protein-protein interactions. Three different LIE models were therefore tested for each of the series of protein-protein complexes included, and the best models for each system turn out to be very similar. The difference in parameterization between small drug-like compounds and protein point mutations is attributed to the preorganization of the binding surface. Our results clearly demonstrate the potential of free energy calculations for probing the effect of point mutations at protein-protein interfaces and for exploring the principles of specificity of hot spots at the interface.

Animals↗

PubMed2006

The molecular basis for ligand specificity in a mouse olfactory receptor: a network of functionally important residues.

Sequence differences between members of the mouse olfac-tory receptor MOR42 subfamily (MOR42-3 and MOR42-1) are likely to be the basis for variation in ligand binding preference among these receptors. We investigated the specificity of MOR42-3 for a variety of dicarboxylic acids. We used site-directed mutagenesis, guided by homology modeling and ligand docking studies, to locate functionally important residues. Receptors were expressed in Xenopus oocytes and assayed using high throughput electrophysiology. The importance of the Val-113 residue, located deep within the receptor, was analyzed in the context of interhelical interactions. We also screened additional residues predicted to be involved in ligand binding site, based on comparison of ortholog/paralog pairs from the mouse and human olfactory receptor genomes (Man, O., Gilad, Y., and Lancet, D. (2004) Protein Sci. 13, 240-254). A network of 8 residues in transmembrane domains III, V, and VI was identified. These residues form part of the ligand binding pocket of MOR42-3. C12 dicarboxylic acid did not activate the receptor in our functional assay, yet our docking simulations predicted its binding site in MOR42-3. Binding without activation implied that C12 dicarboxylic acid might act as an antagonist. In our functional assay, C12 dicarboxylic acid did indeed act as an antagonist of MOR42-3, in agreement with molecular docking studies. Our results demonstrate a powerful approach based on the synergy between computational predictions and physiological assays.

Animals↗

PubMed2005

Structural basis for a broad but selective ligand spectrum of a mouse olfactory receptor: mapping the odorant-binding site.

The olfactory receptor (OR) superfamily provides a basis for the remarkable ability to recognize and discriminate a large number of odorants. In mice, the superfamily includes approximately 1000 members, and they recognize overlapping sets of odorants with distinct affinities and specificities. To address the molecular basis of odor discrimination by the mammalian OR superfamily, we performed functional analysis on a series of site-directed mutants and performed ligand docking simulation studies to define the odorant-binding site of a mouse OR. Our results indicate that several amino acids in the transmembrane domains formed a ligand-binding pocket. Although other G-protein-coupled receptors (GPCRs) recognize biogenic ligands mainly with ionic or hydrogen bonding interactions, ORs recognize odorants mostly via hydrophobic and van der Waals interactions. This accounts for the broad but selective binding by ORs as well as their relatively low ligand-binding affinities. Furthermore, we succeeded in rational receptor design, inserting point mutations in the odorant-binding site that resulted in predicted changes in ligand specificity and antagonist activity. This ability to rationally design the receptor validated the binding site structure that was deduced with our mutational and ligand docking studies. Such broad and specific sensitivity suggests an evolutionary process during which mutations in the active site led to an enormous number of ORs with a wide range of ligand specificity. The current study reveals the molecular environment of the odorant-binding site, and it further advances the understanding of GPCR pharmacology.

Allosteric Regulation↗

PubMed1998

Mapping the active site of factor Xa by selective inhibitors: an NMR and MD study.

The structure of two selective inhibitors, Ac-Tyr-Ile-Arg-Ile-Pro-NH2 and Ac-(4-Amino-Phe)-(Cyclohexyl-Gly)-Arg-NH2, in the active site of the blood clotting enzyme factor Xa was determined by using transferred nuclear Overhauser effect nuclear magnetic resonance (NMR) spectroscopy. They represent a family of peptidic inhibitors obtained by the screening of a vast combinatorial library. Each structure was first calculated by using standard computational procedures (distance geometry, simulated annealing, energy minimization) and then further refined by systematic search of the conformation of the inhibitor docked in the active site and repeating the simulated annealing and energy minimization. The final structure was optimized by molecular dynamics simulations of the inhibitor-complex in water. The NMR restraints were kept throughout the refinement. The inhibitors assume a compact, very well defined conformation, embedded into the substrate binding site not in the same way as a substrate, blocking thus the catalysis. The model allows to explain the mode of action, affinity, and specificity of the peptides and to map the active site.

Amino Acid Sequence↗

PubMed1999

Homology modeling and substrate binding study of human CYP4A11 enzyme.

Although both bacterial CYP102 (P450BM3) and mammalian CYP4A isozymes share a common function as fatty acid hydroxylases, distinctly different preferred sites of oxidation are observed with the CYP102 performing the usual non-terminal hydroxylation or epoxidation and the CYP4A enzymes performing the unusual and enigmatic terminal hydroxylation. The origin of this unique product specificity in human CYP4A11 has been explored in this work, focusing on possible differences in the binding site architecture of the two isozymes as the cause. To this end, 3D model structures of the human CYP4A11 enzyme were built and compared to the X-ray structure of CYP102. The substrate-binding channel identified in CYP4A11 was found to have a much more sterically restricted active site than that in CYP102 that could cause limited access of long-chain fatty acid to the ferryl oxygen leading to the preferred omega-hydroxylation. Results of docking of a common substrate, lauric acid, into the binding site of both CYP4A11 and CYP102 and molecular dynamics simulations provided additional support for this hypothesis. Specifically, in the CYP4A11-lauric acid simulations, the omega hydrogens were closest to the ferryl oxygen most of the time. By contrast, in the CYP102-lauric acid complex, the substrate could penetrate further into the active site providing access of the non-terminal (omega-1, omega-2) positions to the ferryl oxygen. These results, taken together, have elucidated the origin of the unusual product specificity of CYP4A11 and illustrated the central role of binding site architecture in subtle modulation of function.

Amino Acid Sequence↗

PubMed2006

Mapping residues in the ligand-binding domain of the 5-HT(3) receptor onto d-tubocurarine structure.

The serotonin 5-HT(3) receptor (5-HT(3)R) is a member of the cys-loop ligand-gated ion channel family. We have used the combination of site-directed mutagenesis, homology modeling of the 5-HT(3)R extracellular domain, and ligand docking simulations as a way to map the architecture of the 5-HT(3)R ligand binding domain. Mutation of Phe226 in loop C of the binding site to tyrosine (F226Y) has no effect on the apparent affinity of the competitive antagonist d-tubocurarine (dTC) for the receptor. On the other hand, replacement of Asn128 in loop A of the binding site with alanine (N128A) increases the apparent affinity of dTC by approximately 10-fold. Double-mutant cycle analysis employing a panel of dTC analogs with substitutions at various positions to identify specific points of interactions between the dTC analogs and Asn128 suggests that Asn128 makes a direct interaction with the 2'N of dTC. Molecular modeling of the 5-HT(3)R extracellular domain using the antagonist-bound conformation of the Aplysia californica acetylcholine binding protein as a template followed by ligand docking simulations produces two classes of structures of the 5-HT(3)R/dTC complex; only one of these has the 2'N of dTC positioned at Asn128 and thus is consistent with the data from this study and previously published data. The use of the rigid dTC analogs as "molecular rulers" in conjunction with double-mutant cycle analysis of mutant receptors can allow the spatial mapping of the position of various residues in the ligand-binding site.

Binding Sites↗

PubMed1996

The mu- and delta-opioid pharmacophore conformations of cyclic beta-casomorphin analogues indicate docking of the Phe3 residue to different domains of the opioid receptors.

Cyclic beta-casomorphin analogues with a D-configured amino acid residue in position 2, such as Tyr-c[-Xaa-Phe-Pro-Gly-] and Tyr-c[-Xaa-Phe-D-Pro-Gly-] (Xaa = D-A2bu, D-Orn, D-Lys) were found to bind to the mu-opioid receptor as well as to the delta-opioid receptor, whereas the corresponding L-Xaa2 derivatives are nearly inactive at both. Low-energy conformers of both active and nearly inactive derivatives have been determined in a systematic conformational search or by molecular dynamics simulations using the TRIPOS force field. The obtained conformations were compared with regard to a model for mu-selective opiates developed by Brandt et al. [Drug Des. Discov., 10 (1993) 257]. Superpositions as well as electrostatic, lipophilic and hydrogen bounding similarities with the delta-opioid receptor pharmacophore conformation of t-Hpp-JOM-13 proposed by Mosberg et al. [J. Med. Chem., 37 (1994) 4371, 4384] were used to establish the probable delta-pharmacophoric cyclic beta-casomorphin conformations. These conformations were also compared with a delta-opioid agonist (SNC 80) and the highly potent antagonist naltrindole. These investigations led to a prediction of the mu- and delta-pharmacophore structures for the cyclic beta-casomorphins. Interestingly, for the inactive compounds such conformations could not be detected. The comparison between the mu- and delta-pharmacophore conformations of the cyclic beta-casomorphins demonstrates not only differences in spatial orientation of both aromatic groups, but also in the backbone conformations of the ring part. In particular, the differences on phi2 and psi2 (mu approximately 70 degrees, -80 degrees; delta approximately 165 degrees, 55 degrees) cause a completely different spatial arrangement of the cyclized peptide rings when all compounds are matched with regard to maximal spatial overlap of the tyrosine residue. Assuming that both the mu- and delta-pharmacophore conformations bind with the tyrosine residue in a similar orientation at the same transmembrane domain X of their receptors, the side chain of Phe3 as a second binding site has to dock with different domains.

Computer Graphics↗

PubMed2001

Structural basis of binding of high-affinity ligands to protein kinase C: prediction of the binding modes through a new molecular dynamics method and evaluation by site-directed mutagenesis.

The structural basis of protein kinase C (PKC) binding to several classes of high-affinity ligands has been investigated through complementary computational and experimental methods. Employing a recently developed q-jumping molecular dynamics (MD) simulation method, which allows us to consider the flexibility of both the ligands and the receptor in docking studies, we predicted the binding models of phorbol-13-acetate, phorbol-12,13-dibutyrate (PDBu), indolactam V (ILV), ingenol-3-benzoate, and thymeleatoxin to PKC. The "predicted" binding model for phorbol-13-acetate is virtually identical to the experimentally determined binding model for this ligand. The predicted binding model for PDBU is the same as that for phorbol-13-acetate in terms of the hydrogen-bonding network and hydrophobic contacts. The predicted binding model for ILV is the same as that obtained in a previous docking study using a Monte Carlo method and is consistent with the structure-activity relationships for this class of ligands. Together with the X-ray structure of phorbol-13-acetate in complex with PKCdelta C1b, the predicted binding models of PDBu, ILV, ingenol-3-benzoate, and thymeleatoxin in complex with PKC showed that the binding of these ligands to PKC is governed by a combination of several highly specific and optimal hydrogen bonds and hydrophobic contacts. However, the hydrogen-bonding network for each class of ligand is somewhat different and the number of hydrogen bonds formed between PKC and these ligands has no correlation with their binding affinities. To provide a direct and quantitative assessment of the contributions of several conserved residues around the binding site to PKC-ligand binding, we have made 11 mutations and measured the binding affinities of the high-affinity PKC ligands to these mutants. The results obtained through site-directed mutagenic analysis support our predicted binding models for these ligands and provide new insights into PKC-ligand binding. Although all the ligands have high affinity for the wild-type PKCdelta C1b, our site-directed mutagenic results showed that ILV is the ligand most sensitive to structural perturbations of the binding site while ingenol-3-benzoate is the least sensitive among the four classes of ligands examined here. Finally, we have employed conventional MD simulations to investigate the structural perturbations caused by each mutation to further examine the role played by each individual residue in PKC-ligand binding. MD simulations revealed that several mutations, including Pro11 --> Gly, Leu21 --> Gly, Leu24 --> Gly, and Gln27 --> Gly, cause a rather large conformational alteration to the PKC binding site and, in some cases, to the overall structure of the protein. The complete abolishment or the significant reduction in PKC-ligand binding observed for these mutants thus reflects the loss of certain direct contacts between the side chain of the mutated residue in PKC and ligands as well as the large conformational alteration to the binding site caused by the mutation.

Crystallography, X-Ray↗

PubMed2003

Protein-protein docking with a reduced protein model accounting for side-chain flexibility.

A protein-protein docking approach has been developed based on a reduced protein representation with up to three pseudo atoms per amino acid residue. Docking is performed by energy minimization in rotational and translational degrees of freedom. The reduced protein representation allows an efficient search for docking minima on the protein surfaces within. During docking, an effective energy function between pseudo atoms has been used based on amino acid size and physico-chemical character. Energy minimization of protein test complexes in the reduced representation results in geometries close to experiment with backbone root mean square deviations (RMSDs) of approximately 1 to 3 A for the mobile protein partner from the experimental geometry. For most test cases, the energy-minimized experimental structure scores among the top five energy minima in systematic docking studies when using both partners in their bound conformations. To account for side-chain conformational changes in case of using unbound protein conformations, a multicopy approach has been used to select the most favorable side-chain conformation during the docking process. The multicopy approach significantly improves the docking performance, using unbound (apo) binding partners without a significant increase in computer time. For most docking test systems using unbound partners, and without accounting for any information about the known binding geometry, a solution within approximately 2 to 3.5 A RMSD of the full mobile partner from the experimental geometry was found among the 40 top-scoring complexes. The approach could be extended to include protein loop flexibility, and might also be useful for docking of modeled protein structures.

Algorithms↗

PubMed2006

Conformational and orientational guidance of the analgesic dipeptide kyotorphin induced by lipidic membranes: putative correlation toward receptor docking.

The analgesic dipeptide kyotorphin (L-Tyr-L-Arg) and an acylated kyotorphin derivative were studied by a combination of theoretical (molecular dynamics simulation and quantum mechanics methods) and experimental (fluorescence and infrared spectroscopies) approaches both in solution and in model systems of membranes. At biological pH the peptides have a neutral net charge. Nevertheless, their phenolic rings interact with phospholipid molecules (partition coefficient varies from 6 x 10(2) to 2 x 10(4), depending on the lipid and pH used) despite being exposed to the aqueous bulk medium. The lowest energy transition dipole moment is displaced from the normal to the lipid bilayer by 20 degrees on average. The observed extensive interaction, pK(a), precise location, and well-defined orientation in membranes combined with the ability to discriminate rigid raftlike membrane domains suggest that kyotorphin meets the structural constraints needed for receptor-ligand interaction. The acylated kyotorphin derivative mimics kyotorphin properties and represents a promising way for entrapment in a drug carrier and transport across the blood-brain barrier.

Blood-Brain Barrier↗

PubMed2006

DynaPred: a structure and sequence based method for the prediction of MHC class I binding peptide sequences and conformations.

MOTIVATION: The binding of endogenous antigenic peptides to MHC class I molecules is an important step during the immunologic response of a host against a pathogen. Thus, various sequence- and structure-based prediction methods have been proposed for this purpose. The sequence-based methods are computationally efficient, but are hampered by the need of sufficient experimental data and do not provide a structural interpretation of their results. The structural methods are data-independent, but are quite time-consuming and thus not suited for screening of whole genomes. Here, we present a new method, which performs sequence-based prediction by incorporating information obtained from molecular modeling. This allows us to perform large databases screening and to provide structural information of the results. RESULTS: We developed a SVM-trained, quantitative matrix-based method for the prediction of MHC class I binding peptides, in which the features of the scoring matrix are energy terms retrieved from molecular dynamics simulations. At the same time we used the equilibrated structures obtained from the same simulations in a simple and efficient docking procedure. Our method consists of two steps: First, we predict potential binders from sequence data alone and second, we construct protein-peptide complexes for the predicted binders. So far, we tested our approach on the HLA-A0201 allele. We constructed two prediction models, using local, position-dependent (DynaPred(POS)) and global, position-independent (DynaPred) features. The former model outperformed the two sequence-based methods used in our evaluation; the latter shows a much higher generalizability towards other alleles than the position-dependent models. The constructed peptide structures can be refined within seconds to structures with an average backbone RMSD of 1.53 A from the corresponding experimental structures.

Algorithms↗

PubMed2005

Exploration of a binding mode of indole amide analogues as potent histone deacetylase inhibitors and 3D-QSAR analyses.

Docking simulations and three-dimensional quantitative structure-activity relationship (3D-QSAR) analyses were conducted on a series of indole amide analogues as potent histone deacetylase inhibitors. The studies include comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Selected ligands were docked into the active site of human HDAC1. Based on the docking results, a novel binding mode of indole amide analogues in the human HDAC1 catalytic core is presented, and enzyme/inhibitor interactions are discussed. The indole amide group is located in the open pocket, and anchored to the protein through a pair of hydrogen bonds with Asp99 O-atom and amide NH group on ligand. Based on the binding mode, predictive 3D-QSAR models were established, which had conventional r2 and cross-validated coefficient values (r(cv)2) up to 0.982 and 0.601 for CoMFA and 0.954 and 0.598 for CoMSIA, respectively. A comparison of the 3D-QSAR field contributions with the structural features of the binding site showed good correlation between the two analyses. The results of 3D-QSAR and docking studies validate each other and provided insight into the structural requirements for activity of this class of molecules as HDAC inhibitors. The CoMFA and CoMSIA PLS contour maps and MOLCAD-generated active site electrostatic, lipophilicity, and hydrogen-bonding potential surface maps, as well as the docking studies, provided good insights into inhibitor-HDAC interactions at the molecular level. Based on these results, novel molecules with improved activity can be designed.

Amides↗

PubMed2004

Optimization and dynamics of protein-protein complexes using B-splines.

A moving-grid approach for optimization and dynamics of protein-protein complexes is introduced, which utilizes cubic B-spline interpolation for rapid energy and force evaluation. The method allows for the efficient use of full electrostatic potentials joined smoothly to multipoles at long distance so that multiprotein simulation is possible. Using a recently published benchmark of 58 protein complexes, we examine the performance and quality of the grid approximation, refining cocrystallized complexes to within 0.68 A RMSD of interface atoms, close to the optimum 0.63 A produced by the underlying MMFF94 force field. We quantify the theoretical statistical advantage of using minimization in a stochastic search in the case of two rigid bodies, and contrast it with the underlying cost of conjugate gradient minimization using B-splines. The volumes of conjugate gradient minimization basins of attraction in cocrystallized systems are generally orders of magnitude larger than well volumes based on energy thresholds needed to discriminate native from nonnative states; nonetheless, computational cost is significant. Molecular dynamics using B-splines is doubly efficient due to the combined advantages of rapid force evaluation and large simulation step sizes. Large basins localized around the native state and other possible binding sites are identifiable during simulations of protein-protein motion. In addition to providing increased modeling detail, B-splines offer new algorithmic possibilities that should be valuable in refining docking candidates and studying global complex behavior.

Algorithms↗

PubMed2005

Three-dimensional models of histamine H3 receptor antagonist complexes and their pharmacophore.

Molecular modeling was used to analyze the binding mode and activities of histamine H3 receptor antagonists. A model of the H3 receptor was constructed through homology modeling methods based on the crystal structure of bovine rhodopsin. Known H3 antagonists were interactively docked into the putative antagonist binding pocket and the resultant model was subjected to molecular mechanics energy minimization and molecular dynamics simulations which included a continuum model of the lipid bilayer and intra- and extracellular aqueous environments surrounding the transmembrane helices. The transmembrane helices stayed well embedded in the dielectric slab representing the lipid bilayer and the intra- and extracellular loops remain situated in the aqueous solvent region of the model during molecular dynamics simulations of up to 200 ps in duration. A pharmacophore model was calculated by mapping the features common to three active compounds three-dimensionally in space. The 3D pharmacophore model complements our atomistic receptor/ligand modeling. The H3 antagonist pharmacophore consists of two protonation sites (i.e. basic centers) connected by a central aromatic ring or hydrophobic region. These two basic sites can simultaneously interact with Asp 114 (3.32) in helix III and a Glu 206 (5.46) in helix V which are believed to be the key residues that histamine interacts with to stabilize the receptor in the active state. The interaction with Glu 206 is consistent with the enhanced activity resulting from the additional basic site. In addition to these two salt bridging interactions, the central region of these antagonists contains a lipophilic group, usually an aromatic ring, that is found to interact with several nearby hydrophobic side chains. The picture of antagonist binding provided by these models is consistent with earlier pharmacophore models for H3 antagonists with some exceptions.

Amino Acid Sequence↗

PubMed2003

Synthesis, activity, and molecular modeling studies of novel human aldose reductase inhibitors based on a marine natural product.

Aldose reductase (ALR2) has been implicated in the etiology of diabetic complications, including blindness. Because of the limited number of currently available drugs for the prevention of these long-term complications, the discovery of new ALR2 inhibitors appears highly desirable. In this study, a polybrominated diphenyl ether (1) naturally occurring in a marine sponge was found to inhibit recombinant human ALR2 with an IC(50) of 6.4 microM. A series of polyhalogenated analogues that were synthesized and tested in vitro to explore the structure-activity relationships displayed various degrees of inhibitory activity. The most active compounds were also capable of preventing sorbitol accumulation inside human retinal cells. In this cell-based assay, the most potent synthesized analogue (16) showed a 17-fold increase in inhibitory activity compared to that of sorbinil (IC(50) = 0.24 vs 4 microM). A molecular representation of human ALR2 in complex with the natural product was built using homology modeling, automated docking, and energy refinement methods. AMBER parameters for the halogen atoms were derived and calibrated using condensed phase molecular dynamics simulations of fluorobenzene, chlorobenzene, and bromobenzene. Inhibitor binding is proposed to cause a conformational change similar to that recently reported for zenarestat. A free energy perturbation thermodynamic cycle allowed us to assess the importance of a crucial bromine atom that distinguishes the active lead compound from a much less active close natural analogue. Remarkably, the spatial location of this bromine atom is equivalent to that occupied by the only bromine atom present in zenarestat.

Aldehyde Reductase↗

PubMed1999

Modeling helix-turn-helix protein-induced DNA bending with knowledge-based distance restraints.

A crucial element of many gene functions is protein-induced DNA bending. Computer-generated models of such bending have generally been derived by using a presumed bending angle for DNA. Here we describe a knowledge-based docking strategy for modeling the structure of bent DNA recognized by a major groove-inserting alpha-helix of proteins with a helix-turn-helix (HTH) motif. The method encompasses a series of molecular mechanics and dynamics simulations and incorporates two experimentally derived distance restraints: one between the recognition helix and DNA, the other between respective sites of protein and DNA involved in chemical modification-enabled nuclease scissions. During simulation, a DNA initially placed at a distance was "steered" by these restraints to dock with the binding protein and bends. Three prototype systems of dimerized HTH DNA binding were examined: the catabolite gene activator protein (CAP), the phage 434 repressor (Rep), and the factor for inversion stimulation (Fis). For CAP-DNA and Rep-DNA, the root mean square differences between model and x-ray structures in nonhydrogen atoms of the DNA core domain were 2.5 A and 1.6 A, respectively. An experimental structure of Fis-DNA is not yet available, but the predicted asymmetrical bending and the bending angle agree with results from a recent biochemical analysis.

Base Sequence↗