Hubbry Logo
CHARMMCHARMMMain
Open search
CHARMM
Community hub
CHARMM
logo
7 pages, 0 posts
0 subscribers
Be the first to start a discussion here.
Be the first to start a discussion here.
CHARMM
CHARMM
CHARMM
View on Wikipedia
from Wikipedia
CHARMM
DevelopersMartin Karplus, Accelrys
Initial release1983; 42 years ago (1983)
Stable release
c47b1 / 2022; 3 years ago (2022)[1]
Preview release
c48a1 / 2022; 3 years ago (2022)[1]
Written inFORTRAN 77-95, CUDA
Operating systemUnix-like: Linux, macOS, AIX, iOS[2]
Platformx86, ARM, Nvidia GPU; Cray XT4, XT5[2]
Available inEnglish
TypeMolecular dynamics
LicenseProprietary
Websitewww.academiccharmm.org

Chemistry at Harvard Macromolecular Mechanics (CHARMM) is the name of a widely used set of force fields for molecular dynamics, and the name for the molecular dynamics simulation and analysis computer software package associated with them.[3][4][5] The CHARMM Development Project involves a worldwide network of developers working with Martin Karplus and his group at Harvard to develop and maintain the CHARMM program. Licenses for this software are available, for a fee, to people and groups working in academia.

Force fields

[edit]

The CHARMM force fields for proteins include: united-atom (sometimes termed extended atom) CHARMM19,[6] all-atom CHARMM22[7] and its dihedral potential corrected variant CHARMM22/CMAP, as well as later versions CHARMM27 and CHARMM36 and various modifications such as CHARMM36m and CHARMM36IDPSFF.[8] In the CHARMM22 protein force field, the atomic partial charges were derived from quantum chemical calculations of the interactions between model compounds and water. Furthermore, CHARMM22 is parametrized for the TIP3P explicit water model. Nevertheless, it is often used with implicit solvents. In 2006, a special version of CHARMM22/CMAP was reparametrized for consistent use with implicit solvent GBSW.[9]

The CHARMM22 force field has the following potential energy function:[7][10]

The bond, angle, dihedral, and nonbonded terms are similar to those found in other force fields such as AMBER. The CHARMM force field also includes an improper term accounting for out-of-plane bending (which applies to any set of four atoms that are not successively bonded), where is the force constant and is the out-of-plane angle. The Urey-Bradley term is a cross-term that accounts for 1,3 nonbonded interactions not accounted for by the bond and angle terms; is the force constant and is the distance between the 1,3 atoms.

For DNA, RNA, and lipids, CHARMM27[11] is used. Some force fields may be combined, for example CHARMM22 and CHARMM27 for the simulation of protein-DNA binding. Also, parameters for NAD+, sugars, fluorinated compounds, etc., may be downloaded. These force field version numbers refer to the CHARMM version where they first appeared, but may of course be used with subsequent versions of the CHARMM executable program. Likewise, these force fields may be used within other molecular dynamics programs that support them.

In 2009, a general force field for drug-like molecules (CGenFF) was introduced. It "covers a wide range of chemical groups present in biomolecules and drug-like molecules, including a large number of heterocyclic scaffolds".[12] The general force field is designed to cover any combination of chemical groups. This inevitably comes with a decrease in accuracy for representing any particular subclass of molecules. Users are repeatedly warned in Mackerell's website not to use the CGenFF parameters for molecules for which specialized force fields already exist (as mentioned above for proteins, nucleic acids, etc.).

CHARMM also includes polarizable force fields using two approaches. One is based on the fluctuating charge (FQ) model, also termed Charge Equilibration (CHEQ).[13][14] The other is based on the Drude shell or dispersion oscillator model.[15][16]

Parameters for all of these force fields may be downloaded from the Mackerell website for free.[17]

Molecular dynamics program

[edit]

The CHARMM program allows for generating and analysing a wide range of molecular simulations. The most basic kinds of simulation are minimizing a given structure and production runs of a molecular dynamics trajectory. More advanced features include free energy perturbation (FEP), quasi-harmonic entropy estimation, correlation analysis and combined quantum, and quantum mechanics–molecular mechanics (QM/MM) methods.

CHARMM is one of the oldest programs for molecular dynamics. It has accumulated many features, some of which are duplicated under several keywords with slight variants. This is an inevitable result of the many outlooks and groups working on CHARMM worldwide. The changelog file, and CHARMM's source code, are good places to look for the names and affiliations of the main developers. The involvement and coordination by Charles L. Brooks III's group at the University of Michigan is salient.

Software history

[edit]

Around 1969, there was considerable interest in developing potential energy functions for small molecules. CHARMM originated at Martin Karplus's group at Harvard. Karplus and his then graduate student Bruce Gelin decided the time was ripe to develop a program that would make it possible to take a given amino acid sequence and a set of coordinates (e.g., from the X-ray structure) and to use this information to calculate the energy of the system as a function of the atomic positions. Karplus has acknowledged the importance of major inputs in the development of the (at the time nameless) program, including:

  • Schneior Lifson's group at the Weizmann Institute, especially from Arieh Warshel who went to Harvard and brought his consistent force field (CFF) program with him
  • Harold Scheraga's group at Cornell University
  • Awareness of Michael Levitt's pioneering energy calculations for proteins

In the 1980s, finally a paper appeared and CHARMM made its public début. Gelin's program had by then been considerably restructured. For the publication, Bob Bruccoleri came up with the name HARMM (HARvard Macromolecular Mechanics), but it seemed inappropriate. So they added a C for Chemistry. Karplus said: "I sometimes wonder if Bruccoleri's original suggestion would have served as a useful warning to inexperienced scientists working with the program."[18] CHARMM has continued to grow and the latest release of the executable program was made in 2015 as CHARMM40b2.

Running CHARMM under Unix-Linux

[edit]

The general syntax for using the program is:

