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).     

Identification of natural compound inhibitors against PfDXR: A hybrid structure-based molecular modeling approach and molecular dynamics simulation studies

In the present contribution, multicomplex-based pharmacophore studies were carried out on the structural proteome of Plasmodium falciparum 1-deoxy-D -xylulose-5-phosphate reductoisomerase. Among the constructed models, a representative model with complementary features, accountable for the inhibition was used as a primary filter for the screening of database molecules. Auxiliary evaluations of the screened molecules were performed via drug-likeness and molecular docking studies. Subsequently, the stability of the docked inhibitors was envisioned by molecular dynamics simulations, principle component analysis, and molecular mechanics-Poisson-Boltzmann surface area-based free binding energy calculations. The stability assessment of the hits was done by comparing with the reference (beta-substituted fosmidomycin analog, LC5) to prioritize more potent candidates. All the complexes showed stable dynamic behavior while three of them displayed higher binding free energy compared with the reference. The work resulted in the identification of the compounds with diverse scaffolds, which could be used as initial leads for the design of novel PfDXR inhibitors.     

Structural analysis and molecular modeling of two antitrichosanthin IgE clones from phage antibody library

Recently we constructed a murine IgE phage surface display library and screened out two IgE (Fab) clones with specific binding activity to Trichosanthin (TCS).In this work,the Vε and Vκ genes of the two clones were sequenced and their putative germline gene usages were studied.On the basis of the known 3D structure of Trichosanthin and antibody,molecular modeling was carried out to study the antigen-antibody interaction.The possible antigenic determinant sites on the surface of TCS recognized by both the clones were analyzed,and the reaction forces between TCS and two Fab fragments were also analyzed respectively.     

Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling

Despite GPCRs sharing a common seven helix bundle, analysis of the diverse crystallographic structures available reveal specific features that might be relevant for ligand design. Despite the number of crystallographic structures of GPCRs steadily increasing, there are still challenges that hamper the availability of new structures. In the absence of a crystallographic structure, homology modeling remains one of the important techniques for constructing 3D models of proteins. In the present study we investigated the use of molecular dynamics simulations for the refinement of GPCRs models constructed by homology modeling. Specifically, we investigated the relevance of template selection, ligand inclusion as well as the length of the simulation on the quality of the GPCRs models constructed. For this purpose we chose the crystallographic structure of the rat muscarinic M3 receptor as reference and constructed diverse atomistic models by homology modeling, using different templates. Specifically, templates used in the present work include the human muscarinic M2; the more distant human histamine H1 and the even more distant bovine rhodopsin as shown in the GPCRs phylogenetic tree. We also investigated the use or not of a ligand in the refinement process. Hence, we conducted the refinement process of the M3 model using the M2 muscarinic as template with tiotropium or NMS docked in the orthosteric site and compared with the results obtained with a model refined without any ligand bound.     

Binding properties of palmatine to DNA: spectroscopic and molecular modeling investigations

Palmatine, an isoquinoline alkaloid, is an important medicinal herbal extract with diverse pharmacological and biological properties. In this work, spectroscopic and molecular modeling approaches were employed to reveal the interaction between palmatine and DNA isolated from herring sperm. The absorption spectra and iodide quenching results indicated that groove binding was the main binding mode of palmatine to DNA. Fluorescence studies indicated that the binding constant (K) of palmatine and DNA was ~ 10⁴ L·mol^?1. The associated thermodynamic parameters, ΔG, ΔH, and ΔS, indicated that hydrogen bonds and van der Waals forces played major roles in the interaction. The effects of chemical denaturant, thermal denaturation and pH on the interaction were investigated and provided further support for the groove binding mode. In addition to experimental approaches, molecular modeling was conducted to verify binding pattern of palmatine–DNA. Copyright © 2015 John Wiley & Sons, Ltd.     

Pharmacophore modeling,molecular docking and molecular dynamics studies on natural products database to discover novel skeleton as non-purine xanthine oxidase inhibitors

