PubMed Health⌕ Search

Biomedical subjects

A Goldblum

Publications and source records attributed to A Goldblum.

At least 19 recordsLinked to original sources

Computational protein design: a novel path to future protein drugs.

Computational protein design emerges in recent years as a field that could make a substantial impact on the design of protein drugs. It still consists mainly of redesigning parts of a protein sequence for increasing the stability of a given 3-dimensional conformation of a protein, but has already been extended from redesigning core residues to redesigning in all other protein regions, as well as to the design of backbone conformations. More recently, proteins with new binding functions and new enzymes, protein libraries, designs of full folds and of a new protein fold, have been some of the main highlights. The search and the scoring problems are however not fully solved, and many of the design processes should be examined on much larger scales in order to assess their usefulness. We examine some of the basic assumptions in computational protein design, in particular, the separation between sequence and scaffold designs. Among others, we suggest to include more protein residues in computations, to include relevant parts of the backbone, to use appropriate reference states, to produce the proteins and to validate the designs by structural examination of the protein products.

Computational Biology↗

Stochastic algorithm for kinase homology model construction.

A stochastic algorithm for constructing multiple loops in homology modeling of proteins is presented. The algorithm discards variable values in iterations based on a cost function and on statistical analysis of results. Values that remain are used for constructing an ensemble of best solutions. In test cases, the stochastic algorithm retains all the best solutions, compared to an exhaustive scan of the full set of conformations. Individual loops are constructed by adding dipeptide units. Dipeptide conformations are extracted from a database of proteins and their conformations include bond lengths and all angles. Single loops are constructed from both N- and C- terminals to the center, and loop closure is evaluated by a combination of penalties for the peptide closure and Miyazawa-Jernigan (MJ) [1]. residue-residue interactions with the rest of the protein. Large ensembles of each loop are clustered and re-evaluated with a refined [2]. energy term. The reduced, clustered set of each loop is then employed to construct simultaneously all the loops. The algorithm was applied to construct simultaneously six loops in c-Src kinase family proteins, incorporating a total of 37-40 residues. The best RMSD for reconstructing the loops is 1.45 A for Lck (structure 1QPE in the Protein Data Bank) and 2.54 A for human c-Src (structure 1FMK). The multiple loop conformations with lowest energy have higher RMSD values, of 2.06 A and 3.09 A, respectively. The average RMSD values for the first 1000 conformations are 3.00 A and 3.46 A, respectively. Models for the "open" structures of c-Src and of Jak-2 were constructed on the basis of 1QPE. The Jak-2 model is found to be more flexible in the loops region than its c-Src counterpart.

Algorithms↗

Two cysteine residues in the DNA-binding domain of CREB control binding to CRE and CREB-mediated gene expression.

The cAMP-responsive element-binding protein (CREB) has been implicated in the regulation of numerous physiological functions including those of several hypoxia-responding genes. All CREB transcription-regulated genes harbor the eight base-pair cAMP-responsive element (CRE) or the seven base-pair AP-1 sequence. Utilizing mutational analysis and biochemical assays, we found that reduction of two cysteine residues located in the DNA-binding basic domain of CREB, enhances the binding efficiency of CREB to DNA and regulates CRE-mediated gene expression. Substitution of these residues to serine renders insensitivity to reduction, hypoxia and to the sulfhydryl-specific modifying agent, N-ethylmaleimide. These substitutions enhance the binding of CREB to its cognate DNA sites under oxidative conditions, and of the CREB-dependent gene expression during normoxia. These findings are supported by results of molecular modeling of the CREB-CRE interactions. We also found that HTLV-1 Tax enhancement of CREB binding to the cellular and the viral DNA sites and activation of the CRE-dependent gene expression are independent of CREB activation exerted by redox conditions. The genetic biochemical and molecular modeling presented in this work indicate that the two cysteine residues in the bZIP domain of CREB regulate the binding efficiency of CREB to its cognate DNA sites and as a consequence the activation of CREB-mediated gene expression.

Amino Acid Sequence↗

A novel energy-based stochastic method for positioning polar protons in protein structures from X-rays.

