PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Models, Molecular”

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 145 records · Page 8Linked to original sources

A critical assessment of comparative molecular modeling of tertiary structures of proteins.

In spite of the tremendous increase in the rate at which protein structures are being determined, there is still an enormous gap between the numbers of known DNA-derived sequences and the numbers of three-dimensional structures. In order to shed light on the biological functions of the molecules, researchers often resort to comparative molecular modeling. Earlier work has shown that when the sequence alignment is in error, then the comparative model is guaranteed to be wrong. In addition, loops, the sites of insertions and deletions in families of homologous proteins, are exceedingly difficult to model. Thus, many of the current problems in comparative molecular modeling are minor versions of the global protein folding problem. In order to assess objectively the current state of comparative molecular modeling, 13 groups submitted blind predictions of seven different proteins of undisclosed tertiary structure. This assessment shows that where sequence identity between the target and the template structure is high (> 70%), comparative molecular modeling is highly successful. On the other hand, automated modeling techniques and sophisticated energy minimization methods fail to improve upon the starting structures when the sequence identity is low (approximately 30%). Based on these results it appears that insertions and deletions are still major problems. Successfully deducing the correct sequence alignment when the local similarity is low is still difficult. We suggest some minimal testing of submitted coordinates that should be required of authors before papers on comparative molecular modeling are accepted for publication in journals.

Amino Acid Sequence↗

[Methods of computer-assisted molecular modeling].

We present an overview on modern computer methods of molecular modelling. After treating three main steps of drug evaluation, namely target identification, lead identification and lead optimisation, we shortly discuss molecular graphics, molecular mechanics, molecular orbital and molecular dynamics methods. These are suitable for the more-or-less adequate modelling of real molecular processes both at the microscopic and the macroscopic levels. Molecular graphics provides beautiful pictures for the specialist that allow inspection and manipulation of detailed molecular models. An especially useful tool for the visualisation of molecular entities is the display of various properties on the molecular surface that allows rapid recognition of important relationships. Molecular mechanics is able to predict properties (e.g. geometric parameters, conformer stability) of several classes of molecules with an accuracy close to the experimental one, therefore it plays an important role in complementing molecular graphics. The performance of molecular orbital methods increased considerably in the last decade thus we can calculate parameters for isolated or interacting molecules that are not easily amenable to experiment (e.g. structure and energetics of unstable species or activation energies of elementary processes). Computer simulation methods provide a link between gas-phase models of microscopic structures or processes and macroscopic properties or events that may be derived from the former. Thus, it became possible to apply computerised methods for an adequate simulation of important events, like chemical and biochemical reactions, drug-target interactions, drug delivery and the similar that determine drug action. It is stressed that the hardware and software for computer-aided molecular modelling may not be absent from the arsenal of a drug designer.

Computer Graphics↗

Approaches to molecular modeling studies and specific application to serotonin ligands and receptors.

Molecular modeling studies are useful in as much as they may allow us to understand the activity and selectivity of currently existing agents, and, furthermore, may aid in the design of completely novel therapeutic agents. There are two basic modeling strategies: the ligand-ligand approach and the ligand-receptor approach. Both approaches possess certain inherent advantages and disadvantages and, in addition, make certain assumptions about the agents and/or receptors being investigated. Keeping with the spirit of this minisymposium, we describe these two approaches, their general usefulness, and their limitations. Using serotonin (5-HT) receptors as a focal point, we review and provide novel examples of molecular modeling studies involving both strategies. Presented for the first time are examples of ligand-receptor models to account for the binding of serotonergic agents at 5-HT2 and 5-HT1C receptors.

Amino Acid Sequence↗

Three-dimensional molecular modeling of bovine caseins: a refined, energy-minimized kappa-casein structure.

A refined three-dimensional molecular model of kappa-casein has been produced using energy minimization techniques and a Kollman force field on a previously reported predicted three-dimensional structure. This initial model was constructed via molecular modeling techniques from sequence-based secondary structural prediction algorithms. Both the initial and refined structures agreed with global secondary structure analysis from vibration spectroscopy. The refined structure contained many of the features of the initial model, including two sets of antiparallel beta-sheet structures containing predominantly hydrophobic side chains, which could form interaction sites with alpha s1-casein. Two types of energy-minimized dimer and tetramer models are presented: 1) using Cys as potential intermolecular disulfide binding sites and 2) using the two sheets as possible hydrophobic self-association sites, without Cys interactions. All structures yielded good stabilization energies and are in agreement with chemical, biochemical, and physical chemical results obtained for kappa-casein.

Amino Acid Sequence↗

Identity and geometry of a base triple in 16S rRNA determined by comparative sequence analysis and molecular modeling.