Gout is a common inflammatory arthritis caused by the deposition of urate crystals within joints. It is increasingly in prevalence during the past few decades as shown by the epidemiological survey results. Xanthine oxidase (XO) is a key enzyme to transfer hypoxanthine and xanthine to uric acid, whose overproduction leads to gout. Therefore, inhibiting the activity of xanthine oxidase is an important way to reduce the production of urate. In the study, in order to identify the potential natural products targeting XO, pharmacophore modeling was employed to filter databases. Here, two methods, pharmacophore based on ligand and pharmacophore based on receptor-ligand, were constructed by Discovery Studio. Then GOLD was used to refine the potential compounds with higher fitness scores. Finally, molecular docking and dynamics simulations were employed to analyze the interactions between compounds and protein. The best hypothesis was set as a 3D query to screen database, returning 785 and 297 compounds respectively. A merged set of the above 1082 molecules was subjected to molecular docking, which returned 144 hits with high-fitness scores. These molecules were clustered in four main kinds depending on different backbones. What is more, molecular docking showed that the representative compounds established key interactions with the amino acid residues in the protein, and the RMSD and RMSF of molecular dynamics results showed that these compounds can stabilize the protein. The information represented in the study confirmed previous reports. And it may assist to discover and design new backbones as potential XO inhibitors based on natural products.     

An eco-morphological explanation of individual variability in the shape of the fish otolith: comparison of the otolith of Hoplostethus atlanticus with other species by depth

Variation in otolith shape (otolith polymorphism) in Hoplostethus atlanticus, Hoplostethus mediterraneus, Paratrichichthys trailli, Pagrus major and Trachurus murphyi , quantified using principal components analysis based on Fourier transform decompositions of outlines of otolith shape, was particularly high among the samples of H. atlanticus from the North Atlantic, New Zealand, Australia and Namibia. The scatter was uniform, however, and did not show any significant differences among regions. The implication drawn from the high variability in otolith shape of H. atlanticus was that otolith shape polymorphism was maintained by some form of balancing selection across many small local environments which may result in k-selection with consequent poor response by H. atlanticus to maximum sustainable yield harvesting strategies. The variation in otolith shape defined by the otolith morphospace of the five species that were measured, showed a decreasing trend in scatter (i.e. decreasing complexity of shape) proceeding from the species with the deepest habitat ( H. atlanticus ) to the most shallow ( T. murphyi ).     

Prediction of antigenic determinants of trichosanthin by molecular modeling

The antigenic determinants of trichosanthin were predicted by molecular modeling.First,the threedimensional structure model of the antigen-binding fragment of anti-trichosanthin immunoglobulin E was built on the basis of its amino-acid sequence and the known three-dimensional structure of an antibody with similar sequence.Secondly,the preferable antigen-antibody interactions were obtained based on the known three-dimensional structure of trichosanthin and of the hypervariable regions of anti-trichosanthin immunoglobulin E.Two regions in the molecular surface of trichosanthin were found to form extensive interactions with the hypervariable regions of the antibody and have been predicted to be the possible antigenic determinants:one is composed of two polypeptide segments,Ile201-Glu210 and Ile225-Asp229,which are close to each other in the three-dimensional structure;and the other is the segment Lys173-Thr178.The former region seems to be the more reasonable antigenic determinant than the latter one.     

Homology modeling and molecular dynamics study on Schwanniomyces occidentalis alpha-amylase

With consumers growing increasingly aware of environmental issues, industries find enzymes as a reasonable alternative over physical conditions and chemical catalysts. Amylases are important hydrolase enzymes, which have been widely used in variety of industrial process such as pharmaceutical, food, and fermentation industries. Among amylases α-Amylase is in maximum demand due to its wide range of applications. The homology modeling study on Schwanniomyces occidentalis amylase (AMY1, UniProt identifier number: P19269) was performed by Modeller using Aspergillus oryzae (6TAA) as the template. The resulting structure was analyzed for validity and subjected to 14 ns of molecular dynamics (MD) simulation trough GROMACS. The validity of obtained model may represent that utilized OPLS force field is suitable for calcium-containing enzymes. DSSP secondary structure and contact map analysis represent the conservation of domain A TIM barrel feature together with calcium ion coordination sphere. Investigating the covariance matrix followed by principle component analyses for the first five eigenvectors of both trajectories indicate a little more flexibility for AMY1 structure. The electrostatic calculation for the final structures shows similar isoelectric point and superimposed buffering zone in the 5–8 pH range.     

