Proteomics and systems biology: application in drug discovery and development.

Studies of complex biological systems aimed at understanding their functions at a global level are the goals of systems biology. Proteomics, generally regarded as the comprehensive study of the expression of all the proteins at a particular time in different organs, tissues, and cell types is a key enabling technology for the systems biology approach. Rapid advances in this regard have been made following the success of the human genome project as well as those of various animals and microorganisms. Possibly, one of the most promising outcomes from studies on the human genome and proteome is the identification of potential new drugs for the treatment of different diseases and tailoring the drugs for individualized patient therapy. Following the identification of a new drug candidate, knowledge on organ and system-level responses helps prioritize the drug targets and design clinical trials based on their efficacy and safety. Toxicoproteomics is playing an important role in that respect. In essence, over the past decade, proteomics has played a major role in drug discovery and development. In this review article, we explain systems biology, discuss the current proteomic technologies, and highlight some important applications of proteomics and systems biology approaches in drug discovery and development.

Comparison of SYPRO Ruby and Flamingo fluorescent stains for application in proteomic research.

Fluorescent dyes are widely used for the detection and quantitation of proteins separated by polyacrylamide gel electrophoresis. SYPRO Ruby is one such fluorescent dye widely used for this purpose. More recently, another fluorescent dye, Flamingo, is available for expression proteomic research. Using a standard ultraviolet (UV) transilluminator and a charge-coupled device (CCD)-based imaging system, the relative sensitivity of these two different fluorescent stains with regard to detection of protein spots separated by two-dimensional gel electrophoresis (2D-GE) and identification by liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) were compared. Using mouse kidney and liver homogenates as well as Escherichia coli extract, we detected a greater number of protein spots using Flamingo compared with SYPRO Ruby. In addition, when we compared the number of matched peptides and the percentage of amino acid residues identified for 22 different protein spots of mouse kidney proteome, we observed a higher number of matched peptides and a higher percentage of amino acid residues for the majority of the proteins using Flamingo compared with SYPRO Ruby. Also, we were able to characterize a protein spot that can be detected by Flamingo only. Therefore, we recommend Flamingo over SYPRO Ruby to be used for studies on expression proteomics.

Structural study of Ac-Phe-[Orn-Pro-dCha-Trp-Arg], a potent C5a receptor antagonist, by NMR.

The cyclic hexapeptide Ac(0)-Phe(1)-[Orn(2)-Pro(3)-dCha(4)-Trp(5)-Arg(6)] (the square brackets denote cyclization) is a potent antagonist against C5a (the a-fragment of complement protein C5) binding to C5a receptor (C5aR) and an excellent candidate to become a therapeutic agent against diseases that involve unregulated activation of the complement system. We present the solution structure determination of this cyclic C5aR peptide antagonist (cC5aR-pa), using nuclear magnetic resonance (NMR) data and restrained molecular dynamics-based simulated annealing in torsion angle space with NMR-derived distance and torsion angle restraints. The calculated NMR ensemble of structures demonstrates the presence of a predominant conformation of a distorted type II' beta-turn in the segment Pro(3)-dCha(4)-Trp(5)-Arg(6). We critically examine the calculated structure with measured NMR parameters, such as nuclear Overhauser enhancement (NOE) connectivity patterns and intensities characteristic of specific structures, (3)J(H(N)-H(alpha)) scalar coupling constants, temperature coefficients for NH groups, and differences between observed chemical shifts and their random coil values. The raw NMR data are consistent with the presence of the type II' beta-turn, but also indicate the presence of conformational inter-conversion. The calculated three-dimensional coordinates for cC5aR-pa will form the basis for further computational studies and for the development of pharmacophore models.

Proteomic profiling of aging in the mouse heart: Altered expression of mitochondrial proteins.

