Automated multiple analysis of protein structures: Application to homology modeling of cytochromes P450

A computational strategy for homology modeling, using several protein structures comparison, is described. This strategy implies a formalized definition of structural blocks common to several protein structures, a new program to compare these structures simultaneously, and the use of consensus matrices to improve sequence alignment between the structurally known and target proteins. Applying this method to cytochromes P450 led to the definition of 15 substructures common to P450_cam, P450_BM3, and P450_terp, and to proposing a 3D model of P450_eryF. Proteins 28:388–404, 1997 © 1997 Wiley-Liss, Inc.     

Exploring the structure of opioid receptors with homology modeling based on single and multiple templates and subsequent docking: A comparative study

Opioid receptors are the principal targets for opioids, which have been used as analgesics for centuries. Opioid receptors belong to the rhodopsin family of G-protein coupled receptors (GPCRs). In the absence of crystal structures of opioid receptors, 3D homology models have been reported with bovine rhodopsin as a template, though the sequence homology is low. Recently, it has been reported that use of multiple templates results in a better model for a target having low sequence identity with a single template. With the objective of carrying out a comparative study on the structural quality of the 3D models based on single and multiple templates, the homology models for opioid receptors (mu, delta and kappa) were generated using bovine rhodopsin as single template and the recently deposited crystal structures of squid rhodopsin, turkey β-1 and human β-2 adrenoreceptors along with bovine rhodopsin as multiple templates. In this paper we report the results of comparison between the refined 3D models based on multiple sequence alignment (MSA) and models built with bovine rhodopsin as template, using validation programs PROCHECK, PROSA, Verify 3D, Molprobity and docking studies. The results indicate that homology models of mu and kappa with multiple templates are better than those built with only bovine rhodopsin as template, whereas, in many aspects, the homology model of delta opioid receptor with single template is better with respect to the model based on multiple templates. Three nonselective ligands were docked to both the models of mu, delta and kappa opioid receptors using GOLD 3.1. The results of docking complied well with the pharamacophore, reported for nonspecific opioid ligands. The comparison of docking results for models with multiple templates and those with single template have been discussed in detail. Three selective ligands for each receptor were also docked. As the crystallographic structures are not yet known, this comparison will help in choosing better homology models of opioid receptors for studying ligand receptor interactions to design new potent opioid antagonists.     

Computational investigations of structural changes resulting from point mutations in a collagen-like peptide

The results of 0.5-1.0 ns molecular dynamics simulations of the collagen-like peptides [(POG)4(POA)(POG)4]3 and [(POG)9]3 (POG: proline-hydroxyproline-glycine) are presented. All simulations were performed using the AMBER-94 molecular mechanical force field with a shell of TIP3P waters surrounding the peptides. The initial geometries for the collagen-like peptides included an x-ray crystallographic structure, a computer-generated structure, a [(POG)9]3 structure modeled from the x-ray structure, and the x-ray structure with crystallographic waters replaced with a shell of modeled TIP3P waters. We examined the molecular dynamics peptide residue rms deviation fluctuations, dihedral angles, molecular and chain end-to-end distances, helical parameters, and peptide-peptide and peptide-solvent hydrogen-bonding patterns. Our molecular dynamics simulations of [(POG)4(POA)(POG)4]3 show average structures and internal coordinates similar to the x-ray crystallographic structure. Our results demonstrate that molecular dynamics can be used to reproduce the experimental structures of collagen-like peptides. We have demonstrated the feasibility of using the AMBER-94 molecular mechanical force field, which was parameterized to model nucleic acids and globular proteins, for fibril proteins. We provide a new interpretation of peptide-solvent hydrogen bonding and a peptide-peptide hydrogen bonding pattern not previously reported in x-ray studies. Last, we report on the differences; in particular with respect to main-chain dihedral angles and hydrogen bonding, between the native and mutant collagen-like peptides.     

Retardation of the unfolding process by single N-glycosylation of ribonuclease A based on molecular dynamics simulations

