Mechanism of androgen receptor antagonism by bicalutamide in the treatment of prostate cancer.

The androgen receptor (AR) plays a key role in regulating gene expression in a variety of tissues, including the prostate. In that role, it is one of the primary targets in the development of new chemotherapeutics for treatment of prostate cancer and the target of the most widely prescribed current drug, bicalutamide (Bcu), for this disease. In view of its importance, and the absence of a crystal structure for any antagonist--AR complex, we have conducted a series of molecular dynamics-based simulations of the AR--Bcu complex and quantum mechanical (QM) calculations of Bcu, to elucidate the structural basis for antagonism of this key target. The structures that emerge show that bicalutamide antagonizes AR by accessing an additional binding pocket (B-site) adjacent to the hormone binding site (HBS), induced by displacing helix 12. This distorts the coactivator binding site and results in the inactivation of transcription. An alternative equienergetic conformational state of bicalutamide was found to bind in an expanded hormone pocket without materially perturbing either helix 12 or the coactivator binding site. Thus, both the structural basis of antagonism and the mechanism underlying agonist properties displayed by bicalutamide in different environments may be rationalized in terms of these structures. In addition, the antagonist structure and especially the induced second site (B-site) provide a structural framework for the design of novel antiandrogens.

Ab initio protein folding.

Ab initio protein folding methods have been developing rapidly over the past few years and, at the last Critical assessment of methods of protein structure prediction (CASP) meeting, it was shown that important progress has been made in generating structure from sequence. Both methods based on statistical potentials and methods using physics-based potentials have shown improvements. Most current methods use statistics-based potentials and the development of these is ongoing. Additionally, the inclusion of multiple sequence data in the algorithms in order to aid in finding the native structure is a common theme. The use of physics-based potentials is less developed, which means that less progress has been made in understanding why a sequence forms a structure.

Improved ab initio predictions with a simplified, flexible geometry model.

Structure predictions for nine targets from the CASP3 meeting are presented and compared with the experimental structure. These predictions are made using the simplified flexible geometry representation of protein structure and potentials which mimic the physical forces involved in protein folding with no help from multiple sequences. The major differences from the CASP2 potentials are identified, and the prediction successes and failures are related to the underlying potentials. Target T0065 was successfully folded from an extended chain to 3.8 A CA RMS of the native structure. The quality of the secondary structure component of all predictions was significantly improved, averaging a Q3 of 64%, which was spread evenly across the fold types, all alpha, alpha/beta mixed, and mainly beta. A number of the other predictions had a spatial arrangement of the secondary structure segments close to native, although the major problem with most predictions was the over-extension of secondary structure segments which often coalesced independent segments together. For target T0056 the ab initio prediction was the closest to the native structure of all methods including threading according to the Hubbard RMS coverage plots.

Analysis of the predicted structures of domain 1 of protein g3 (T0030) and NK-lysin (T0042).

Structure predictions for two targets from the CASP2 meeting are presented and compared with the experimental structure. These predictions were made using a novel simplified flexible geometry representation of protein structure. The method uses potentials which mimic the physical forces involved in protein folding in a simplified representation of protein structure, and not by directly using data derived from statistical analyses of known protein structures. Additionally, the method is designed to work with a single protein sequence. The method was successful in generating reasonable protein-like structures, with mainly buried hydrophobic residues, exposed charges, and a good fraction of secondary structure. Specific details of the structure were remarkably close to the experimental structure. However, the overall fold in both cases was totally wrong. Some specific causes of this incorrect folding are suggested and a major improvement to the algorithm proposed.

Molecular dynamics: deciphering the data.

The dynamic behaviour of molecules is important in determining their activity. Molecular dynamics (MD) simulations give a detailed description of motion, from small fluctuations to conformational transitions, and can include solvent effects. However, extracting useful information about conformational motion from a trajectory is not trivial. We have used digital signal-processing techniques to characterise the motion in MD simulations, including: calculating the frequency distribution, applying filtering functions, and extraction of vectors defining the characteristic motion for each frequency in an MD simulation. We describe here some typical results obtained for peptides and proteins. The nature of the low-frequency modes of motion, as obtained from MD and normal mode (NM) analysis, of Ace-(Ala)31-Nma and of a proline mutant is discussed. Low-frequency modes extracted from the MD trajectories of Rop protein and phospholipase A2 reveal characteristic motions of secondary structure elements, as well as concerned motions that are of significance to the protein's biological activity. MD simulations are also used frequently as a tool for conformational searches and for investigating protein folding/unfolding. We have developed a novel method that uses time-domain filtering to channel energy into conformational motion and thus enhance conformational transitions. The selectively enhanced molecular dynamics method is tested on the small molecule hexane.