Using two-dimensional gel electrophoresis and liquid chromatography-tandem mass spectrometry, we have used a systems biology approach to study the molecular basis of aging of the mouse heart. We have identified 8 protein spots whose expression is up-regulated due to aging and 36 protein spots whose expression is down-regulated due to aging (p0.05 as judged by Wilcoxon Rank Sum test). Among the up-regulated proteins, we have characterized 5 protein spots and 2 of them, containing 3 different enzymes, are mitochondrial proteins. Among the down-regulated proteins, we have characterized 27 protein spots and 16 of them are mitochondrial proteins. Mitochondrial damage is believed to be a key factor in the aging process. Our current study provides molecular evidence at the level of the proteome for the alteration of structural and functional parameters of the mitochondria that contribute to impaired activity of the mouse heart due to aging.

pH dependence of stability of the 10th human fibronectin type III domain: a computational study.

We present detailed computational studies based on electrostatic calculations to evaluate the origins of pKa values and the pH dependence of stability for the 10th type III domain of human fibronectin (FNfn10). One of our goals is to validate the calculation protocols by comparison to experimental data (Koide, A.; Jordan, M. R.; Horner, S.; Batori, V.; Koide, S. Biochemistry 2001, 40, 10326-10333). Another goal is to evaluate the sensitivity of the calculated ionization free energies and apparent pKa values on local structural fluctuations, which do not alter the structural convergence to a particular architecture, by using a complete ensemble of solution NMR structures and the NMR average minimized structure of FNfn10 (Main, A. L.; Harvey, T. S.; Baron, M.; Boyd, J.; Campbell, I. D. Cell 1992, 71, 671-678). Our calculations demonstrate that, at high ionic strength, FNfn10 is more stable at low pH compared to neutral pH, in overall agreement with experimental data. This behavior is attributed to contributions from unfavorable Coulombic interactions in a surface patch for the pairs Asp7-Glu9 and Asp7-Asp23. The unfavorable interactions are decreased at low pH, where the acidic residues become neutral, and are further decreased at high ionic strength because of increased screening by salt ions. Elimination of the unfavorable interactions in the theoretical mutants Asp7Asn (D7N) and Asp7Lys (D7K) produce higher calculated stabilities at neutral pH and any ionic strength compared to the wild-type, in agreement with the experimental data. We also discuss subtleties in the calculated apparent pKa values and ionization free energies, which are not in agreement with the experimental data. This work demonstrates that comparative electrostatic calculations can provide rapid predictions of pH-dependent properties of proteins and can be significant aids in guiding the design of proteins with tailored properties.

Computational studies of CXCR1, the receptor of IL-8/CXCL8, using molecular dynamics and electrostatics.

The three-dimensional structure of IL-8/CXCL8 has been previously determined using NMR spectroscopy and X-ray crystallography, but the structure of the receptors for this chemokine has not been determined experimentally. We present here the development of a model for the structure of the IL-8/CXCL8 receptor CXCR1, using a combination of homology modeling and a molecular dynamics simulation. Based on this model, we discuss the analysis of structural, dynamic, and physicochemical properties of CXCR1. We focused on the role of pairwise ionic interactions in local structural stability of CXCR1 and the role of electrostatic potentials in recognition of CXCR1 with IL-8/CXCL8. We have performed theoretical mutations of six charged amino acids in CXCR1, which abolish binding as suggested by earlier experimental data, to shed light on the effect of charge on association ability. We propose that the observed loss of binding in the six CXCR1 mutants is owed to loss of local structural stability, rather than hindrance of the recognition process because of changes in the overall electrostatic properties of the receptor. Based on further structural analysis, we propose some mutations of charged residues involving ion pairs in different elements of transmembrane helices and extracellular loops, which are expected to alter the local structure and possibly affect binding.

Immunophysical exploration of C3d-CR2(CCP1-2) interaction using molecular dynamics and electrostatics.

