Biomedical subjects

Cezary Czaplewski

Publications and source records attributed to Cezary Czaplewski.

14 recordsLinked to original sources

PubMed2007

Modification and optimization of the united-residue (UNRES) potential energy function for canonical simulations. I. Temperature dependence of the effective energy function and tests of the optimization method with single training proteins.

We report the modification and parametrization of the united-residue (UNRES) force field for energy-based protein structure prediction and protein folding simulations. We tested the approach on three training proteins separately: 1E0L (beta), 1GAB (alpha), and 1E0G (alpha + beta). Heretofore, the UNRES force field had been designed and parametrized to locate native-like structures of proteins as global minima of their effective potential energy surfaces, which largely neglected the conformational entropy because decoys composed of only lowest-energy conformations were used to optimize the force field. Recently, we developed a mesoscopic dynamics procedure for UNRES and applied it with success to simulate protein folding pathways. However, the force field turned out to be largely biased toward -helical structures in canonical simulations because the conformational entropy had been neglected in the parametrization. We applied the hierarchical optimization method, developed in our earlier work, to optimize the force field; in this method, the conformational space of a training protein is divided into levels, each corresponding to a certain degree of native-likeness. The levels are ordered according to increasing native-likeness; level 0 corresponds to structures with no native-like elements, and the highest level corresponds to the fully native-like structures. The aim of optimization is to achieve the order of the free energies of levels, decreasing as their native-likeness increases. The procedure is iterative, and decoys of the training protein(s) generated with the energy function parameters of the preceding iteration are used to optimize the force field in a current iteration. We applied the multiplexing replica-exchange molecular dynamics (MREMD) method, recently implemented in UNRES, to generate decoys; with this modification, conformational entropy is taken into account. Moreover, we optimized the free-energy gaps between levels at temperatures corresponding to a predominance of folded or unfolded structures, as well as to structures at the putative folding-transition temperature, changing the sign of the gaps at the transition temperature. This enabled us to obtain force fields characterized by a single peak in the heat capacity at the transition temperature. Furthermore, we introduced temperature dependence to the UNRES force field; this is consistent with the fact that it is a free-energy and not a potential energy function. beta

Biophysics↗

PubMed2006

Molecular dynamics study of amyloid formation of two Abl-SH3 domain peptides.

Molecular dynamics (MD) simulations were carried out for two-strand and ten-strand beta-sheets constructed from two peptides corresponding to the diverging turn of two homologous Abl-SH3 domains, DLSFMKGE (MK; from Drosophila) and DLSFKKGE (KK; from man), in explicit water at the temperatures of 30, 170/190 and 300 K. It was found that the 2 x MK beta-sheet is more stable than the 2 x KK beta-sheet, and that the 10 x MK beta-sheet is more stable than the 10 x KK beta-sheet; this suggests that the MK systems are fibril-creating and the KK systems are not. These results might explain why most SH3 domains possess two conserved basic residues at the diverging turn, which may act as gatekeepers in order to avoid aggregation.

Amino Acid Sequence↗

PubMed2005

Protein structure prediction with the UNRES force-field using Replica-Exchange Monte Carlo-with-Minimization; Comparison with MCM, CSA, and CFMC.

Two current methods of global optimization are coupled to produce the Replica-Exchange method together with Monte Carlo-with-Minimization (REMCM). Its performance is compared with each separate component and with other global optimization techniques. REMCM was applied to search the conformational space of coarse grain protein systems described by the UNRES force field. The method consists of several noninteracting copies of Monte Carlo simulation, and minimization was used after every perturbation to enhance the sampling of low-energy conformations. REMCM was applied to five proteins of different topology, and the results were compared to those from other optimization methods, namely Monte Carlo-with-Minimization (MCM), Conformational Space Annealing (CSA), and Conformational Family Monte Carlo (CFMC). REMCM located global minima for four proteins faster and more consistently than either MCM or CFMC, and it converged faster than CSA for three of the five proteins tested. A performance comparison was also carried out between REMCM and the traditional Replica Exchange method (REM) for one protein, with REMCM showing a significant improvement. Moreover, because of its simplicity, REMCM was easy to implement, thereby offering an alternative to other global optimization methods used in protein structure prediction.

