Protein-Protein Docking with Large-Scale Backbone Flexibility Using Coarse-Grained Monte-Carlo Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Kurcinski M
PROVIDER: S-EPMC8306105 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A simple coarse-grained model of mucus structure and dynamics is proposed and evaluated. The model is based on simple cubic, face-centered lattice representation. Mucins are simulated as lattice chains in which each bead of the model chains represents a mucin domain, equivalent to its Kuhn segment. The remaining lattice sites are considered to be occupied by the solvent. Model mucins consist of three types of domains: polar (glycosylated central segments), hydrophobic, and cysteine-rich, located at the terminal part of the mucin chains. The sequence of these domains mimics the sequence of real mucins. Static and dynamic properties of the system were studied by means of Monte Carlo dynamics. It was shown that the model system undergoes sol-gel transition and that the interactions between hydrophobic domains are responsible for the transition and characteristic properties of the dynamic network in the gel phase. Cysteine-rich domains are essential for frictional properties of the system. Structural and dynamic properties of the model mucus observed in simulations are in qualitative agreement with known experimental facts and provide mechanistic explanation of complex properties of real mucus.
| S-EPMC2998598 | biostudies-literature
Project description:Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this minireview, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions, and highly simplified structure-based elastic network models of protein flexibility. In contrast to classical all-atom molecular dynamics, the modeling strategies discussed here allow the quite accurate modeling of much larger systems and longer-time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.
| S-EPMC6274762 | biostudies-literature
Project description:Simulating the nature of voltage-activated systems is a problem of major current interest, ranging from the action of voltage-gated ion channels to energy storage batteries. However, fully microscopic converging molecular simulations of external voltage effects present a major challenge, and macroscopic models are associated with major uncertainties about the dielectric treatment and the underlying physical basis. Recently we developed a coarse-grained (CG) model that represents explicitly the electrodes, the electrolytes, and the membrane/protein system. The CG model provides a semimacroscopic way of capturing the microscopic physics of voltage-activated systems. Our method was originally validated by reproducing macroscopic and analytical results for key test cases and then used in modeling voltage-activated ion channels and related problems. In this work, we further establish the reliability of the CG voltage model by comparing it to the results of Monte Carlo (MC) simulations with a microscopic electrolyte model. The comparison explores different aspects of membrane, electrolyte, and electrode systems ranging from the Gouy-Chapman model to the determination of the electrolyte charge distribution in the solution between two electrodes (without and with a separating membrane), as well as the evaluation of gating charges. Overall the agreement is very impressive. This provides confidence in the CG model and also shows that the MC model can be used in realistic simulation of voltage activation of membrane proteins with sufficient computer time.
| S-EPMC6513005 | biostudies-literature
Project description:BackgroundAn adaptive coarse-grained (kinetic) Monte Carlo (ACGMC) simulation framework is applied to reaction and diffusion dynamics in inhomogeneous domains. The presented model is relevant to the diffusion and dimerization dynamics of epidermal growth factor receptor (EGFR) in the presence of plasma membrane heterogeneity and specifically receptor clustering. We perform simulations representing EGFR cluster dissipation in heterogeneous plasma membranes consisting of higher density clusters of receptors surrounded by low population areas using the ACGMC method. We further investigate the effect of key parameters on the cluster lifetime.ResultsCoarse-graining of dimerization, rather than of diffusion, may lead to computational error. It is shown that the ACGMC method is an effective technique to minimize error in diffusion-reaction processes and is superior to the microscopic kinetic Monte Carlo simulation in terms of computational cost while retaining accuracy. The low computational cost enables sensitivity analysis calculations. Sensitivity analysis indicates that it may be possible to retain clusters of receptors over the time scale of minutes under suitable conditions and the cluster lifetime may depend on both receptor density and cluster size.ConclusionsThe ACGMC method is an ideal platform to resolve large length and time scales in heterogeneous biological systems well beyond the plasma membrane and the EGFR system studied here. Our results demonstrate that cluster size must be considered in conjunction with receptor density, as they synergistically affect EGFR cluster lifetime. Further, the cluster lifetime being of the order of several seconds suggests that any mechanisms responsible for EGFR aggregation must operate on shorter timescales (at most a fraction of a second), to overcome dissipation and produce stable clusters observed experimentally.
| S-EPMC2868014 | biostudies-literature
Project description:Structures of hitherto unknown protein complexes can be predicted by docking the solved protein monomers. Here, we present a method to refine initial docking estimates of protein complex structures by a Monte Carlo approach including rigid-body moves and side-chain optimization. The energy function used is comprised of van der Waals, Coulomb, and atomic contact energy terms. During the simulation, we gradually shift from a novel smoothed van der Waals potential, which prevents trapping in local energy minima, to the standard Lennard-Jones potential. Following the simulation, the conformations are clustered to obtain the final predictions. Using only the first 100 decoys generated by a fast Fourier transform (FFT)-based rigid-body docking method, our refinement procedure is able to generate near-native structures (interface RMSD <2.5 A) as first model in 14 of 59 cases in a benchmark set. In most cases, clear binding funnels around the native structure can be observed. The results show the potential of Monte Carlo refinement methods and emphasize their applicability for protein-protein docking.
| S-EPMC2222816 | biostudies-literature
Project description:We present a novel multi-conformation Monte Carlo simulation method that enables the modeling of protein-protein interactions and aggregation in crowded protein solutions. This approach is relevant to a molecular-scale description of realistic biological environments, including the cytoplasm and the extracellular matrix, which are characterized by high concentrations of biomolecular solutes (e.g., 300-400 mg/mL for proteins and nucleic acids in the cytoplasm of Escherichia coli). Simulation of such environments necessitates the inclusion of a large number of protein molecules. Therefore, computationally inexpensive methods, such as rigid-body Brownian dynamics (BD) or Monte Carlo simulations, can be particularly useful. However, as we demonstrate herein, the rigid-body representation typically employed in simulations of many-protein systems gives rise to certain artifacts in protein-protein interactions. Our approach allows us to incorporate molecular flexibility in Monte Carlo simulations at low computational cost, thereby eliminating ambiguities arising from structure selection in rigid-body simulations. We benchmark and validate the methodology using simulations of hen egg white lysozyme in solution, a well-studied system for which extensive experimental data, including osmotic second virial coefficients, small-angle scattering structure factors, and multiple structures determined by X-ray and neutron crystallography and solution NMR, as well as rigid-body BD simulation results, are available for comparison.
| S-EPMC7797190 | biostudies-literature
Project description:Loop flexibility is often crucial to protein biological function in solution. We report a new Monte Carlo method for generating conformational ensembles for protein loops and cyclic peptides. The approach incorporates the triaxial loop closure method which addresses the inverse kinematic problem for generating backbone move sets that do not break the loop. Sidechains are sampled together with the backbone in a hierarchical way, making it possible to make large moves that cross energy barriers. As an initial application, we apply the method to the flexible loop in triosephosphate isomerase that caps the active site, and demonstrate that the resulting loop ensembles agree well with key observations from previous structural studies. We also demonstrate, with 3 other test cases, the ability to distinguish relatively flexible and rigid loops within the same protein.
| S-EPMC3129859 | biostudies-literature
Project description:Using the endosomal sorting complex required for transport (ESCRT)-III membrane remodeling complex as an example, we analyze three popular coarse-grained models (the regular MARTINI, polarizable MARTINI (POL-MARTINI), and big multipole water MARTINI (BMW-MARTINI)) for the description of membrane curvature sensing and generation activities of peripheral proteins. Although the three variants of the MARTINI model provide consistent descriptions for the protein-protein interface in a linear filament model of ESCRT-III, they differ considerably in terms of protein-membrane interface and therefore membrane curvature sensing and generation behaviors. In particular, BMW-MARTINI provides the most consistent description of the protein-membrane interface as compared to all-atom simulations, whereas the regular MARTINI is most consistent with atomistic simulations in terms of the qualitative sign of membrane curvature sensing and generation. With POL-MARTINI, the ESCRT-III model interacts weakly with the membrane and therefore does not exhibit any curvature-sensitive activities. Analysis suggests that the incorrect membrane curvature activities predicted by BMW-MARTINI are due to overestimated insertion depth of an amphipathic helix and incorrect sign for the spontaneous curvature of anionic lipids. These results not only point to ways that coarse-grained models can be improved but also explicitly highlight local lipid composition and insertion depth of protein motifs as essential regulatory factors for membrane curvature sensing and generation.
| S-EPMC8390925 | biostudies-literature
Project description:Many commonly used coarse-grained models for proteins are based on simplified interaction sites and consequently may suffer from significant limitations, such as the inability to properly model protein secondary structure without the addition of restraints. Recent work on a benzene fluid (Lettieri S.; Zuckerman D. M.J. Comput. Chem.2012, 33, 268-275) suggested an alternative strategy of tabulating and smoothing fully atomistic orientation-dependent interactions among rigid molecules or fragments. Here we report our initial efforts to apply this approach to the polar and covalent interactions intrinsic to polypeptides. We divide proteins into nearly rigid fragments, construct distance and orientation-dependent tables of the atomistic interaction energies between those fragments, and apply potential energy smoothing techniques to those tables. The amount of smoothing can be adjusted to give coarse-grained models that range from the underlying atomistic force field all the way to a bead-like coarse-grained model. For a moderate amount of smoothing, the method is able to preserve about 70-90% of the α-helical structure while providing a factor of 3-10 improvement in sampling per unit computation time (depending on how sampling is measured). For a greater amount of smoothing, multiple folding-unfolding transitions of the peptide were observed, along with a factor of 10-100 improvement in sampling per unit computation time, although the time spent in the unfolded state was increased compared with less smoothed simulations. For a β hairpin, secondary structure is also preserved, albeit for a narrower range of the smoothing parameter and, consequently, for a more modest improvement in sampling. We have also applied the new method in a "resolution exchange" setting, in which each replica runs a Monte Carlo simulation with a different degree of smoothing. We obtain exchange rates that compare favorably to our previous efforts at resolution exchange (Lyman E.; Zuckerman D. M.J. Chem. Theory Comput.2006, 2, 656-666).
| S-EPMC4230378 | biostudies-literature
Project description:Cells are complex, dynamic systems that actively adapt to various stimuli including mechanical alterations. Central to understanding cellular response to mechanical stimulation is the organization of the cytoskeleton and its actin filament network. In this manuscript, we present a minimalistic network Monte Carlo based approach to model actin filament organization under cyclic stretching. Utilizing a coarse-grained model, a filament network is prescribed within a two-dimensional circular space through nodal connections. When cyclically stretched, the model demonstrates that a perpendicular alignment of the filaments to the direction of stretch emerges in response to nodal repositioning to minimize net nodal forces from filament stress states. In addition, the filaments in the network rearrange and redistribute themselves to reduce the overall stress by decreasing their individual stresses. In parallel, we cyclically stretch NIH 3T3 fibroblasts and find a similar cytoskeletal response. With this work, we test the hypothesis that a first-principles mechanical model of filament assembly in a confined space is by itself capable of yielding the remodeling behavior observed experimentally. Identifying minimal mechanisms sufficient to reproduce mechanical influences on cellular structure has important implications in a diversity of fields, including biology, physics, medicine, computer science, and engineering.
| S-EPMC3501734 | biostudies-literature