QCell: Comprehensive Quantum-Mechanical Dataset Spanning Diverse Biomolecular Fragments

Basic Information

• Paper ID: 2510.09939
• Title: QCell: Comprehensive Quantum-Mechanical Dataset Spanning Diverse Biomolecular Fragments
• Authors: Adil Kabylda, Sergio Suárez-Dou, Nils Davoine, Florian N. Brünig, Alexandre Tkatchenko
• Classification: physics.chem-ph
• Publication Date: October 11, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.09939

Abstract

Recent advances in machine learning force fields (MLFFs) are revolutionizing molecular simulation by bridging the gap between quantum mechanical accuracy and computational efficiency of classical potentials. However, the development of reliable MLFFs for biomolecular systems remains constrained by the scarcity of high-quality, chemically diverse quantum mechanical datasets that must encompass all major biomolecular categories expressed in living cells. Critically, such comprehensive datasets must be computed using non-empirical or minimally empirical Schrödinger equation solving approximations. To address these limitations, the authors introduce the QCell dataset—a curated collection of 525,000 new quantum mechanical calculations spanning biomolecular fragments of carbohydrates, nucleic acids, lipids, dimers, and ionic clusters. QCell complements existing datasets, bringing the total number of available data points to 41 million molecular systems, all computed using hybrid density functional theory with non-local many-body dispersion interactions, captured at the PBE0+MBD(-NL) quantum mechanical level.

Research Background and Motivation

Problem Definition

1. Core Issue: Existing quantum mechanical datasets primarily cover small molecules and proteins, with significant gaps in three major biomolecular categories—nucleic acids, lipids, and carbohydrates—which collectively constitute approximately 40% of cellular biomass.
2. Significance:
   • Biomolecular chemical space possesses unique characteristics, with complexity arising primarily from the conformational space of relatively limited repeated chemical building blocks
   • Accurate modeling of biomolecular interactions is essential for computational chemistry and biophysics
   • MLFFs require diverse and high-quality QM datasets to faithfully represent the chemical space encountered in biomolecular systems
3. Limitations of Existing Approaches:
   • Traditional QM methods offer high accuracy but low computational efficiency
   • Empirical force fields provide high efficiency but limited accuracy
   • Existing datasets such as GEMS, QCML, and OMol25, while representing progress, still exhibit significant gaps in the three major biomolecular categories
4. Research Motivation:
   • Fill gaps in biomolecular datasets
   • Employ a consistent non-empirical quantum mechanical theoretical level
   • Provide comprehensive training resources for next-generation MLFFs

Core Contributions

1. Constructed the QCell Dataset: Contains 525,881 new QM calculations of biomolecular fragments, spanning nucleic acids, lipids, carbohydrates, ions/water, and non-covalent dimers
2. Extended Data Coverage: Combined with existing datasets, total data points reach 41 million molecular systems, encompassing 82 chemical elements
3. Unified Theoretical Level: All calculations employ the PBE0+MBD(-NL) level, ensuring data consistency
4. Deep Conformational Sampling: Focuses on conformational diversity in biologically relevant chemical environments
5. Technical Validation: Dataset quality verified through structural analysis and machine learning force field training

Methodology Details

Dataset Construction Workflow

The QCell dataset employs a five-step workflow:

1. Building Block Library Management and Initial 3D Structure Generation
2. Extensive Conformational Sampling (molecular dynamics or dedicated conformational generation tools)
3. Representative Fragment Selection
4. DFTB+MBD Method Pre-optimization
5. High-Quality PBE0+MBD(-NL) Quantum Mechanical Calculations

Specific Methodologies for Each Molecular Category

Nucleic Acids

• Utilized Nucleic Acid Builder to construct solvated DNA heptamer double helices (A-, B-, Z-DNA forms)
• Performed molecular dynamics simulations using the OL21 force field
• Extracted central double-stranded trinucleotide fragments from heptamer trajectories
• Included DNA base pair dimers and gas-phase RNA fragments

Lipids

• Generated phospholipid membrane structures using CHARMM-GUI Membrane Builder
• Covered POPC, POPE, POPG, POPS phospholipids and cholesterol
• Conducted 500 ns production simulations using the Lipid21 force field
• Selected fatty acid monomers, dimers, and trimers based on geometric proximity

Carbohydrates

• Constructed a library of 52 common monosaccharides, including α/β anomeric configurations of pentoses and hexoses
• Built disaccharides and sugar-peptide linkages using PyMOL
• Generated conformations using the CREST program with 12 kcal/mol maximum energy threshold
• Clustered by linking dihedral angles and selected representative conformations

Ions and Water

• Prepared solvated ionic systems with ions centered in water boxes
• Applied MBpol force field for monovalent ions and AMBER force field for divalent ions
• Captured solvation effects at varying hydration levels (1-100 water molecules)

Quantum Mechanical Calculation Details

• Theoretical Level: PBE0+MBD(-NL) - non-empirical hybrid functional with many-body dispersion treatment
• Software: FHI-aims code
• Basis Set: "Tight" basis set for small molecules, "intermediate" basis set for molecules >350 atoms
• Convergence Criteria: Total energy 10^-5 eV, eigenvalue sum 10^-3 eV, charge density 10^-5 electrons/ų, forces 10^-4 eV/Å