The formation of the complex between the d-fragment of the complement component C3 (C3d) and the modular complement receptor-2 (CR2) is important for cross-linking foreign antigens with surface-bound antibodies and C3d on the surface of B cells. The first two modules of CR2, complement control protein modules (CCPs), participate in non-bonded interactions with C3d. We have used computational methods to analyze the dynamic and electrostatic properties of the C3d-CR2(CCP1-2) complex. The interaction between C3d and CR2 is known to depend on pH and ionic strength. Also, the intermodular mobility of the CR2 modules has been questioned before. We performed a 10 ns molecular dynamics simulation to generate a relaxed structure from crystal packing effects for the C3d-CR2(CCP1-2) complex and to study the energetics of the C3d-CR2(CCP1-2) association. The MD simulation suggests a tendency for intermodular twisting in CR2(CCP1-2). We propose a two-step model for recognition and binding of C3d with CR2(CCP1-2), driven by long and short/medium-range electrostatic interactions. We have calculated the matrix of specific short/medium-range pairwise electrostatic free energies of interaction involved in binding and in intermodular communications. Electrostatic interactions may mediate allosteric effects important for C3d-CR2(CCP1-2) association. We present calculations for the pH and ionic strength-dependence of C3d-CR2(CCP1-2) ionization free energies, which are in overall agreement with experimental binding data. We show how comparison of the calculated and experimental data allows for the decomposition of the contributions of electrostatic from other effects in association. We critically compare predicted stabilities for several mutants of the C3d-CR2(CCP1-2) complex with the available experimental data for binding ability. Finally, we propose that CR2(CCP1-2) is capable of assuming a large array of intermodular topologies, ranging from closed V-shaped to open linear states, with similar recognition properties for C3d, but we cannot exclude an additional contact site with C3d.

Development of a quasi-dynamic pharmacophore model for anti-complement peptide analogues.

Three quasi-dynamic pharmacophore models have been constructed for the complement inhibitor peptide compstatin, using first principles. Uniform sampling along 5-ns molecular dynamics trajectories provided dynamic conformers that are thought to represent the entire conformational space for nine training set molecules, compstatin, four active analogues, and four inactive analogues. The pharmacophore models were built using mixed physicochemical and structural properties of residues indispensable for structural stability and activity. Owing to the size and flexibility of compstatin, one-dimensional probability distributions of intrapharmacophore point distances, angles, and dihedral angles of different analogues spread over wide and overlapping ranges. More robust two-dimensional distance-angle probability distributions for two pharmacophore models discriminated individual analogues in terms of specific distance-angle pairs, but overall failed to identify the active and the inactive analogues as two distinct groups. Two-dimensional distance-dihedral angle probability distributions in a third pharmacophore model allowed discrimination of the groups of active and inactive analogues more effectively, with the highest-activity analogue having distinct behavior. The present study indicates that more stringent structural constraints should be used for a set of structurally similar but flexible peptides, as opposed to organic molecules, to convert dynamic conformers into pharmacophore models. Flexibility is a general aspect of the structure and function of peptides and should be taken into account in ligand-based pharmacophore design. However, the discrimination of activity using multidimensional probability surfaces depends on the peptide system, the selection of the training set, the molecular dynamics protocol, and the selection of the type and number of pharmacophore points.

Design and NMR characterization of active analogues of compstatin containing non-natural amino acids.

We present new findings in our drug discovery effort to develop an anticomplement therapeutic. We have designed several active analogues of compstatin by altering its amino acid composition at positions 4 and 9. The most effective analogues have tryptophan or fused-ring non-natural amino acids at position 4 and alanine or an unbranched single-methyl amino acid at position 9. Twenty-one of these analogues have 2-99-fold higher activities compared to the parent peptide compstatin. The analogue Ac-V4(2Nal)/H9A-NH(2) has the highest inhibitory activity with IC(50) 500 nM. NMR data, through NOE and chemical shift analysis, suggest the presence of interconverting conformers spanning the extended and helical regions of the Ramachandran plot, and they detect a predominant averaged conformer with coil structure and at least one flexible beta-turn, of type I. The fused-ring non-natural amino acids at position 4 contribute to the formation of the hydrophobic cluster of compstatin, which has been previously proposed, together with the beta-turn and a disulfide bridge, to be essential for binding to the target of compstatin, complement component C3. We propose that additional mechanisms may contribute to the structural stability of the analogues and to binding to C3, involving intra- and intermolecular electrostatic interactions of the pi-electron system of side chain aromatic rings. The presence of pi-pi interactions for Trp4-Trp7 was confirmed with a molecular dynamics simulation for the most active analogue with natural amino acids, Ac-V4W/H9A-NH(2). Alanine or aminobutyric acid at position 9 contribute to the weak propensity for helical structure of the residue segment 4-10 of the analogues, which may also play a role in increased activity.