A novel automated method for the optimal placement of polar hydrogens in a protein structure is presented. The algorithm adds initially, to a protein data bank file of the protein, nonrotatable hydrogens such as peptide backbone hydrogens according to geometric considerations. Then, water protons and polar side chain protons of lysine, serine, threonine, tyrosine, aspartic acid, glutamic acid, and the C and N termini of a protein are added according to energy considerations. A unique stochastic approach has been developed to overcome a combinatorial explosion in the search for the lowest energy structure. First, the system is divided into ensembles. Each ensemble is treated separately: N conformations are sampled at random, their energies computed, whereas common components of high-energy combinations are gathered on one hand, and low-energy combinations on the other. Components that yield only high-energy conformations and do not contribute to any low energies are excluded. This is reiterated while the total amount of combinations is decreased along the iterative process. When the total number of combinations is lower than a user defined threshold, all remaining combinations are evaluated by exhaustive search. Energy evaluations use nonbonding energy expressions alone. The program was tested on five high-resolution crystal structures: bovine pancreatic trypsin inhibitor (Brookhaven Protein Data Bank file 5PTI), RNase-A (5RSA), trypsin (1NTP), and carbon monoxymyoglobin (2MB5), for which neutron diffraction structures are available, as well as phosphate binding protein (1IXH) for which very high resolution X-ray crystallography was used. The low RMS values prove the efficiency of this algorithm as a tool for positioning protons in proteins. It may be used for other biological structures.

Algorithms↗

A novel computational method for predicting the transmembrane structure of G-protein coupled receptors: application to human C5aR and C3aR.

A novel algorithm was applied to the sequences of bacteriorhodopsin (BRh), of rhodopsin (Rh), and of the two human anaphylatoxin receptors, C5a-receptor (hC5aR) and C3a-receptor (hC3aR), that predicts their transmembrane domains (TMD) according to energy criteria alone, on the basis of their sequences and a template structure for each. Two consecutive criteria were applied for the predictions: the first is hydrophobicity of a sequence of residues, which determines the candidate stretches of residues that form one of the transmembrane helices. The second criterion is an energy function composed of inter residue contact energies, of hydrophobic contributions due to membrane exposure and of the interactions of a few residues with the phospholipid head groups. The sequence of candidate residues for each helix is longer than that of the template, and is finally determined by threading each of the candidate stretches on each of the template helices and evaluating the energy for all possible configurations. Contact energies between residues were taken from a database (Miyazawa S and Jernigan RL (1996) J Mol Biol 256 623-44). The algorithm predicts well the TMD structure of BRh based on its own template, and the TMD structure of Rh conforms well with the model of Baldwin et al (Baldwin JM Schertler GFX and Unger VM (1997) J Biol Chem 272 144-64). Results for the construction of the TMD of hC5aR and hC3aR were compared, employing the template structure of Rh. Most of the results for these receptors are in accord with alignments and with mutation experiments on hC5aR and hC3aR. The predictions may serve as a basis for future mutagenesis experiments of these receptors.

Algorithms↗

Point mutation in avian sarcoma leukaemia virus protease which increases its activity but impairs infectious virus production.

The retrovirus protease (PR), an aspartic PR, is composed of two identical subunits, each containing a conserved tripeptide sequence present at the active site of the enzyme. Asp-Ser-Gly is found in avian sarcoma leukaemia viruses (ASLV) and Asp-Thr-Gly in mammalian oncoretroviruses. We have mutated the conserved sequence at the active site of ASLV PR by converting the Ser and Gly residues to Thr and Ala, respectively. Replacement of Gly with Ala yielded an ASLV PR devoid of proteolytic activity. The Ser to Thr conversion did not alter the substrate specificity of the enzyme. Both wild-type and mutated PRs correctly cleaved viral precursors expressed in bacterial cells, as well as synthetic peptides homologous to ASLV and human immunodeficiency virus type 1 cleavage sites. Bacterially produced ASLV PR with Thr instead of Ser had increased enzymatic activity, as shown by hydrolysis of synthetic peptides. However, this mutation reduced the production of reverse transcriptase-containing particles and infectious virus following transfection of permissive cells with virus DNA.

Alpharetrovirus↗

A molecular model for an anionic opiate mu-receptor: affinity and activation of morphine conformers.

