Search PubMed⌕ Search

Biomedical subjects

D G Covell

Publications and source records attributed to D G Covell.

At least 37 records · Page 2Linked to original sources

Analysis of protein-protein interactions and the effects of amino acid mutations on their energetics. The importance of water molecules in the binding epitope.

A modeling analysis has been conducted to assess the determinants of binding strength and specificity for three crystal complexes; the anti-hen egg white lysozyme antibody D1.3 complexed with hen egg white lysozyme (HEL), the D1.3 antibody complexed with the anti-lysozyme antibody E5.2, and barnase complexed with barstar. The strengths of individual binding components within these interfaces are evaluated using a model of binding free energy that is based on pairwise surface preferences. In all cases the energetics of binding are dominated by a relatively small number of interfacial residues that define the binding epitope. A precise geometric arrangement of these residues was not found; they were either localized to one region, or distributed throughout the binding interface. Surprisingly, interfacial crystal water molecules were calculated to contribute around 25% of the total calculated binding strength. Theoretical alanine mutations were completed by atomic deletions of the wild-type complexes. Strong correlations were observed between the calculated changes in binding free energy (deltadeltaG(calculated)) and the experimental values (deltadeltaG(observed)) for all but three of the 30 single residue mutations in the D1.3-HEL, D1.3-E5.2 and barnase-barstar systems and for all of the double mutations in the barnase-barstar system. This analysis finds that the observed differences in binding strength are consistent with a model that accounts for the changes in binding energy from the direct contacts between each member of the complex and indirect changes due to released crystallographic water molecules that are near the mutation site. The observed energy changes for double mutations in the barnase-barstar system is fully accounted for by considering water molecules bound jointly by each member of the complex.

Amino Acids↗

Differential contact of disparate class I/peptide complexes as the basis for epitope cross-recognition by a single T cell receptor.

In an effort to better understand functional recognition of structurally dissimilar ligands by a single TCR, a model system for studying cross-recognition of disparate peptide/class I complexes was developed using the murine (H-2b) CTL clone AHIII12.2, which is reactive to a human self-peptide (p1049) bound to an HLA-A2.1 molecule. We identified a second complex comprised of a synthetic peptide, designated p1058, bound to H-2Db that is recognized by clone AHIII12.2. In cytolysis assays, dose-response profiles for peptides p1049 and p1058 pulsed onto the appropriate target cells were comparable, suggesting that p1049/A2.1 and p1058/Db form functionally equivalent epitopes. To probe the interaction between each complex and the TCR of AHIII12.2, singly substituted analogues of each peptide were tested for their activity in lysis assays. Differences were observed between the two epitopes with respect to permissible residue substitutions at each peptide position from P3 to P8; marked differences were evident at P3 and at P8. The results obtained suggest that this TCR forms critical contacts with atoms at peptide positions P3 and P5 of p1049/A2.1 and at P5 and P8 of p1058/Db, and that TCR cross-recognition of these ligands is a function of both shared and complex-specific contacts made with each epitope. These findings further highlight the versatile reactivity that may be shown by a single TCR and suggest a basis for the recognition of peptide ligands sharing only a limited set of structural features.

Animals↗

Identification of cooperative folding units in a set of native proteins.

Cooperative unfolding penalties are calculated by statistically evaluating an ensemble of denatured states derived from native structures. The ensemble of denatured states is determined by dividing the native protein into short contiguous segments and defining all possible combinations of native, i.e., interacting, and non-native, i.e., non-interacting, segments. We use a novel knowledge-based scoring function, derived from a set of non-homologous proteins in the Protein Data Bank, to describe the interactions among residues. This procedure is used for the structural identification of cooperative folding cores for four globular proteins: bovine pancreatic trypsin inhibitor, horse heart cytochrome c, French bean plastocyanin, and staphylococcal nuclease. The theoretical folding units are shown to correspond to regions that exhibit enhanced stability against denaturation as determined from experimental hydrogen exchange protection factors. Using a sequence similarity score for related sequences, we show that, in addition to residues necessary for enzymatic function, those amino acids comprising structurally important folding cores are also preferentially conserved during evolution. This implies that the identified folding cores may be part of an array of fundamental structural folding units.