Algorithms↗

PubMed2005

Molecular origin of anticooperativity in hydrophobic association.

The structuring of water molecules in the vicinity of nonpolar solutes is responsible for hydrophobic hydration and association thermodynamics in aqueous solutions. Here, we studied the potential of mean force (PMF) for the formation of a dimer and trimers of methane molecules in three specific configurations in explicit water to explain multibody effects in hydrophobic association on a molecular level. We analyzed the packing and orientation of water molecules in the vicinity of the solute to explain the effect of ordering of the water around nonpolar solutes on many-body interactions. Consistent with previous theoretical studies, we observed cooperativity, manifested as a reduction of the height of the desolvation barrier for the trimer in an isosceles triangle geometry, but for linear trimers, we observed only anticooperativity. A simple mechanistic picture of hydrophobic association is drawn. The free energy of hydrophobic association depends primarily on the difference in the number of water molecules in the first solvation shell of a cluster and that in the monomers of a cluster; this can be approximated by the molecular surface area. However, there are unfavorable electrostatic interactions between the water molecules from different parts of the solvation shell of a trimer because of their increased orientation induced by the nonpolar solute. These electrostatic interactions make an anticooperative contribution to the PMF, which is clearly manifested for the linear trimer where the multibody contribution due to changes in the molecular surface area is equal to zero. The information theory model of hydrophobic interactions of Hummer et al. also explains the anticooperativity of hydrophobic association of the linear trimers; however, it predicts anticooperativity with a qualitatively identical distance dependence for nonlinear trimers, which disagrees with the results of simulations.

Hydrophobic and Hydrophilic Interactions↗

PubMed2005

An efficient molecular docking using conformational space annealing.

Molecular docking falls into the general category of global optimization problems because its main purpose is to find the most stable complex consisting of a receptor and its ligand. Conformational space annealing (CSA), a powerful global optimization method, is incorporated with the Tinker molecular modeling package to perform molecular docking simulations of six receptor-ligand complexes (3PTB, 1ULB, 2CPP, 1STP, 3CPA, and 1PPH) from the Protein Data Bank. In parallel, Monte Carlo with the minimization (MCM) method is also incorporated into the Tinker package for comparison. The energy function, consisting of electrostatic interactions, van der Waals interactions, and torsional energy terms, is calculated using the AMBER94 all-atom empirical force field. Rigid docking simulations for all six complexes and flexible docking simulations for three complexes (1STP, 3CPA, and 1PPH) are carried out using the CSA and the MCM methods. The simulation results show that the docking procedures using the CSA method generally find the most stable complexes as well as the native-like complexes more efficiently and accurately than those using the MCM, demonstrating that CSA is a promising search method for molecular docking problems.

Algorithms↗

PubMed2005

Charge loss in gas-phase multiply negatively charged oligonucleotides.

In an attempt to shed light on the mechanism by which gaseous samples of negatively charged oligonucleotides undergo extremely slow (i.e., over 1-1000 s) charge loss, we have carried out molecular dynamics simulations on an oligonucleotide anion, T(5)(3-), containing five thymine, deoxyribose, and phosphate units in which the first, third, and fifth phosphates are negatively charged. The study is aimed at determining the rate at which an electron is detached from such a trianion by way of an internal Coulomb repulsion induced event. In this process, the intrinsic 5.0-5.1 eV electron binding strength of each phosphate site is reduced by the repulsive Coulomb potentials of the other two negative sites. As geometrical fluctuations cause the distances among the three negative phosphate sites to change, this causes the Coulomb repulsion energies at these sites to fluctuate. Once the Coulomb potential at any phosphate site exceeds ca. 5 eV, the electron on that site is able to undergo autodetachment. Although such an electron must tunnel through a barrier to escape, it is shown that the tunneling rate is not the rate-limiting step in electron loss; instead, it is the rate at which geometrical fluctuations cause the Coulomb potentials to exceed 5 eV that determines the rate of electron loss. Because these rates are extremely slow, special techniques had to be introduced to allow results of dynamics simulations on more flexible models of T(5)(3-) to be extrapolated to predict the behavior of the actual T(5)(3-).