Comparative sequence analysis complements experimental methods for the determination of RNA three-dimensional structure. This approach is based on the concept that different sequences within the same gene family form similar higher-order structures. The large number of rRNA sequences with sufficient variation, along with improved covariation algorithms, are providing us with the opportunity to identify new base triples in 16S rRNA. The three-dimensional conformations for one of our strongest candidates involving U121 (C124:G237) and/or U121 (U125:A236) (Escherichia coli sequence and numbering) are analyzed here with different molecular modeling tools. Molecular modeling shows that U121 interacts with C124 in the U121 (C124:G237) base triple. This arrangement maintains isomorphic structures for the three most frequent sequence motifs (approximately 93% of known bacterial and archaeal sequences), is consistent with chemical reactivity of U121 in E. coli ribosomes, and is geometrically favorable. Further, the restricted set of observed canonical (GU, AU, GC) base-pair types at positions 124:237 and 125:236 is consistent with the fact that the canonical base-pair sets (for both base pairs) that are not observed in nature prevent the formation of the 121 (124:237) base triple. The analysis described here serves as a general scheme for the prediction of specific secondary and tertiary structure base pairing where there is a network of correlated base changes.

Base Pairing↗

Selection of putative immunogenic peptides by molecular modelling of the urease of Helicobacter pylori.

Computational methods such as molecular modelling are becoming an increasingly useful means of rationalising experimental data and creating a hypothesis that can suggest new experiments. In this report we discuss the application of molecular modelling methods to aid the selection of feasible peptide epitopes of the urease enzyme from Helicobacter pylori, an important vaccine candidate. Surface exposure was chosen as a criterion for the selection of three peptides which each had different levels of accessibility according to the 3D model. Antibodies raised against these peptides were analysed for their immunoreactivity with the holo enzyme. Only one anti-peptide antibody showed good reactivity with the urease. Our findings emphasise that surface exposure of peptide is not the only important criterion for the selection of immunogenic peptides.

Enzyme-Linked Immunosorbent Assay↗

Refinement of a molecular model for lamprey hemoglobin from Petromyzon marinus.

A molecular model for the protein and ambient solvent of the complex of cyanide with methemoglobin V from the sea lamprey Petromyzon marinus yields an R-factor of 0.142 against X-ray diffraction data to 2.0 A resolution. The root-mean-square discrepancies from ideal bond length and angle are, respectively, 0.014 A and 1.5 degrees. Atoms that belong to planar groups deviate by 0.012 A from planes determined by a least-squares procedure. The average standard deviation for chiral volumes, peptide torsion angle and torsion angles of side-chains are 0.150 A3, 2.0 degrees and 19.4 degrees, respectively. The root-mean-square variation in the thermal parameters of bonded atoms of the polypeptide backbone is 1.21 A2; the variation in thermal parameters for side-chain atoms is 2.13 A2. The model includes multiple conformations for 11 side-chains of the 149 amino acid residues of the protein. We identify 231 locations as sites of water molecules in full or partial occupancy. The sum of occupancy factors for these sites is approximately 154, representing 28% of the 550 molecules of water within the crystallographic asymmetric unit. The environment of the heme in the cyanide complex of lamprey methemoglobin resembles the deoxy state of the mammalian tetramer. In particular, the bond between atom NE2 of the proximal histidine and the Fe lies 5.1 degrees from the normal of the heme plane. In deoxy- and carbonmonoxyhemoglobins, the deviations from the normal to the heme plane are 7 to 8 degrees and 1 degree, respectively. Furthermore, the inequality in the distance of atom CD2 of the proximal histidine from the pyrrole nitrogen of ring-C of the heme (distance = 3.29 A) and CE1 from the pyrrole nitrogen of ring-A (distance = 3.06 A) is characteristic of deoxyhemoglobin, not carbonmonoxyhemoglobin, where these distances are equal. Finally, a hydrogen bond exists between carbonyl 111 and the hydroxyl of tyrosine 149. The corresponding hydrogen link in the mammalian tetramer is central to the T to R state transition and is present in deoxyhemoglobin but absent in carbonmonoxyhemoglobin. We suggest that the low affinity of oxygen for lamprey hemoglobin may be a consequence of these T-state geometries.

Amino Acid Sequence↗

Abasic sites in duplex DNA: molecular modeling of sequence-dependent effects on conformation.