FOCUS: a program for analyzing molecular dynamics simulations, featuring digital signal-processing techniques.

FOCUS is a program for analyzing molecular dynamics simulations. It enables the researcher to monitor structural and energetic properties during the trajectory, and to calculate the corresponding statistical averages, correlation functions and Fourier transforms. In addition to these conventional methods, the program also utilizes novel methods based on digital signal-processing techniques to characterize the various motions. The characteristic frequencies in the system are revealed by the frequency distribution function g(v), which is calculated from the Fourier transform of the atomic coordinates. A filtering technique is employed to remove uninteresting motion (e.g., high-frequency bond stretching) while retaining and focusing on important motion (e.g., low-frequency conformational motion). The filtering technique enables fast display of slow events without getting a blurry or jittery picture due to the high-frequency motions. Another new way for analyzing the motion is by extracting "characteristic modes" and associated frequencies. This yields a pictorial description of the oscillatory motions in a manner analogous to normal mode analysis.

Modeling of substrate and inhibitor binding to phospholipase A2.

Molecular graphics and molecular mechanics techniques have been used to study the mode of ligand binding and mechanism of action of the enzyme phospholipase A2. A substrate-enzyme complex was constructed based on the crystal structure of the apoenzyme. The complex was minimized to relieve initial strain, and the structural and energetic features of the resultant complex analyzed in detail, at the molecular and residue level. The minimized complex was then used as a basis for examining the action of the enzyme on modified substrates, binding of inhibitors to the enzyme, and possible reaction intermediate complexes. The model is compatible with the suggested mechanism of hydrolysis and with experimental data about stereoselectivity, efficiency of hydrolysis of modified substrates, and inhibitor potency. In conclusion, the model can be used as a tool in evaluating new ligands as possible substrates and in the rational design of inhibitors, for the therapeutic treatment of diseases such as rheumatoid arthritis, atherosclerosis, and asthma.

BIOSITE: a program for the interactive comparison of aligned homologous protein sequences.

A program, BIOSITE, providing for the interactive visual comparison of aligned homologous amino-acid sequences is presented, including an example of its application. The program allows for two types of comparison sequence to be generated: an 'identity' sequence and a 'difference' sequence. These may be used on subsets of sequences and in further comparisons to identify candidate sites involved in a distinct functional property. The program should prove a useful tool for biologists engaged in understanding sequence--function relationships.

Relative-residue surface-accessibility patterns reveal myoglobin and catalase similarity.

A novel sliding-window search method using relative-residue surface-accessibility patterns identified extensive, but unsuspected, structural similarity over a 3-helix region in the C-terminus of the evolutionarily unrelated proteins sperm-whale myoglobin and beef liver catalase. This clear example of structural similarity between non-homologous proteins highlights the importance of relative-residue surface-accessibility patterns in understanding the local folded structure in proteins.

Conformational analysis of peptide surrogates. Reduced and retro-amide links in blocked alanine and in secondary structures.

We have investigated the conformational effects of modifying the amide link of peptides. We studied a reverse amide bond psi [NHCO], a reduced amide bond psi [CH2NH] and a retro-reduced bond psi [NHCH2] as surrogates for the amide link [CONH] in native peptides. A complete search of the conformational space available to residues with these modified links was carried out. The local minima and the rotational barriers were described and compared to the minima of the native residue. The results are compatible with the available observed structural data. These modified links have been incorporated in secondary structure units such as beta turns, alpha helices, and parallel and anti-parallel beta sheets. It was found that a reduced amide link can lead to stabilised beta turns, while the retro modification can be incorporated in stable beta sheets. A significant reduction in the stability of alpha helices is caused by the retro links, while a reduced amide link results in only a small destabilisation.

Structure-activity studies with fragments and analogous of salmonid melanin-concentrating hormone.

