Researcher profile

Maksim Kouza

Maksim Kouza contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
9works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2015arXiv

Modeling of protein-peptide interactions using the CABS-dock web server for binding site search and flexible docking

Protein-peptide interactions play essential functional roles in living organisms and their structural characterization is a hot subject of current experimental and theoretical research. Computational modeling of the structure of protein-peptide interactions is usually divided into two stages: prediction of the binding site at a protein receptor surface, and then docking (and modeling) the peptide structure into the known binding site. This paper presents a comprehensive CABS-dock method for the simultaneous search of binding sites and flexible protein-peptide docking, available as a users friendly web server. We present example CABS-dock results obtained in the default CABS-dock mode and using its advanced options that enable the user to increase the range of flexibility for chosen receptor fragments or to exclude user-selected binding modes from docking search. Furthermore, we demonstrate a strategy to improve CABS-dock performance by assessing the quality of models with classical molecular dynamics. Finally, we discuss the promising extensions and applications of the CABS-dock method and provide a tutorial appendix for the convenient analysis and visualization of CABS-dock results. The CABS-dock web server is freely available at http://biocomp.chem.uw.edu.pl/CABSdock/

preprint2013arXiv

A structure-based model fails to probe the mechanical unfolding pathways of the titin I27 domain

We discuss the use of a structure based C$α$-Go model and Langevin dynamics to study in detail the mechanical properties and unfolding pathway of the titin I27 domain. We show that a simple Go-model does detect correctly the origin of the mechanical stability of this domain. The unfolding free energy landscape parameters $x_u$ and $ΔG^{\ddagger}$, extracted from dependencies of unfolding forces on pulling speeds, are found to agree reasonably well with experiments. We predict that above $v=10^4$ nm/s the additional force-induced intermediate state is populated at an end-to-end extension of about $75 \mathring{A}$. The force-induced switch in the unfolding pathway occurs at the critical pulling speed $v_{crit} \approx 10^6-10^7$ nm/s. We argue that this critical pulling speed is an upper limit of the interval where Bell's theory works. However, our results suggest that the Go-model fails to reproduce the experimentally observed mechanical unfolding pathway properly, yielding an incomplete picture of the free energy landscape. Surprisingly, the experimentally observed intermediate state with the A strand detached is not populated in Go-model simulations over a wide range of pulling speeds. The discrepancy between simulation and experiment is clearly seen from the early stage of the unfolding process which shows the limitation of the Go model in reproducing unfolding pathways and deciphering the complete picture of the free energy landscape.

preprint2013arXiv

Combining Coarse-Grained Protein Models with Replica-Exchange All-Atom Molecular Dynamics

We describe a combination of all-atom simulations with CABS, a well-established coarse-grained protein modeling tool, into a single multiscale protocol. The simulation method has been tested on the C-terminal beta hairpin of protein G, a model system of protein folding. After reconstructing atomistic details, conformations derived from the CABS simulation were subjected to replica-exchange molecular dynamics simulations with OPLS-AA and AMBER99sb force fields in explicit solvent. Such a combination accelerates system convergence several times in comparison with all-atom simulations starting from the extended chain conformation, demonstrated by the analysis of melting curves, the number of native-like conformations as a function of time and secondary structure propagation. The results strongly suggest that the proposed multiscale method could be an efficient and accurate tool for high-resolution studies of protein folding dynamics in larger systems.

preprint2013arXiv

Numerical Simulation of Folding and Unfolding of Proteins

The thesis examines in detail the folding and unfolding processes of a number of proteins including hbSBD, DDLNF4, single and multi Ubiquitin. Using simplified coarse-grained off-lattice Go model and CD experiments we have shown the two-state behavior of hbSBD protein. It was shown that refolding pathways of single Ubiquitin depend on what end is anchored to the surface. Namely, the fixation of the N-terminal changes refolding pathways but anchoring the C-terminal leaves them unchanged. Interestingly, the end fixation has no effect on multi-domain Ubiquitin. Using the Go modeling and all-atom models with explicit water, we have studied the mechanical unfolding mechanism of DDFLN4 in detail. We predict that, contrary to the AFM experiments, an additional unfolding peak would occur at the end-to-end $ΔR \approx 1.5 $nm in the force-extension curve. Our study reveals the important role of non-native interactions which are responsible for a peak located at $ΔR \approx 22 $nm. This peak can not be encountered by the Go models in which the non-native interactions are neglected. Our finding may stimulate further experimental and theoretical studies on this protein. The cooperativity of denaturation transition was investigated using lattice and off-latice models. Our studies reveal that the sharpness of the transition enhances as the number of amino acids grows. The corresponding scaling behavior is governed by an universal critical exponent. A new force replica exchange (RE) method was developed for efficient configuration sampling of biomolecules pulled by an external mechanical force. Contrary to the standard temperature RE, the exchange is carried out between different forces (replicas). Our method was successfully applied to study thermodynamics of a three-domain Ubiquitin.