Computer Simulation↗

PubMed2004

Towards gelsolin amyloid formation.

Amyloid diseases result from protein misfolding and aggregation into fibrils. Some features of gelsolin amyloidogenic fragments comprised of residues 173-243 (G173-243) and residues 173-202 (G173-202) were investigated by the method of molecular dynamics (MD). The alpha-helical structure of G173-243 present in the whole protein unwinds during the course of MD simulation of the fragment G173-243, suggesting that the G173-243 structure is not stable and could unfold before becoming involved in gelsolin amyloid fibril formation. Twelve fragments of G173-202 were used to build a possible beta-fibril. During the course of the simulation, G173-202 fragments formed hydrogen bonds and tended to turn by an angle of 10 degrees -20 degrees towards each other.

Amyloid↗

PubMed2004

Prediction of the structures of proteins with the UNRES force field, including dynamic formation and breaking of disulfide bonds.

The presence of disulfide bonds is essential for maintaining the structure and function of many proteins. The disulfide bonds are usually formed dynamically during folding. This process is not accounted for in present algorithms for protein-structure prediction, which either deduce the possible positions of disulfide bonds only after the structure is formed or assume fixed disulfide bonds during the course of simulated folding. In this work, the conformational space annealing (CSA) method and the UNRES united-residue force field were extended to treat dynamic formation of disulfide bonds. A harmonic potential is imposed on the distance between disulfide-bonded cysteine side-chain centroids to describe the energetics of bond distortion and an energy gain of 5.5 kcal/mol is added for disulfide-bond formation. Formation, breaking and rearrangement of disulfide bonds are included in the CSA search by introducing appropriate operations; the search can also be carried out with a fixed disulfide-bond arrangement. The algorithm was applied to four proteins: 1EI0 (alpha), 1NKL (alpha), 1L1I (beta-helix) and 1ED0 (alpha + beta). For 1EI0, a low-energy structure with correct fold was obtained both in the runs without and with disulfide bonds; however, it was obtained as the lowest in energy only with the native disulfide-bond arrangement. For the other proteins studied, structures with the correct fold were obtained as the lowest (1NKL and 1L1I) or low-energy structures (1ED0) only in runs with disulfide bonds, although the final disulfide-bond arrangement was non-native. The results demonstrate that, by including the possibility of formation of disulfide bonds, the predictive power of the UNRES force field is enhanced, even though the disulfide-bond potential introduced here rarely produces disulfide bonds in native positions. To the best of our knowledge, this is the first algorithm for energy-based prediction of the structure of disulfide-bonded proteins without any assumption as to the positions of native disulfides or human intervention. Directions for improving the potentials and the search method are suggested.

Algorithms↗

PubMed2004

Use of NMR and fluorescence spectroscopy as well as theoretical conformational analysis in conformation-activity studies of cyclic enkephalin analogues.