In this study, the minimal neglect of differential overlap, hydrogen bonding corrected (MNDO/H) method was used to construct an explicit three component mu-opioid receptor binding site composed of formate (Glu-, Asp-), H2O (Ser, Thr) and NH4+ (Lys+, Arg+) which has optimum interactions with the protonated amine and polar oxygen regions of morphine, respectively. These moieties are common to most classes of opioids and are thought to be involved in key interactions of the protonated form with the receptor leading to recognition and activation. Two predominant conformers of morphineH+, N-Me (equatorial) and N-Me (axial) were used as templates for construction of the binding site. In these studies, a plausible recognition mechanism involving electrostatic interaction between the protonated amine and an anionic receptor site, together with a proton-donating phenolic group and proton-accepting polar oxygens was characterized. The role of hydration of both the receptor site and morphine in determining affinity was also explicitly considered. The results strongly indicate that the high-affinity binding of morphineH+ to an 'anionic' receptor site is due mainly to a large entropy term resulting from expulsion of H2O from the receptor site upon introduction of morphine. A possible mechanism of receptor activation was also explored involving proton transfer from morphine to the receptor. Two results obtained support the plausibility of this mechanism: the barrier to proton transfer is reduced by receptor interaction with the polar oxygens and a conformational change occurs in the model receptor during this process.

Ammonia↗

Modulation of the affinity of aspartic proteases by the mutated residues in active site models.

The active sites of 3 types of aspartic proteases are modeled, based on crystallographic coordinates of endothiapepsin and of a model of HIV-1 protease. The enthalpies of deprotonation from neutral to mono-anion and to dianion are calculated with semiempirical minimal neglect of differential overlap, hydrogen bonding corrected (MNDO/H). This quantum mechanical study of models for the active sites of pepsins, human renin and retroviral aspartic proteases demonstrates that the replacements of Thr-218 from pepsins by Ala in human renin and of both Ser-35 and Thr-218 by alanines in retroviral proteases increases the proton affinity and modulates the charge distribution of those active sites compared to the pepsins.

Amino Acid Sequence↗

Quantum mechanical modeling of aspartic proteinase interactions: difference in binding of diastereomeric statine models.

Quantum mechanical calculations were carried out for the interaction of two diastereomeric model inhibitors of aspartic proteinases with a model for the active site, based on crystallographic coordinates of endothiapepsin. The model inhibitor is formamide N-(2-hydroxy 3-methyl propane) and the active site is represented by the full backbone and most of the side chains of the two partial sequences D32-T33-G34-S35 and D215-T216-G217-T218. Those calculations demonstrate that the best binding mode for this short inhibitor is consistent with the X-ray experiments and somewhat stronger with the inhibitor in a 2(S) configuration, compared to 2(R). Another binding mode is possible for this model inhibitor only in the 2(S)-configuration, and is weaker than the first.

Amino Acids↗

Theoretical calculations on the acidity of the active site in aspartic proteinases.

Semiempirical minimal neglect of differential overlap-self-consistent field calculations, corrected and modified for multiple hydrogen-bonding interactions, were applied to models of the active site of aspartic proteinases (AP). The propensities of the two active-site aspartates to ionize were compared under the influence of various neighboring residues and of water molecules. Asp-32 and Asp-215 in three aspartic proteinases (endothiapepsin, Rhizopus pepsin, and penicillopepsin) are found to be basically asymmetric, Asp-32 being preferentially (by 2-3 kcal) ionized with respect to Asp-215. In penicillopepsin, this asymmetry is compensated by effects of surrounding residues. In our largest model for the active site, which includes such other residues, near equality is found for the ionizing tendency of Asp-32 and Asp-215. The pK difference is rationalized in terms of first and second ionizations of the full active-site model. Its ionization enthalpies correlate well with those of other small organic diacids. This "gas-phase" approach to AP active-site interactions represents the main possible contributions to the acidity of the active site.

Aspartic Acid Endopeptidases↗

Geometrical analysis of Cys-Cys bridges in proteins and their prediction from incomplete structural information.

Analysis of C-alpha atom positions from cysteines involved in disulphide bridges in protein crystals shows that their geometric characteristics are unique with respect to other Cys-Cys, non-bridging pairs. They may be used for predicting disulphide connections in incompletely determined protein structures, such as low resolution crystallography or theoretical folding experiments. The basic unit for analysis and prediction is the 3 x 3 distance matrix for Cx positions of residues (i - 1), Cys(i), (i +1) with (j - 1), Cys(j), (j + 1). In each of its columns, row and diagonal vector--outer distances are larger than the central distance. This analysis is compared with some analytical models.

Algorithms↗

