CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed.

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published.

CHARMM-GUI Multicomponent Assembler for modeling and simulation of complex multicomponent systems.

Atomic-scale molecular modeling and simulation are powerful tools for computational biology. However, constructing models with large, densely packed molecules, non-water solvents, or with combinations of multiple biomembranes, polymers, and nanomaterials remains challenging and requires significant time and expertise. Furthermore, existing tools do not support such assemblies under the periodic boundary conditions (PBC) necessary for molecular simulation.

Pectin Nanoporous Structures Prepared via Salt-Induced Phase Separation and Ambient Azeotropic Evaporation Processes.

Polysaccharide nanoporous structures are suitable for various applications, ranging from biomedical scaffolds to adsorption materials, owing to their biocompatibility and large surface areas. Pectin, in particular, can create 3D nanoporous structures in aqueous solutions by binding with calcium cations and creating nanopores by phase separation; this process involves forming hydrogen bonds between alcohols and pectin chains in water and alcohol mixtures and the resulting penetration of alcohols into calcium-bound pectin gels. However, owing to the dehydration and condensation of polysaccharide chains during drying, it has proven to be challenging to maintain the 3D nanoporous structure without using a freeze-drying process or supercritical fluid.

Anillin-related Mid1 as an adaptive and multimodal contractile ring anchoring protein: A simulation study.

Cytokinesis of animal and fungi cells depends crucially on the anillin scaffold proteins. Fission yeast anillin-related Mid1 anchors cytokinetic ring precursor nodes to the membrane. However, it is unclear if both of its Pleckstrin Homology (PH) and C2 C-terminal domains bind to the membrane as monomers or dimers, and if one domain plays a dominant role.

Accurate Force Fields for Atomistic Simulations of Oxides, Hydroxides, and Organic Hybrid Materials up to the Micrometer Scale.

The simulation of metals, oxides, and hydroxides can accelerate the design of therapeutics, alloys, catalysts, cement-based materials, ceramics, bioinspired composites, and glasses. Here we introduce the INTERFACE force field (IFF) and surface models for α-AlO, α-CrO, α-FeO, NiO, CaO, MgO, β-Ca(OH), β-Mg(OH), and β-Ni(OH). The force field parameters are nonbonded, including atomic charges for Coulomb interactions, Lennard-Jones (LJ) potentials for van der Waals interactions with 12-6 and 9-6 options, and harmonic bond stretching for hydroxide ions.

CHARMM-GUI Multicomponent Assembler for Modeling and Simulation of Complex Multicomponent Systems.

Atomic-scale molecular modeling and simulation are powerful tools for computational biology. However, constructing models with large, densely packed molecules, non-water solvents, or with combinations of multiple biomembranes, polymers, and nanomaterials remains challenging and requires significant time and expertise. Furthermore, existing tools do not support such assemblies under the periodic boundary conditions (PBC) necessary for molecular simulation.

CHARMM-GUI : Past, Current, and Future Developments and Applications.

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes.

Anillin Related Mid1 as an Adaptive and Multimodal Contractile Ring Anchoring Protein: A Simulation Study.

The organization of the cytokinetic ring at the cell equator of dividing animal and fungi cells depends crucially on the anillin scaffold proteins. In fission yeast, anillin related Mid1 binds to the plasma membrane and helps anchor and organize a medial broad band of cytokinetic nodes, which are the precursors of the contractile ring. Similar to other anillins, Mid1 contains a C terminal globular domain with two potential regions for membrane binding, the Pleckstrin Homology (PH) and C2 domains, and an N terminal intrinsically disordered region that is strongly regulated by phosphorylation.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane.

() is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance.

Suppression of Lung Cancer Malignancy by Micellized siRNA through Cell Cycle Arrest.

UBA6-specific E2 conjugation enzyme 1 (USE1) is frequently overexpressed in lung cancer patients. Moreover, the critical role of USE1 in the progression of human lung cancer is also indicated. As the next step, the authors aim to develop USE1-targeted therapeutic agents based on RNA interference (RNAi).

Structural basis for Y2 receptor-mediated neuropeptide Y and peptide YY signaling.

Neuropeptide Y (NPY) and its receptors are expressed in various human tissues including the brain where they regulate appetite and emotion. Upon NPY stimulation, the neuropeptide Y1 and Y2 receptors (YR and YR, respectively) activate G signaling, but their physiological responses to food intake are different. In addition, deletion of the two N-terminal amino acids of peptide YY (PYY(3-36)), the endogenous form found in circulation, can stimulate YR but not YR, suggesting that YR and YR may have distinct ligand-binding modes.

Structural basis of neuropeptide Y signaling through Y1 receptor.

Neuropeptide Y (NPY) is highly abundant in the brain and involved in various physiological processes related to food intake and anxiety, as well as human diseases such as obesity and cancer. However, the molecular details of the interactions between NPY and its receptors are poorly understood. Here, we report a cryo-electron microscopy structure of the NPY-bound neuropeptide Y1 receptor (YR) in complex with G protein.

CHARMM-GUI Nanomaterial Modeler for Modeling and Simulation of Nanomaterial Systems.

Molecular modeling and simulation are invaluable tools for nanoscience that predict mechanical, physicochemical, and thermodynamic properties of nanomaterials and provide molecular-level insight into underlying mechanisms. However, building nanomaterial-containing systems remains challenging due to the lack of reliable and integrated cyberinfrastructures. Here we present Nanomaterial Modeler in CHARMM-GUI, a web-based cyberinfrastructure that provides an automated process to generate various nanomaterial models, associated topologies, and configuration files to perform state-of-the-art molecular dynamics simulations using most simulation packages.