charmm -i filename.inp -o filename.out

  • charmm – The name of the program (or script which runs the program) on the computer system being used.
  • filename.inp – A text file which contains the CHARMM commands. It starts by loading the molecular topologies (top) and force field (par). Then one loads the molecular structures' Cartesian coordinates (e.g. from PDB files). One can then modify the molecules (adding hydrogens, changing secondary structure). The calculation section can include energy minimization, dynamics production, and analysis tools such as motion and energy correlations.
  • filename.out – The log file for the CHARMM run, containing echoed commands, and various amounts of command output. The output print level may be increased or decreased in general, and procedures such as minimization and dynamics have printout frequency specifications. The values for temperature, energy pressure, etc. are output at that frequency.

Volunteer computing

[edit]

Docking@Home, hosted by University of Delaware, one of the projects which use an open-source platform for the distributed computing, BOINC, used CHARMM to analyze the atomic details of protein-ligand interactions in terms of molecular dynamics (MD) simulations and minimizations.

World Community Grid, sponsored by IBM, ran a project named The Clean Energy Project[19] which also used CHARMM in its first phase, which has been completed.

See also

[edit]

References

[edit]
[edit]
Revisions and contributorsEdit on WikipediaRead on Wikipedia

CHARMM

View on Grokipedia
from Grokipedia
CHARMM (Chemistry at HARvard Macromolecular Mechanics) is a versatile molecular simulation program designed for atomic-level modeling of biomolecular systems, including proteins, nucleic acids, , carbohydrates, and ligands in environments such as solutions, crystals, and membranes. It employs a comprehensive set of empirical force fields, along with hybrid methods, to perform simulations that elucidate structure, dynamics, and thermodynamics. Developed initially in the late 1970s at by Bruce Gelin and for studies of macromolecules like , CHARMM has evolved into a flexible, extensible tool supporting techniques such as (MD), energy minimization, free energy calculations, analysis, and advanced sampling methods like replica-exchange MD. The program's history traces back to its first formal description in 1983, marking the transition from precursor efforts to a robust framework for classical and semiempirical simulations. Over the subsequent decades, CHARMM has been continuously enhanced by a global developer community under the long-term leadership of , with periodic releases managed through version control systems since 1994, incorporating features like , Ewald summation for electrostatics, parallel computing support, and interfaces to quantum chemistry packages such as GAMESS, Gaussian, and Q-Chem. The latest version, c49b2 (as of 2024), includes enhancements in accessibility, functionality, and community tools. Key innovations include multi-scale modeling (e.g., MM/coarse-grained hybrids), implicit and explicit solvent representations, and tools for model building and analysis, enabling high-performance computations on clusters and GPUs. Academic users can access CHARMM freely upon registration, while it remains commercially available through . CHARMM's applications span , , and , facilitating investigations into , , ligand binding, , and large-scale complexes like the system. It integrates with experimental data from and NMR for atomic-resolution structure refinement and supports conformational sampling, path integrals, and methods to probe phenomena inaccessible to direct observation. Associated resources, such as CHARMM-GUI for system preparation and forums for user support, further enhance its utility in research and education.

Overview

Definition and Purpose

CHARMM, an acronym for Chemistry at HARvard Macromolecular Mechanics, is a versatile molecular simulation program designed for modeling biomolecular systems, including proteins, nucleic acids, , and carbohydrates, using approaches. It enables detailed investigations into the structure, dynamics, and interactions of these systems at atomic resolution, supporting applications in computational and . The core of CHARMM consists of empirical force fields that define functions for biomolecular interactions and a computational program that implements algorithms for energy minimization, (MD) simulations, and calculations. These components allow users to perform energy evaluations and manipulations essential for simulating conformational changes, ligand binding, and thermodynamic properties in complex macromolecular environments. As one of the first comprehensive biomolecular simulation packages, CHARMM has facilitated pioneering studies of biomolecular behavior since its inception, providing a foundational tool for atomic-level modeling that integrates empirical potentials with advanced simulation techniques. The general form of the CHARMM potential energy function, UU, captures these interactions through additive terms: U=∑bondskb(r−r0)2+∑angleskθ(θ−θ0)2+∑dihedralskϕ(1+cos⁡(nϕ−δ))+∑nonbonded(qiqjrij+Aijrij12−Bijrij6)U = \sum_{\text{bonds}} k_b (r - r_0)^2 + \sum_{\text{angles}} k_\theta (\theta - \theta_0)^2 + \sum_{\text{dihedrals}} k_\phi (1 + \cos(n\phi - \delta)) + \sum_{\text{nonbonded}} \left( \frac{q_i q_j}{r_{ij}} + \frac{A_{ij}}{r_{ij}^{12}} - \frac{B_{ij}}{r_{ij}^6} \right) Here, the first three sums represent bonded interactions—harmonic potentials for bond lengths (rr, equilibrium r0r_0, force constant kbk_b), bond angles (θ\theta, equilibrium θ0\theta_0, kθk_\theta), and periodic dihedral angles (ϕ\phi, multiplicity nn, phase δ\delta, kϕk_\phi)—while the nonbonded sum includes Coulombic (qi,qjq_i, q_j charges, distance rijr_{ij}) and Lennard-Jones van der Waals terms (AijA_{ij} repulsive, BijB_{ij} attractive parameters).

Licensing and Availability

CHARMM is a molecular simulation program originally developed at and commercially licensed through (formerly Accelrys). The academic version, known as CHARMM, became freely available to academic, government, and non-profit users starting in , distributed via the official site academiccharmm.org without any licensing fees for eligible institutions. In contrast, for-profit entities must acquire commercial s for the CHARMm variant directly from , ensuring controlled access to its full capabilities in industrial applications. The software remains overall, with no open-source release, though the academic distribution includes comprehensive access to its features for non-commercial research. Academic users gain access by registering at brooks.chem.lsa.umich.edu/register, after which they can download the complete release package containing , , test cases, topology and parameter files, and pre-built binaries for select platforms. Commercial access involves contacting for tailored licensing agreements, often integrated into broader software suites like . Building from source requires a 95-compliant compiler, such as GCC gfortran (version 4.4 or later, excluding 4.5.1), ifort (11.1 or later), or PGI pgf95 (11.1 or later), along with MPI and for parallel execution. The package unpacks into a directory like ~/c50b1 for version c50b1, with installation handled via configure scripts and make commands. CHARMM primarily supports Unix/ environments, with confirmed compatibility for platforms including em64t, gnu , osx (macOS), and GPU-accelerated systems via interfaces like DOMDEC-GPU and OpenMM. Binaries are available for macOS and certain distributions, while Windows users typically compile from source or use compatibility tools like , as native binaries are not standard. follow a cXX naming convention, such as c48a1 in 2022 or the current c50b1 as of 2025, with major releases occurring annually; detailed changelogs outlining enhancements and fixes are hosted on the site. Community resources at academiccharmm.org include extensive covering installation, usage, and advanced features, along with tutorials for setup on various platforms. User support is facilitated through dedicated forums at forums-academiccharmm.org, where researchers discuss installation issues, share best practices, and access developer guides. This , enhanced by the 2022 shift to free academic access, has broadened CHARMM's reach within the .

