Biomedical subjects

Marek Cieplak

Publications and source records attributed to Marek Cieplak.

At least 19 recordsLinked to original sources

PubMed2006

Using the principle of entropy maximization to infer genetic interaction networks from gene expression patterns.

We describe a method based on the principle of entropy maximization to identify the gene interaction network with the highest probability of giving rise to experimentally observed transcript profiles. In its simplest form, the method yields the pairwise gene interaction network, but it can also be extended to deduce higher-order interactions. Analysis of microarray data from genes in Saccharomyces cerevisiae chemostat cultures exhibiting energy metabolic oscillations identifies a gene interaction network that reflects the intracellular communication pathways that adjust cellular metabolic activity and cell division to the limiting nutrient conditions that trigger metabolic oscillations. The success of the present approach in extracting meaningful genetic connections suggests that the maximum entropy principle is a useful concept for understanding living systems, as it is for other complex, nonequilibrium systems.

Computational Biology↗

PubMed2006

Stretching of proteins in a uniform flow.

Stretching of a protein by a fluid flow is compared to that in a force-clamp apparatus. The comparison is made within a simple topology-based dynamical model of a protein in which the effects of the flow are implemented using Langevin dynamics. We demonstrate that unfolding induced by a uniform flow shows a richer behavior than that in the force clamp. The dynamics of unfolding is found to depend strongly on the selection of the amino acid, usually one of the termini, which is anchored. These features offer potentially wider diagnostic tools to investigate structure of proteins compared to experiments based on the atomic force microscopy.

Models, Biological↗

PubMed2006

Protein folding in a force clamp.

Kinetics of folding of a protein held in a force clamp are compared to an unconstrained folding. The comparison is made within a simple topology-based dynamical model of ubiquitin. We demonstrate that the experimentally observed variations in the end-to-end distance reflect microscopic events during folding. However, the folding scenarios in and out of the force clamp are distinct.

Protein Folding↗

PubMed2006

Pulling single bacteriorhodopsin out of a membrane: Comparison of simulation and experiment.

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.

Bacterial Proteins↗

PubMed2006

Nanoscale fluid flows in the vicinity of patterned surfaces.

Molecular dynamics simulations of dense and rarefied fluids comprising small chain molecules in chemically patterned nanochannels predict a novel switching from Poiseuille to plug flow along the channel. We also demonstrate behavior akin to the lotus effect for a nanodrop on a chemically patterned substrate. Our results show that one can control and exploit the behavior of fluids at the nanoscale using chemical patterning.

Journal Article↗

PubMed2006

Geometry of proteins: hydrogen bonding, sterics, and marginally compact tubes.

The functionality of proteins is governed by their structure in the native state. Protein structures are made up of emergent building blocks of helices and almost planar sheets. A simple coarse-grained geometrical model of a flexible tube barely subject to compaction provides a unified framework for understanding the common character of globular proteins. We argue that a recent critique of the tube idea is not well founded.

Computer Simulation↗

PubMed2005

Mechanical unfolding of ubiquitin molecules.

Mechanical stretching of ubiquitin and of its several repeats are studied through molecular-dynamics simulations. A Go-type model [H. Abe and N. Go, Biopolymers 20, 1013 (1981)] with a realistic contact map and with Lennard-Jones contact interactions is used. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on polyubiquitins stretched by various linkages. The terminal-to-terminal stretching of polyubiquitin results in peak forces similar to those measured for titin-based polyproteins and of a magnitude that matches measurements. Consistent with the experimental measurements, the simulated peak forces depend on the pulling speed logarithmically when thermal fluctuations are explicitly introduced. These results validate the application of topology-based models in the study of the mechanical stretching of proteins.

Amino Acids↗

PubMed2005

Thermal unfolding of proteins.

Thermal unfolding of proteins is compared to folding and mechanical stretching in a simple topology-based dynamical model. We define the unfolding time and demonstrate its low-temperature divergence. Below a characteristic temperature, contacts break at separate time scales and unfolding proceeds approximately in a way reverse to folding. Features in these scenarios agree with experiments and atomic simulations on titin.

Chemistry, Physical↗

PubMed2005

Mechanical properties of the domains of titin in a Go-like model.

Comparison of properties of three domains of titin, I1, I27, and I28, in a simple geometry-based model shows that despite a high structural homology between their native states different domains show similar but distinguishable mechanical properties. Folding properties of the separate domains are predicted to be diversified which reflects sensitivity of the kinetics to the details of native structures. The Go-like model corresponding to the experimentally resolved native structure of the I1 domain is found to provide the biggest thermodynamic and mechanical stability compared to the other domains studied here. We analyze elastic, thermodynamic, and kinetic properties of several structures corresponding to the I28 domain as obtained through homology-based modeling. We discuss the ability of the models of the I28 domain to reproduce experimental results qualitatively. A strengthening of contacts that involve hydrophobic amino acids does not affect theoretical comparisons of the domains. Tandem linkages of up to five identical or different domains unravel in a serial fashion at low temperatures. We study the nature of the intermediate state that arises in the early stages of the serial unraveling and find it to qualitatively agree with the results of Marszalek et al.

Journal Article↗

PubMed2004

Lattice tube model of proteins.

We present a new lattice model for proteins that incorporates a tubelike anisotropy by introducing a preference for mutually parallel alignments in the conformations. The model is demonstrated to capture many aspects of real proteins.

Anisotropy↗

PubMed2004

