Three-dimensional structure prediction of bovine AP lyase, BAP1: prediction of interaction with DNA and alterations as a result of Arg176-->Ala, Asp282-->Ala, and His308-->Asn mutations

BAP1 is an apurinic/apyrimidinic lyase (AP lyase) that plays an important role in the repair of DNA damage. The present study deals with the prediction of the 3D structure of bovine AP lyase based on its sequence homology with human AP lyase. The predicted 3D model of bovine AP1 shows remarkable similarity with human endonuclease in the overall 3D fold. However, significant differences in the model and the X-ray structure were located at some of the important sites. We have analyzed the active center of the enzyme and other sites that are involved in DNA repair. A number of amino acids bind the bases located in the major/minor grooves of DNA. An insertion of Arg176 in the major groove and Met270 in the minor groove caps the DNA bound enzyme's active site, stabilizing the extra helical AP site conformation and effectively locking the protein onto the AP-DNA. Three BAP1 mutants were also modeled and analyzed as regards the changes in the structure. Substitution of Arg176-->Ala leads to the loss of DNA binding whereas mutation of Asp282-->Ala and His308-->Asn leads to a decrease in the enzymatic activity.     

In silico determination of intracellular glycosylation and phosphorylation sites in human selectins: implications for biological function

Post-translational modifications provide the proteins with the possibility to perform functions in addition to those determined by their primary sequence. However, analysis of multifunctional protein structures in the environment of cells and body fluids is made especially difficult by the presence of other interacting proteins. Bioinformatics tools are therefore helpful to predict protein multifunctionality through the identification of serine and threonine residues wherein the hydroxyl group is likely to become modified by phosphorylation or glycosylation. Moreover, serines and threonines where both modifications are likely to occur can also be predicted (YinYang sites), to suggest further functional versatility. Structural modifications of hydroxyl groups of P-, E-, and L-selectins have been predicted and possible functions resulting from such modifications are proposed. Functional changes of the three selectins are based on the assumption that transitory and reversible protein modifications by phosphate and O-GlcNAc cause specific conformational changes and generate binding sites for other proteins. The computer-assisted prediction of glycosylation and phosphorylation sites in selectins should be helpful to assess the contribution of dynamic protein modifications in selectin-mediated inflammatory responses and cell-cell adhesion processes that are difficult to determine experimentally.     

SPINFAST: using structure alignment profiles to enhance the accuracy and assess the reliability of protein side-chain modeling

We present a novel, knowledge-based method for the side-chain addition step in protein structure modeling. The foundation of the method is a conditional probability equation, which specifies the probability that a side-chain will occupy a specific rotamer state, given a set of evidence about the rotamer states adopted by the side-chains at aligned positions in structurally homologous crystal structures. We demonstrate that our method increases the accuracy of homology model side-chain addition when compared with the widely employed practice of preserving the side-chain conformation from the homology template to the target at conserved residue positions. Furthermore, we demonstrate that our method accurately estimates the probability that the correct rotamer state has been selected. This interesting result implies that our method can be used to understand the reliability of each and every side-chain in a protein homology model.     

Optimization of the catalytic properties of Aspergillus fumigatus phytase based on the three-dimensional structure

Previously, we determined the DNA and amino acid sequences as well as biochemical and biophysical properties of a series of fungal phytases. The amino acid sequences displayed 49-68% identity between species, and the catalytic properties differed widely in terms of specific activity, substrate specificity, and pH optima. With the ultimate goal to combine the most favorable properties of all phytases in a single protein, we attempted, in the present investigation, to increase the specific activity of Aspergillus fumigatus phytase. The crystal structure of Aspergillus niger NRRL 3135 phytase known at 2.5 A resolution served to specify all active site residues. A multiple amino acid sequence alignment was then used to identify nonconserved active site residues that might correlate with a given favorable property of interest. Using this approach, Gln27 of A. fumigatus phytase (amino acid numbering according to A. niger phytase) was identified as likely to be involved in substrate binding and/or release and, possibly, to be responsible for the considerably lower specific activity (26.5 vs. 196 U x [mg protein](-1) at pH 5.0) of A. fumigatus phytase when compared to Aspergillus terreus phytase, which has a Leu at the equivalent position. Site-directed mutagenesis of Gln27 of A. fumigatus phytase to Leu in fact increased the specific activity to 92.1 U x (mg protein)(-1), and this and other mutations at position 27 yielded an interesting array of pH activity profiles and substrate specificities. Analysis of computer models of enzyme-substrate complexes suggested that Gln27 of wild-type A. fumigatus phytase forms a hydrogen bond with the 6-phosphate group of myo-inositol hexakisphosphate, which is weakened or lost with the amino acid substitutions tested. If this hydrogen bond were indeed responsible for the differences in specific activity, this would suggest product release as the rate-limiting step of the A. fumigatus wild-type phytase reaction.     

A coupling of homology modeling with multiple molecular dynamics simulation for identifying representative conformation of GPCR structures: a case study on human bombesin receptor subtype-3