CHARMM-GUI Membrane Builder for Lipid Nanoparticles with Ionizable Cationic Lipids and PEGylated Lipids.

A lipid nanoparticle (LNP) formulation is a state-of-the-art delivery system for genetic drugs such as DNA, messenger RNA, and small interfering RNA, which is successfully applied to COVID-19 vaccines and gains tremendous interest in therapeutic applications. Despite its importance, a molecular-level understanding of the LNP structures and dynamics is still lacking, which makes rational LNP design almost impossible. In this work, we present an extension of CHARMM-GUI to model and simulate all-atom LNPs with various (ionizable) cationic lipids and PEGylated lipids (PEG-lipids).

Dynamic Interactions of Fully Glycosylated SARS-CoV-2 Spike Protein with Various Antibodies.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work and performed all-atom molecular dynamics simulations to gain insight into the structural dynamics and interactions between S protein and antibodies.

Carrier-free micellar CpG interacting with cell membrane for enhanced immunological treatment of HIV-1.

Unmethylated CpG motifs activate toll-like receptor 9 (TLR9), leading to sequence- and species-specific immune stimulation. Here, we engineered a CpG oligodeoxyribonucleotide (ODN) with multiple hydrophobic moieties, so-called lipid-modified uracil, which resulted in a facile micelle formation of the stimulant. The self-assembled CpG nanostructure (U4CpG) containing the ODN 2216 sequence was characterized by various spectroscopic and microscopic methods together with molecular dynamics simulations.

CHARMM-GUI Membrane Builder for Lipid Nanoparticles with Ionizable Cationic Lipids and PEGylated Lipids.

A lipid nanoparticle (LNP) formulation is a state-of-the-art delivery system for genetic drugs such as DNA, mRNA, and siRNA, which is successfully applied to COVID-19 vaccines and gains tremendous interest in therapeutic applications. Despite its importance, a molecular-level understanding of the LNP structures and dynamics is still lacking, which makes a rational LNP design almost impossible. In this work, we present an extension of CHARMM-GUI to model and simulate all-atom LNPs with various (ionizable) cationic lipids and PEGylated lipids (PEG-lipids).

Dynamic Interactions of Fully Glycosylated SARS-CoV-2 Spike Protein with Various Antibodies.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 (ACE2) binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work, and performed all-atom molecular dynamics simulations to get insight into the structural dynamics and interactions between S protein and antibodies.

CHARMM-GUI Polymer Builder for Modeling and Simulation of Synthetic Polymers.

Molecular modeling and simulations are invaluable tools for polymer science and engineering, which predict physicochemical properties of polymers and provide molecular-level insight into the underlying mechanisms. However, building realistic polymer systems is challenging and requires considerable experience because of great variations in structures as well as length and time scales. This work describes in CHARMM-GUI (http://www.

Structure, Dynamics, Receptor Binding, and Antibody Binding of the Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein in a Viral Membrane.

The spike (S) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mediates host cell entry by binding to angiotensin-converting enzyme 2 (ACE2) and is considered the major target for drug and vaccine development. We previously built fully glycosylated full-length SARS-CoV-2 S protein models in a viral membrane including both open and closed conformations of the receptor-binding domain (RBD) and different templates for the stalk region. In this work, multiple μs-long all-atom molecular dynamics simulations were performed to provide deeper insights into the structure and dynamics of S protein and glycan functions.

Structure, Dynamics, and Interactions of GPI-Anchored Human Glypican-1 with Heparan Sulfates in a Membrane.

Glypican-1 and its heparan sulfate (HS) chains play important roles in modulating many biological processes including growth factor signaling. Glypican-1 is bound to a membrane surface via a glycosylphosphatidylinositol (GPI)-anchor. In this study, we used all-atom molecular modeling and simulation to explore the structure, dynamics, and interactions of GPI-anchored glypican-1, three HS chains, membranes, and ions.

Developing a Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein Model in a Viral Membrane.

This technical study describes all-atom modeling and simulation of a fully glycosylated full-length SARS-CoV-2 spike (S) protein in a viral membrane. First, starting from PDB: 6VSB and 6VXX, full-length S protein structures were modeled using template-based modeling, de-novo protein structure prediction, and loop modeling techniques in GALAXY modeling suite. Then, using the recently determined most occupied glycoforms, 22 N-glycans and 1 O-glycan of each monomer were modeled using Glycan Reader & Modeler in CHARMM-GUI.

Developing a Fully-glycosylated Full-length SARS-CoV-2 Spike Protein Model in a Viral Membrane.

This technical study describes all-atom modeling and simulation of a fully-glycosylated full-length SARS-CoV-2 spike (S) protein in a viral membrane. First, starting from PDB:6VSB and 6VXX, full-length S protein structures were modeled using template-based modeling, de-novo protein structure prediction, and loop modeling techniques in GALAXY modeling suite. Then, using the recently-determined most occupied glycoforms, 22 N-glycans and 1 O-glycan of each monomer were modeled using Glycan Reader & Modeler in CHARMM-GUI.

Formation of nanopores in DiynePC-DPPC complex lipid bilayers triggered by on-demand photo-polymerization.

Vesicles have unique characteristics that enable the release of drugs as well as encapsulation while maintaining biocompatibility. A photo-polymerizable liposome composed of 1,2-bis(10,12-tricosadiynoyl)--3-phosphocholine (23:2 DiynePC) has been investigated as vehicles for triggered delivery of drugs to cells. In this study, we confirmed for the first time that supported lipid bilayers (SLBs) prepared with a 1,2-dipalmitoyl--3-phosphocholine (DPPC)/DiynePC mixture generated pores 100-300 nm in size on the membrane after UV polymerization.