What can one learn from experiments about the elusive transition state?

We present the results of an exact analysis of a model energy landscape of a protein to clarify the idea of the transition state and the physical meaning of the phi values determined in protein engineering experiments. We benchmark our findings to various theoretical approaches proposed in the literature for the identification and characterization of the transition state.

Models, Molecular↗

PubMed2004

Thermal effects in stretching of Go-like models of titin and secondary structures.

The effect of temperature on mechanical unfolding of proteins is studied using a Go-like model with a realistic contact map and Lennard-Jones contact interactions. The behavior of the I27 domain of titin and its serial repeats is contrasted to that of simple secondary structures. In all cases, thermal fluctuations accelerate the unraveling process, decreasing the unfolding force nearly linearly at low temperatures. However, differences in bonding geometry lead to different sensitivity to temperature and different changes in the unfolding pattern. Due to its special native-state geometry, titin is much more thermally and elastically stable than the secondary structures. At low temperatures, serial repeats of titin show a parallel unfolding of all domains to an intermediate state, followed by serial unfolding of the domains. At high temperatures, all domains unfold simultaneously, and the unfolding distance decreases monotonically with the contact order, that is, the sequence distance between the amino acids that form the native contact.

Computer Simulation↗

PubMed2004

Stretching of homopolymers and contact order.

Mechanical stretching of self-interacting homopolymers is studied through molecular dynamics simulations in which the polymers are pulled with constant speed. At temperatures below the compactification temperature, the force-extension curves show a plateau that corresponds to the situation in which the polymer adopts "ball-string" configurations. The dependence of rupture distances on contact order and the effects of temperature are similar to those found in the case of model proteins. The dependence of behavior on the pulling speed is logarithmic. In the entropic limit, above the compactification temperature, the rupturing of contacts shows a monotonic decrease of extension with the contact order. The attainment of this limit depends on the system size and the pulling speed.

Biopolymers↗

PubMed2004

Cooperativity and contact order in protein folding.

The effects of cooperativity are studied within Go-Lennard-Jones models of proteins by making the contact interactions dependent on the proximity to the native conformation. The kinetic universality classes are found to remain the same as in the absence of cooperativity. For a fixed native geometry, small changes in the effective contact map may affect the folding times in a chance way, and, to an extent that is comparable to the shift in the folding times due to cooperativity. The contact order controls folding scenarios: the average times necessary to bring pairs of amino acids into their near native separations depend on the sequential distances within the pairs. This dependence is largely monotonic, regardless of the cooperativity, and the dominant trend could be described by a single parameter like the average contact order. However, it is the deviations from the trend which are usually found to set the net folding times.

Amino Acids↗

PubMed2004

Stretching of proteins in the entropic limit.

Mechanical stretching of six proteins is studied through molecular dynamics simulations. The model is Go-like, with Lennard-Jones interactions at native contacts. Low-temperature unfolding scenarios are remarkably complex and sensitive to small structural changes. Thermal fluctuations reduce the peak forces and the number of metastable states during unfolding. The unfolding pathways also simplify as temperature rises. In the entropic limit, all proteins show a monotonic decrease of the extension where bonds rupture with their separation along the backbone (contact order).

Computer Simulation↗

PubMed2003

Assembly of protein tertiary structures from secondary structures using optimized potentials.

We present a simulated annealing-based method for the prediction of the tertiary structures of proteins given knowledge of the secondary structure associated with each amino acid in the sequence. The backbone is represented in a detailed fashion whereas the sidechains and pairwise interactions are modeled in a simplified way, following the LINUS model of Srinivasan and Rose. A perceptron-based technique is used to optimize the interaction potentials for a training set of three proteins. For these proteins, the procedure is able to reproduce the tertiary structures to below 3 A in root mean square deviation (rmsd) from the PDB targets. We present the results of tests on twelve other proteins. For half of these, the lowest energy decoy has a rmsd from the native state below 6 A and, in 9 out of 12 cases, we obtain decoys whose rmsd from the native states are also well below 5 A.

Amino Acids↗

PubMed2003

Universality classes in folding times of proteins.

Molecular dynamics simulations in simplified models allow one to study the scaling properties of folding times for many proteins together under a controlled setting. We consider three variants of the Go models with different contact potentials and demonstrate scaling described by power laws and no correlation with the relative contact order parameter. We demonstrate existence of at least three kinetic universality classes that are correlated with the types of structure: the alpha-, alpha-beta-, and beta- proteins have the scaling exponents of approximately 1.7, 2.5, and 3.2, respectively. The three classes merge into one when the contact range is truncated at a reasonable value. We elucidate the role of the potential associated with the chirality of a protein.

Computer Simulation↗

PubMed2002

Molecular dynamics simulations of crystallization of hard spheres.

We have carried out molecular dynamics simulations of the crystallization of hard spheres modeling colloidal systems that are studied in conventional and space-based experiments. We use microscopic probes to investigate the effects of gravitational forces, polydispersity, and of bounding walls on the phase structure. The simulations employed an extensive exclusive particle grid method and the type and degree of crystalline order was studied in two independent ways: by the structure factor, as in experiments, and through local rotational invariants. We present quantitative comparisons of the nucleation rates of monodisperse and polydisperse hard-sphere systems and benchmark them against experimental results. We show how the presence of bounding walls leads to wall-induced nucleation and rapid crystallization, and discuss the role of gravity on the dynamics of crystallization.

Journal Article↗