An improved approach to the analysis of drug-protein binding by distance geometry.

The calculation of side chain centers of coordinates and the subsequent generation of side chain-side chain and side chain-backbone distance matrices is suggested as an improved method for viewing interactions inside proteins and for the comparison of protein structures. The use of side chain distance matrices is demonstrated with free PTI, and the use of difference distance matrices for side chains is shown for free and trypsin-bound PTI as well as for the X-ray structures of trypsin complexes with PTI and with benzamidine. It is found that conformational variations are reflected in the side chain distance matrices much more than in the standard C-C distance representations.

Benzamidines↗

Beta-adrenergic activity and conformation of the antihypertensive specific alpha 2-agonist drug, guanabenz.

In recent research a new series of specific drugs, one of which is guanabenz (GBZ, 2,6(dichlorobenzyliden)-aminoguanidine) has been introduced into the clinical treatment of centrally mediated hypertension. Guanabenz (GBZ) is considered to be among the most specific alpha 2-adrenergic agonists, acting similarly to clonidine by decreasing the sympathetic outflow from the brain to the peripheral circulatory system. In the present report we show that GBZ displays a significant affinity for beta-adrenoceptors. In displacement studies of the iodinated beta-antagonist [125I]cyanopindolol (CYP) from turkey erythrocyte membranes, the dissociation constant of GBZ was 3.8 microM. Inhibition of the (-) epinephrine induced adenylate cyclase activity by GBZ is competitive, with an apparent dissociation constant of 30 microM. A similar value was obtained by studies of GBZ's effect on the (-) epinephrine-induced [3H]cAMP accumulation in intact turkey erythrocytes. In view of its unexpected affinity for beta-adrenoceptors, we examined the three-dimensional structure of crystalline GBZ. In these studies substantial differences between clonidine and GBZ were observed, despite their strong structural resemblance. These dissimilarities (angle of rotation phi = 39.7 degrees as compared to 76 degrees in clonidine, and the rotational restriction of clonidine as compared to the greater mobility in rotation of GBZ) could explain the difference of specificity between these two compounds.

Adenylyl Cyclases↗

Metabolic activation and toxicity of acetaminophen and related analogs. A theoretical study.

Reaction thermodynamics have been calculated for an oxene model for cytochrome P-450 oxidations of four related arylamines: aniline, p-hydroxyaniline, acetanilide, and acetaminophen, by both radical and nonradical mechanisms, using a semiempirical molecular orbital method (modified neglect of differential overlap). The results indicate that for both p-hydroxyaniline and acetaminophen, a recently proposed peroxidase-like mechanism leading directly to p-benzoquinoneimines via radical intermediates is thermodynamically favored over N-hydroxylamine formation by H abstraction or addition rearrangement. These studies also provide a detailed characterization of three candidate species for the toxic reactive intermediate of acetaminophen: 1) p-benzoquinoneimines, 2) the radical intermediate formed by H abstraction from the nitrogen, and 3) the radical intermediate formed by H abstraction from the phenol. Calculated electron and spin densities indicate that the radical formed by H abstraction from the phenol oxygen does not remain localized on the oxygen, but is primarily a semiquinone aryl radical with significant unpaired spin density on the ring carbon atoms, particularly on C-3 and C-5. This result is consistent with the hyperfine splitting pattern observed for a transient radical species in a hydroxyl radical-mediated chemical oxidation of acetaminophen. The radical formed by H abstraction from the nitrogen also delocalizes on the ring carbons, but to a lesser extent and at the 2- and 4-positions. A closed shell mechanism of N oxidation of arylamines appears to lead directly to the hydroxylamines with less likelihood of precursor reactive intermediates. Toxic species could then be formed by loss of H2O from the hydroxylamines.

Acetaminophen↗

Structure-activity relationships of cholinesterase inhibitors. I. Quantum mechanical study of affinities of phenyl N-methyl carbamates.

Ab initio self-consistent field-molecular orbitals-linear combination of atomic orbitals calculations of electrostatic potentials were performed in a few regions about some phenyl N-methyl carbamate derivatives. A similar active-site conformation was assumed and a low-energy common conformation was obtained from semiempirical perturbative configuration interaction over localized orbitals calculations. The results indicate that affinity constants of the carbamate derivatives used for this study may be an outcome of phenyl ring electrostatic interactions with the active site.

Binding Sites↗