System theoretical investigation of human epidermal growth factor receptor-mediated signalling.

The partitioning of biological networks into coupled-functional modules is being increasingly applied for developing predictive models of biological systems. This approach has the advantage that predicting a system-level response does not require a mechanistic description of the internal dynamics of each module. Identification of the input-output characteristics of the network modules and the connectivity between the modules provide the necessary quantitative representation of system dynamics. However, the determination of the input-output relationships of the modules is not trivial; it requires the controlled perturbation of module inputs and systematic analysis of experimental data. In this report, the authors apply a system theoretical analysis approach to derive the time-dependent input-output relationships of the functional module for the human epidermal growth factor receptor (HER) mediated Erk and Akt signalling pathways. Using a library of cell lines expressing endogenous levels of epidermal growth factor receptor (EGFR) and varying levels of HER2, the authors show that a transfer function-based representation can be successfully applied to quantitatively characterise information transfer in this system.

A model of cytokine shedding induced by low doses of gamma radiation.

A model for sustained shedding of epidermal growth factor (EGF) in response to low doses of gamma radiation was developed based on a time delay in the feedback from mitogen-activated protein kinase (MAPK) activation to metalloprotease activity in an autocrine signaling process. We determined the kinetic parameters of our model using the data available for MAPK activation by gamma irradiation in the 1-2-Gy dose range and then showed that predictions of the model were consistent with experimental results for the kinetics of EGF shedding into the growth medium after exposure of human mammary epithelial cells to 1-5 cGy of gamma radiation in the presence of antibodies that block ligand binding to EGF receptors. The model allowed us to estimate the rate of radiation-induced cytokine release per cell from measurements of EGF concentration in the growth medium and to assess the effectiveness of EGF shedding and subsequent diffusion through the medium as a mechanism for signal transmission between hit cells and bystanders.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP.

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25-30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras-p120GAP-GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH(2) group of the Gln-61 and the gamma-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for PO bond cleavage.

Calculating the local solvent chemical potential in crystal hydrates.

Determining solvation patterns in biological systems is crucial in investigating the functional role water may play in structural stabilization and molecular recognition. Determining whether a particular position would be occupied by a solvent molecule requires the local thermodynamics to be known. In this work we introduce a simple and inexpensive approach based on grand canonical molecular simulations to determine the occupancy factors of the cavities. The method is applied to the test case of the sodium salt of hyaluronic acid. The results agree very well with experimental results and demonstrate the success of the method.

Conformational properties of amphotericin B amide derivatives--impact on selective toxicity.

Even though it is highly toxic, Amphotericin B (AmB), an amphipathic polyene macrolide antibiotic, is used in the treatment of severe systemic fungal infections as a life-saving drug. To examine the influence of conformational factors on selective toxicity of these compounds, we have investigated the conformational properties of five AmB amide derivatives. It was found that the extended conformation with torsional angles (phi,psi)=(290 degrees,180 degrees) is a common minimum of the potential energy surfaces (PES) of unsubstituted AmB and its amide derivatives. The extended conformation of the studied compounds allows for the formation of an intermolecular hydrogen bond network between adjacent antibiotic molecules in the open channel configuration. Therefore, the extended conformation is expected to be the dominant conformer in an open AmB (or its amide derivatives) membrane channel. The derivative compounds for calculations were chosen according to their selective toxicity compared to AmB and they had a wide range of selective toxicity. Except for two AmB derivatives, the PES maps of the derivatives reveal that the molecules can coexist in more than one conformer. Taking into account the cumulative conclusions drawn from the earlier MD simulation studies of AmB membrane channel, the results of the potential energy surface maps, and the physical considerations of the molecular structures, we hypothesize a new model of structure-selective toxicity of AmB derivatives. In this proposed model the presence of the extended conformation as the only well defined global conformer for AmB derivatives is taken as the indicator of their higher selective toxicity. This model successfully explains our results. To further test our model, we also investigated an AmB derivative whose selective toxicity has not been experimentally measured before. Our prediction for the selective toxicity of this compound can be tested in experiments to validate or invalidate the proposed model.

Solvation studies of DMP323 and A76928 bound to HIV protease: analysis of water sites using grand canonical Monte Carlo simulations.

