Porphobilinogen deaminase and uroporphyrinogen III synthase: Structure,molecular biology,and mechanism

Porphobilinogen deaminase (hydroxymethylbilane synthase) and uroporphyrinogen III synthase (uroporphyrinogen III cosynthase) catalyze the transformation of four molecules of porphobilinogen, via the 1-hydroxymethylbilane, preuroporphyrinogen, into uroporphyrinogen III. A combination of studies involving protein chemistry, molecular biology, site-directed mutagenesis, and the use of chemically synthesized substrate analogs and inhibitors is helping to unravel the complex mechanisms by which the two enzymes function. The determination of the X-ray structure ofE. coli porphobilinogen deaminase at 1.76 Å resolution has provided the springboard for the design of further experiments to elucidate the precise mechanism for the assembly of both the dipyrromethane cofactor and the tetrapyrrole chain. The human deaminase structure has been modeled from theE. coli structure and has led to a molecular explanation for the disease acute intermittent porphyria. Molecular modeling has also been employed to simulate the spiro-mechanism of uroporphyrinogen III synthase.     

EzMol: A Web Server Wizard for the Rapid Visualization and Image Production of Protein and Nucleic Acid Structures

EzMol is a molecular visualization Web server in the form of a software wizard, located at http://www.sbg.bio.ic.ac.uk/ezmol/. It is designed for easy and rapid image manipulation and display of protein molecules, and is intended for users who need to quickly produce high-resolution images of protein molecules but do not have the time or inclination to use a software molecular visualization system. EzMol allows the upload of molecular structure files in PDB format to generate a Web page including a representation of the structure that the user can manipulate. EzMol provides intuitive options for chain display, adjusting the color/transparency of residues, side chains and protein surfaces, and for adding labels to residues. The final adjusted protein image can then be downloaded as a high-resolution image. There are a range of applications for rapid protein display, including the illustration of specific areas of a protein structure and the rapid prototyping of images.     

Software extensions to UCSF chimera for interactive visualization of large molecular assemblies

Many structures of large molecular assemblies such as virus capsids and ribosomes have been experimentally determined to atomic resolution. We consider four software problems that arise in interactive visualization and analysis of large assemblies: how to represent multimers efficiently, how to make cartoon representations, how to calculate contacts efficiently, and how to select subassemblies. We describe techniques and algorithms we have developed and give examples of their use. Existing molecular visualization programs work well for single protein and nucleic acid molecules and for small complexes. The methods presented here are proposed as features to add to existing programs or include in next-generation visualization software to allow easy exploration of assemblies containing tens to thousands of macromolecules. Our approach is pragmatic, emphasizing simplicity of code, reliability, and speed. The methods described have been distributed as the Multiscale extension of the UCSF Chimera (www.cgl.ucsf.edu/chimera) molecular graphics program.     

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.     

SIGNATURE: a single-particle selection system for molecular electron microscopy

SIGNATURE is a particle selection system for molecular electron microscopy. It applies a hierarchical screening procedure to identify molecular particles in EM micrographs. The user interface of the program provides versatile functions to facilitate image data visualization, particle annotation and particle quality inspection. The system design emphasizes both functionality and usability. This software has been released to the EM community and has been successfully applied to macromolecular structural analyses.     

Insight from the structural molecular model of cytidylate kinase from Mycobacterium tuberculosis

Mycobacterium tuberculosis is a gram-positive bacterium causes tuberculosis in human. H37Rv strain is a pathogenic strain utilized for tuberculosis research. The cytidylate mono-phosphate (CMP) kinase of Mycobacterium tuberculosis belongs to the family nucleoside mono-phosphate kinase (NMK), this enzyme is required for the bacterial growth. Therefore, it is important to study the structural and functional features of this enzyme in the control of the disease. Hence, we developed the structural molecular model of the CMP kinase protein from Mycobacterium tuberculosis by homology modeling using the software MODELLER (9v10). Based on sequence similarity with protein of known structure (template) of Mycobacterium smegmatis (PDB ID: 3R20) was chosen from protein databank (PDB) by using BLASTp. The energy of constructed models was minimized and the qualities of the models were evaluated by PROCHECK and VERRIFY-3D. Resulted Ramachandran plot analysis showed that conformations for 100.00% of amino acids residues are within the most favored regions. A possible homologous deep cleft active site was identified in the Model using CASTp program. Amino acid composition and polarity of that protein was observed by CLC-Protein Workbench tool. Expasy''s Prot-param server and CYC_REC tool were used for physiochemical and functional characterization of the protein. Studied of secondary structure of that protein was carried out by computational program, ProFunc. The structure is finally submitted in Protein Model Database. The predicted model permits initial inferences about the unexplored 3D structure of the CMP kinase and may be promote in relational designing of molecules for structure-function studies.     