Amino Acid Sequence↗

Azodicarbonamide inhibits HIV-1 replication by targeting the nucleocapsid protein.

Nucleocapsid p7 (NCp7) proteins of human immunodeficiency virus type 1 (HIV-1) contain two zinc binding domains of the sequence Cys-(X)2-Cys-(X)4-His-(X)4-Cys (CCHC). The spacing pattern and metal-chelating residues (3 Cys, 1 His) of these nucleocapside CCHC zinc fingers are highly conserved among retroviruses. These CCHC domains are required during both the early and late phases of retroviral replication, making them attractive targets for antiviral agents. toward that end, we have identified a number of antiviral chemotypes that electrophilically attack the sulfur atoms of the zinc-coordinating cysteine residues of the domains. Such nucleocapside inhibitors were directly virucidal by preventing the initiation of reverse transcription and blocked formation of infectious virus from cells through modification of CCHC domains within Gag precursors. Herein we report that azodicarbonamide (ADA) represents a new compound that inhibits HIV-1 and a broad range of retroviruses by targeting the the nucleocapsid CCHC domains. Vandevelde et al. also recently disclosed that ADA inhibits HIV-1 infection via an unidentified mechanism and that ADA was introduced into Phase I/II clinical trials in Europe for advanced AIDS. These studies distinguish ADA as the first known nucleocapsid inhibitor to progress to human trials and provide a lead compound for drug optimization.

Anti-HIV Agents↗

Evaluation of selected chemotypes in coupled cellular and molecular target-based screens identifies novel HIV-1 zinc finger inhibitors.

Conservation of the Cys-Xaa2-Cys-Xaa4-His-Xaa4-Cys retroviral zinc finger sequences and their absolute requirement in both the early and late phases of retroviral replication make these chemically reactive structures prime antiviral targets. We recently reported that select 2,2'-dithiobisbenzamides (DIBAs) chemically modify the zinc finger Cys residues, resulting in release of zinc from the fingers and inhibition of HIV replication. In the current study we surveyed 21 categories of disulfide-based compounds from the chemical repository of the National Cancer Institute for their capacity to act as retroviral zinc finger inhibitors. Aromatic disulfides that exerted anti-HIV activity tended to cluster in the substituted aminobenzene, benzoate, and benzenesulfonamide disulfide subclasses. Only one thiuram derivative exerted moderate anti-HIV activity, while a number of nonaromatic thiosulfones and miscellaneous disulfide congeners were moderately antiviral. Two compounds (NSC 20625 and NSC 4493) demonstrated anti-cultures. The two compounds chemically modified the p7NC zinc fingers in two separate in vitro assays, and interatomic surface molecular modeling docked the compounds efficiently but differentially into the zinc finger domains. The combined efforts of rational drug selection, cell-based screening, and molecular target-based screening led to the identification of zinc finger inhibitors that can now be optimized by medicinal chemistry for the development of biopharmaceutically useful anti-HIV agents.

Anti-HIV Agents↗

Docking enzyme-inhibitor complexes using a preference-based free-energy surface.

We present a docking scheme that utilizes both a surface complementarity screen as well as an energetic criterion based on surface area burial. Twenty rigid enzyme/inhibitor complexes with known coordinate sets are arbitrarily separated and reassembled to an average all-atom rms (root mean square) deviation of 1.0 A from the native complexes. Docking is accomplished by a hierarchical search of geometrically compatible triplets of surface normals on each molecule. A pruned tree of possible bound configurations is built up using successive consideration of larger and larger triplets. The best scoring configurations are then passed through a free-energy screen where the lowest energy member is selected as the predicted native state. The free energy approximation is derived from observations of surface burial by atom pairs across the interface of known enzyme/inhibitor complexes. The occurrence of specific atom-atom surface burial, for a set of complexes with well-defined secondary structure both in the bound and unbound states, is parameterized to mimic the free energy of binding. The docking procedure guides the inhibitor into its native state using orientation and distance-dependent functions that reproduce the ideal model of free energies with an average rms deviation of 0.9 kcal/mol. For all systems studied, this docking procedure identifies a single, unique minimum energy configuration that is highly compatible with the native state.

Enzyme Inhibitors↗