A number of cyclic and linear fragments and analogues of MCH were synthesized and their biological potencies tested using the isolated carp scale melanophore assay. In this system the cyclic portion MCH(5-14) exhibited only 0.1% bioactivity, which was markedly enhanced by the addition of the exocyclic sequences MCH(15-17) and MCH(1-4). The exocyclic sequence itself, MCH(1-4,15-17), had minimal activity, however. Substitution of Tyr11 with phenylalanine reduced the potency of the ring structure MCH(5-14) by about 4-fold. Substitution of Gly8 with D-alanine reduced the potency of MCH(5-14) 16-fold, while both substitutions together caused a still more marked reduction (200-fold) in bioactivity. Linearized fragments of MCH, extending from MCH(15-17) to [Cys(Acm)5,14]MCH(1-17), showed a progressive increase in potency. The linearized forms of MCH, MCH(5-17) and MCH(5-14), were approximately 100-fold or less potent than their cyclic forms. The significant increases in bioactivity produced by the addition of the C- and N-terminal exocyclic sequence even to these linearized forms further emphasizes the importance of these regions for interaction at the receptor site.

Extraction of the energetics of selected types of motion from molecular dynamics trajectories by filtering.

A novel method for analyzing molecular dynamics trajectories has been developed which enables the study of selected motions and the corresponding energetics. In particular, it is possible to filter out the high-frequency motions and focus on the structural and energetic features of low-frequency collective motions. The trajectories of the properties of interest are Fourier transformed to the frequency domain, a filtering function is applied, and then an inverse transformation back to the time domain yields the filtered trajectory. The method is demonstrated for harmonic fluctuations and conformational transitions of acetamide and N-acetylalanine N-methylamide, as models for peptides and proteins.

The conformational preferences of gamma-lactam and its role in constraining peptide structure.

The conformational constraints imposed by gamma-lactams in peptides have been studied using valence force field energy calculations and flexible geometry maps. It has been found that while cyclisation restrains the psi of the lactam, non-bonded interactions contribute to the constraints on psi of the lactam. The gamma-lactam also affects the (psi, psi) of the residue after it in a peptide sequence. For an L-lactam, the ring geometry restricts psi to about -120 degrees, and psi has two minima, the lowest energy around -140 degrees and a higher minimum (5 kcal/mol higher) at 60 degrees, making an L-gamma-lactam more favourably accommodated in a near extended conformation than in position 2 of a type II' beta-turn. The energy of the psi approximately +60 degrees minimum can be lowered substantially until it is more favoured than the -140 degrees minimum by progressive substitution of bulkier groups on the amide N of the L-gamma-lactam. The (psi, psi) maps of the residue succeeding a gamma-lactam show subtle differences from those of standard N-methylated residues. The dependence of the constraints on the chirality of gamma-lactams and N-substituted gamma-lactams, in terms of the formation of secondary structures like beta-turns is discussed and the comparison of the theoretical conformations with experimental results is highlighted.

Melanin-concentrating hormone: a structural and conformational study based on synthesis, biological activity, high-field NMR, and molecular modeling techniques.

A series of Melanin-concentrating hormone (MCH) fragments have been synthesized and their biological activities compared with the parent peptide. The substructural units, 5-14 linear and 5-14 cyclic, have been used as models for MCH-- H-Asp1-Thr-Met-Arg-Cys-Met-Val-Gly-Arg HO-Val17-Glu-Trp-Cys-Pro-Arg-Tyr-Val in 1H-nmr conformational studies. Conformational features predicted by molecular dynamics analyses find support in the nmr experiments.

Accessible conformations of melanin-concentrating hormone: a molecular dynamics approach.

Molecular dynamics simulations have been used to search for the accessible conformations of the melanin-concentrating hormone (MCH). The studies have been performed on native MCH and two of its peptide fragments, a cyclic MCH(5-14) fragment and a linear MCH(5-14) fragment. An analysis of the molecular dynamics trajectories of the three peptides indicates that two regions of the peptide have characteristic conformational properties that may be important for the biological activity. One is a region around Gly8, which is conformationally mobile, and the other is around Pro13, which shows unusual rigidity. The molecular dynamics simulation results are discussed in terms of backbone structural features like beta turns, side-chain interactions, and orientations of the disulfide bridge. The results of this analysis are used to suggest new analogues that will modify the conformational features of the peptide and further define the conformational requirements for activity. Finally, the results are related to nmr studies of the peptide and reveal agreements between the experimental nuclear Overhauser effect constraints and some of the accessible conformations obtained from the simulation.

Modelling of binding sites of the nicotinic acetylcholine receptor and their relation to models of the whole receptor.