Molecular modeling studies of atorvastatin analogues as HMGR inhibitors using 3D-QSAR,molecular docking and molecular dynamics simulations

3-Hydroxy-3-methylglutaryl coenzyme-A reductase (HMGR) is generally regarded as targets for the treatment of hypercholesterolemia. HMGR inhibitors (more commonly known as statins) are discovered as plasma cholesterol lowering molecules. In this work, 120 atorvastatin analogues were studied using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulation. The results show that the best CoMFA (comparative molecular field analysis) model has q² = 0.558 and r² = 0.977, and the best CoMSIA (comparative molecular similarity indices analysis) model has q² = 0.582 and r² = 0.919. Molecular docking and MD simulation explored the binding relationship of the ligand and the receptor protein. The calculation results indicated that the hydrophobic and electrostatic fields play key roles in QSAR model. After MD simulation, we found four vital residues (Lys735, Arg590, Asp690 and Asn686) and three hydrophobic regions in HMGR binding site. The calculation results show that atorvastatin analogues obtained by introduction of F atoms or gem-difluoro groups could obviously improve the inhibitory activity. The new HMGR inhibitor analogues design in this Letter had been submitted which is being currently synthesized by our laboratories.     

3D-QSAR (CoMFA,CoMSIA), molecular docking and molecular dynamics simulations study of 6-aryl-5-cyano-pyrimidine derivatives to explore the structure requirements of LSD1 inhibitors

Recently, Histone Lysine Specific Demethylase 1 (LSD1) was regarded as a promising anticancer target for the novel drug discovery. And several small molecules as LSD1 inhibitors in different structures have been reported. In this work, we carried out a molecular modeling study on the 6-aryl-5-cyano-pyrimidine fragment LSD1 inhibitors using three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. The results show that the best CoMFA model has q² = 0.802, r²_ncv = 0.979, and the best CoMSIA model has q² = 0.799, r²_ncv = 0.982. The electrostatic, hydrophobic and H-bond donor fields play important roles in the models. Molecular docking studies predict the binding mode and the interactions between the ligand and the receptor protein. Molecular dynamics simulations results reveal that the complex of the ligand and the receptor protein are stable at 300 K. All the results can provide us more useful information for our further drug design.     

Homology modeling, molecular dynamics and QM/MM study of the regulatory protein PhoP from Corynebacterium pseudotuberculosis

Corynebacterium pseudotuberculosis is a facultatively intracellular Gram-positive bacterium that causes caseous lymphadenitis, principally in sheep and goats, though sometimes in other species of animals, leading to considerable economic losses. This pathogen has a TCS known as PhoPR, which consists of a sensory histidine kinase protein (PhoR) and an intracellular response regulator protein (PhoP). This system is involved in the regulation of proteins present in various processes, including virulence. The regulation is activated by PhoP protein phosphorylation, an event that requires a magnesium (Mg²⁺) ion. Here we describe the 3D structure of the regulatory response protein (PhoP) of C. pseudotuberculosis through molecular modeling by homology. The model generated provides the first structural information on a full-length member of the OmpR/PhoP subfamily. Classical molecular dynamics was used to investigate the stability of the model. In addition, we used quantum mechanical/molecular mechanical techniques to perform (internal, potential) energy optimizations to determine the interaction energy between the Mg²⁺ ion and the structure of the PhoP protein. Analysis of the interaction energy residue by residue shows that Asp-16 and Asp-59 play an important role in the protein–Mg²⁺ ion interactions. These results may be useful for the future development of a new vaccine against tuberculosis based on genetic attenuation via a point mutation that results in the polar residue Asp-16 and/or Asp-59 being replaced with a nonpolar residue in the DNA-binding domain of PhoP of C. pseudotuberculosis.     

Synthesis,molecular docking,and antioxidant activity of new fluorescent tetrafluoroterphenyl analogues