Identification of protein functions from a molecular surface database,eF-site

A bioinformatics method was developed to identify the protein surface around the functional site and to estimate the biochemical function, using a newly constructed molecular surface database named the eF-site (electrostatic surface of Functional site. Molecular surfaces of protein molecules were computed based on the atom coordinates, and the eF-site database was prepared by adding the physical properties on the constructed molecular surfaces. The electrostatic potential on each molecular surface was individually calculated solving the Poisson–Boltzmann equation numerically for the precise continuum model, and the hydrophobicity information of each residue was also included. The eF-site database is accessed by the internet (http://pi.protein.osaka-u.ac.jp/eF-site/). We have prepared four different databases, eF-site/antibody, eF-site/prosite, eF-site/P-site, and eF-site/ActiveSite, corresponding to the antigen binding sites of antibodies with the same orientations, the molecular surfaces for the individual motifs in PROSITE database, the phosphate binding sites, and the active site surfaces for the representatives of the individual protein family, respectively. An algorithm using the clique detection method as an applied graph theory was developed to search of the eF-site database, so as to recognize and discriminate the characteristic molecular surfaces of the proteins. The method identifies the active site having the similar function to those of the known proteins.     

Binding of ibuprofen to human hemoglobin: elucidation of their molecular recognition by spectroscopy,calorimetry, and molecular modeling techniques

Ibuprofen, used for the treatment of acute and chronic pain, osteoarthritis, rheumatoid arthritis, and related conditions has ample affinity to globular proteins. Here we have explored this fundamental study pertaining to the interaction of ibuprofen with human hemoglobin (HHb), using multispectroscopic, calorimetric, and molecular modeling techniques to gain insights into molecular aspects of binding mechanism. Ibuprofen-induced graded decrease in absorption spectra indicates protein disruption along with sedimentation of HHb particle. Red shifting of absorption peak at 195 nm indicates alteration in the secondary structure of HHb upon interaction with ibuprofen. Flouremetric and isothermal titration calorimetric (ITC) studies suggested one binding site in HHb for ibuprofen at 298.15 K. However, with increase in temperature, ITC revealed increasing number of binding sites. The negative values of Gibbs energy change (ΔG⁰) and enthalpy change (ΔH⁰) along with positive value of entropy change (ΔS⁰) strongly suggest that it is entropy-driven spontaneous exothermic reaction. Moreover, hydrophobic interaction, hydrogen bonding, and π–π interaction play major role in this binding process as evidenced from ANS (8-anilino-1-napthalenesulphonic acid), sucrose binding, and molecular modeling studies. The interaction impacts on structural integrity and functional aspects of HHb as confirmed by CD spectroscopy, increased free iron release, increased rate of co-oxidation and decreased rate of esterase activity. These findings suggest us to conclude that ibuprofen upon interaction perturbs both structural and functional aspects of HHb.     

Identification of protein biochemical functions by similarity search using the molecular surface database eF-site

The identification of protein biochemical functions based on their three-dimensional structures is strongly required in the post-genome-sequencing era. We have developed a new method to identify and predict protein biochemical functions using the similarity information of molecular surface geometries and electrostatic potentials on the surfaces. Our prediction system consists of a similarity search method based on a clique search algorithm and the molecular surface database eF-site (electrostatic surface of functional-site in proteins). Using this system, functional sites similar to those of phosphoenoylpyruvate carboxy kinase were detected in several mononucleotide-binding proteins, which have different folds. We also applied our method to a hypothetical protein, MJ0226 from Methanococcus jannaschii, and detected the mononucleotide binding site from the similarity to other proteins having different folds.     

Binding of polarity-sensitive hydrophobic ligands to erythroid and nonerythroid spectrin: fluorescence and molecular modeling studies

We have used three polarity-sensitive fluorescence probes, 6-propionyl 2-(N,N-dimethyl-amino) naphthalene (Prodan), pyrene and 8-anilino 1-naphthalene sulphonic acid, to study their binding with erythroid and nonerythroid spectrin, using fluorescence spectroscopy. We have found that both bind to prodan and pyrene with high affinities with apparent dissociation constants (K_d) of .50 and .17?μM, for prodan, and .04 and .02?μM, for pyrene, respectively. The most striking aspect of these bindings have been that the binding stoichiometry have been equal to 1 in erythroid spectrin, both in dimeric and tetrameric form, and in tetrameric nonerythroid spectrin. From an estimate of apparent dielectric constants, the polarity of the binding site in both erythroid and nonerythroid forms have been found to be extremely hydrophobic. Thermodynamic parameters associated with such binding revealed that the binding is favored by positive change in entropy. Molecular docking studies alone indicate that both prodan and pyrene bind to the four major structural domains, following the order in the strength of binding to the Ankyrin binding domain?>?SH3 domain?>?Self-association domain?>?N-terminal domain of α-spectrin of both forms of spectrin. The binding experiments, particularly with the tetrameric nonerythroid spectrin, however, indicate more toward the self association domain in offering the unique binding site, since the binding stoichiometry have been 1 in all forms of dimeric and tetrameric spectrin, so far studied by us. Further studies are needed to characterize the hydrophobic binding sites in both forms of spectrin.     

A molecular modeling analysis of diethylstilbestrol conformations and their similarity to estradiol-17β

Crystallographic and computer modeling studies throughout the last 25 years have shown the structure of diethylstilbestrol (DES) to exist in two symmetrical or one asymmetrical conformation. As a result of specific comparisons to estradiol-17β (E₂), the asymmetrical DES conformer has been suggested as the geometry possessing estrogenic activity. In the present study, a more complete set of DES conformations has been elucidated through the use of computer modeling. All previously defined DES geometries were found within this new set of ten structural forms. Differences between the molecular mechanics heat of formation energies of the ten conformers, as well as the transition energies separating them from each other, were found to be less than I kcal/mol. Additionally, a computer-based molecular alignment method was employed to quantitatively compare the steric and electrostatic molecular features of each DES conformer relative to E₂. All ten DES structures were found to have shape relationships similar to E₂. Thus, a model for the estrogen action of DES is presented whereby this stilbene can favorably interact with the estrogen receptor regardless of the conformation or orientation of the initial ligand-receptor association.     

米曲霉木糖醇脱氢酶的分子模建与对接

[目的]研究米曲霉木糖醇脱氢酶基因的结构与功能.[方法]克隆测序来源于米曲霉的木糖醇脱氢酶(XDH)基因,利用Swiss-MODEL和Modeller对XDH进行三级结构模建,通过PROCHECK和Prosa2003对得到的4个目标模型进行评价,从中得到一个最佳模型.在同源建模的基础上,通过分子对接软件MolsoftICM-Pro,对辅因子进行对接,预测了XDH与NAD+、Zn2+作用的相关残基.寻找底物木糖醇与XDH结合的可能活性口袋,用Molsoft模拟XDH与木糖醇的对接,预测了酶与底物作用的关键氨基酸残基.[结果]结构分析显示,米曲霉XDH含有醇脱氢酶家族锌指纹结构和典型醇脱氢酶Rossmann折叠的辅酶结合域,属于Medium-chain脱氢酶(MDR)家族.通过对接研究,预测了XDH与NAD+之间形成氢键的氨基酸有Asp206、Arg211、Ser255、Ser301和Arg303,这些氨基酸位于结合域,与Zn2+形成氢键的氨基酸有His72和Glu73,位于催化域,与天然底物木糖醇形成氢键的氨基酸有Ile46、Ile349、Lys350和Thr351,位于催化域.[结论]所得信息对XDH分子定向改造、拓展米曲霉工业应用范围有重要意义.     

Inhibition of trypsin by cowpea thionin: characterization,molecular modeling,and docking

Higher plants produce several families of proteins with toxic properties, which act as defense compounds against pests and pathogens. The thionin family represents one family and comprises low molecular mass cysteine-rich proteins, usually basic and distributed in different plant tissues. Here, we report the purification and characterization of a new thionin from cowpea (Vigna unguiculata) with proteinase inhibitory activity. Cowpea thionin inhibits trypsin, but not chymotrypsin, binding with a stoichiometry of 1:1 as shown with the use of mass spectrometry. Previous annotations of thionins as proteinase inhibitors were based on their erroneous identification as homologues of Bowman-Birk family inhibitors. Molecular modeling experiments were used to propose a mode of docking of cowpea thionin with trypsin. Consideration of the dynamic properties of the cowpea thionin was essential to arrive at a model with favorable interface characteristics comparable with structures of trypsin-inhibitor complexes determined by X-ray crystallography. In the final model, Lys11 occupies the S1 specificity pocket of trypsin as part of a canonical style interaction.     

Study on the interaction of taiwaniaquinoids with FTO by spectroscopy and molecular modeling

In this work, an attempt has been made to study the interaction of four taiwaniaquinoids with fat mass and obesity-associated protein (FTO) by UV–vis absorption, fluorescence spectroscopy, and molecular docking techniques. The results indicated that taiwaniaquinoids effectively quenched the intrinsic fluorescence of FTO via static quenching. According to the binding constants and thermodynamic parameters at three different temperatures, the hydrophobic force and electrostatic interactions appeared be the predominant intermolecular forces in stabilizing the complex. Results revealed that W-4 was the strongest quencher and W-3 was the weakest. The results of synchronous and three-dimensional fluorescence spectra showed that the conformation of FTO was changed. In addition, the influence of molecular structure on the quenching effect has been investigated.     

Investigation of the molecular interactions of triclocarban with human serum albumin using multispectroscopies and molecular modeling

Triclocarban (TCC), as a broad spectrum antibacterial agent widely used in personal care products, has recently been recognized as environmental pollutant with the potential of adversely affecting wildlife and human health. However, the behavior of TCC in blood circulatory system and the potential toxicity of TCC at the molecular level have been poorly investigated. In this study, the effect of TCC on human serum albumin (HSA) and the binding mechanism of TCC to HSA were examined using spectroscopic techniques and molecular modeling methods. The fluorescence results suggested that the fluorescence of HSA was quenched by TCC through a static quenching mechanism and nonradiation energy transfer, and TCC was bound to HSA with moderately strong binding affinity via hydrophobic interaction based on the analysis of the thermodynamic parameters. The site marker competitive experiments revealed that TCC bound into subdomain IIA (site I) of HSA. In addition, the results obtained from the circular dichroism, Fourier transform infrared (FT-IR), 8-anilino-1-naphthalenesulfonic acid fluorescence, synchronous fluorescence, three-dimensional fluorescence spectra and dynamic light scattering suggested the change in the microenvironment and conformation of HSA during the binding reaction. Finally, the best binding mode of TCC and specific interaction of TCC with amino acid residues were determined using molecular docking and molecular dynamics simulations. In a word, the present studies can provide a way to help us well understand the transport, distribution and toxicity effect of TCC when it diffused in the human body.

Communicated by Ramaswamy H. Sarma     

Conformational analysis of the nonapeptide leuprorelin using NMR and molecular modeling

Summary The nonapeptide Leuprorelin, one of the LHRH agonists, was studied by means of 2D nuclear magnetic resonance spectroscopy and molecular modeling. NOESY spectra in aqueous/deuterated methanol solution (50% H₂O/CD₃OD) at low temperature (268 K) revealed short-range nOe connectivities (i, i+1), characteristic of flexibility of the molecule. The H ^N -H ^N sequential connectivities observed provide evidence that the sequence has the propensity to form a bend involving residues 5 and 6 and the N-terminal segment. The α-proton chemical shifts compared to random coil and additional data from the amide proton temperature coefficients support this assumption. One long-range nOe cross peak between H ₂ ^α -H ^NEth is indicative of proximity between C- and N-termini.     

Conformational analysis of the nonapeptide leuprorelin using NMR and molecular modeling

The nonapeptide Leuprorelin, one of the LHRH agonists, was studied by means of 2D nuclear magnetic resonance spectroscopy and molecular modeling. NOESY spectra in aqueous/deuterated methanol solution (50%H₂O/CD₃OD) at low temperature (268 K) revealed short-range nOe connectivities (i, i+1), characteristic of flexibility of the molecule. The H^N–H^N sequential connectivities observed provide evidence that the sequence has the propensity to form a bend involving residues 5 and 6 and the N-terminal segment. The -proton chemical shifts compared to random coil and additional data from the amide proton temperature coefficients support this assumption. One long-range nOe cross peak between H₂ –H^NEth is indicative of proximity between C- and N-termini.