preprint2012arXiv

Universal geometrical factor of protein conformations as a consequence of energy minimization

The biological activity and functional specificity of proteins depend on their native three-dimensional structures determined by inter- and intra-molecular interactions. In this paper, we investigate the geometrical factor of protein conformation as a consequence of energy minimization in protein folding. Folding simulations of 10 polypeptides with chain length ranging from 183 to 548 residues manifest that the dimensionless ratio (V/(A<r>)) of the van der Waals volume V to the surface area A and average atomic radius <r> of the folded structures, calculated with atomic radii setting used in SMMP [Eisenmenger F., et. al., Comput. Phys. Commun., 138 (2001) 192], approach 0.49 quickly during the course of energy minimization. A large scale analysis of protein structures show that the ratio for real and well-designed proteins is universal and equal to 0.491\pm0.005. The fractional composition of hydrophobic and hydrophilic residues does not affect the ratio substantially. The ratio also holds for intrinsically disordered proteins, while it ceases to be universal for polypeptides with bad folding properties.

preprint2009arXiv

Dependence of protein mechanical unfolding pathways on pulling speeds

Mechanical unfolding of the fourth domain of Distyostelium discoideum filamin (DDFLN4) was studied in detail using the C$_α$-Go model. We show that unfolding pathways of this protein depend on the pulling speed. The agreement between theoretical and experimental results on the sequencing of unfolding events is achieved at low loading rates. The unfolding free energy landscape is also constructed using dependencies of unfolding forces on pulling speeds.

preprint2009arXiv

Protein mechanical unfolding: importance of non-native interactions

Mechanical unfolding of the fourth domain of Distyostelium discoideum filamin (DDFLN4) was studied by all-atom molecular dynamics simulations, using the GROMOS96 force field 43a1 and the simple point charge explicit water solvent. Our study reveals an important role of non-native interactions in the unfolding process. Namely, the existence of a peak centered at the end-to-end extension 22 nm in the force-extension curve, is associated with breaking of non-native hydrogen bonds. Such a peak has been observed in experiments but not in Go models, where non-native interactions are neglected. We predict that an additional peak occurs at 2 nm using not only GROMOS96 force field 43a1 but also Amber 94 and OPLS force fields. This result would stimulate further experimental studies on elastic properties of DDFLN4.

preprint2006arXiv

Effect of finite size on cooperativity and rates of protein folding

We analyze the dependence of cooperativity of the thermal denaturation transition and folding rates of globular proteins on the number of amino acid residues, $N$, using lattice models with side chains,off-lattice Go models and the available experimental data. A dimensionless measure of cooperativity, $Ω_c$ ($0 < Ω_c < \infty$), scales as $Ω_c \sim N^ζ$. The results of simulations and the analysis of experimental data further confirm the earlier prediction that $ζ$ is universal with $ζ= 1 +γ$, where exponent $γ$ characterizes the susceptibility of a self-avoiding walk. This finding suggests that the structural characteristics in the denaturated state are manifested in the folding cooperativity at the transition temperature. The folding rates $k_F$ for the Go models and a dataset of 69 proteins can be fit using $k_F = k_F^0 \exp(-cN^β)$. Both $β= 1/2$ and 2/3 provide a good fit of the data. We find that $k_F = k_F^0 \exp(-cN^{1/2})$, with the average (over the dataset of proteins) $k_F^0 \approx (0.2μs)^{-1}$ and $c \approx 1.1$, can be used to estimate folding rates to within an order of magnitude in most cases. The minimal models give identical $N$ dependence with $c \approx 1$. The prefactor for off-lattice Go models is nearly four orders of magnitude larger than the experimental value.

preprint2006arXiv

Folding of the Protein Domain hbSBD

The folding of the alpha-helice domain hbSBD of the mammalian mitochondrial branched-chain alpha-ketoacid dehydrogenase (BCKD) complex is studied by the circular dichroism technique in absence of urea. Thermal denaturation is used to evaluate various thermodynamic parameters defining the equilibrium unfolding, which is well described by the two-state model with the folding temperature T_f = 317.8 K and the enthalpy change Delta H_g = 19.67 kcal/mol. The folding is also studied numerically using the off-lattice coarse-grained Go model and the Langevin dynamics. The obtained results, including the population of the native basin, the free energy landscape as a function of the number of native contacts and the folding kinetics, also suggest that the hbSBD domain is a two-state folder. These results are consistent with the biological function of hbSBD in BCKD.