Nucleophilic aromatic substitution (S_NAr) chemistry has been applied to develop many functionalized pentafluorobenzene derivatives. Those compounds are highly specific at the para position of the fluorinated ring. Therefore, they are typical adducts for the preparation of antioxidant molecular systems. In this context, we report the use of S_NAr chemistry as a suitable and simple approach for the synthesis of fluorescent antioxidant perfluorinated materials bearing ether bonds in various para-substituted alkoxy chains and with high purity and excellent yields. The fluoroterphenyl core was prepared via alkylation, Cu(I)-assisted decarboxylation, and cross-coupling using the potassium salt of fluorobenzoate, followed by the reaction with different alcohols. The structures of the synthesized fluoroterphenyl adducts were investigated using FT-IR, ¹H NMR, ¹³C NMR, and ¹⁹F NMR spectroscopy. The emission spectra and absorption spectra showed solvatochromism. The newly prepared tetrafluoroterphenyl analogues were investigated by antioxidant examination using the 2,2-diphenyl-1-picrylhydrazyl assay. Results were compared with ascorbic acid and butylated hydroxytoluene as references, and revealed that the tetrafluoroterphenyl analogues containing a decyl chain had the highest activity, with an IC₅₀ value of 22.36 ± 0.19 g/ml. The produced tetrafluoroterphenyl analogues were used in molecular docking strategies with a Protein Data Bank protein ID 5IKQ. The antioxidant investigations and docking results were convergent.     

PSPP: A Protein Structure Prediction Pipeline for Computing Clusters

Background

Protein structures are critical for understanding the mechanisms of biological systems and, subsequently, for drug and vaccine design. Unfortunately, protein sequence data exceed structural data by a factor of more than 200 to 1. This gap can be partially filled by using computational protein structure prediction. While structure prediction Web servers are a notable option, they often restrict the number of sequence queries and/or provide a limited set of prediction methodologies. Therefore, we present a standalone protein structure prediction software package suitable for high-throughput structural genomic applications that performs all three classes of prediction methodologies: comparative modeling, fold recognition, and ab initio. This software can be deployed on a user''s own high-performance computing cluster.

Methodology/Principal Findings

The pipeline consists of a Perl core that integrates more than 20 individual software packages and databases, most of which are freely available from other research laboratories. The query protein sequences are first divided into domains either by domain boundary recognition or Bayesian statistics. The structures of the individual domains are then predicted using template-based modeling or ab initio modeling. The predicted models are scored with a statistical potential and an all-atom force field. The top-scoring ab initio models are annotated by structural comparison against the Structural Classification of Proteins (SCOP) fold database. Furthermore, secondary structure, solvent accessibility, transmembrane helices, and structural disorder are predicted. The results are generated in text, tab-delimited, and hypertext markup language (HTML) formats. So far, the pipeline has been used to study viral and bacterial proteomes.

Conclusions

The standalone pipeline that we introduce here, unlike protein structure prediction Web servers, allows users to devote their own computing assets to process a potentially unlimited number of queries as well as perform resource-intensive ab initio structure prediction.     

Discovery of novel inhibitors of Mycobacterium tuberculosis MurG: homology modelling,structure based pharmacophore,molecular docking,and molecular dynamics simulations

MurG (Rv2153c) is a key player in the biosynthesis of the peptidoglycan layer in Mycobacterium tuberculosis (Mtb). This work is an attempt to highlight the structural and functional relationship of Mtb MurG, the three-dimensional (3D) structure of protein was constructed by homology modelling using Discovery Studio 3.5 software. The quality and consistency of generated model was assessed by PROCHECK, ProSA and ERRAT. Later, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with substrate Uridine-diphosphate-N-acetylglucosamine (UD1) facilitated us to employ structure-based virtual screening approach to obtain new hits from Asinex database using energy-optimized pharmacophore modelling (e-pharmacophore). The pharmacophore model was validated using enrichment calculations, and finally, validated model was employed for high-throughput virtual screening and molecular docking to identify novel Mtb MurG inhibitors. This study led to the identification of 10 potential compounds with good fitness, docking score, which make important interactions with the protein active site. The 25 ns MD simulations of three potential lead compounds with protein confirmed that the structure was stable and make several non-bonding interactions with amino acids, such as Leu290, Met310 and Asn167. Hence, we concluded that the identified compounds may act as new leads for the design of Mtb MurG inhibitors.     

Studies of H4R antagonists using 3D-QSAR,molecular docking and molecular dynamics