The conformational characteristics of glycosylated- and unglycosylated bovine pancreatic ribonuclease A (RNaseA) were traced with unfolding molecular dynamics simulations using CHARMM program at 470 K. The glycosylated RNase (Glc_RNase) possesses nearly identical protein structure with RNaseA, differing only by presence of a single acetylglucosamine residue N-linked to Asn34 in the RNaseA. Attaching of monomeric N-acetylglucosamine residue to the Asn34 in RNaseA resulted in a change of denaturing process of Glc_RNase. Simulations showed that the unfolding of RNaseA involved significant weakening of nonlocal interactions whereas the glycosylation led Glc_RNase to preserve the nonlocal interactions even in its denatured form. Even in simulations over 8 ns at 470 K, Glc_RNase remained relatively stable as a less denatured conformation. However, conformation of RNaseA was changed to a fully unfolded state before 3 ns of the simulations at 470 K. This difference was due to fact that formation of hydrogen bond bridges and nonlocal contacts induced by the attached N-acetylglucosamine of Glc_RNase showing in the unfolding simulations. These high-temperature unfolding MD simulations provided a theoretical basis for the previous experimental work in which Glc_RNase showed slower unfolding kinetics compared with unglycosylated RNaseA, suggesting that single N-glycosylation induced retardation of unfolding process of the ribonuclease protein.     

Integrated Homology Modelling and X-Ray Study of Herpes Simplex Virus I Thymidine Kinase: A Case Study

Abstract

Knowledge-based homology modelling together with site-directed mutagenesis, epitope and conformational mapping is an approach to predict the structures of proteins and for the rational design of new drugs. In this study we present how this procedure has been applied to model the structure of herpes simplex virus type 1 thymidine kinase (HSV1 TK, HSV1 ATP-thymidine-5′-phosphotransferase, EC 2.7.1.21). We have used, and evaluated, several secondary structure prediction methods, such as the classical one based on Chou and Fastman algorithm, neural networks using the Kabsch and Sander classification, and the PRISM method. We have validated the algorithms by applying them to the porcine adenylate kinase (ADK), whose three-dimensional structure is known and that has been used for the alignment of the TKs as well. The resulting first model of HSV1-TK consisted of the first β-strand connected to the phosphate binding loop and its subsequent α-helix, the fourth β-strand connected to the conserved FDRH sequence and two α-helix with basic amino acids. The 3D structure was built using the X-ray structure of ADK as template and following the general procedure for homology modelling. We extended the model by means of COMPOSER, an automatic process for protein modelling. Site-directed mutagenesis was used to experimentally verify the predicted active-site model of HSV1-TK. The data measured in our lab and by others support the suggestion that the FDRH motif is part of the active site and plays an important role in the phosphorylation of substrates. The structure of HSV1 TK, recently solved in collaboration with Prof. G. Schulz at 2.7 Å resolution, includes 284 of 343 residues of the N-terminal truncated TK. The secondary structures could be clearly assigned and fitted to the density. The comparison between crystallographically determined structure and the model shows that nearly 70% of the HSV1 TK structure has been correctly modelled by the described integrated approach to knowledge based ligand protein complex structure prediction. This indicate that computer assisted methods, combined with manual” correction both for alignment and 3D construction are useful and can be successful.     

Simulations of a G protein-coupled receptor homology model predict dynamic features and a ligand binding site

A computational approach to predict structures of rhodopsin-like G protein-coupled receptors (GPCRs) is presented and evaluated by comparison to the X-ray structural models. By combining sequence alignment, the rhodopsin crystal structure, and point mutation data on the beta2 adrenoreceptor (b2ar), we predict a (-)-epinephrine-bound computational model of the beta2 adrenoreceptor. The model is evaluated by molecular dynamics simulations and by comparison with the recent X-ray structures of b2ar. The overall correspondence between the predicted and the X-ray structural model is high. Especially the prediction of the ligand binding site is accurate. This shows that the proposed dynamic homology modelling approach can be used to create reasonable models for the understanding of structure and dynamics of other rhodopsin-like GPCRs.     