History

Origins and Early Development

CHARMM, or Chemistry at HARvard Macromolecular Mechanics, was initiated by in the early 1970s at as a computational tool initially designed for simulating protein structures and dynamics. The program's inception stemmed from Karplus's visit to Schneior Lifson's group at the Weizmann Institute in 1969, where there was growing interest in developing empirical potential energy functions to model the conformations of small molecules and extend these approaches to larger biomolecules. At the time, quantum mechanical calculations were computationally prohibitive for systems as complex as proteins, necessitating the use of classical empirical potentials to approximate intramolecular interactions and enable studies of structural perturbations, such as those induced by ligand binding in . Early development of CHARMM was driven by the need to bridge the gap between static data and dynamic behavior in biological macromolecules, with initial efforts focusing on energy minimization and normal mode analysis for proteins. Key collaborators included graduate students Bruce Gelin, who contributed significantly to the program's coding and implementation, and J. Andrew McCammon, who helped pioneer its application to . What began as scripts for specific calculations evolved into a more structured software package, emphasizing for handling atomic coordinates, force field parameters, and simulation algorithms. The initial scope was narrow, targeting proteins using simple empirical force fields that parameterized bonded and non-bonded interactions based on available experimental data. The program's first major milestone came in 1977 with the publication of the inaugural simulation of a protein, the bovine pancreatic inhibitor (BPTI), conducted using an early version of CHARMM. This simulation, spanning just 9.2 picoseconds, demonstrated the feasibility of capturing atomic fluctuations in a vacuum environment and revealed dynamic elements like hydrogen bonding networks that were invisible in static structures. Running on mainframe computers such as the , these early computations were severely limited by hardware constraints, including slow processing speeds and modest , restricting simulations to short timescales and small systems of a few hundred atoms. CHARMM remained an in-house tool at Harvard for research purposes until its public debut in 1983 as version c19, marking the transition to a distributable package for broader scientific use.

Key Milestones and Versions

The development of CHARMM began in the late , with the first formal releases occurring in the under versions c20 through c25, which introduced core capabilities for energy minimization and simulations of proteins, nucleic acids, and crystalline solids. These early versions, such as c20, laid the foundation for biomolecular modeling by supporting isolated molecules, solutions, and solids, with initial force fields like PARAM19 providing polar hydrogen representations for proteins and nucleic acids. In the 1990s, CHARMM advanced through versions c26 to c30, incorporating lipid parameters to enable simulations of membrane systems and enhancing nucleic acid support with the CHARMM27 force field in 1998, which improved accuracy for DNA and RNA structures. Key releases included c26 in 1998 and c27 in 2000, alongside the introduction of targeted molecular dynamics in 1993 for studying conformational transitions. The 2000s saw versions c31 to c36, marked by the addition of cross-term map (CMAP) corrections in 2004 via c30a1 to better capture protein backbone dihedral interactions, significantly enhancing simulation fidelity for folded states. This period also initiated a shift toward polarizable force fields, with the oscillator model prototyped by 2007 in c34b1 for inducible dipoles in biomolecules, and support for systems scaling to 10^10 atoms in c31b1 by 2003. Lipid force fields were refined in 2005, building on parameters for bilayers. During the 2010s, versions progressed from c37 to c41, with CHARMM36 released in 2012 featuring optimized CMAP terms for proteins, , and nucleic acids, improving agreement with NMR data and membrane properties. Polarizable models expanded with Drude-2013 for proteins, and academic licensing began broadening access. In , received the , shared with Michael Levitt and , for techniques that underpinned CHARMM's foundational simulations of chemical reactions in proteins. The 2020s brought versions c42 to c50, including developmental builds up to c50a1 in 2024 and releases like c49b1, integrating GPU acceleration through the CHARMM/ introduced in c37b1 and advanced with domain decomposition in 2014 for faster . CHARMM became freely available for academic and non-profit use starting in August 2022, expanding accessibility via platforms like academiccharmm.org. Polarizable force fields continued evolving, with Drude-2023 for lipids and bilayers. , the longtime leader of CHARMM development, passed away on December 28, 2024. As of November 2025, CHARMM has received minor patches for compatibility with emerging hardware like advanced GPUs, without a major force field overhaul, maintaining stability across c50 series builds.

Force Fields

Additive Force Fields