In this study, a computational pipeline was therefore devised to overcome homology modeling (HM) bottlenecks. The coupling of HM with molecular dynamics (MD) simulation is useful in that it tackles the sampling deficiency of dynamics simulations by providing good-quality initial guesses for the native structure. Indeed, HM also relaxes the severe requirement of force fields to explore the huge conformational space of protein structures. In this study, the interaction between the human bombesin receptor subtype-3 and MK-5046 was investigated integrating HM, molecular docking, and MD simulations. To improve conformational sampling in typical MD simulations of GPCRs, as in other biomolecules, multiple trajectories with different initial conditions can be employed rather than a single long trajectory. Multiple MD simulations of human bombesin receptor subtype-3 with different initial atomic velocities are applied to sample conformations in the vicinity of the structure generated by HM. The backbone atom conformational space distribution of replicates is analyzed employing principal components analysis. As a result, the averages of structural and dynamic properties over the twenty-one trajectories differ significantly from those obtained from individual trajectories.     

S-adenosylhomocysteine hydrolase from the archaeon Pyrococcus furiosus: biochemical characterization and analysis of protein structure by comparative molecular modeling

S-adenosylhomocysteine hydrolase (AdoHcyHD) is an ubiquitous enzyme that catalyzes the breakdown of S-adenosylhomocysteine, a powerful inhibitor of most transmethylation reactions, to adenosine and L-homocysteine. AdoHcyHD from the hyperthermophilic archaeon Pyrococcus furiosus (PfAdoHcyHD) was cloned, expressed in Escherichia coli, and purified. The enzyme is thermoactive with an optimum temperature of 95 degrees C, and thermostable retaining 100% residual activity after 1 h at 90 degrees C and showing an apparent melting temperature of 98 degrees C. The enzyme is a homotetramer of 190 kDa and contains four cysteine residues per subunit. Thiol groups are not involved in the catalytic process whereas disulfide bond(s) could be present since incubation with 0.8 M dithiothreitol reduces enzyme activity. Multiple sequence alignment of hyperthermophilic AdoHcyHD reveals the presence of two cysteine residues in the N-terminus of the enzyme conserved only in members of Pyrococcus species, and shows that hyperthermophilic AdoHcyHD lack eight C-terminal residues, thought to be important for structural and functional properties of the eukaryotic enzyme. The homology-modeled structure of PfAdoHcyHD shows that Trp220, Tyr181, Tyr184, and Leu185 of each subunit and Ile244 from a different subunit form a network of hydrophobic and aromatic interactions in the central channel formed at the subunits interface. These contacts partially replace the interactions of the C-terminal tail of the eukaryotic enzyme required for tetramer stability. Moreover, Cys221 and Lys245 substitute for Thr430 and Lys426, respectively, of the human enzyme in NAD-binding. Interestingly, all these residues are fairly well conserved in hyperthermophilic AdoHcyHDs but not in mesophilic ones, thus suggesting a common adaptation mechanism at high temperatures.     

Generalized comparative modeling (GENECOMP): a combination of sequence comparison, threading, and lattice modeling for protein structure prediction and refinement

An improved generalized comparative modeling method, GENECOMP, for the refinement of threading models is developed and validated on the Fischer database of 68 probe-template pairs, a standard benchmark used to evaluate threading approaches. The basic idea is to perform ab initio folding using a lattice protein model, SICHO, near the template provided by the new threading algorithm PROSPECTOR. PROSPECTOR also provides predicted contacts and secondary structure for the template-aligned regions, and possibly for the unaligned regions by garnering additional information from other top-scoring threaded structures. Since the lowest-energy structure generated by the simulations is not necessarily the best structure, we employed two structure-selection protocols: distance geometry and clustering. In general, clustering is found to generate somewhat better quality structures in 38 of 68 cases. When applied to the Fischer database, the protocol does no harm and in a significant number of cases improves upon the initial threading model, sometimes dramatically. The procedure is readily automated and can be implemented on a genomic scale.     

Homology modeling and molecular dynamics study of West Nile virus NS3 protease: a molecular basis for the catalytic activity increased by the NS2B cofactor

The West Nile virus (WNV) NS3 serine protease, which plays an important role in assembly of infective virion, is an attractive target for anti-WNV drug development. Cofactors NS2B and NS4A increase the catalytic activity of NS3 in dengue virus and Hepatitis C virus, respectively. Recent studies on the WNV-NS3 characterize the catalytically active form of NS3 by tethering the 40-residue cofactor NS2B. It is suggested that NS2B is essential for the NS3 activity in WNV, while there is no information of the WNV-NS3-related crystal structure. To understand the role of NS2B/substrate in the NS3 catalytic activity, we built a series of models: WNV-NS3 and WNV-NS3-NS2B and WNV-NS3-NS2B-substrate using homology modeling and molecular modeling techniques. Molecular dynamics (MD) simulations were performed for 2.75 ns on each model, to investigate the structural stabilization and catalytic triad motion of the WNV NS3 protease with and without NS2B/substrate. The simulations show that the NS3 rearrangement occurs upon the NS2B binding, resulting in the stable D75-OD1...H51-NH hydrogen bonding. After the substrate binds to the NS3-NS2B active site, the NS3 protease becomes more stable, and the catalytic triad is formed. These results provide a structural basis for the activation and stabilization of the enzyme by its cofactor and substrate.     