In this review the conformational studies of natural enkephalins (H-Tyr-Gly-Gly-Phe-Met-OH; the [Met(5)]enkephalin and H-Tyr-Gly-Gly-Phe-Leu-OH; the [Leu(5)]enkephalin), their acyclic and cyclic analogues, including those carried out in our laboratory, performed by experimental and theoretical methods and their combination, are described. Emphasis is given on the role of conformational constraints introduced by cyclization on activity at the micro and delta opioid receptors. Comparison of the conformations of cyclic enkephalin analogues with high delta-receptor activity with those of potent rigid non-peptide delta-receptor agonists indicates that the proximity of the aromatic side chains in positions 1 and 4 as well as the N-terminal amino group is desirable for the activity at the delta opioid receptors; early conformational studies also suggested that spatial separation of the aromatic side chains and rigidity of the cyclic backbone is desirable for micro-receptor activity. The results of our recent conformational studies performed with the use of fluorescence and NMR spectroscopy as well as theoretical calculations indicate, however, that these structural features are not necessary for activity at the micro opioid receptors. Methods applied to the determination of the conformation of flexible peptides, such as Nuclear Magnetic Resonance (NMR), fluorescence spectroscopy, and theoretical conformational analysis are also discussed briefly.

Enkephalins↗

PubMed2003

Molecular dynamics study of a gelsolin-derived peptide binding to a lipid bilayer containing phosphatidylinositol 4,5-bisphosphate.

Gelsolin is an actin-severing protein whose action is initiated by Ca(2+) and inhibited by binding to phosphorylated inositol lipid or phosphoinositides. The regions of gelsolin responsible for phosphoinositide binding are comprised of residues 150-169 (G150-169) and 135-142 (G135-142). The corresponding peptides possess similar binding potency as native gelsolin. Their common feature is the presence of arginine and lysine residues that can bind to negatively charged phosphate groups of phosphoinositides. In this work the binding of the G150-169 peptide to a phosphatidylinositol 4,5-bisphosphate (PIP2) cluster in a lipid membrane model was investigated by molecular dynamics calculations (MD) with the AMBER 4.1 force field, taking into account explicit solvent molecules. Initially the structure of G150-169 was simulated by using the electrostatically driven Monte Carlo (EDMC) and MD methods, and the resulting structure agreed within 3.7 A backbone-atom root mean square deviation with the corresponding experimentally derived structure (PDB code: 1SOL). Using this model for the peptide, a subsequent MD simulation of G150-169 in a periodic box containing a model of dimyristoyl-phosphatidylcholine (DMPC) lipids with a cluster of four PIP2 molecules was carried out. During the simulation G150-169 interacted strongly with PIP2 molecules, initially by formation of salt bridges between its N-terminal basic groups and the phosphate groups of PIP2, followed by formation of hydrophobic bonds between the hydrophobic side chains of the peptide and the fatty acid tail of the lipid. As a result of the formation of hydrophobic bonds, the PIP2 molecules were pulled out from the lipid bilayer. This mode of binding differs from those of other PIP2-binding protein motifs such as PH domains that interact solely with the hydrophilic head group of PIP2. These results suggest that dissociation of gelsolin from actin by PIP2 lipids may involve entering of the PIP2 molecules to the gelsolin-actin interface, thereby weakening the interactions between these proteins.

Binding Sites↗

PubMed2003

Molecular simulation study of cooperativity in hydrophobic association: clusters of four hydrophobic particles.

The multibody contribution to the potential of mean force (PMF) of hydrophobic association of four methane molecules in water was investigated by means of umbrella-sampling molecular dynamics. Two systems were considered: (i). a trigonal pyramid with three methane molecules at contact distance forming a fixed base, the fourth molecule being placed on the top with variable distance from the base; and (ii). a regular uniformly expanding tetrahedron. Methane-methane distances as far as 12.5 A, i.e. beyond the second solvent-separated minimum of the PMF, were considered to address the baseline problem. In contrast to the small effect in the three-body case studied previously (Protein Sci 9 (2000) 1235), the multibody contribution was found to amount to approximately 0.2 kcal/mol per methane-methane pair, or approximately 25% of the depth of the contact minimum in the PMF. The main effect of the multibody contribution to the PMF is a reduction of the height of the barrier between the contact and solvent separated minima and a narrowing of the region of its maximum, while the region of the contact minimum is affected only weakly. The reduction of the barrier is due to four-body contributions. The cooperative contributions to the PMF agree very well with those computed from the molecular surface of the systems under consideration, which further supports earlier observations that the molecular surface can be used with good accuracy to describe the energetics of hydrophobic association.