On the accuracy of homology modeling and sequence alignment methods applied to membrane proteins

In this study, we investigate the extent to which techniques for homology modeling that were developed for water-soluble proteins are appropriate for membrane proteins as well. To this end we present an assessment of current strategies for homology modeling of membrane proteins and introduce a benchmark data set of homologous membrane protein structures, called HOMEP. First, we use HOMEP to reveal the relationship between sequence identity and structural similarity in membrane proteins. This analysis indicates that homology modeling is at least as applicable to membrane proteins as it is to water-soluble proteins and that acceptable models (with C alpha-RMSD values to the native of 2 A or less in the transmembrane regions) may be obtained for template sequence identities of 30% or higher if an accurate alignment of the sequences is used. Second, we show that secondary-structure prediction algorithms that were developed for water-soluble proteins perform approximately as well for membrane proteins. Third, we provide a comparison of a set of commonly used sequence alignment algorithms as applied to membrane proteins. We find that high-accuracy alignments of membrane protein sequences can be obtained using state-of-the-art profile-to-profile methods that were developed for water-soluble proteins. Improvements are observed when weights derived from the secondary structure of the query and the template are used in the scoring of the alignment, a result which relies on the accuracy of the secondary-structure prediction of the query sequence. The most accurate alignments were obtained using template profiles constructed with the aid of structural alignments. In contrast, a simple sequence-to-sequence alignment algorithm, using a membrane protein-specific substitution matrix, shows no improvement in alignment accuracy. We suggest that profile-to-profile alignment methods should be adopted to maximize the accuracy of homology models of membrane proteins.     

A comparative study of available software for high-accuracy homology modeling: from sequence alignments to structural models

An open question in protein homology modeling is, how well do current modeling packages satisfy the dual criteria of quality of results and practical ease of use? To address this question objectively, we examined homology‐built models of a variety of therapeutically relevant proteins. The sequence identities across these proteins range from 19% to 76%. A novel metric, the difference alignment index (DAI), is developed to aid in quantifying the quality of local sequence alignments. The DAI is also used to construct the relative sequence alignment (RSA), a new representation of global sequence alignment that facilitates comparison of sequence alignments from different methods. Comparisons of the sequence alignments in terms of the RSA and alignment methodologies are made to better understand the advantages and caveats of each method. All sequence alignments and corresponding 3D models are compared to their respective structure‐based alignments and crystal structures. A variety of protein modeling software was used. We find that at sequence identities >40%, all packages give similar (and satisfactory) results; at lower sequence identities (<25%), the sequence alignments generated by Profit and Prime, which incorporate structural information in their sequence alignment, stand out from the rest. Moreover, the model generated by Prime in this low sequence identity region is noted to be superior to the rest. Additionally, we note that DSModeler and MOE, which generate reasonable models for sequence identities >25%, are significantly more functional and easier to use when compared with the other structure‐building software.     

Structure alignment of membrane proteins: Accuracy of available tools and a consensus strategy

Protein structure alignment methods are used for the detection of evolutionary and functionally related positions in proteins. A wide array of different methods are available, but the choice of the best method is often not apparent to the user. Several studies have assessed the alignment accuracy and consistency of structure alignment methods, but none of these explicitly considered membrane proteins, which are important targets for drug development and have distinct structural features. Here, we compared 13 widely used pairwise structural alignment methods on a test set of homologous membrane protein structures (called HOMEP3). Each pair of structures was aligned and the corresponding sequence alignment was used to construct homology models. The model accuracy compared to the known structures was assessed using scoring functions not incorporated in the tested structural alignment methods. The analysis shows that fragment‐based approaches such as FR‐TM‐align are the most useful for aligning structures of membrane proteins. Moreover, fragment‐based approaches are more suitable for comparison of protein structures that have undergone large conformational changes. Nevertheless, no method was clearly superior to all other methods. Additionally, all methods lack a measure to rate the reliability of a position within a structure alignment. To solve both of these problems, we propose a consensus‐type approach, combining alignments from four different methods, namely FR‐TM‐align, DaliLite, MATT, and FATCAT. Agreement between the methods is used to assign confidence values to each position of the alignment. Overall, we conclude that there remains scope for the improvement of structural alignment methods for membrane proteins. Proteins 2015; 83:1720–1732. © 2015 Wiley Periodicals, Inc.     