Molecular modeling calculations using JUnction Minimization of Nucleic Acids (JUMNA) have been used to study sequence effects on the conformation of abasic sites within duplex DNA. We have considered lesions leading to all possible unpaired bases (X), adenine, guanine, cytosine, or thymine contained within two distinct sequence contexts, CXC and GXG. Calculations were carried out on DNA 11-mers using extensive conformational search techniques to locate the most stable abasic conformations and using Poisson-Boltzmann corrected electrostatics to account for solvation effects. The results, which are in very good agreement with available experimental data, point to strong sequence effects on both the position of the unpaired base (intra or extrahelical) and on the overall curvature induced by the abasic lesion. For CXC, unpaired purines are found to lie within the helix, while unpaired pyrimidines are either extrahelical or in equilibrium between the intra and extrahelical forms. For GXG, all unpaired bases lead to intrahelical forms, but with marked, sequence-dependent differences in induced curvature.

Base Pairing↗

BALLView: a tool for research and education in molecular modeling.

We present BALLView, a molecular viewer and modeling tool. It combines state-of-the-art visualization capabilities with powerful modeling functionality including implementations of force field methods and continuum electrostatics models. BALLView is a versatile and extensible tool for research in structural bioinformatics and molecular modeling. Furthermore, the convenient and intuitive graphical user interface offers novice users direct access to the full functionality, rendering it ideal for teaching. Through an interface to the object-oriented scripting language Python it is easily extensible.

Computer Graphics↗

Review of the molecular modelling studies of the cytochrome P-450 estrogen synthetase enzyme, aromatase.

The biosynthesis of the family of female hormones, the estrogens, is mediated by the cytochrome P-450 based enzyme Aromatase (AR), an enzyme which has been the focus of extensive research for some time and has resulted in the synthesis of numerous inhibitory compounds with varying structural properties. The crystal structure of this membrane-bound enzyme has not been determined, although that of the bacterial P-450 camphor hydroxylase enzyme has, and the latter has been used for homology modelling of AR and other related enzymes. Several workers have attempted to obtain an approximate model of the active site of AR through the use of various molecular modelling techniques using inhibitors of this enzyme, as well as a combination of experiments involving active site mutagenesis of specific areas of the P-450 backbone together with molecular modelling. In this report these studies are reviewed. In summary, although the crystal structure of AR has yet to be determined, the use of molecular modelling studies has provided useful information in the design and synthesis of its inhibitors.

Aromatase↗

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics.

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.

Amino Acid Sequence↗

Molecular modeling of mammalian cytochrome P450s.

The cytochrome P450 enzymes collectively metabolize a wide range of xenobiotic and endogenous compounds. The broad substrate specificity of this superfamily derives from the multiplicity of P450s whose unique substrate specificity profiles reflect underlying differences in primary sequence. Experimental structures of P450s, where available, have provided great insight into the basis of substrate recognition. However, for those mammalian P450s whose structures have not been determined, homology modeling has become an increasingly important tool for understanding substrate specificity and mechanism. P450 modeling is often a challenging task, owing to the rather low sequence identity between target and template proteins. Although mammalian P450 models have previously been based on bacterial P450 structures, the recent advent of mammalian P450 structures holds great potential for generating more accurate homology models. Consequently, the substrate recognition properties of several mammalian P450s have been rationalized using the predicted substrate binding site of recently developed models. This review summarizes the major concepts and current approaches of molecular modeling of P450s.

Algorithms↗

Molecular mechanism of DApd/DXG against zidovudine- and lamivudine- drug resistant mutants: a molecular modelling approach.

In order to understand molecular mechanism of antiviral drug resistance of HIV-1 reverse transcriptase (RT) as well as potent antiviral activity of 2,6-diaminopurine dioxolane (DAPD) [prodrug of (-)-beta-D-dioxolane guanine (DXG)] against drug-resistant RTs, molecular modelling studies of three structurally distinct nucleoside RT inhibitor (NRTI)-triphosphates (TP) [zidovudine (AZT)-TP, lamivudine (3TC)-TP and DXG-TP] complexed with the wild-type (WT) and mutated RT were conducted. The computational analyses indicated that the antiviral activity and the calculated relative binding energy of the RT inhibitor triphosphates can be correlated, and the minimized structures gave information on the molecular mechanism of drug resistance conferred by mutations. The interactions between the NRTI-TP and adjacent amino acid residues (Lys65, Lys70, Arg72, Tyr115 and/or Gln151) played important roles in stabilizing the enzyme-inhibitor complex. Particularly, Arg72 was found to stabilize the dioxolane and oxathiolane sugar moiety through hydrogen bonding, which was responsible for favourable binding affinity of DXG-TP to AZT- as well as 3TC-resistant mutants. The conformational changes in these amino acid residues caused by mutation always affected the changes in the tertiary structures of enzyme-inhibitor complexes through either closing or opening the gap between the fingers and palm domains. The enzyme-inhibitor complexes with good binding affinity showed tight binding modes by closing the gap between the two domains, whereas weak inhibitors gave open and loose complexes.

Anti-HIV Agents↗

Using surface envelopes for discrimination of molecular models.

