Biomedical subjects

P Beroza

Publications and source records attributed to P Beroza.

10 recordsLinked to original sources

PubMed2000

A rapid computational method for lead evolution: description and application to alpha(1)-adrenergic antagonists.

The high failure rate of drugs in the development phase requires a strategy to reduce risks by generating lead candidates from different chemical classes. We describe a new three-dimensional computational approach for lead evolution, based on multiple pharmacophore hypotheses. Using full conformational models for both active and inactive compounds, a large number of pharmacophore hypotheses are analyzed to select the set or "ensemble" of hypotheses that, when combined, is most able to discriminate between active and inactive molecules. The ensemble hypothesis is then used to search virtual chemical libraries to identify compounds for synthesis. This method is very rapid, allowing very large virtual libraries on the order of a million compounds to be filtered efficiently. In applying this method to alpha(1)-adrenergic receptor ligands, we have demonstrated lead evolution from heterocyclic alpha(1)-adrenergic receptor ligands to highly dissimilar active N-substituted glycine compounds. Our results also show that the active N-substituted glycines are part of our smaller filtered library and thus could have been identified by synthesizing only a portion of the N-substituted glycine library.

Adrenergic alpha-Antagonists↗

PubMed2000

Designing chiral libraries for drug discovery.

Identification of new biological targets combined with combinatorial medicinal chemistry has enormous potential to facilitate the discovery of new therapeutics. However, these new targets might require the use of different or more complex screening libraries. This complexity might be achieved through the use of new templates or the identification of new chemistries. One approach is the use of chiral libraries that have been designed using three-dimensional pharmacophore and shape descriptors to provide maximal structure-activity information after screening. The advantages of using a designed chiral library in primary screening and the information one should obtain are discussed.

Journal Article↗

PubMed1998

Calculations of proton-binding thermodynamics in proteins.

Computational models of proton binding can range from the chemically complex and statistically simple (as in the quantum calculations) to the chemically simple and statistically complex. Much progress has been made in the multiple-site titration problem. Calculations have improved with the inclusion of more flexibility in regard to both the geometry of the proton binding and the larger scale protein motions associated with titration. This article concentrated on the principles of current calculations, but did not attempt to survey their quantitative performance. This is (1) because such comparisons are given in the cited papers and (2) because continued developments in understanding conformational flexibility and interaction energies will be needed to develop robust methods with strong predictive power. Nevertheless, the advances achieved over the past few years should not be underestimated: serious calculations of protonation behavior and its coupling to conformational change can now be confidently pursued against a backdrop of increasing understanding of the strengths and limitations of such models. It is hoped that such theoretical advances will also spur renewed experimental interest in measuring both overall titration curves and individual pKa values or pKa shifts. Exploration of the shapes of individual titration curves (as measured by Hill coefficients and other parameters) would also be useful in assessing the accuracy of computations and in drawing connections to functional behavior.

Hemoglobins↗

PubMed1997

Computational, pulse-radiolytic, and structural investigations of lysine-136 and its role in the electrostatic triad of human Cu,Zn superoxide dismutase.

Key charged residues in Cu,Zn superoxide dismutase (Cu,Zn SOD) promote electrostatic steering of the superoxide substrate to the active site Cu ion, resulting in dismutation of superoxide to oxygen and hydrogen peroxide, Lys-136, along with the adjacent residues Glu-132 and Glu-133, forms a proposed electrostatic triad contributing to substrate recognition. Human Cu,Zn SODs with single-site replacements of Lys-136 by Arg,Ala, Gln, or Glu or with a triple-site substitution (Glu-132 and Glu-133 to Gln and Lys-136 to Ala) were made to test hypotheses regarding contributions of these residues to Cu,Zn SOD activity. The structural effects of these mutations were modeled computationally and validated by the X-ray crystallographic structure determination of Cu,Zn SOD having the Lys-136-to-Glu replacement. Brownian dynamics simulations and multiple-site titration calculations predicted mutant reaction rates as well as ionic strength and pH effects measured by pulse-radiolytic experiments. Lys-136-to-Glu charge reversal decreased dismutation activity 50% from 2.2 x 10(9) to 1.2 x 10(9) M-1 s-1 due to repulsion of negatively charged superoxide, whereas charge-neutralizing substitutions (Lys-136 to Gln or Ala) had a less dramatic influence. In contrast, the triple-mutant Cu,Zn SOD (all three charges in the electrostatic triad neutralized) surprisingly doubled the reaction rate compared with wild-type enzyme but introduced phosphate inhibition. Computational and experimental reaction rates decreased with increasing ionic strength in all of the Lys-136 mutants, with charge reversal having a more pronounced effect than charge neutralization, implying that local electrostatic effects still govern the dismutation rates. Multiple-site titration analysis showed that deprotonation events throughout the enzyme are likely responsible for the gradual decrease in SOD activity above pH 9.5 and predicted a pKa value of 11.7 for Lys-136. Overall, Lys-136 and Glu-132 make comparable contributions to substrate recognition but are less critical to enzyme function than Arg-143, which is both mechanistically and electrostatically essential. Thus, the sequence-conserved residues of this electrostatic triad are evidently important solely for their electrostatic properties, which maintain the high catalytic rate and turnover of Cu,Zn SOD while simultaneously providing specificity by selecting against binding by other anions.