Experimental Setup

Dataset Composition

Evaluation Metrics

• Structural geometric descriptor validation
• Mean absolute error (MAE) of forces in machine learning force fields
• Radial distribution function comparison with experimental reference values

Machine Learning Validation

Trained MLFFs using SO3LR architecture to assess dataset quality:

• Three model sizes: small, medium, large
• Combined loss function: forces, dipole moments, Hirshfeld ratios, energies (weights 100:10:10:1)
• 10 Å long-range cutoff, trained on A100 GPU for 180 hours

Experimental Results

Structural Validation Results

1. Nucleic Acids: Phosphate-phosphate distances and backbone bending angle distributions of DNA fragments reproduced expected values for A-, B-, Z-DNA
2. Lipids: Radius of gyration distributions of fatty acid fragments reasonably reflected chain extension and packing
3. Carbohydrates: N/O-glycosidic dihedral angles spanned the full conformational space, reproducing all major rotamers
4. Ions/Water: Radial distribution functions matched experimental hydration distances, with accurate peak positions for monovalent ion-oxygen and O-O

Machine Learning Performance

Force MAE results for different dataset subsets:

• Nucleic Acids: ~0.8 kcal/mol/Å (large model)
• Lipids: ~0.6 kcal/mol/Å (large model)
• Carbohydrates: ~0.5 kcal/mol/Å (large model)
• Ions/Water: ~0.7 kcal/mol/Å (large model)
• DES370k: ~0.8 kcal/mol/Å (large model)

Errors systematically decreased with model capacity, with most subsets achieving below 1 kcal/mol/Å, demonstrating dataset internal consistency and the generalization capability of modern MLFFs across chemically diverse systems.

Related Work

Existing QM Datasets

• QM7-X: Small organic molecules, 4.19 million data points
• MD22: Molecular dynamics trajectories
• GEMS: Hierarchical protein fragmentation strategy
• SPICE: Drug-like molecules and peptides
• QCML: Systematic mapping of small molecule chemical space
• OMol25: Chemically heterogeneous ensemble

Advantages of This Work

• First systematic coverage of three major biomolecular categories: nucleic acids, lipids, and carbohydrates
• Unified non-empirical theoretical level ensures data consistency
• Deep conformational sampling focuses on biologically relevant chemical environments
• Perfect compatibility with existing datasets for unified training

Conclusions and Discussion

Main Conclusions

1. The QCell dataset successfully fills important gaps in biomolecular QM data
2. The unified PBE0+MBD(-NL) theoretical level ensures compatibility with existing datasets
3. Structural validation confirms the chemical reasonableness and diversity of the dataset
4. Machine learning validation demonstrates excellent predictive performance

Limitations

1. Radial distribution functions for divalent ions show slight deviations from experimental values
2. Fragment size limited to 402 atoms
3. Primarily focuses on biologically relevant elements with relatively limited elemental diversity
4. Balance between gas-phase and solution-phase environments requires further optimization

Future Directions

1. Extension to larger biomolecular fragments
2. Inclusion of additional solvent effects and environmental conditions
3. Further validation and calibration against experimental data
4. Development of new MLFF architectures specifically designed for biomolecules

In-Depth Evaluation

Strengths

1. Fills Important Gaps: First systematic solution to data insufficiency in nucleic acids, lipids, and carbohydrates
2. Rigorous Methodology: Employs non-empirical quantum mechanical methods with solid theoretical foundations
3. High Data Quality: Multiple validations ensure structural and energetic reasonableness
4. Large Practical Value: Compatible with existing datasets, directly applicable to MLFF training
5. Open Access: Publicly available dataset promotes field development

Limitations

1. Computational Cost: PBE0+MBD(-NL) calculations are computationally expensive, limiting dataset expansion
2. Fragment Constraints: Maximum 402-atom limitation may inadequately capture long-range interactions
3. Environmental Simplification: Primarily considers gas-phase and simple solvation, insufficient for complex biological environments
4. Limited Validation: Lacks direct comparison with high-accuracy methods (e.g., CCSD(T))

Impact

1. Academic Contribution: Provides important data foundation for biomolecular MLFF development
2. Practical Value: Directly applicable to drug design, biomolecular simulation, and related fields
3. Reproducibility: Detailed methodology and open data ensure reproducibility
4. Advancement Promotion: May catalyze development of new biomolecular modeling methods

Application Scenarios

1. Biomolecular MLFF Training: Direct application to training universal force fields covering multiple biomolecular types
2. Drug Design: Provides data for protein-ligand and DNA-drug interaction modeling
3. Membrane Biology: Lipid data applicable to membrane protein and membrane interaction research
4. Glycobiology: Carbohydrate data supports glycoprotein and glycolipid research
5. Method Development: Provides benchmark testing data for new quantum chemistry methods and MLFF architectures

References

This paper cites 58 important references, encompassing key works in quantum chemistry methods, machine learning force fields, biomolecular simulation, and related datasets, providing solid theoretical foundations and technical support for the research.