3 Nsec molecular dynamics simulation of the protein ubiquitin and comparison with X-ray crystal and solution NMR structures.

Mainly due to computational limitations, past protein molecular dynamics simulations have rarely been extended to 300 psec; we are not aware of any published results beyond 350 psec. The present work compares a 3000 psec simulation of the protein ubiquitin with the available x-ray crystallographic and solution NMR structures. Aside from experimental structure availability, ubiquitin was studied because of its relatively small size (76 amino acids) and lack of disulfide bridges. An implicit solvent model was used except for explicit treatment of waters of crystallization. We found that the simulated average structure retains most of the character of the starting x-ray crystal structure. In two highly surface accessible regions, the simulation was not in agreement with the x-ray structure. In addition, there are six backbone-backbone hydrogen bonds that are in conflict between the solution NMR and x-ray crystallographic structures; two are bonds that the NMR does not locate, and four are ones that the two methods disagree upon the donor. Concerning these six backbone-backbone hydrogen bonds, the present simulation agrees with the solution NMR structure in five out-of-the six cases, in that if a hydrogen bond is present in the x-ray structure and not in the NMR structure, the bond breaks within 700 psec. Of the two hydrogen bonds that are found in the NMR structure and not in the x-ray structure, one forms at 1400 psec and the other forms rarely. The present results suggest that relatively long molecular dynamics simulations, that use protein x-ray crystal coordinates for the starting structure and a computationally efficient solvent representation, may be used to gain an understanding of conformational and dynamic differences between the solid-crystal and dilute-solution states.     

Revision of the Classical Dopamine D2 Agonist Pharmacophore Based on an Integrated Medicinal Chemistry,Homology Modelling and Computational Docking Approach

The scientific advances during the 1970ies and 1980ies within the field of dopaminergic neurotransmission enabled the development of a pharmacophore that became the template for design and synthesis of dopamine D₂ agonists during the following four decades. A major drawback, however, is that this model fails to accommodate certain classes of restrained dopamine D₂ agonists including ergoline structures. To accommodate these, a revision of the original model was required. The present study has addressed this by an extension of the original model without compromising its obvious qualities. The revised pharmacophore contains an additional hydrogen bond donor feature, which is required for it to accommodate ergoline structures in a low energy conformation and in accordance with the steric restrictions dictated by the original model. The additional pharmacophore feature suggests ambiguity in the binding mode for certain compounds, including a series of ergoline analogues, which was reported recently. The ambiguity was confirmed by docking to a homology model of the D₂ receptor as well as by pharmacological characterization of individual enantiomers of one of the analogues. The present research also addresses the potential of designing ligands that interact with the receptor in a large, distal cavity of the dopamine D₂ receptor that has not previously been studied systematically. The pharmacological data indicate that this area may be a major determinant for both the dopamine D₂ affinity and efficacy, which remains to be explored in future studies.     

Secondary structure model for the ITS-2 precursor rRNA of strongyloid nematodes of equids: implications for phylogenetic inference

In order to maximise the positional homology in the primary sequence alignment of the second internal transcribed spacer for 30 species of equine strongyloid nematodes, the secondary structures of the precursor ribosomal RNA were predicted using an approach combining an energy minimisation method and comparative sequence analysis. The results indicated that a common secondary structure model of the second internal transcribed spacer of these nematodes was maintained, despite significant interspecific differences (2–56%) in primary sequences. The secondary structure model was then used to refine the primary second internal transcribed spacer sequence alignment. The “manual” and “structure” alignments were both subjected to phylogenetic analysis using three different tree-building methods to compare the effect of using different sequence alignments on phylogenetic inference. The topologies of the phylogenetic trees inferred from the manual second internal transcribed spacer alignment were usually different to those derived from the structure second internal transcribed spacer alignment. The results suggested that the positional homology in the second internal transcribed spacer primary sequence alignment was maximised when the secondary structure model was taken into consideration.     