The additive force fields in CHARMM represent the standard non-polarizable models, utilizing fixed atomic partial charges and Lennard-Jones parameters to describe electrostatic and van der Waals interactions, respectively, without accounting for inducible polarization effects. These force fields form the core of CHARMM's empirical potential energy function, enabling efficient simulations of biomolecular systems by balancing computational cost with accuracy in reproducing structural and thermodynamic properties. For proteins, the CHARMM22 force field, released in 2002, marked a significant advancement in all-atom modeling, with the subsequent addition of the Cross-term map (CMAP) correction in to better capture backbone dihedral energetics and improve secondary structure stability, such as alpha-helices and beta-sheets. Building on this, the CHARMM36m force field, introduced in 2017, refines protein parameters through targeted adjustments to dihedral and non-bonded terms, enhancing performance for both folded domains and intrinsically disordered regions by achieving closer agreement with experimental NMR chemical shifts, residual dipolar couplings, and profiles. Nucleic acid simulations rely on the CHARMM27 force field, released in 2004, which provides optimized parameters for DNA and RNA, including glycosidic torsion potentials that stabilize helical conformations and base stacking interactions. For lipids, the CHARMM36 force field, developed in 2012, incorporates refined aliphatic chain parameters and headgroup interactions to accurately reproduce phase transition temperatures, bilayer thicknesses, and area per lipid in simulations of phosphatidylcholine and other membrane lipids. The CHARMM General Force Field (CGenFF), introduced in 2009, extends the additive framework to drug-like small molecules and organic ligands, covering a broad range of functional groups compatible with biomolecular parameters, and supports automated parameterization through the CGenFF server for rapid generation. The update, CGenFF version 5.0 (published 2025), expands the training set by adding 1,390 new molecules to the previous approximately 930, resulting in over 2,300 molecules total, improving charge assignment and bonded terms for better prediction of intramolecular geometries and non-covalent binding affinities. Validation of these additive force fields emphasizes quantitative comparisons with experimental data, including NMR-derived order parameters and J-couplings for proteins, diffraction-derived densities for lipid bilayers, and thermodynamic quantities like free energies of for small molecules, where CHARMM36m and CGenFF achieve root-mean-square deviations of approximately 2 kcal/mol for solvation free energies and similar accuracy for other key observables. Early limitations in monovalent parameters, such as overestimation of Na⁺ hydration free energies, have been mitigated in updates through quantum mechanical refinements and experimental calibration against osmotic pressures and ion-DNA binding constants.

Polarizable Force Fields

CHARMM incorporates polarizable force fields to account for induced electronic polarization, which allows for more accurate modeling of environmental effects on molecular interactions compared to fixed-charge additive models. These force fields dynamically adjust electrostatic properties in response to the local , improving simulations of complex systems such as biomolecular interfaces and ionic environments. The primary polarizable model in CHARMM is the Drude oscillator approach, where atomic is represented by attaching a positively charged "Drude particle" to each non-hydrogen atom via a virtual harmonic spring; this particle oscillates in response to external , mimicking the displacement of clouds. The force field includes additional terms for induced interactions between these oscillators, screened using Thole's damping to prevent polarization catastrophe. The polarization energy contribution is given by Upol=∑i12kd(rd−r0)2+∑i,jqiqj′rij,U_{\text{pol}} = \sum_i \frac{1}{2} k_d (r_d - r_0)^2 + \sum_{i,j} \frac{q_i q_j'}{r_{ij}}, where kdk_d is the spring constant, rdr_d and r0r_0 are the Drude particle position and equilibrium distance, and qj′q_j' denotes charges including the induced Drude charges. An alternative polarizable model in CHARMM is the fluctuating charge (FQ) approach, which allows partial atomic charges to vary dynamically based on electronegativity equalization principles, enabling charge transfer and polarization effects without additional particles. This model derives from density functional theory-inspired charge responses and has been parameterized for proteins and organic liquids. Key implementations include the -2013 force field, developed for proteins and models like SWM4-NDP, which explicitly treats for and nucleic acids. Extensions to emerged in the 2020s, with Drude polarizable parameters for phospholipids like DPPC, enabling simulations of biomembranes with explicit long-range . These polarizable models incur approximately 2-3 times the computational cost of additive force fields due to the extra and extended electrostatic calculations. Polarizable force fields in CHARMM offer advantages in capturing electronic effects at protein-ion interfaces, lipid-water boundaries, and even in excited states through integrations, providing superior accuracy over additive models in these regimes. Validation studies demonstrate close agreement with quantum mechanical calculations for dipole moments, solvation free energies, and interaction energies, such as ion-protein binding affinities and responses.

Parameterization and Validation

Parameter derivation in CHARMM force fields primarily relies on quantum mechanical (QM) calculations to determine bonded parameters such as bond and angle force constants, which are fitted to surfaces obtained from high-level methods like MP2/6-31G(d). These QM targets ensure accurate representation of intramolecular interactions, with geometries optimized and vibrational frequencies scaled to match experimental spectra where available. Empirical fitting complements this by adjusting nonbonded parameters, such as Lennard-Jones terms, to reproduce experimental observables including liquid densities and heats of vaporization from pure solvent simulations. For example, in the development of the CHARMM General Force Field (CGenFF), partial charges are derived from QM electrostatic potentials and refined against experimental thermodynamic to enhance compatibility with biomolecular simulations. Tools like FFParam facilitate this process by automating the optimization of electrostatic and bonded parameters for both additive and polarizable models, integrating QM target data for geometry and energy scans alongside empirical condensed-phase properties such as free energies. The CGenFF server provides an accessible platform for parameterizing small molecules, employing QM calculations for charges and conformational energies while targeting experimental densities and vibrational spectra to generate transferable parameters compatible with CHARMM biomolecular force fields. Validation of CHARMM parameters involves direct comparison to experimental observables, such as radii of gyration from (SAXS) for disordered proteins and helix propensities assessed via NMR chemical shifts and J-couplings, ensuring structural accuracy across folded and unfolded states. Benchmarking against other force fields, like ff99SB-ILDN, reveals CHARMM36m's competitive performance in reproducing experimental order parameters and secondary structure distributions, though AMBER variants sometimes show lower deviations in gyration radii for . Key metrics include root-mean-square error (RMSE) for hydration free energies, typically around 2.04 kcal/mol, and Pearson correlation coefficients exceeding 0.88 for structural alignments, indicating robust predictive power. Early challenges in CHARMM lipid force fields, such as overestimation of chain ordering in saturated leading to gel-like bilayers in versions like C27r, were addressed through targeted refinements in C36, including adjustments to torsional and nonbonded parameters based on QM and experimental bilayer data, resulting in surface areas within 2% of experiment. The 2025 release of CGenFF v5.0 further improves small-molecule transferability by expanding the training set by adding 1,390 new compounds to the previous approximately 930, resulting in over 2,300 compounds total, enhancing agreement with QM geometries, vibrations, and dipole moments while maintaining low errors in solvent properties. Ongoing refinements incorporate community feedback through the MacKerell lab's parameter repository, iteratively updating parameters to resolve discrepancies in diverse chemical spaces.

Software Features

Molecular Dynamics Capabilities