Implicit solvent simulations of DPC micelle formation.

The formation of micelles by dodecylphosphocholine (DPC) is modeled by treating the surfactants in atomic detail and the solvent implicitly, in the spirit of the EEF1 solvation model for proteins. The solvation parameters of the DPC atoms are carried over from those of similar atoms in proteins. A slight adjustment of the parameters for the headgroup was found necessary for obtaining an aggregation number consistent with experiment. Molecular dynamics simulations of 960 DPC molecules at different concentrations are used to obtain the aggregation number, the micelle size distribution, and the CMC. At 20 mM concentration we obtain an aggregation number of 53-56 and a CMC of 1.25 mM, values close to the experimental ones. At 100 mM the aggregation number increases to 90. Simulations of individual micelles of varying size show that the effective energy per surfactant molecule is initially a decreasing function of aggregation number but stabilizes at about 60 molecules. The van der Waals term and the desolvation of nonpolar groups contribute to micellization, whereas the desolvation of polar groups opposes it. From the difference between the effective energy and the free energy (calculated from the CMC), the translational and rotational entropy contributions to the free energy are estimated at about 7 kcal/mol per monomer. The micelles obtained here are more irregular than those obtained in explicit water simulations. This modeling approach allows the study of larger surfactant aggregates for longer times and the extraction of thermodynamic in addition to structural information.

Conformational interconversion in compstatin probed with molecular dynamics simulations.

Compstatin is a 13-residue cyclic peptide that has the potential to become a therapeutic agent against unregulated complement activation. In our effort to understand the structural and dynamic characteristics of compstatin that form the basis for rational and combinatorial optimization of structure and activity, we performed 1-ns molecular dynamics (MD) simulations. We used as input in the MD simulations the ensemble of 21 lowest energy NMR structures, the average minimized structure, and a global optimization structure. At the end of the MD simulations we identified five conformations, with populations ranging between 9% and 44%. These conformations are as follows: 1) coil with alphaR-alphaR beta-turn, as was the conformation of the initial ensemble of NMR structures; 2) beta-hairpin with epsilon-alphaR beta-turn; 3) beta-hairpin with alphaR-alphaR beta-turn; 4) beta-hairpin with alphaR-beta beta-turn; and 5) alpha-helical. Conformational switch was possible with small amplitude backbone motions of the order of 0.1-0.4 A and free energy barrier crossing of 2-11 kcal/mol. All of the 21 MD structures corresponding to the NMR ensemble possessed a beta-turn, with 14 structures retaining the alphaR-alphaR beta-turn type, but the average minimized structure and the global optimization structures were converted to alpha-helical conformations. Overall, the MD simulations have aided to gain insight into the conformational space sampled by compstatin and have provided a measure of conformational interconversion. The calculated conformers will be useful as structural and possibly dynamic templates for optimization in the design of compstatin using structure-activity relations (SAR) or dynamics-activity relations (DAR).

Distance and exposure dependent effective dielectric function.

In an effort to develop a dielectric screening function for molecular dynamics simulations of biomolecules in implicit solvent, effective dielectric constants (D(eff)) for a large number of atom pairs in a typical globular protein are calculated by continuum electrostatics. Plots of D(eff) versus the intercharge distance are in general sigmoidal with the characteristics of the curve depending on the distance of the two charges from the dielectric boundary and, secondarily, on the extent to which the area surrounding each charge is occupied by solvent (the "exposure"). The D(eff) values were fitted to an empirical, analytical function of these parameters that reproduces the data reasonably well, although considerable scatter exists in the range of D(eff) from 30 to 80. In the system used for parameterization, the mean square deviation of electrostatic interaction energies with this function is 0.48 kcal/mol, compared to 1.45 for an analytical Generalized Born model and 1.52 for the linear distance-dependent dielectric model. When tested in other proteins of varying size and compactness, the present function is superior to both of the above models, except for a fully unfolded polypeptide chain, where the Generalized Born model is superior.