A revised model of the active site of alternative oxidase

The plant mitochondrial protein alternative oxidase catalyses dioxygen dependent ubiquinol oxidation to yield ubiquinone and water. A structure of this protein has previously been proposed based on an assumed structural homology to the di-iron carboxylate family of proteins. However, these authors suggested the protein has a very different topology than the known structures of di-iron carboxylate proteins. We have re-examined this model and based on comparison of recent sequences and structural data on di-iron carboxylate proteins we present a new model of the alternative oxidase which allows prediction of active site residues and a possible membrane binding motif.     

Protein structure alignment

A new method of comparing protein structures is described, based on distance plot analysis. It is relatively insensitive to insertions and deletions in sequence and is tolerant of the displacement of equivalent substructures between the two molecules being compared. When presented with the co-ordinate sets of two structures, the method will produce automatically an alignment of their sequences based on structural criteria. The method uses the dynamic programming optimization technique, which is widely used in the comparison of protein sequences and thus unifies the techniques of protein structure and sequence comparison. Typical structure comparison problems were examined and the results of the new method compared to the published results obtained using conventional methods. In most examples, the new method produced a result that was equivalent, and in some cases superior, to those reported in the literature.     

Integrating multi-scale data on homologous recombination into a new recognition mechanism based on simulations of the RecA-ssDNA/dsDNA structure

RecA protein is the prototypical recombinase. Members of the recombinase family can accurately repair double strand breaks in DNA. They also provide crucial links between pairs of sister chromatids in eukaryotic meiosis. A very broad outline of how these proteins align homologous sequences and promote DNA strand exchange has long been known, as are the crystal structures of the RecA-DNA pre- and postsynaptic complexes; however, little is known about the homology searching conformations and the details of how DNA in bacterial genomes is rapidly searched until homologous alignment is achieved. By integrating a physical model of recognition to new modeling work based on docking exploration and molecular dynamics simulation, we present a detailed structure/function model of homology recognition that reconciles extremely quick searching with the efficient and stringent formation of stable strand exchange products and which is consistent with a vast body of previously unexplained experimental results.     

An in vitro binding assay for angiogenin using placental ribonuclease inhibitor

A convenient in vitro assay for angiogenin has been developed which greatly facilitates its routine detection and quantitation. The assay is based on the capacity of angiogenin to bind placental ribonuclease inhibitor (PRI); it is less tedious and more versatile than existing procedures that measure blood vessel growth or cleavage of rRNA. The test sample is added to a reaction mixture containing a known quantity of PRI, which complexes any angiogenin present in the sample. A slight excess of RNase A, relative to PRI, is then added, and the amount of RNase A which remains unbound is determined by measuring the generation of acid-soluble fragments from yeast RNA. The assay is sensitive to 30 fmol of angiogenin and is linear over a 17-fold concentration range. Use of the binding assay in parallel with a conventional RNase A assay provides a means of detecting angiogenin in chromatographic fractions and differentiating it from RNases. This procedure makes possible the isolation of angiogenin from new sources, such as nonhuman sera. It may also be applicable to other biologically active proteins with sequence homology to RNase A, e.g., eosinophil cationic protein or eosinophil derived neurotoxin.     

A comparison of the predicted and X-ray structures of angiogenin. Implications for further studies of model building of homologous proteins

The three-dimensional structure of human angiogenin has been determined by X-ray crystallography and is compared here with an earlier model which predicted its structure, based on the homology of angiogenin with bovine pancreatic ribonuclease A. Comparison of the predicted model and crystal structure shows that the active-site histidine residues and the core of the angiogenin molecule, including most of the-strands and-helices, were predicted reasonably well. However, the structure of the surface loop regions and residues near the truncated C-terminus differs significantly. The C-terminal segment includes the active-site residues Asp-116, Gln-117, and Ser-118; Gln-117 in particular has been shown to be important in affecting the ribonucleolytic activity of angiogenin. Also, the orientation of one helix in the model differed from the orientation observed experimentally by about 20°, resulting in a large displacement of this chain segment. The difficulty encountered in predicting the surface loop regions has led to a new algorithm [Palmer and Scheraga (1991),J. Comput. Chem.,12, 505–526; (1992),J. Comput. Chem.,13, 329–350] for predicting the conformations of surface loops.     