On the origins of the hydrophobic effect: observations from simulations of n-dodecane in model solvents.

The importance of the small size of a water molecule as contributing to the hydrophobic effect is examined from simulations of n-dodecane in different solvents. The earlier observations of the origin of hydrophobicity, derived from cavity formations by Pratt and Pohorille (1992, Proc. Natl. Acad. Sci. USA. 89:2995-2999) and Madan and Lee (1994, Biophys. Chem, 51:279-289), are shown to be largely consistent for a hydrocarbon-induced water pocket. In effect, the small size of a water molecule limits the probability (and hence free energy) of finding an appropriate void in the fluid that will accommodate a solute. In this work a simulated collapse of an n-dodecane molecule in H2O, CCl4, and a water-like Lennard-Jones solvent indicates that the induced entropy and enthalpy changes are qualitatively similar for hydrogen-bonded and Lennard-Jones water solvents. These results suggest that a large part of the hydrophobic response of solutes in aqueous solutions is due to the small size of the solvent. Important quantitative differences between the studied water solvents indicate that the hydrogen-bonded properties for water are still needed to determine the overall hydrophobic response.

Alkanes↗

A preference-based free-energy parameterization of enzyme-inhibitor binding. Applications to HIV-1-protease inhibitor design.

The interface between protein receptor-ligand complexes has been studied with respect to their binary interatomic interactions. Crystal structure data have been used to catalogue surfaces buried by atoms from each member of a bound complex and determine a statistical preference for pairs of amino-acid atoms. A simple free energy model of the receptor-ligand system is constructed from these atom-atom preferences and used to assess the energetic importance of interfacial interactions. The free energy approximation of binding strength in this model has a reliability of about +/- 1.5 kcal/mol, despite limited knowledge of the unbound states. The main utility of such a scheme lies in the identification of important stabilizing atomic interactions across the receptor-ligand interface. Thus, apart from an overall hydrophobic attraction (Young L, Jernigan RL, Covell DG, 1994, Protein Sci 3:717-729), a rich variety of specific interactions is observed. An analysis of 10 HIV-1 protease inhibitor complexes is presented that reveals a common binding motif comprised of energetically important contacts with a rather limited set of atoms. Design improvements to existing HIV-1 protease inhibitors are explored based on a detailed analysis of this binding motif.

Alcohols↗

Lattice model simulations of polypeptide chain folding.

Simulated annealing methods are applied to simple cubic lattice C alpha models of eight small monomeric globular proteins and their transition from a random chain to a low energy compact state is examined. The lowest energy structures are compared to their crystal forms using coordinate distance deviations, dRMS and RMS, and by distance contact maps. Analysis of the transition region indicates that, for this model, collapse begins with a rapid decline in radius of gyration followed continuously by chain repackings that lead to progressively lower values of chain energy. Chain repackings represent a highly cooperative interplay between the formation of local and non-local interactions. The components of this transition are characterized by rapid relaxation of shorter chain segments to form local contacts and slower relaxations of longer chain segments to form non-local contacts. Final structures obtained with this procedure contain many of the gross topologies of their native structures.

Computer Simulation↗

A role for surface hydrophobicity in protein-protein recognition.

The role of hydrophobicity as a determinant of protein-protein interactions is examined. Surfaces of apo-protein targets comprising 9 classes of enzymes, 7 antibody fragments, hirudin, growth hormone, and retinol-binding protein, and their associated ligands with available X-ray structures for their complexed forms, are scanned to determine clusters of surface-accessible amino acids. Clusters of surface residues are ranked on the basis of the hydrophobicity of their constituent amino acids. The results indicate that the location of the co-crystallized ligand is commonly found to correspond with one of the strongest hydrophobic clusters on the surface of the target molecule. In 25 of 38 cases, the correspondence is exact, with the position of the most hydrophobic cluster coinciding with more than one-third of the surface buried by the bound ligand. The remaining 13 cases demonstrate this correspondence within the top 6 hydrophobic clusters. These results suggest that surface hydrophobicity can be used to identify regions of a protein's surface most likely to interact with a binding ligand. This fast and simple procedure may be useful for identifying small sets of well-defined loci for possible ligand attachment.