CHARMM employs the Verlet/leap-frog as its primary algorithm for propagating trajectories, enabling the simulation of atomic motions under Newtonian mechanics. This , specified via the DYNAmics command with the LEAP keyword, updates positions and velocities in a staggered manner, offering stability and suitable for biomolecular systems. For energy minimization prior to dynamics, CHARMM supports the steepest descent (SD) method, which rapidly reduces high-energy configurations by following the negative gradient of the , and the conjugate gradient (CONJ) technique, which converges more efficiently for refined optimizations by incorporating curvature information. These minimization algorithms are invoked through the MINImize command and are essential for preparing stable starting structures. Advanced simulation methods in CHARMM extend beyond standard dynamics to address complex thermodynamic and reactive processes. Free energy perturbation (FEP) calculations, implemented via the PERTurb command, allow estimation of free energy differences by scaling interactions between perturbed states, often used for alchemical transformations like ligand binding. Umbrella sampling, facilitated by the UMBRel command, applies biasing potentials along a reaction coordinate to enhance sampling of rare events, enabling the reconstruction of potential of mean force profiles. For regions involving chemical reactivity, CHARMM integrates quantum mechanics/molecular mechanics (QM/MM) hybrid approaches through the QMMM module, treating active sites quantum mechanically (e.g., via semiempirical methods like PM6) while the surrounding environment uses classical force fields. Boundary conditions in CHARMM simulations accommodate diverse system sizes and environments. (PBC), defined using the command, replicate the simulation cell to mimic bulk phases, with long-range electrostatics handled by invoked via the EWALD keyword in nonbonded options for accurate treatment of charged systems. For solvated biomolecules, stochastic boundary molecular dynamics (SBMD) confines dynamics to a reaction region with Langevin friction and random forces at the boundary, reducing computational cost while maintaining realistic solvation effects. On modern hardware, CHARMM supports molecular dynamics simulations spanning nanosecond (ns) to microsecond (μs) timescales, particularly for systems up to tens of thousands of atoms, leveraging optimized integrators and parallelization. Implicit solvent models, such as generalized Born (GB) with solvent-accessible surface area (SA) nonpolar terms, are available via the GBNP command, approximating solvation without explicit water molecules to accelerate longer runs. CHARMM simulations are scripted using stream files with the .inp extension, which define , coordinates, parameters, and execution steps in a command-based syntax. A basic run typically begins with reading (READ RTFs) and parameter (READ PARAmeters) files, followed by generating structure (GENERate), assigning coordinates (READ COORdinates), minimizing energy (MINImize), and initiating dynamics (DYNAmics) with specified timestep, steps, and output frequencies, concluding with coordinate writes (WRITE COORdinates). For example:

* Basic MD Example READ RTFS CARD TOP_ALL36_PROT.RTF READ PARA CARD PAR_ALL36_PROT.PAR GENER SEGID PROT RESI 1 100 READ COOR CARD COORDS.PDB MINI SD NSTEP 1000 DYNA LEAP NSTEP 10000 TIMESTEP 0.002 \ IPRFRQ 1000 IUNCRD 20 NTWF 1000 \ NTWE 1000 WRITE COOR CARD DCD OUT.DCD STOP

* Basic MD Example READ RTFS CARD TOP_ALL36_PROT.RTF READ PARA CARD PAR_ALL36_PROT.PAR GENER SEGID PROT RESI 1 100 READ COOR CARD COORDS.PDB MINI SD NSTEP 1000 DYNA LEAP NSTEP 10000 TIMESTEP 0.002 \ IPRFRQ 1000 IUNCRD 20 NTWF 1000 \ NTWE 1000 WRITE COOR CARD DCD OUT.DCD STOP

This structure ensures reproducible, modular workflows for dynamics propagation.

Analysis and Utility Tools

CHARMM provides a suite of built-in tools for analyzing molecular dynamics (MD) trajectories, enabling researchers to extract structural and dynamic insights from simulation outputs. The COOR module facilitates root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) calculations, which quantify structural deviations and atomic fluctuations relative to a reference structure. For instance, the coor orient rms command aligns selected atoms, such as alpha carbons, and computes RMSD values across trajectory frames, while RMSF is derived by averaging deviations over time for each residue. Hydrogen bonding analysis is supported via the coor hbond command, which identifies donor-acceptor pairs based on geometric criteria (e.g., distance < 2.4 Å and angle > 120°) and outputs statistics like average bond counts and lifetimes for intra- or intermolecular interactions. Secondary structure assignment employs DSSP-like algorithms through the coor secs command, classifying residues into helices, sheets, or coils based on hydrogen bonding patterns and dihedral angles, with options to track temporal evolution in trajectories. Energy decomposition tools in CHARMM allow dissection of the potential energy into contributions from specific residues or atom groups, aiding in the identification of stabilizing interactions. The INTEraction command computes pairwise interaction energies (e.g., van der Waals and electrostatic) between selected subsets, such as a ligand and protein residues, while the ENERGY module extends this to per-residue breakdowns by summing intra- and intermolecular terms for each residue. Correlation functions for dynamics are handled by the CORREL module, which processes time series data from trajectories to compute autocorrelation functions for quantities like dihedral angles or energies, revealing timescales of motions (e.g., via exponential fitting). These tools support quasi-harmonic analysis through the VIBRAN facility, which derives covariance matrices from trajectory fluctuations to estimate entropic contributions and low-frequency modes. Utility functions in CHARMM streamline preprocessing and postprocessing tasks through its internal , which supports conditional statements, loops, variable substitution, and subroutine calls for automating workflows. PDB file manipulation is achieved with READ and WRITE COOR PDB commands, allowing atom selection, renumbering, and formatting adjustments, while the IC (internal coordinates) module enables mutations by parameterizing new residue topologies and refining geometries via energy minimization. Solvation box generation uses the SOLV command to add molecules within a defined spherical or cubic boundary around the solute, followed by ion placement via the IONize command to neutralize charge. These scripts can chain operations, such as building solvated systems from initial coordinates. Visualization integration is inherent in CHARMM's output formats, with trajectory data saved in DCD binary files compatible with external tools like VMD and PyMOL for interactive rendering of dynamics, hydrogen bonds, and secondary structures. Built-in plotting capabilities via the CORREL and GRAPHX modules generate graphs for energies, forces, and RMSD, outputting to text or files for further analysis. The GRAPHX facility supports basic 3D visualization with features like atom coloring and bond rendering, though it is often supplemented by external software. Recent additions since 2023 enhance CHARMM's extensibility through the pyCHARMM Python interface, which embeds core functionality into Python scripts for custom trajectory analyses, such as integrating for advanced statistical processing of RMSD/RMSF data. This interface facilitates hooks, exemplified by the MLPot module, which couples CHARMM force fields with potentials like PhysNet for enhanced sampling in free energy calculations, enabling on-the-fly potential corrections during . These developments, including support for Regression in QM/MM simulations via delta-ML potentials, broaden utility for complex workflows while maintaining compatibility with existing tools. As of 2024, further enhancements include apoCHARMM for GPU-accelerated simulations, the MIST approach for third-order conformational entropy calculations, and the COOR command for hydration maps, as detailed in version c50b1.