Protein multiple alignment incorporating primary and secondary structure information.

Identifying common local segments, also called motifs, in multiple protein sequences plays an important role for establishing homology between proteins. Homology is easy to establish when sequences are similar (sharing an identity > 25%). However, for distant proteins, it is much more difficult to align motifs that are not similar in sequences but still share common structures or functions. This paper is a first attempt to align multiple protein sequences using both primary and secondary structure information. A new sequence model is proposed so that the model assigns high probabilities not only to motifs that contain conserved amino acids but also to motifs that present common secondary structures. The proposed method is tested in a structural alignment database BAliBASE. We show that information brought by the predicted secondary structures greatly improves motif identification. A website of this program is available at www.stat.purdue.edu/~junxie/2ndmodel/sov.html.     

Role of two consecutive α,β-dehydrophenylalanines in peptide structure: Crystal and molecular structure of Boc-Leu-ΔPhe-ΔPhe-Ala-Phe-NHMe

An N^α-protected model pentapeptide containing two consecutive ΔPhe residues, Boc-Leu-ΔPhe-ΔPhe-Ala-Phe-NHMe, has been synthesized by solution methods and fully characterized. ¹H-nmr studies provided evidence for the occurrence of a significant population of a conformer having three consecutive, intramolecularly H-bonded β-bends in solution. The solid state structure has been determined by x-ray diffraction methods. The crystals grown from aqueous methanol are orthorhombic, space group P2₁2₁2₁, a = 11.503(2), b = 16.554(2), c = 22.107(3) Å, V = 4209(1) Å,³ and Z = 4. The x-ray data were collected on a CAD4 diffractometer using CuK_a radiation (λ = 1.5418 Å). The structure was determined using direct methods and refined by full-matrix least-squares procedure. The R factor is 5.3%. The molecule is characterized by a right handed 3₁₀-helical conformation (〈ϕ〉 = −68.2°, 〈ψ〉 = −26.3°), which is made up of two consecutive type III β-bends and one type I β-bend. In the solid state the helical molecules are aligned head-to-tail, thus forming long rod like structures. A comparison with other peptide structures containing consecutive ΔPhe residues is also provided. The present study confirms that the -ΔPhe-ΔPhe-sequence can be accommodated in helical structures. © 1997 John Wiley & Sons, Inc. Biopoly 42: 373–382, 1997     

Computer modeling and molecular dynamics simulations of ligand bound complexes of bovine angiogenin: dinucleotide topology at the active site of RNase a family proteins

We have undertaken the modeling of substrate-bound structures of angiogenin. In our recent study, we modeled the dinucleotide ligand binding to human angiogenin. In the present study, the substrates CpG, UpG, and CpA were docked onto bovine angiogenin. This was achieved by overcoming the problem of an obstruction to the B1 site by the C-terminus and identifying residues that bind to the second base. The modeled complexes retain biochemically important interactions. The docked models were subjected to 1 ns of molecular dynamics, and structures from the simulation were refined by using simulated annealing. Our models explained the enzyme's specificity for both B1 and B2 bases as observed experimentally. The nature of binding of the dinucleotide substrate was compared with that of the mononucleotide product. The models of these complexes were also compared with those obtained earlier with human angiogenin. On the basis of the simulations and annealed structures, we came up with a consensus topology of dinucleotide ligands that binds to human and bovine angiogenins. This dinucleotide conformation can serve as a starting model for ligand-bound complex structures for RNase A family of proteins. We demonstrated this capability by generating the complex structure of CpA bound to eosinophil-derived neurotoxin (EDN) by fitting the consensus topology of CpA to the crystal structure of native EDN.