Three-dimensional quantitative structure–activity relationship studies were performed on a series of 88 histamine receptor 4 (H4R) antagonists in an attempt to elucidate the 3D structural features required for activity. Several in silico modeling approaches, including comparative molecular field analysis (CoMFA), comparative similarity indices analysis (CoMSIA), molecular docking, and molecular dynamics (MD), were carried out. The results show that both the ligand-based CoMFA model (Q ² = 0.548, R _ncv² = 0.870, R _pre² = 0.879, SEE = 0.410, SEP = 0.386) and the CoMSIA model (Q ² = 0.526, R _ncv² =0.866, R _pre² = 0.848, SEE = 0.416, SEP = 0.413) are acceptable, as they show good predictive capabilities. Furthermore, a combined analysis incorporating CoMFA, CoMSIA contour maps and MD results shows that (1) compounds with bulky or hydrophobic substituents at positions 4–6 in ring A (R2 substituent), positively charged or hydrogen-bonding (HB) donor groups in the R1 substituent, and hydrophilic or HB acceptor groups in ring C show enhanced biological activities, and (2) the key amino acids in the binding pocket are TRP67, LEU71, ASP94, TYR95, PHE263 and GLN266. To our best knowledge, this work is the first to report the 3D-QSAR modeling of these H4R antagonists. The conclusions of this work may lead to a better understanding of the mechanism of antagonism and aid in the design of new, more potent H4R antagonists.     

Computational Design of a PAK1 Binding Protein

We describe a computational protocol, called DDMI, for redesigning scaffold proteins to bind to a specified region on a target protein. The DDMI protocol is implemented within the Rosetta molecular modeling program and uses rigid-body docking, sequence design, and gradient-based minimization of backbone and side-chain torsion angles to design low-energy interfaces between the scaffold and target protein. Iterative rounds of sequence design and conformational optimization were needed to produce models that have calculated binding energies that are similar to binding energies calculated for native complexes. We also show that additional conformation sampling with molecular dynamics can be iterated with sequence design to further lower the computed energy of the designed complexes. To experimentally test the DDMI protocol, we redesigned the human hyperplastic discs protein to bind to the kinase domain of p21-activated kinase 1 (PAK1). Six designs were experimentally characterized. Two of the designs aggregated and were not characterized further. Of the remaining four designs, three bound to the PAK1 with affinities tighter than 350 μM. The tightest binding design, named Spider Roll, bound with an affinity of 100 μM. NMR-based structure prediction of Spider Roll based on backbone and ¹³C^β chemical shifts using the program CS-ROSETTA indicated that the architecture of human hyperplastic discs protein is preserved. Mutagenesis studies confirmed that Spider Roll binds the target patch on PAK1. Additionally, Spider Roll binds to full-length PAK1 in its activated state but does not bind PAK1 when it forms an auto-inhibited conformation that blocks the Spider Roll target site. Subsequent NMR characterization of the binding of Spider Roll to PAK1 revealed a comparably small binding ‘on-rate’ constant (? 10⁵ M^− 1 s^− 1). The ability to rationally design the site of novel protein-protein interactions is an important step towards creating new proteins that are useful as therapeutics or molecular probes.     

Protein flexibility and conformational states of Leishmania antigen eIF-4A: identification of a novel plausible protein adjuvant using comparative genomics and molecular modeling

Recent homology modeling studies have identified specific residues (epitope) of the Leishmania RNA helicase protein (LmeIF) that stimulates production of IL-12 cytokine. However, question remains concerning how LmeIF’s N-terminal moiety initiates adjuvant effects. Extensive molecular modeling combining the normal mode analysis (NMA) and molecular dynamics simulations, in the present study, has demonstrated that the LmeIF structure may exist in two different forms corresponding to the extended and collapsed (closed) states of the entire structure. The computational results showed that the two domains of the LmeIF structure tend to undergo large fluctuations in a concerted fashion and have strong effect on the solvent accessible surface of the epitope situated on the N-terminal structure. The conformational freedom of the C-terminal domains may explain why the entire LmeIF protein is not as active as the N-terminal moiety. Thereafter, a comparative genome analysis with subsequent homology modeling and molecular electrostatic potential (MEP) techniques allowed us to predict a novel and plausible RNA helicase (LI-helicase) from the Listeria source with adjuvant property as observed for the Leishmania eIF-4A protein. The structural folding and MEP maps revealed similar topologies of the epitope of both LmeIF and LI-helicase proteins and striking identity in the local disposition of the charged groups.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:7