Amino Acid Sequence↗

Analysis of hydrophobicity in the alpha and beta chemokine families and its relevance to dimerization.

The chemokine family of chemotactic cytokines plays a key role in orchestrating the immune response. The family has been divided into 2 subfamilies, alpha and beta, based on the spacing of the first 2 cysteine residues, function, and chromosomal location. Members within each subfamily have 25-70% sequence identity, whereas the amino acid identity between members of the 2 subfamilies ranges from 20 to 40%. A quantitative analysis of the hydrophobic properties of 11 alpha and 9 beta chemokine sequences, based on the coordinates of the prototypic alpha and beta chemokines, interleukin-8 (IL-8), and human macrophage inflammatory protein-1 beta (hMIP-1 beta), respectively, is presented. The monomers of the alpha and beta chemokines have their strongest core hydrophobic cluster at equivalent positions, consistent with their similar tertiary structures. In contrast, the pattern of monomer surface hydrophobicity between the alpha and beta chemokines differs in a manner that is fully consistent with the observed differences in quaternary structure. The most hydrophobic surface clusters on the monomer subunits are located in very different regions of the alpha and beta chemokines and comprise in each case the amino acids that are buried at the interface of their respective dimers. The theoretical analysis of hydrophobicity strongly supports the hypothesis that the distinct dimers observed for IL-8 and hMIP-1 beta are preserved for all the alpha and beta chemokines, respectively. This provides a rational explanation for the lack of receptor crossbinding and reactivity between the alpha and beta chemokine subfamilies.

Amino Acid Sequence↗

Conformations of t-RNA: base pairing and stacking.

The Phe t-RNA structure can be fit with one point per nucleotide to lattice models, and a fit for the 76 points to a face-centered cubic lattice is achieved with an RMS of 1.76 A. There are 32 chain folds possible upon these points. Because it is impossible to calculate directly all combinations of potential base pairs for these cases, an alternative is to determine low energy secondary structures and subsequently the tertiary pairs. For each lattice fold, the low energy secondary structures are generated from a list of proximal bases. From the lists of remaining possible tertiary pairs, all combinations are generated, and these include 2,365,440 allowed conformers. Among the possible types of non-native conformational variations observed is slip pairing, accompanied by a bulge, at the end of a stem. Small changes in secondary structure can result in different tertiary pairs. Other calculations, not constrained to the t-RNA shape, are presented that involve the packing of rigid stems on a flexible internal loop. For a simple cubic lattice there are 36,484,128 lattice folds for the sixteen bases enclosing the internal loop. By attaching rigid stems and accounting for their excluded volume these are reduced to only 258,979 possible configurations. The most common stacking arrangements involve the usual two pairs of stacked stems indicated in the crystal structure. The present enumerations suggest that a completely thorough exploration of three dimensional RNA structures is feasible only with prior specification of restrictions on conformational freedom, such as those given by secondary structures.

Base Composition↗

Folding protein alpha-carbon chains into compact forms by Monte Carlo methods.

A method is presented for generating folded chains of specific amino acid sequences on a simple cubic lattice. Monte Carlo simulations are used to transform extended geometries of simplified alpha-carbon chains for eight small monomeric globular proteins into folded states. Permitted chain transitions are limited to a few types of moves, all restricted to occur on the lattice. Crude residue-residue potentials derived from statistical structure data are used to describe the energies for each conformer. The low resolution structures obtained by this procedure contain many of the correct gross features of the native folded architectures with respect to average residue energy per nonbonded contact, segment density, and location of surface loops and disulfide pairs. Rms deviations between these and the native X-ray structures and percentage of native long-range contacts found in these final folded structures are 7.6 +/- 0.7 A and 48 +/- 3%, respectively. This procedure can be useful for predicting approximate tertiary interactions from amino acid sequence.

Models, Molecular↗

Conformations of folded proteins in restricted spaces.

A new method is presented to examine the complete range of folded topologies accessible in the compact state of globular proteins. The procedure is to generate all conformations, with volume exclusion, upon a lattice in a space restricted to the individual protein's known compact conformational space. Using one lattice point per residue, we find 10(2)-10(4) possible compact conformations for the five small globular proteins studied. Subsequently, these conformations are evaluated in terms of residue-specific, pairwise contact energies that favor nonbonded, hydrophobic interactions. Native structures for the five proteins are always found within the best 2% of all conformers generated. This novel method is simple and general and can be used to determine a small group of most favorable overall arrangements for the folding of specific amino acid sequences within a restricted space.

