No Arabic abstract
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.
Nearly a quarter of genomic sequences and almost half of all receptors that are likely to be targets for drug design are integral membrane proteins. Understanding the detailed mechanisms of the folding of membrane proteins is a largely unsolved, key problem in structural biology. Here, we introduce a general model and use computer simulations to study the equilibrium properties and the folding kinetics of a $C_{alpha}$-based two helix bundle fragment (comprised of 66 amino-acids) of Bacteriorhodopsin. Various intermediates are identified and their free energy are calculated toghether with the free energy barrier between them. In 40% of folding trajectories, the folding rate is considerably increased by the presence of non-obligatory intermediates acting as traps. In all cases, a substantial portion of the helices is rapidly formed. This initial stage is followed by a long period of consolidation of the helices accompanied by their correct packing within the membrane. Our results provide the framework for understanding the variety of folding pathways of helical transmembrane proteins.
External forces acting on a microswimmer can feed back on its self-propulsion mechanism. We discuss this load response for a generic microswimmer that swims by cyclic shape changes. We show that the change in cycle frequency is proportional to the Lighthill efficiency of self-propulsion. As a specific example, we consider Najafis three-sphere swimmer. The force-velocity relation of a microswimmer implies a correction for a formal superposition principle for active and passive motion.
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.
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.