Shape information about macromolecules is increasingly available but is difficult to use in modeling efforts. We demonstrate that shape information alone can often distinguish structural models of biological macromolecules. By using a data structure called a surface envelope (SE) to represent the shape of the molecule, we propose a method that generates a fitness score for the shape of a particular molecular model. This score correlates well with root mean squared deviation (RMSD) of the model to the known test structures and can be used to filter models in decoy sets. The scoring method requires both alignment of the model to the SE in three-dimensional space and assessment of the degree to which atoms in the model fill the SE. Alignment combines a hybrid algorithm using principal components and a previously published iterated closest point algorithm. We test our method against models generated from random atom perturbation from crystal structures, published decoy sets used in structure prediction, and models created from the trajectories of atoms in molecular modeling runs. We also test our alignment algorithm against experimental electron microscopic data from rice dwarf virus. The alignment performance is reliable, and we show a high correlation between model RMSD and score function. This correlation is stronger for molecular models with greater oblong character (as measured by the ratio of largest to smallest principal component).

Algorithms↗

Molecular modeling and dynamics of the sodium channel inactivation gate.

The intracellular linker L(III-IV) of voltage-gated sodium channels is known to be involved in their mechanism of inactivation. Its primary sequence is well conserved in sodium channels from different tissues and species. However, the role of charged residues in this region, first thought to play an important role in inactivation, has not been well identified, whereas the IFM triad (I1488-M1490) has been characterized as the crucial element for inactivation. In this work, we constructed theoretical models and performed molecular dynamics simulations, exploring the role of L(III-IV)-charged residues in the presence of a polar/nonpolar planar interface represented by a dielectric discontinuity. From structural predictions, two alpha-helical segments are proposed. Moreover, from dynamics simulations, a time-conserved motif is detected and shown to play a relevant role in guiding the inactivation particle toward its receptor site.

Amino Acid Sequence↗

Strategic approaches to drug design. I. An integrated software framework for molecular modelling.

An integrated molecular graphics and computational chemistry framework is described which has been designed primarily to handle small molecules of up to 300 atoms. The system provides a means of integrating software from any source into a single framework. It is split into two functional subsystems. The first subsystem, called COSMIC, runs on low-cost, serial-linked colour graphics terminals and allows the user to prepare and examine structural data and to submit them for extensive computational chemistry. Links also allow access to databases, other modelling systems and user-written modules. Much of the output from COSMIC cannot be examined with low level graphics. A second subsystem, called ASTRAL, has been developed for the high-resolution Evans & Sutherland PS300 colour graphics terminal and is designed to manipulate complex display structures. The COSMIC minimisers, geometry investigators, molecular orbital displays, electrostatic isopotential generators and various interfaces and utilities are described.

Computer Graphics↗

Molecular modeling of calmodulin: a comparison with crystallographic data.

Two methods of side-chain placement on a modeled protein have been examined. Two molecular models of calmodulin were constructed that differ in the treatment of side chains prior to optimization of the molecule. A virtual bond analysis program developed by Purisima and Scheraga was used to determine the backbone conformation based on 2.2 angstroms resolution C alpha coordinates for the molecules. In the first model, side chains were initially constructed in an extended conformation. In the second model, a conformational grid search technique was employed. Calcium ions were treated explicitly during energy optimization using CHARMM. The models are compared to a recently published refined crystal structure of calmodulin. The results indicate that the initial choices for side-chains, but also significant effects on the main-chain conformation and supersecondary structure. The conformational differences are discussed. Analysis of these and other methods makes possible the formulation of a methodology for more appropriate side-chain placement in modeled proteins.

Amino Acid Sequence↗

Conformation of receptor-associated PGI2: an investigation by molecular modeling.

To elucidate the conformation of receptor-associated prostacyclin (PGI2), we first performed structure-activity correlation analysis of over 200 PGI2 analogues and derived from this analysis several crucial features pertaining to structural requirements for PGI2 activity [Ah-lim Tsai and Kenneth K. Wu, Eicosanoids, 2 (1989) 131-143]. These structural features proved to be useful guidelines for selecting 'model molecules' for further investigations by molecular mechanics. By properly selecting four analogues with either rigid or uniquely oriented alpha-side chain structure for geometric fitting, we succeeded in maximally minimizing the degree of freedom of the carboxylate terminus of PGI2. We were able to define the spatial relationship among the four critical functional groups, i.e., C1-COOH, C6a-O, C11-OH and C15-OH. More information is needed, however, to define the geometry of the omega-side chain, particularly for the moiety beyond C15. Nevertheless, results from structure-activity correlation analysis and molecular modeling provide useful information regarding the conformation of receptor-associated PGI2, which assumes an 'elongated' conformation instead of the traditional 'hairpin' structure.

Binding Sites↗