Binding Sites↗

PubMed1996

NMR solution structure of Cu(I) rusticyanin from Thiobacillus ferrooxidans: structural basis for the extreme acid stability and redox potential.

The solution structure of the Cu(I) form of the rusticyanin from Thiobacillus ferrooxidans has been calculated from a total of 1979 distance and dihedral angle constraints derived from 1H, 13C and 15N NMR spectra. The structures reveal two beta-sheets, one of six strands and one of seven strands that are tightly packed in a beta-barrel or beta-sandwich arrangement, and a short helix that extends on the outside of one of the sheets to form a second hydrophobic core. The copper coordination sphere is composed of the standard type I ligands (His2CysMet) in a distorted tetrahedral arrangement. The copper-binding site is located within a hydrophobic region at one end of the molecule, surrounded by a number of aromatic rings and hydrophobic residues. This configuration probably contributes to the acid stability of the copper site, since close association of the aromatic rings with the histidine ligands would sterically hinder their dissociation from the copper. An electrostatic analysis based on a comparison of the structures of rusticyanin and French bean plastocyanin shows that factors determining the high redox potential of rusticyanin include contributions from charged side-chains and from the disposition of backbone peptide dipoles, particularly in the 81 to 86 region of the sequence and the ligand cysteine residue. These interactions should also contribute to the acid stability by inhibiting protonation of His143.

Azurin↗

PubMed1996

Co-crystallization and characterization of the photosynthetic reaction center-cytochrome c2 complex from Rhodobacter sphaeroides.

The photosynthetic reaction center (RC) of Rhodobacter sphaeroides and cytochrome c2 (cyt c2), its physiological secondary electron donor, have been co-crystallized. The molar ratio of RC/cyt c2 was found by SDS-PAGE and optical absorbance changes in the co-crystals to be 4. The crystals diffracted X-rays to 3.5 angstroms. However, the resolution degraded during data collection. A data set, 82.5% complete, was collected to 4.5 angstroms. The crystals belong to the tetragonal space group P4(3)2(1)2, with unit cell dimensions of a = b = 142.7 angstroms and c = 254.8 angstroms. The positions of the RCs in the unit cell were determined by molecular replacement. A comparable search for the cyt c2 by this method was unsuccessful because of the small contribution of the cytochrome to the total scattering and because of its low occupancy. The cyt c2 was positioned manually into patches of difference electron density, adjacent to the periplasmic surface of the M polypeptide subunit of the RC. The difference electron density was not sufficient for precise positioning of the cyt c2, and its orientation was modeled by placing the exposed edge of the heme toward the primary donor of the reaction center D and by forming pairs for electrostatically interacting RC and cyt c2 amino acid residues. The RC-cyt c2 structure derived from the co-crystal data was supported by use of omit maps and structure refinement analyses. Cyt c2 reduces the photooxidized primary donor D+ in 0.9 +/- 0.1 micros in the co-crystals, which is the same as the fast electron transfer rate in vivo and in solution. This result provides strong evidence that the structure of the complex in the co-crystal is the same as in solution. Two additional methods were used to investigate the structure of the RC-cyt c2 complex: (i) Docking calculations based on interprotein electrostatic interactions identified possible binding positions of the cyt c2 on the RC. The cyt c2 position with the lowest electrostatic energy is very similar to that of the cyt c2 in the proposed co-crystal structure. (ii) Site-directed mutagenesis was used to modify two aspartic acid residues (M184 and L155) on the periplasmic surface of the RC. Cyt c2 binding affinity to these RCs and electron transfer rates to D+ in these RCs support the co-crystal structure of th RC-cyt c2 complex.

Binding Sites↗

PubMed1995

Pathways, pathway tubes, pathway docking, and propagators in electron transfer proteins.