Models for the acetylcholine (ACh)-binding site of the nicotinic acetylcholine receptor (nAChR) are proposed. These models have been developed by using the concept of the ligand-gated ion-channel (LGIC) superfamily of receptors that have evolved from a common ancestor. An initial component of the binding site was identified as a highly conserved 15-residue stretch of primary structure in the N-terminal extracellular region of all known LGIC subunits, based on aligned sequence data of LGICs. This subregion, termed the Cys-loop, was modelled as an amphiphilic beta-hairpin and we propose that it forms a major determinant of the binding cleft for agonists. This initial, partial binding-site model has been extended to include residues biochemically identified as spatially adjacent to the binding cleft. A recently developed technique for rapidly scanning the known protein structural database for 'non-homologous similarity' using just sequence information identified the known structure of the enzyme pyrophosphatase (PPase) as a candidate scaffold for the N-terminal domain of the nAChR. This similarity was investigated further using sequence alignments. A framework model of the full N-terminal domain in which the position of the Cys-loop and other binding-site determinants, as well as the main immunogenic region (MIR), have been mapped on to the PPase structure.

Modeling of agonist binding to the ligand-gated ion channel superfamily of receptors.

A generalized model is presented of agonist binding to ligand-gated ion channels (LGICs). Broad similarity in the structure of agonists suggests that the binding sites of LGICs may have evolved from a protobinding site. Aligned sequence data identified as a candidate for such a site a highly conserved 15 residue stretch of primary structure in the N-terminal extracellular region of all known LGIC subunits. We modeled this subregion, termed the cys-loop, as a rigid, amphiphilic beta-hairpin and propose that it may form a major determinant of a conserved structural binding cleft. In the model of the binding complex (1) an invariant aspartate residue at position 11 of the cys-loop is the anionic site interacting with the positively charged amine group of agonists, (2) a local dipole within the pi-electron system of agonists is favorably oriented in the electrostatic field of the invariant aspartate, (3) the epsilon ring-proton of a conserved aromatic residue at the turn of the cys-loop interacts orthogonally with the agonist pi-electron density at its electronegative center, and (4) selective recognition is partly a result of the type of amino acid residue at position 6 of the cys-loop. Additionally, formation of a hydrogen bond between the electronegative atom of the pi-electron system of agonist and a complementary group in the receptor may be important in the high-affinity binding of agonists.

Filtering molecular dynamics trajectories to reveal low-frequency collective motions: phospholipase A2.

A novel method for analysing molecular dynamics trajectories has been developed, which filters out high frequencies using digital signal processing techniques and facilitates focusing on the low-frequency collective motions of proteins. These motions involve low energy slow motions, which lead to important biological phenomena such as domain closure and allosteric effects in enzymes. The filtering method treats each of the atomic trajectories obtained from the molecular dynamics simulation as a "signal". The trajectories of each of the atoms in the system (or any subset of interest) are Fourier transformed to the frequency domain, a filtering function is applied and then an inverse transformation back to the time domain yields the filtered trajectory. The filtering method has been used to study the dynamics of the enzyme phospholipase A2. In the filtered trajectory, all the high frequency bond and valence angle vibrations were eliminated, leaving only low-frequency motion, mainly fluctuations in torsions and conformational transitions. Analysis of this trajectory revealed interesting motions of the protein, including concerted movements of helices, and changes in shape of the active site cavity. Unlike normal mode analysis, which has been used to study the motion of proteins, this method does not require converged minimizations or diagonalization of a matrix of second derivatives. In addition, anharmonicity, multiple minima and conformational transitions are treated explicitly. Thus, the filtering method avoids most of the approximations implicit in other investigations of the dynamic behaviour of large systems.

A novel beta-turn location in an LHRH antagonist: a combined conformational search and molecular dynamics study.

A 50 pico-second molecular dynamics simulation on a cyclic LHRH antagonist analogue Ac-D-Phe1-D-Phe2-D-Trp3-Ser4-Glu5-D-Arg6-Leu7-Lys8+ ++-Pro9-D-Ala10-NH2 (where the cyclisation is via an amide linkage between the Glu5 and Lys8 side chains), reveals some hitherto unseen conformational features. The LHRH analogue is found to adopt a near beta-sheet type of conformation with the reversal in the chain being brought about by a D-Trp3-Ser4-Glu5-D-Arg6 beta turn. The N- and C-terminal ends of the peptide come close together and interact through a network of hydrogen bonds. Additional hydrogen bonds expected of a sheet type of conformation stabilise the lowest energy minima. A conformational search of all possible cyclic structures of a model system c(Glu-D-Ala-Ala-Lys) which was used to determine the starting structure for the simulation studies of the cyclic LHRH antagonist analogue is also highlighted. The influence of the cyclic part on the conformation of this LHRH analogue is discussed.