Implementation

Running CHARMM on Unix/Linux

CHARMM installation on Unix/Linux systems begins with downloading the source package from the official academic distribution site, academiccharmm.org, which provides access to the latest release, such as c50b1, including source files, , test cases, and / files. Unpack the tarball into a , typically ~/c50b1 or similar, ensuring sufficient disk space for compilation and libraries. Compilation requires a ; recommended options include gfortran version 4.4 or later (excluding 4.5.1) or ifort version 11.1 or later, with icc for C components if needed. To build, navigate to the unpacked directory and execute the configuration script, such as ./configure --with-gcc for gfortran or --with-intel for ifort, followed by make -jN -C build/cmake install using for modern builds, where N is the number of parallel jobs. Optional switches during configuration enable features like FFTW support via --enable-fftw or via --with-netcdf=/path/to/netcdf. The resulting executable, named charmm, is placed in the bin subdirectory, such as ~/c50b1/bin/charmm. Environment variables facilitate execution and customization. Set CHARMMEXEC to the full path of the compiled executable (e.g., export CHARMMEXEC=~/c50b1/bin/charmm) to simplify invocation from scripts or other tools. Additionally, include the and paths in PATH, and for optional libraries, define FFTW_HOME or NETCDF_DIR pointing to their installation directories (e.g., /usr/local/netcdf). These variables ensure CHARMM locates dependencies during runtime, particularly for I/O formats like coordinates. Basic execution of CHARMM on Unix/ uses command-line redirection for input and output files. The standard syntax is charmm < input.inp > output.out, where input.inp contains the sequence of CHARMM commands (starting with a * title line) and output.out captures the log and results. For interactive sessions, omit redirection and enter commands directly at the CHARMM prompt. Graphics output, if enabled via the OPEN GRAPH command in the input, requires X11 forwarding (e.g., ssh -X). CHARMM relies on specific file structures for molecular systems. Topology files, typically in Residue TOPology (.rtf) or extended .top format, define atom types, bonds, angles, and dihedrals for residues. Parameter files (.prm) provide force field constants like bond lengths and angles, loaded via READ PARA CARD or similar commands. Coordinate files specify atomic positions, commonly in (.pdb) format for initial structures or binary Coordinate (.crd) for dynamics trajectories. A typical loads these sequentially: READ RTF CARD topology.rtf, READ PARA CARD parameters.prm, READ COOR PDBATOMS coord.pdb. For batch scripting, wrap executions in a , such as:

#!/bin/bash export CHARMMEXEC=~/c50b1/bin/charmm $CHARMMEXEC < my_simulation.inp > my_simulation.out

#!/bin/bash export CHARMMEXEC=~/c50b1/bin/charmm $CHARMMEXEC < my_simulation.inp > my_simulation.out

This example runs a full simulation non-interactively, suitable for job schedulers like SLURM on Linux clusters. Troubleshooting common issues enhances reliability. Missing libraries often cause linking errors during compilation; for NetCDF, install via package managers (e.g., sudo apt install libnetcdf-dev on Ubuntu) and specify the path in configuration, as it supports advanced I/O for large trajectories. Similarly, FFTW is required for fast Fourier transforms in simulations; install with sudo yum install fftw-devel on CentOS/Rocky Linux and enable the switch to avoid "undefined reference" errors. Performance optimization involves compiler flags like -O3 -march=native passed via FFLAGS or FCFLAGS environment variables (e.g., export FFLAGS="-O3 -funroll-loops" before configure), which can accelerate builds by 20-30% on modern x86_64 hardware without altering correctness. If the executable fails to produce (e.g., due to mismatched MPI modules), clean the build directory with make clean and verify compiler consistency. Platform specifics vary across Linux distributions. On (e.g., 24.04 LTS), use apt for dependencies like gfortran, libfftw3-dev, and libnetcdf-dev, with configuration targeting gnu machine type for seamless integration. CentOS or its successors like 8/9 require dnf for packages such as gcc-gfortran and fftw-devel, often with compilers preferred for HPC environments due to better vectorization. For , as recommended in 2024 documentation, use Apptainer (successor to Singularity) over Docker for security in shared clusters; build from a base image like apptainer build charmm.sif image.def, binding data directories, to encapsulate CHARMM and dependencies portably across distributions. This approach avoids system conflicts and supports reproducible runs on or CentOS-based nodes.

Parallel and Distributed Computing