The simplest views of long-range electron transfer utilize flat one-dimensional barrier tunneling models, neglecting structural details of the protein medium. The pathway model of protein electron transfer reintroduces structure by distinguishing between covalent bonds, hydrogen bonds, and van der Waals contacts. These three kinds of interactions in a tunneling pathway each have distinctive decay factors associated with them. The distribution and arrangement of these bonded and nonbonded contacts in a folded protein varies tremendously between structures, adding a richness to the tunneling problem that is absent in simpler views. We review the pathway model and the predictions that it makes for protein electron transfer rates in small proteins, docked proteins, and the photosynthetic reactions center. We also review the formulation of the protein electron transfer problem as an effective two-level system. New multi-pathway approaches and improved electronic Hamiltonians are described briefly as well.

Electron Transport↗

PubMed1995

Electrostatic calculations of amino acid titration and electron transfer, Q-AQB-->QAQ-B, in the reaction center.

The titration of amino acids and the energetics of electron transfer from the primary electron acceptor (QA) to the secondary electron acceptor (QB) in the photosynthetic reaction center of Rhodobacter sphaeroides are calculated using a continuum electrostatic model. Strong electrostatic interactions between titrating sites give rise to complex titration curves. Glu L212 is calculated to have an anomalously broad titration curve, which explains the seemingly contradictory experimental results concerning its pKa. The electrostatic field following electron transfer shifts the average protonation of amino acids near the quinones. The pH dependence of the free energy between Q-AQB and QAQ-B calculated from these shifts is in good agreement with experiment. However, the calculated absolute free energy difference is in severe disagreement (by approximately 230 meV) with the observed experimental value, i.e., electron transfer from Q-A to QB is calculated to be unfavorable. The large stabilization energy of the Q-A state arises from the predominantly positively charged residues in the vicinity of QA in contrast to the predominantly negatively charged residues near QB. The discrepancy between calculated and experimental values for delta G(Q-AQB-->QAQ-B) points to limitations of the continuum electrostatic model. Inclusion of other contributions to the energetics (e.g., protein motion following quinone reduction) that may improve the agreement between theory and experiment are discussed.

Amino Acids↗

PubMed1991

Protonation of interacting residues in a protein by a Monte Carlo method: application to lysozyme and the photosynthetic reaction center of Rhodobacter sphaeroides.

We used Monte Carlo methods to treat statistical problem of electrostatic interactions among many titrating amino acids and applied these methods to lysozyme and the photosynthetic reaction center of Rhodobacter sphaeroides, including all titrating sites. We computed the average protonation of residues as a function of pH from an equilibrium distribution of states generated by random sampling. Electrostatic energies were calculated from a finite difference solution to the linearized Poisson-Boltzmann equation using the coordinates from solved protein structures. For most calculations we used the Metropolis algorithm to sample protonation states; for strongly coupled sites, we substantially reduced sampling errors by using a modified algorithm that allows multiple site transitions. The Monte Carlo method agreed with calculations for a small test system, lysozyme, for which the complete partition function was calculated. We also calculated the pH dependence of the free energy change associated with electron transfer from the primary to the secondary quinone in the photosynthetic reaction center. The shape of the resulting curve agreed fairly well with experiment, but the proton uptake from which the free energy was calculated agreed only to within a factor of two with the observed values. We believe that this discrepancy resulted from errors in the individual electrostatic energy calculations rather than from errors in the Monte Carlo sampling.

Electron Transport↗

PubMedDate not supplied

Applications of random sampling to virtual screening of combinatorial libraries.

We describe statistical techniques for effective evaluation of large virtual combinatorial libraries (> 10(10) potential compounds). The methods described are used for computationally evaluating templates (prioritization of candidate libraries for synthesis and screening) and for the design of individual combinatorial libraries (e.g., for a given diversity site, reagents can be selected based on the estimated frequency with which they appear in products that pass a computational filter). These statistical methods are powerful because they provide a simple way to estimate the properties of the overall library without explicitly enumerating all of the possible products. In addition, they are fast and simple, and the amount of sampling required to achieve a desired precision is calculable. In this article, we discuss the computational methods that allow random product selection from a combinatorial library and the statistics involved in estimating errors from quantities obtained from such samples. We then describe three examples: (1) an estimate of average molecular weight for the several billion possible products in a four-component Ugi reaction, a quantity that can be calculated exactly for comparison; (2) the prioritization of four templates for combinatorial synthesis using a computational filter based on four-point pharmacophores; and (3) selection of reagents for the four-component Ugi reaction based on their frequency of occurrence in products that pass a pharmacophore filter.

Combinatorial Chemistry Techniques↗