CHARMM‐GUI,
Nanodiscs are discoidal protein–lipid complexes that have wide applications in membrane protein studies. Modeling and simulation of nanodiscs are challenging due to the absence of structures of many membrane scaffold proteins (MSPs) that wrap around the membrane bilayer. We have developed CHARMM‐GUI
- NSF-PAR ID:
- 10462911
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Journal of Computational Chemistry
- Volume:
- 40
- Issue:
- 7
- ISSN:
- 0192-8651
- Page Range / eLocation ID:
- p. 893-899
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
http://www.charmm-gui.org , is a web‐based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM‐GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM‐GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1)PDB Reader & Manipulator ,Glycan Reader , andLigand Reader & Modeler for reading and modifying molecules; (2)Quick MD Simulator ,Membrane Builder ,Nanodisc Builder ,HMMM Builder ,Monolayer Builder ,Micelle Builder , andHex Phase Builder for building all‐atom simulation systems in various environments; (3)PACE CG Builder andMartini Maker for building coarse‐grained simulation systems; (4)DEER Facilitator andMDFF/xMDFF Utilizer for experimentally guided simulations; (5)Implicit Solvent Modeler ,PBEQ‐Solver , andGCMC/BD Ion Simulator for implicit solvent related calculations; (6)Ligand Binder for ligand solvation and binding free energy simulations; and (7)Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM‐GUI, such asGlycolipid Modeler for generation of various glycolipid structures, andLPS Modeler for generation of lipopolysaccharide structures from various Gram‐negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM‐GUI development project. © 2016 Wiley Periodicals, Inc. -
null (Ed.)We developed coarse-grained models of spike proteins in SARS-CoV-2 coronavirus and angiotensin-converting enzyme 2 (ACE2) receptor proteins to study the endocytosis of a whole coronavirus under physiologically relevant spatial and temporal scales. We first conducted all-atom explicit-solvent molecular dynamics simulations of the recently characterized structures of spike and ACE2 proteins. We then established coarse-grained models using the shape-based coarse-graining approach based on the protein crystal structures and extracted the force field parameters from the all-atom simulation trajectories. To further analyze the coarse-grained models, we carried out normal mode analysis of the coarse-grained models to refine the force field parameters by matching the fluctuations of the internal coordinates with the original all-atom simulations. Finally, we demonstrated the capability of these coarse-grained models by simulating the endocytosis of a whole coronavirus through the host cell membrane. We embedded the coarse-grained models of spikes on the surface of the virus envelope and anchored ACE2 receptors on the host cell membrane, which is modeled using a one-particle-thick lipid bilayer model. The coarse-grained simulations show the spike proteins adopt bent configurations due to their unique flexibility during their interaction with the ACE2 receptors, which makes it easier for them to attach to the host cell membrane than rigid spikes.more » « less
-
Computational methodologies are increasingly addressing modeling of the whole cell at the molecular level. Proteins and their interactions are the key component of cellular processes. Techniques for modeling protein interactions, thus far, have included protein docking and molecular simulation. The latter approaches account for the dynamics of the interactions but are relatively slow, if carried out at all-atom resolution, or are significantly coarse grained. Protein docking algorithms are far more efficient in sampling spatial coordinates. However, they do not account for the kinetics of the association (i.e., they do not involve the time coordinate). Our proof-of-concept study bridges the two modeling approaches, developing an approach that can reach unprecedented simulation timescales at all-atom resolution. The global intermolecular energy landscape of a large system of proteins was mapped by the pairwise fast Fourier transform docking and sampled in space and time by Monte Carlo simulations. The simulation protocol was parametrized on existing data and validated on a number of observations from experiments and molecular dynamics simulations. The simulation protocol performed consistently across very different systems of proteins at different protein concentrations. It recapitulated data on the previously observed protein diffusion rates and aggregation. The speed of calculation allows reaching second-long trajectories of protein systems that approach the size of the cells, at atomic resolution.
-
Protein–peptide interactions play a crucial role in a variety of cellular processes. The protein–peptide complex structure is a key to understand the mechanisms underlying protein–peptide interactions and is critical for peptide therapeutic development. We present a user‐friendly protein–peptide docking server, MDockPeP. Starting from a peptide sequence and a protein receptor structure, the MDockPeP Server globally docks the all‐atom, flexible peptide to the protein receptor. The produced modes are then evaluated with a statistical potential‐based scoring function, ITScorePeP. This method was systematically validated using the peptiDB benchmarking database. At least one near‐native peptide binding mode was ranked among top 10 (or top 500) in 59% (85%) of the bound cases, and in 40.6% (71.9%) of the challenging unbound cases. The server can be used for both protein–peptide complex structure prediction and initial‐stage sampling of the protein–peptide binding modes for other docking or simulation methods. MDockPeP Server is freely available at
http://zougrouptoolkit.missouri.edu/mdockpep . © 2018 Wiley Periodicals, Inc. -
Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.more » « less