Predicting the three-dimensional structure of human P-glycoprotein in absence of ATP by computational techniques embodying crosslinking data: insight into the mechanism of ligand migration and binding sites

P-glycoprotein is a membrane protein involved in the phenomenon of multidrug resistance. Its activity and transport function have been largely characterized by various biochemical studies and a low-resolution image has been obtained by electron microscopy. Obtaining a high-resolution structure is, however, still remote due to the inherent difficulties in the experimental determination of membrane protein structures. We present here a three-dimensional (3D) atomic model of P-glycoprotein in absence of ATP. This model was obtained using a combination of computational techniques including comparative modeling and rigid body dynamics simulations that embody all available cysteine disulfide crosslinking data characterizing the whole protein in absence of ATP. The model features rather well most of the experimental interresidue distances derived both in the transmembrane domains and in the nucleotide binding domains. The model is also in good agreement with electron microscopy data, particularly in terms of size and topology. It features a large cavity detected in the protein core into which seven ligands were successfully docked. Their predicted affinity correlates well with experimental values. Locations of docked ligands compare favorably with those suggested by cysteine-scanning data. The finding of different positions both for a single ligand and for different ligands corroborates the experimental evidence indicating the existence of multiple drug binding sites. The interactions identified between P-glycoprotein and the docked ligands reveal that different types of interactions such as H-bonds, pi-pi and cation-pi interactions occur in agreement with a recently proposed pharmacophore model of P-glycoprotein ligands. Furthermore, the model also displays a lateral opening located in the transmembrane domains connecting the lipid bilayer to the central cavity. This feature supports rather well the commonly admitted mechanism of substrate uptake from the lipid bilayer. We propose that this 3D model may be an important tool to understand the structure-function relationship of P-glycoprotein.     

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

Protein conformational transitions coupled to binding in molecular recognition of unstructured proteins: deciphering the effect of intermolecular interactions on computational structure prediction of the p27Kip1 protein bound to the cyclin A-cyclin-dependent kinase 2 complex

The relationship between folding mechanism coupled to binding and structure prediction of the tertiary complexes is studied for the p27(Kip) (1) protein which has an intrinsically disordered unbound form and undergoes a functional folding transition during complex formation with the phosphorylated cyclin A-cyclin-dependent kinase 2 (Cdk2) binary complex. Hierarchy of p27(Kip1) structural loss determined in our earlier studies from temperature-induced Monte Carlo simulations and subsequent characterization of the transition state ensemble (TSE) for the folding reaction have shown that simultaneous ordering of the p27(Kip1) native intermolecular interface for the beta-hairpin and beta-strand secondary structure elements is critical for nucleating a rapid kinetic transition to the native tertiary complex. In the present study, we investigate the effect of forming specific intermolecular interactions on structure prediction of the p27(Kip1) tertiary complex. By constraining different secondary structure elements of p27(Kip1) in their native bound conformations and conducting multiple simulated annealing simulations, we analyze differences in the success rate of predicting the native structure of p27(Kip1) in the tertiary complex. In accordance with the nucleation-condensation mechanism, we have found that further stabilization of the native intermolecular interface for the beta-hairpin and beta-strand elements of p27(Kip1), that become ordered in the TSE, but are hardly populated in the unbound state, results in a consistent acquisition of the native bound structure. Conversely, the excessive stablization of the local secondary structure elements, which are rarely detected in the TSE, has a detrimental effect on convergence to the native bound structure.     

Predicting the three-dimensional structure of the human facilitative glucose transporter glut1 by a novel evolutionary homology strategy: insights on the molecular mechanism of substrate migration, and binding sites for glucose and inhibitory molecules

The glucose transporters (GLUT/SLC2A) are members of the major facilitator superfamily. Here, we generated a three-dimensional model for Glut1 using a two-step strategy: 1), GlpT structure as an initial homology template and 2), evolutionary homology using glucose-6-phosphate translocase as a template. The resulting structure (PDB No. 1SUK) exhibits a water-filled passageway communicating the extracellular and intracellular domains, with a funnel-like exofacial vestibule (infundibulum), followed by a 15 A-long x 8 A-wide channel, and a horn-shaped endofacial vestibule. Most residues which, by mutagenesis, are crucial for transport delimit the channel, and putative sugar recognition motifs (QLS, QLG) border both ends of the channel. On the outside of the structure there are two positively charged cavities (one exofacial, one endofacial) delimited by ATP-binding Walker motifs, and an exofacial large side cavity of yet unknown function. Docking sites were found for the glucose substrate and its inhibitors: glucose, forskolin, and phloretin at the exofacial infundibulum; forskolin, and phloretin at an endofacial site next to the channel opening; and cytochalasin B at a positively charged endofacial pocket 3 A away from the channel. Thus, 1SUK accounts for practically all biochemical and mutagenesis evidence, and provides clues for the transport process.