CHARMM supports parallel execution through the (MPI) for distributed-memory systems across multiple nodes and for shared-memory parallelism within nodes, enabling efficient scaling for simulations of large biomolecular systems. The domain decomposition (DOMDEC) module divides the simulation domain into subdomains assigned to processors, facilitating load balancing and communication minimization, which has demonstrated effective utilization on hundreds of CPU cores for systems like protein complexes. This approach allows CHARMM to scale to thousands of cores in environments, though optimal performance requires careful partitioning to handle varying computational loads from bonded and nonbonded interactions. GPU acceleration in CHARMM was introduced with interfaces in version c41 (2016), primarily targeting nonbonded calculations such as and van der Waals forces through the DOMDEC-GPU module. This integration offloads intensive computations to graphics processing units, yielding performance gains of up to 10 times compared to multi-core CPU runs for suitable benchmarks, such as simulations. Additionally, interfaces to external libraries like OpenMM provide further GPU support for broader force field compatibility. For distributed and volunteer computing, CHARMM integrates with the BOINC platform to enable fault-tolerant job distribution across volunteer resources, allowing simulations to resume from checkpoints after interruptions. This was utilized in the Docking@Home project (2008–2012), where CHARMM performed protein-ligand docking calculations on global volunteer networks. Similarly, the Clean Energy Project on World Community Grid employed CHARMM in its first phase to screen organic molecules for solar cell applications, leveraging BOINC's decentralized architecture for massive parallel screening. Recent developments in have enhanced hybrid CPU-GPU workflows, incorporating specialized kernels and adaptor APIs for improved interoperability and speed in heterogeneous environments. CHARMM also demonstrates compatibility with cloud platforms like AWS and Cloud, where users can deploy parallel jobs on virtual clusters for scalable simulations without dedicated hardware. However, polarizable simulations, such as those using the oscillator model, incur a 3–4-fold computational overhead compared to additive force fields, potentially amplifying parallel inefficiencies due to increased inter-processor communications. Best practices for load balancing include activating DOMDEC with appropriate cutoff radii and monitoring domain sizes to minimize migration overhead during runs.

Applications and Extensions

Research and Scientific Applications

CHARMM has been instrumental in advancing the understanding of and dynamics since its early applications. The program's first major demonstration in this area came with the 1977 simulation of bovine pancreatic trypsin inhibitor (BPTI), which revealed atomic fluctuations and dynamic behavior in a folded over a 9.2-picosecond trajectory, marking a pioneering effort in all-atom (MD). This foundational work laid the groundwork for exploring protein conformational changes, with subsequent CHARMM simulations extending to modern studies of aggregation, such as those of Aβ42 oligomers, where CHARMM36m parameters captured membrane-inserted β-sheet edge structures critical to pore formation in models. Furthermore, CHARMM simulations have elucidated allosteric mechanisms in proteins, for instance, by modeling the modulation of dynamics, highlighting how ligand binding propagates conformational changes across distant sites using CHARMM36 force fields. In , CHARMM facilitates through its CGenFF parameterization, enabling accurate modeling of small-molecule ligands for protein targets. For example, CGenFF has been applied in high-throughput docking and to evaluate ligand-protein interactions in cancer-related targets, prioritizing candidates with favorable binding poses and energetics. During the , CHARMM-driven (FEP) calculations assessed binding affinities of inhibitors to the main protease (Mpro), quantifying mutational impacts on nanobody affinity with ΔΔG values around -2 to +1 kcal/mol, aiding the design of resilient therapeutics. CHARMM simulations have significantly contributed to membrane biophysics, particularly in modeling rafts and function. Using the CHARMM36 force field, studies have probed formation in ternary mixtures of sphingomyelin, , and phospholipids, revealing and domain stability through order parameters like tailgroup alignment (S_CH2 ≈ 0.3-0.5). In research, CHARMM36m has simulated gating in the TRAAK , demonstrating how interactions influence conductance and voltage-dependent opening, with root-mean-square fluctuations indicating flexible selectivity filter dynamics during activation. The program's impact is underscored by its role in Nobel-recognized work, such as Martin Karplus's 1980s simulations of dynamics, which used CHARMM to map multiple conformational states on the protein's energy landscape, revealing subnanosecond fluctuations that validated experimental neutron scattering data and advanced paradigms. Recent applications (2023-2025) highlight CHARMM's ongoing relevance in cutting-edge research, including investigations into biomolecular complexes and stability under environmental stress. CHARMM-GUI serves as a prominent web-based for CHARMM, facilitating the construction of complex biomolecular systems and the generation of simulation inputs since its in 2006. It streamlines tasks such as protein , builder for bilayers, and preparation of inputs for multiple simulation engines including NAMD, , , OpenMM, and others. Version 3.8, released in July 2022, introduced enhancements like the Multicomponent Assembler, which automates the assembly of diverse molecular components such as multiple s combined with sheet-like or polymers under . A 2024 extension of this tool further supports modeling of intricate multicomponent systems, enabling efficient setup for advanced simulations. CHARMM integrates with other molecular dynamics packages through dedicated conversion tools, allowing users to leverage CHARMM force fields in alternative environments. The charmm2gmx utility automates the porting of CHARMM additive force fields to , ensuring compatibility for topology and parameter files while validating energy conservation. Similarly, CHARMM-GUI's FF-Converter module supports force fields by generating compatible inputs and converting CHARMM topologies to AMBER formats, accommodating proteins, lipids, and glycans. Python interfaces enhance accessibility, with pyCHARMM providing an embedding framework that exposes CHARMM's core functionalities—such as calculations and dynamics—for scripting in Python, including compatibility with extensions. MDAnalysis, a Python library for trajectory analysis, natively supports reading and writing CHARMM formats, enabling seamless post-simulation processing. Extensions like QwikMD and HTMD expand CHARMM's workflow capabilities for specialized applications. QwikMD, a plugin integrated with VMD, offers a user-friendly interface for plugin-based preparation, execution, and of CHARMM simulations, targeting both novices and experts in biomolecular studies. HTMD provides a Python-based platform for high-throughput , incorporating CHARMM force fields for automated system building, production, and Markov state model to accelerate molecular discovery. Community-driven tools further support CHARMM's ecosystem, particularly for force field parameterization and large-scale applications. ParamChem offers an online service for applying the CGenFF program, automating atom typing, parameter assignment, and charge derivation for small molecules to refine CHARMM topologies. The Clean Energy Project, a initiative, utilizes CHARMM's for of organic materials in solar cells and , processing millions of candidates to identify promising photovoltaic properties. As of 2025, integrations with in workflows involving CHARMM-GUI have supported advanced biomolecular simulations, such as in the design of through diffusion models followed by MD validation.