Animals↗

Pharmacokinetic analysis of blood distribution of intravenously administered 153Gd-labeled Gd(DTPA)2- and 99mTc(DTPA) in rats.

Rat plasma distribution data obtained following IV administration of 99mTc(DTPA) alone or after co-administration of 99mTc(DTPA) and 153Gd-labeled Gd(DTPA)2- at 0.001, 0.1, and 1.0 mmol Gd/kg were evaluated using compartmental modeling techniques. A three-compartment open model was found to fit the data significantly better (P less than 0.01) than a two- or four-compartment open model. This model incorporates and links the plasma and urine data and includes a delay to account for the transit time through the kidneys/ureters. The two nonplasma compartments of the model were assumed to be related to rapidly and slowly equilibrating tissues. Tc(DTPA) and Gd(DTPA)2- had nearly identical pharmacokinetic profiles in plasma and the rate constants were essentially the same. No significant dose dependent pharmacokinetic differences were found for the range of Gd(DTPA)2- doses tested. Simulations of the proposed three-compartment model were used to generate concentration-time curves for each of the three compartments.

Animals↗

Conformational analysis of the tachykinins in solution: substance P and physalaemin.

A determination of the solution conformational behavior of two tachykinins, substance P and physalaemin, is described. Two-dimensional homonuclear Hartmann-Hahn (HOHAHA) and rotating-frame cross relaxation spectroscopy (ROESY) are used to obtain complete proton resonance assignments. Interproton distance restraints obtained from ROESY spectroscopy are used to characterize the conformational behavior. These data show that in solution both substance P and physalaemin exist in a mixture of conformational states, rather than as a single three-dimensional structure. In water both peptides prefer to be in an extended chain structure. In methanol, their behavior is described as a mixture of beta-turn conformations in dynamic equilibrium. Solvent titration data and chemical shift temperature coefficients complement the NMR estimate of interproton distances by locating hydrogen bonds and serving to identify predominant conformational states. The C-terminal tetrapeptide segment has the same conformational behavior for both substance P and physalaemin. In physalaemin, the midsegment of the peptide may also be constrained by formation of a salt bridge. The conformational behavior of substance P and physalaemin is discussed in relation to potency and receptor binding properties.

Magnetic Resonance Spectroscopy↗

Recent studies of human calcium metabolism using stable isotopic tracers.

Stable isotopes of calcium are used safely as tracers for calcium in human populations ranging in age from infants to postmenopausal women. Thermal ionization mass spectrometry is used to measure calcium isotope ratios with relative accuracies of about 1% for natural abundance ratios at precisions of about 1% relative to the mean. Perturbations of natural abundance ratios are determined for the calcium in blood, urine, and feces with a limit of detection of about 2 delta % excess. The mathematical rationale for clinical studies of fractional absorption of dietary calcium and the kinetics of calcium's internal distribution are presented.

Animals↗

A modeling analysis of monoclonal antibody percolation through tumors: a binding-site barrier.

For successful use of radiolabeled monoclonal antibodies (MAbs) for diagnosis and therapy, it is helpful to understand both global and microscopic aspects of antibody biodistribution. In this study, antibody distribution in a tumor is simulated by splicing together information on global pharmacokinetics: transport across the capillary wall, diffusive penetration through the tumor interstitial space, and antigen-antibody interaction. The geometry simulated corresponds to spherical nodules of densely packed tumor cells. This modeling analysis demonstrates that: 1) antigen-antibody binding in tumors can retard antibody percolation; 2) high antibody affinity at a given dose tends to decrease antibody percolation because there are fewer free antibody molecules. The result is a more heterogeneous distribution; 3) the average antibody concentration in the tumor does not increase linearly with affinity; and 4) increasing antibody dose leads to better percolation and more uniform distribution. This mathematical model and the general principles developed here can be applied as well to other biologic ligands.

Antibodies, Monoclonal↗