Role of protein interactions in defining HIV-1 viral capsid shape and stability: a coarse-grained analysis.
Ontology highlight
ABSTRACT:
SUBMITTER: Krishna V
PROVIDER: S-EPMC2800964 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The maturation of HIV-1 viral particles is essential for viral infectivity. During maturation, many copies of the capsid protein (CA) self-assemble into a capsid shell to enclose the viral RNA. The mechanistic details of the initiation and early stages of capsid assembly remain to be delineated. We present coarse-grained simulations of capsid assembly under various conditions, considering not only capsid lattice self-assembly but also the potential disassembly of capsid upon delivery to the cytoplasm of a target cell. The effects of CA concentration, molecular crowding, and the conformational variability of CA are described, with results indicating that capsid nucleation and growth is a multi-stage process requiring well-defined metastable intermediates. Generation of the mature capsid lattice is sensitive to local conditions, with relatively subtle changes in CA concentration and molecular crowding influencing self-assembly and the ensemble of structural morphologies.
| S-EPMC4869257 | biostudies-literature
Project description:Coarse-grained (CG) models of biomolecules have recently attracted considerable interest because they enable the simulation of complex biological systems on length-scales and timescales that are inaccessible for atomistic molecular dynamics simulation. A CG model is defined by a map that transforms an atomically detailed configuration into a CG configuration. For CG models of relatively small biomolecules or in cases that the CG and all-atom models have similar resolution, the construction of this map is relatively straightforward and can be guided by chemical intuition. However, it is more challenging to construct a CG map when large and complex domains of biomolecules have to be represented by relatively few CG sites. This work introduces a new and systematic methodology called essential dynamics coarse-graining (ED-CG). This approach constructs a CG map of the primary sequence at a chosen resolution for an arbitrarily complex biomolecule. In particular, the resulting ED-CG method variationally determines the CG sites that reflect the essential dynamics characterized by principal component analysis of an atomistic molecular dynamics trajectory. Numerical calculations illustrate this approach for the HIV-1 CA protein dimer and ATP-bound G-actin. Importantly, since the CG sites are constructed from the primary sequence of the biomolecule, the resulting ED-CG model may be better suited to appropriately explore protein conformational space than those from other CG methods at the same degree of resolution.
| S-EPMC2586547 | biostudies-literature
Project description:The interaction of ?-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an ?-helix within a membrane, or that change the probability that the ?-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an ?-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.
| S-EPMC3077683 | biostudies-literature
Project description:The environment of a living cell is vastly different from that of an in vitro reaction system, an issue that presents great challenges to the use of in vitro models, or computer simulations based on them, for understanding biochemistry in vivo. Virus capsids make an excellent model system for such questions because they typically have few distinct components, making them amenable to in vitro and modeling studies, yet their assembly can involve complex networks of possible reactions that cannot be resolved in detail by any current experimental technology. We previously fit kinetic simulation parameters to bulk in vitro assembly data to yield a close match between simulated and real data, and then used the simulations to study features of assembly that cannot be monitored experimentally. The present work seeks to project how assembly in these simulations fit to in vitro data would be altered by computationally adding features of the cellular environment to the system, specifically the presence of nucleic acid about which many capsids assemble. The major challenge of such work is computational: simulating fine-scale assembly pathways on the scale and in the parameter domains of real viruses is far too computationally costly to allow for explicit models of nucleic acid interaction. We bypass that limitation by applying analytical models of nucleic acid effects to adjust kinetic rate parameters learned from in vitro data to see how these adjustments, singly or in combination, might affect fine-scale assembly progress. The resulting simulations exhibit surprising behavioral complexity, with distinct effects often acting synergistically to drive efficient assembly and alter pathways relative to the in vitro model. The work demonstrates how computer simulations can help us understand how assembly might differ between the in vitro and in vivo environments and what features of the cellular environment account for these differences.
| S-EPMC4887116 | biostudies-literature
Project description:Understanding the fundamental biophysics behind protein-nanoparticle (NP) interactions is essential for the design and engineering bio-NP systems. The authors describe the development of a coarse-grained protein-NP model that utilizes a structure centric protein model. A key feature of the protein-NP model is the quantitative inclusion of the hydrophobic character of residues in the protein and their interactions with the NP surface. In addition, the curvature of the NP is taken into account, capturing the protein behavior on NPs of different size. The authors evaluate this model by comparison with experimental results for structure and adsorption of a model protein interacting with an NP. It is demonstrated that the simulation results recapitulate the structure of the small ?/? protein GB1 on the NP for data from circular dichroism and fluorescence spectroscopy. In addition, the calculated protein adsorption free energy agrees well with the experimental value. The authors predict the dependence of protein folding on the NP size, surface chemistry, and temperature. The model has the potential to guide NP design efforts by predicting protein behavior on NP surfaces with various chemical properties and curvatures.
| S-EPMC5551056 | biostudies-literature
Project description:The major protective coat of most viruses is a highly symmetric protein capsid that forms spontaneously from many copies of identical proteins. Structural and mechanical properties of such capsids, as well as their self-assembly process, have been studied experimentally and theoretically, including modeling efforts by computer simulations on various scales. Atomistic models include specific details of local protein binding but are limited in system size and accessible time, while coarse grained (CG) models do get access to longer time and length scales but often lack the specific local interactions. Multi-scale models aim at bridging this gap by systematically connecting different levels of resolution. Here, a CG model for CCMV (Cowpea Chlorotic Mottle Virus), a virus with an icosahedral shell of 180 identical protein monomers, is developed, where parameters are derived from atomistic simulations of capsid protein dimers in aqueous solution. In particular, a new method is introduced to combine the MARTINI CG model with a supportive elastic network based on structural fluctuations of individual monomers. In the parametrization process, both network connectivity and strength are optimized. This elastic-network optimized CG model, which solely relies on atomistic data of small units (dimers), is able to correctly predict inter-protein conformational flexibility and properties of larger capsid fragments of 20 and more subunits. Furthermore, it is shown that this CG model reproduces experimental (Atomic Force Microscopy) indentation measurements of the entire viral capsid. Thus it is shown that one obvious goal for hierarchical modeling, namely predicting mechanical properties of larger protein complexes from models that are carefully parametrized on elastic properties of smaller units, is achievable.
| S-EPMC3628857 | biostudies-literature
Project description:Coarse-grained (CG) models of large biomolecular complexes enable simulations of these systems over long timescales that are not accessible for atomistic molecular dynamics (MD) simulations. A systematic methodology, called essential dynamics coarse-graining (ED-CG), has been developed for defining coarse-grained sites in a large biomolecule. The method variationally determines the CG sites so that key dynamic domains in the protein are preserved in the CG representation. The original ED-CG method relies on a principal component analysis (PCA) of a MD trajectory. However, for many large proteins and multi-protein complexes such an analysis may not converge or even be possible. This work develops a new ED-CG scheme using an elastic network model (ENM) of the protein structure. In this procedure, the low-frequency normal modes obtained by ENM are used to define dynamic domains and to define the CG representation accordingly. The method is then applied to several proteins, such as the HIV-1 CA protein dimer, ATP-bound G-actin, and the Arp2/3 complex. Numerical results show that ED-CG with ENM (ENM-ED-CG) is much faster than ED-CG with PCA because no MD is necessary. The ENM-ED-CG models also capture functional essential dynamics of the proteins almost as well as those using full MD with PCA. Therefore, the ENM-ED-CG method may be better suited to coarse-grain a very large biomolecule or biomolecular complex that is too computationally expensive to be simulated by conventional MD, or when a high resolution atomic structure is not even available.
| S-EPMC2764101 | biostudies-literature
Project description:Kalata B1 has been demonstrated to have bioactivity relating to membrane disruption. In this study, we conducted coarse-grained molecular dynamics simulations to gain further insight into kB1 bioactivity. The simulations were performed at various concentrations of kB1 to capture the overall progression of its activity. Two configurations of kB1 oligomers, termed tower-like and wall-like clusters, were detected. The conjugation between the wall-like oligomers resulted in the formation of a ring-like hollow in the kB1 cluster on the membrane surface. Our results indicated that the molecules of kB1 were trapped at the membrane-water interface. The interfacial membrane binding of kB1 induced a positive membrane curvature, and the lipids were eventually extracted from the membrane through the kB1 ring-like hollow into the space inside the kB1 cluster. These findings provide an alternative view of the mechanism of kB1 bioactivity that corresponds with the concept of an interfacial bioactivity model.
| S-EPMC3910381 | biostudies-literature
Project description:Binding of HIV-1 capsid (CA) to cleavage and polyadenylation specificity factor 6 (CPSF6) is hypothesized to provide a significant fitness advantage to in vivo viral replication, explaining why CA-CPSF6 interactions are strictly conserved in primate lentiviruses. We recently identified a Q4R mutation in CA after propagation of an interferon (IFN)-β-hypersensitive CA mutant, RGDA/Q112D (H87R, A88G, P90D, P93A and Q112D) virus, in IFN-β-treated cells. The Q4R substitution conferred significant IFN-β resistance to the RGDA/Q112D virus by affecting several properties of the virus, including the sensitivity to myxovirus resistance protein B (MxB), the kinetics of reverse transcription, and the initiation of uncoating. Notably, the Q4R substitution restored the CPSF6 interaction of the RGDA/Q112D virus. To better understand how the Q4R substitution modulated the CA-CPSF6 interaction, we generated a series of CA mutants harboring substitutions at the 4th and 112th residues. In contrast to the effect in the RGDA/Q112D background, the Q4R substitution diminished CA-CPSF6 interaction in an otherwise wild-type virus. Our genetic and structural analyses revealed that while either the Q4R or Q112D substitution impaired CA-CPSF6 interaction, the combination of these substitutions restored this interaction. These results suggest that the 4th and 112th residues in HIV-1 CA cooperatively modulate CA-CPSF6 interactions, further highlighting the tremendous levels of plasticity in primate lentivirus CA, which is one of the barriers to antiretroviral therapy in HIV-1-infected individuals.
| S-EPMC7262650 | biostudies-literature
Project description:Physics-based potentials have been developed for the interactions between proteins and DNA for simulations with the UNRES + NARES-2P force field. The mean-field interactions between a protein and a DNA molecule can be divided into eight categories: (1) nonpolar side chain–DNA base, (2) polar uncharged side chain–DNA base, (3) charged side chain–DNA base, (4) peptide group–phosphate group, (5) peptide group–DNA base, (6) nonpolar side chain–phosphate group, (7) polar uncharged side chain–phosphate group, and (8) charged side chain–phosphate group. Umbrella-sampling molecular dynamics simulations in explicit TIP3P water using the AMBER force field were carried out to determine the potentials of mean force (PMF) for all 105 pairs of interacting components. Approximate analytical expressions for the mean-field interaction energy of each pair of the different kinds of interacting molecules were then fitted to the PMFs to obtain the parameters of the analytical expressions. These analytical expressions can reproduce satisfactorily the PMF curves corresponding to different orientations of the interacting molecules. The results suggest that the physics-based mean-field potentials of amino acid–nucleotide interactions presented here can be used in coarse-grained simulation of protein–DNA interactions.
| S-EPMC4455907 | biostudies-literature