Hydrophobic and Hydrophilic Interactions↗

PubMed2002

Conformational studies of cyclic enkephalin analogues with L- or D-proline in position 3.

The conformation of a series of cyclic enkephalin analogues of a general formula X(1)-cyclo[Y(2)-Z(3)-Nal(4)-Leu(5)] (Nal: beta-(2-naphthyl)alanine), where X = Tyr, Phe, or Phe(NO(2)), Y = D-Dab or L-Dab (Dab: 2,4-diaminobutyric acid), and Z = D-Pro or L-Pro, was studied by means of NMR spectroscopy and theoretical conformational analysis with the Empirical Conformational Energy Program for Peptides and Proteins force field plus solvation. The NMR measurements were performed in dimethyl sulfoxide solution. The nuclear Overhauser effect intensities and coupling constants were used to compute the statistical weights of the conformations of the ensemble generated in global conformational searches. The purpose of this study was to determine whether introducing the D- or L-proline residue in position 3 can produce peptides with both rigid backbone and significant separation of the pharmacophore groups in position 1 and 4 (as required for high affinity for the mu-type opioid receptors). It was found that the analogues with D-Dab in position 2 and D-Pro in position 3 possess a stable type II' beta-turn at positions 3 and 4, which rigidifies the cyclic backbone; this finding was confirmed by independent measurements of the temperature coefficients of the amide protons, which indicated very significant screening of the Leu(5) amide proton from the solvent. However, these analogues were found to possess a short interchromophore distance. The analogues containing both Dab and Pro in the L-configuration are characterized by a larger interchromophore distance; however, they do not possess a stable beta-turn and have therefore a higher conformational flexibility. The modifications proposed in this work are therefore not likely to lead to enkephalin analogues with a high affinity for the mu-receptors.

Amino Acid Sequence↗

PubMed2002

A method for optimizing potential-energy functions by a hierarchical design of the potential-energy landscape: application to the UNRES force field.

A method for optimizing potential-energy functions of proteins is proposed. The method assumes a hierarchical structure of the energy landscape, which means that the energy decreases as the number of native-like elements in a structure increases, being lowest for structures from the native family and highest for structures with no native-like element. A level of the hierarchy is defined as a family of structures with the same number of native-like elements (or degree of native likeness). Optimization of a potential-energy function is aimed at achieving such a hierarchical structure of the energy landscape by forcing appropriate free-energy gaps between hierarchy levels to place their energies in ascending order. This procedure is different from methods developed thus far, in which the energy gap and/or the Z score between the native structure and all non-native structures are maximized, regardless of the degree of native likeness of the non-native structures. The advantage of this approach lies in reducing the number of structures with decreasing energy, which should ensure the searchability of the potential. The method was tested on two proteins, PDB ID codes and, with an off-lattice united-residue force field. For, the search of the conformational space with the use of the conformational space annealing method and the newly optimized potential-energy function found the native structure very quickly, as opposed to the potential-energy functions obtained by former optimization methods. After even incomplete optimization, the force field obtained by using located the native-like structures of two peptides, and betanova (a designed three-stranded beta-sheet peptide), as the lowest-energy conformations, whereas for the 46-residue N-terminal fragment of staphylococcal protein A, the native-like conformation was the second-lowest-energy conformation and had an energy 2 kcal/mol above that of the lowest-energy structure.

Algorithms↗

PubMed2002

Evolution of physics-based methodology for exploring the conformational energy landscape of proteins.

The evolution of our physics-based computational methods for determining protein conformation without the introduction of secondary-structure predictions, homology modeling, threading, or fragment coupling is described. Initial use of a hard-sphere potential captured much of the structural properties of polypeptide chains, and subsequent more refined force fields, together with efficient methods of global optimization provide indications that progress is being made toward an understanding of the interresidue interactions that underlie protein folding.

Computational Biology↗