References

  1. https://pmc.ncbi.nlm.nih.gov/articles/PMC2810661/
  2. https://doi.org/10.1073/pnas.80.21.6571
  3. https://academiccharmm.org/
  4. https://pubs.acs.org/doi/10.1021/acs.jpcb.4c04100
  5. https://charmm-gui.org/?doc=license
  6. https://www.charmm-gui.org/
  7. https://pmc.ncbi.nlm.nih.gov/articles/PMC4334745/
  8. https://academiccharmm.org/news/free-charmm
  9. https://academiccharmm.org/program
  10. https://charmm.chemistry.harvard.edu/charmm_lite.php
  11. https://academiccharmm.org/documentation/installation
  12. https://forums-academiccharmm.org/viewtopic.php?t=9590
  13. https://academiccharmm.org/documentation
  14. https://academiccharmm.org/program/versions
  15. https://forums-academiccharmm.org/
  16. https://www.nobelprize.org/prizes/chemistry/2013/karplus/biographical/
  17. https://pubs.acs.org/doi/10.1021/ar020082r
  18. https://www.nobelprize.org/prizes/chemistry/2013/summary/
  19. https://cen.acs.org/people/obituaries/Theoretical-chemist-Martin-Karplus-dies/103/web/2025/01
  20. https://www.sciencedirect.com/science/article/pii/S0304416514002736
  21. https://onlinelibrary.wiley.com/doi/10.1002/jcc.20065
  22. https://onlinelibrary.wiley.com/doi/10.1002/jcc.21367
  23. https://chemrxiv.org/engage/chemrxiv/article-details/6782d5fc6dde43c90825183f
  24. https://pubs.acs.org/doi/10.1021/acs.jctc.4c01234
  25. https://pubs.acs.org/doi/10.1021/acs.chemrev.5b00505
  26. https://pubmed.ncbi.nlm.nih.gov/15224394/
  27. https://pubmed.ncbi.nlm.nih.gov/14634989/
  28. https://pubmed.ncbi.nlm.nih.gov/24459460/
  29. https://pubs.acs.org/doi/10.1021/acs.jctc.0c00057
  30. https://pmc.ncbi.nlm.nih.gov/articles/PMC10404126/
  31. https://pmc.ncbi.nlm.nih.gov/articles/PMC3752912/
  32. https://pubs.acs.org/doi/10.1021/jp510560k
  33. https://pmc.ncbi.nlm.nih.gov/articles/PMC5975207/
  34. https://annualreviews.org/doi/10.1146/annurev-biophys-070317-033349
  35. https://pmc.ncbi.nlm.nih.gov/articles/PMC2888302/
  36. https://pubs.acs.org/doi/10.1021/acs.jpcb.4c01314
  37. https://pmc.ncbi.nlm.nih.gov/articles/PMC6003505/
  38. https://pmc.ncbi.nlm.nih.gov/articles/PMC11275216/
  39. https://pmc.ncbi.nlm.nih.gov/articles/PMC2922408/
  40. https://mackerell.umaryland.edu/charmm_ff_params.html
  41. https://www.charmm-gui.org/charmmdoc/dynamc.html
  42. https://www.charmm-gui.org/charmmdoc/minimiz.html
  43. https://academiccharmm.org/documentation/version/c48b2/minimiz
  44. https://www.charmm-gui.org/charmmdoc/pert.html
  45. https://academiccharmm.org/documentation/version/c48b2/umbrel
  46. https://academiccharmm.org/documentation/version/c44b1/qmmm
  47. https://www.charmm-gui.org/charmmdoc/ewald.html
  48. https://www.charmm-gui.org/charmmdoc/gbim.html
  49. https://www.charmm-gui.org/charmmdoc/usage.html
  50. https://people.chem.umass.edu/jchenlab/teaching/mdlectures/charmm-analysis.pdf
  51. https://doi.org/10.1021/acs.jctc.3c00364
  52. https://academiccharmm.org/documentation/version/c50b1/
  53. https://academiccharmm.org/documentation/version/c50b1/install
  54. https://academiccharmm.org/documentation/version/c43b1/install
  55. https://academiccharmm.org/documentation/version/c50b1/containers
  56. https://pmc.ncbi.nlm.nih.gov/articles/PMC3966901/
  57. https://academiccharmm.org/documentation/version/c41b2/domdec
  58. https://www.researchgate.net/publication/220951864_FENZI_GPU-enabled_molecular_dynamics_simulations_of_large_membrane_regions_based_on_the_CHARMM_force_field_and_PME
  59. https://academiccharmm.org/documentation/version/c45b2/openmm
  60. http://docking.cis.udel.edu/about/project/news.html
  61. https://www.worldcommunitygrid.org/research/cep2/overview.s
  62. https://hpc.nih.gov/apps/charmm/
  63. https://pmc.ncbi.nlm.nih.gov/articles/PMC8741736/
  64. https://pmc.ncbi.nlm.nih.gov/articles/PMC10033322/
  65. https://www.nature.com/articles/s41598-019-41980-x
  66. https://www.nature.com/articles/s41598-022-22992-6
  67. https://pubs.acs.org/doi/10.1021/acs.jcim.3c01933
  68. https://pubs.rsc.org/en/content/articlehtml/2022/cp/d2cp00431c
  69. https://pubs.acs.org/doi/10.1021/acs.jcim.0c01179
  70. https://www.science.org/doi/10.1126/science.3798113
  71. https://charmm-gui.org/?doc=log&version=3.8
  72. https://www.nature.com/articles/s41467-024-49700-4
  73. https://pubs.acs.org/doi/10.1021/acs.jcim.3c00860
  74. https://pubs.aip.org/aip/jcp/article/153/3/035103/1062771/CHARMM-GUI-supports-the-Amber-force-fields
  75. https://pmc.ncbi.nlm.nih.gov/articles/PMC10504603/
  76. https://docs.mdanalysis.org/
  77. https://www.nature.com/articles/srep26536
  78. https://pubs.acs.org/doi/10.1021/ci300363c
  79. https://worldcommunitygrid.org/research/cep1/overview.do
  80. https://www.science.org/doi/10.1126/sciadv.adp7171
Add your contribution
Related Hubs
User Avatar
No comments yet.