Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCombined force-torque spectroscopy of proteins by means of multiscale molecular simulation
78
0
0.0
(
0
)
Added by
Thijs van der Heijden
Publication date
2019
fields
Physics
and research's language is
English
Authors
Thijs W. G. van der Heijden
Ask ChatGPT about the research
No Arabic abstract
Assessing the structural properties of large proteins is important to gain an understanding of their function in, e.g., biological systems or biomedical applications. We propose a method to examine the mechanical properties of proteins subject to applied forces by means of multiscale simulation. We consider both stretching and torsional forces, which can be applied independently of each other. We apply torsional forces to a coarse-grained continuum model of the antibody protein immunoglobulin G (IgG) using Fluctuating Finite Element Analysis and identify the area of strongest deformation. This region is essential to the torsional properties of the molecule as a whole, as it represents the softest, most deformable domain. We subject this part of the molecule to torques and stretching forces on an atomistic level, using molecular dynamics simulations, in order to investigate its torsional properties. We calculate the torsional resistance as a function of the rotation of the domain, while subjecting it to various stretching forces. We learn how these obtained torsion profiles evolve with increasing stretching force and show that they exhibit torsion stiffening, which is in qualitative agreement with experimental findings. We argue that combining the torsion profiles for various stretching forces effectively results in a combined force-torque spectroscopy analysis, which may serve as a mechanical signature for the examined molecule.
rate research
Read More
We use coarse grained molecular dynamics simulations to investigate diffusion properties of sheared lipid membranes with embedded transmembrane proteins. In membranes without proteins, we find normal in-plane diffusion of lipids in all flow conditions. Protein embedded membranes behave quite differently: by imposing a simple shear flow and sliding the monolayers of the membrane over each other, the motion of protein clusters becomes strongly superdiffusive in the shear direction. In such a circumstance, subdiffusion regime is predominant perpendicular to the flow. We show that superdiffusion is a result of accelerated chaotic motions of protein--lipid complexes within the membrane voids, which are generated by hydrophobic mismatch or the transport of lipids by proteins.
Cytoskeletal motor proteins are involved in major intracellular transport processes which are vital for maintaining appropriate cellular function. The motor exhibits distinct states of motility: active motion along filaments, and effectively stationary phase in which it detaches from the filaments and performs passive diffusion in the vicinity of the detachment point due to cytoplasmic crowding. The transition rates between motion and pause phases are asymmetric in general, and considerably affected by changes in environmental conditions which influences the efficiency of cargo delivery to specific targets. By considering the motion of molecular motor on a single filament as well as a dynamic filamentous network, we present an analytical model for the dynamics of self-propelled particles which undergo frequent pause phases. The interplay between motor processivity, structural properties of filamentous network, and transition rates between the two states of motility drastically changes the dynamics: multiple transitions between different types of anomalous diffusive dynamics occur and the crossover time to the asymptotic diffusive or ballistic motion varies by several orders of magnitude. We map out the phase diagrams in the space of transition rates, and address the role of initial conditions of motion on the resulting dynamics.
Protein aggregation in cell membrane is vital for the majority of biological functions. Recent experimental results suggest that transmembrane domains of proteins such as $alpha$-helices and $beta$-sheets have different structural rigidities. We use molecular dynamics simulation of a coarse-grained model of protein-embedded lipid membranes to investigate the mechanisms of protein clustering. For a variety of protein concentrations, our simulations under thermal equilibrium conditions reveal that the structural rigidity of transmembrane domains dramatically affects interactions and changes the shape of the cluster. We have observed stable large aggregates even in the absence of hydrophobic mismatch which has been previously proposed as the mechanism of protein aggregation. According to our results, semi-flexible proteins aggregate to form two-dimensional clusters while rigid proteins, by contrast, form one-dimensional string-like structures. By assuming two probable scenarios for the formation of a two-dimensional triangular structure, we calculate the lipid density around protein clusters and find that the difference in lipid distribution around rigid and semiflexible proteins determines the one- or two-dimensional nature of aggregates. It is found that lipids move faster around semiflexible proteins than rigid ones. The aggregation mechanism suggested in this paper can be tested by current state-of-the-art experimental facilities.
We study the space of all compact structures on a two-dimensional square lattice of size $N=6times6$. Each structure is mapped onto a vector in $N$-dimensions according to a hydrophobic model. Previous work has shown that the designabilities of structures are closely related to the distribution of the structure vectors in the $N$-dimensional space, with highly designable structures predominantly found in low density regions. We use principal component analysis to probe and characterize the distribution of structure vectors, and find a non-uniform density with a single peak. Interestingly, the principal axes of this peak are almost aligned with Fourier eigenvectors, and the corresponding Fourier eigenvalues go to zero continuously at the wave-number for alternating patterns ($q=pi$). These observations provide a stepping stone for an analytic description of the distribution of structural points, and open the possibility of estimating designabilities of realistic structures by simply Fourier transforming the hydrophobicities of the corresponding sequences.
The recently introduced method of excess collisions (MEC) is modified to estimate diffusion-controlled reaction times inside systems of arbitrary size. The resulting MEC-E equations contain a set of empirical parameters, which have to be calibrated in numerical simulations inside a test system of moderate size. Once this is done, reaction times of systems of arbitrary dimensions are derived by extrapolation, with an accuracy of 10 to 15 percent. The achieved speed up, when compared to explicit simulations of the reaction process, is increasing proportional to the extrapolated volume of the cell.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