We examine the water solvation of the complex of the inhibitors DMP323 and A76928 bound to HIV-1 protease through grand canonical Monte Carlo simulations, and demonstrate the ability of this method to reproduce crystal waters and effectively predict water positions not seen in the DMP323 or A76928 structures. The simulation method is useful for identifying structurally important waters that may not be resolved in the crystal structures. It can also be used to identify water positions around a putative drug candidate docked into a binding pocket. Knowledge of these water positions may be useful in designing drugs to utilize them as bridging groups or displace them in the binding pocket. In addition, the method should be useful in finding water sites in homology models of enzymes for which crystal structures are unavailable.

Molecular properties of amphotericin B membrane channel: a molecular dynamics simulation.

Amphotericin B is a powerful but toxic antifungal antibiotic that is used to treat systemic infections. It forms ionic membrane channels in fungal cells. These antibiotic/sterol channels are responsible for the leakage of ions, which causes cell destruction. The detailed molecular properties and structure of amphotericin B channels are still unknown. In the current study, two molecular dynamic simulations were performed of a particular model of amphotericin B/cholesterol channel. The water and phospholipid environment were included in our simulations, and the results obtained were compared with available experimental data. It was found that it is mainly the hydrogen bonding interactions that keep the channel stable in its open form. Our study also revealed the important role of the intermolecular interactions among the hydroxyl, amino, and carboxyl groups of the channel-forming molecules; in particular, some hydroxyl groups stand out as new "hot spots" that are potentially useful for chemotherapeutic investigations. Our results also help to clarify why certain antibiotic derivatives, with a blocked amino group, are less active. We present a hypothesis for the role of membrane lipids and cholesterol in the channel.

Enzyme-inhibitor association thermodynamics: explicit and continuum solvent studies.

Studying the thermodynamics of biochemical association reactions at the microscopic level requires efficient sampling of the configurations of the reactants and solvent as a function of the reaction pathways. In most cases, the associating ligand and receptor have complementary interlocking shapes. Upon association, loosely connected or disconnected solvent cavities at and around the binding site are formed. Disconnected solvent regions lead to severe statistical sampling problems when simulations are performed with explicit solvent. It was recently proposed that, when such limitations are encountered, they might be overcome by the use of the grand canonical ensemble. Here we investigate one such case and report the association free energy profile (potential of mean force) between trypsin and benzamidine along a chosen reaction coordinate as calculated using the grand canonical Monte Carlo method. The free energy profile is also calculated for a continuum solvent model using the Poisson equation, and the results are compared to the explicit water simulations. The comparison shows that the continuum solvent approach is surprisingly successful in reproducing the explicit solvent simulation results. The Monte Carlo results are analyzed in detail with respect to solvation structure. In the binding site channel there are waters bridging the carbonyl oxygen groups of Asp189 with the NH2 groups of benzamidine, which are displaced upon inhibitor binding. A similar solvent-bridging configuration has been seen in the crystal structure of trypsin complexed with bovine pancreatic trypsin inhibitor. The predicted locations of other internal waters are in very good agreement with the positions found in the crystal structures, which supports the accuracy of the simulations.

The sensitivity of conformational free energies of the alanine dipeptide to atomic site charges.

Different atomic point charge sets are obtained for the alpha R and C7,eq conformations of the alanine dipeptide by fitting the charges of each conformation to the respective ab initio electrostatic potential surfaces both individually and simultaneously, in both the united atom and the all-atom representations. Using these charge sets, the sensitivity of the relative conformational aqueous free energies to the atomic site charges is investigated. For this particular system, we find that the solute-water contributions to the conformational free energy differences have a rather weak dependence on site charges; the calculated intramolecular contributions, however, show a rather strong dependence on the atomic site charges. It is suggested that the calculated results for the alanine dipeptide using a single, simultaneously fit set of charges for both conformations are in better agreement with experiments than the calculations carried out with charges determined individually for each conformation.

Grand canonical ensemble Monte Carlo simulation of the dCpG/proflavine crystal hydrate.

The grand canonical ensemble Monte Carlo molecular simulation method is used to investigate hydration patterns in the crystal hydrate structure of the dCpG/proflavine intercalated complex. The objective of this study is to show by example that the recently advocated grand canonical ensemble simulation is a computationally efficient method for determining the positions of the hydrating water molecules in protein and nucleic acid structures. A detailed molecular simulation convergence analysis and an analogous comparison of the theoretical results with experiments clearly show that the grand ensemble simulations can be far more advantageous than the comparable canonical ensemble simulations.