Computational Modeling
NREL uses computational modeling to increase the efficiency of biomass conversion by rational design using multiscale modeling, applying theoretical approaches, and testing scientific hypotheses.
![model of enzymes wrapping on cellulose; colorful circular structures entwined through blue strands](assets/images/biomass-feedstocks_model.jpg)
Cellulosomes are complexes of protein scaffolds and enzymes that are highly effective in decomposing biomass. This is a snapshot of a coarse-grain model of complex cellulosome binding to a coarse grain model of cellulose fibrils. The complex cellulosome is constructed of seven cellulosomes (green, blue, cyan, pink, purple, majenta, and yellow) bound to a secondary scaffold (red) by cohesin-dockerin pairs. The secondary scaffold is very flexible and allows for motion of the primary cellulosomes within a range. The primary cellulosomes are made of nine enzymes each bound to a primary scaffold via cohesin-dockerin pairing with many flexible linkers. All primary scaffolds and some enzymes have cellulose binding domains for attaching to cellulose.
Featured Publications
![](assets/images/computation-modeling_heating-simulations.jpg)
![Illustration of a protein represented as a transparent solvent-accessible surface with the cartoon protein backbone inside. The hydrophobic surface patch is shown as orange on the surface.](assets/images/computation-modeling_cartoon-surface-rep.jpg)
Predicting Enzyme Adsorption to Lignin Films by Calculating Enzyme Surface Hydrophobicity, The Journal of Biological Chemistry (2014)
![Illustration of multi-modular cellobiohydrolase enzyme T. reesei Cel7A, represented as quot;surfacesquot; or quot;blobs,quot; on the surface of cellulose, represented as green and red multi-jointed strands. One of these strands is pulled out of the surface and threaded through the middle of the enzyme's catalytic domain. Shown below this are three panels showing the distinct stages of the cellulose chain's progress across the active site of the enzyme (representing the quot;processive cyclequot; of the enzyme), along with the free energy barriers and free energy differences for each step.](assets/images/carbohydrate-protein-interactions.jpg)
Carbohydrate–Protein Interactions That Drive Processive Polysaccharide Translocation in Enzymes Revealed from a Computational Study of Cellobiohydrolase Processivity, Journal of the American Chemical Society (2014)
![Illustration showing identified minimum-energy pathway (solid curves) of H reduction of the ideal NiO(100) surface. State 1, the initial state, is the ideal surface with two H2 molecules adsorbed on two Ni sites. In state 2, one of the two adsorbed H2 molecules dissociates and an H-H pair is formed. The atomic structure of state 3 is shown in the lower inset indicated by an arrow from the number 3. The dashed curve shows a direct reaction from 1 to 3 without pre-dissociation of the first H2 molecule. The structure and partial charge density (at the highest occupied band) of the transition state for the 1 to 3 reaction step are shown in the upper inset.](assets/images/initial-reduction-nio.jpg)
Initial reduction of the NiO(100) surface in hydrogen, The Journal of Chemical Physics (2013)
View all NREL bioenergy computational modeling publications .
Capabilities
![Illustration of a carbohydrate binding domain (white cartoon with colored licorice-shaped sugars decorating it) on a cellulose surface (depicted as green solvent-accessible surface](assets/images/computation-modeling_h2_capabilities.jpg)
Computer Enzyme Design
NREL's molecular modeling provides insight into enzyme mechanisms and enables higher performance and new functionality for industrial enzymes. For example, the study depicted by the illustration to the left determined the structure and binding orientation of several glycosylation patterns on the binding domain, which we showed in previous studies to greatly enhance the binding to cellulose, a property that has been shown to improve the enzyme activity in breaking down cellulose in biomass.
![Structure of cellulose nanofibrils shown as sticks, showing side view with layers easily seen (top), top view where the sugar rings are shown (middle), and the twisting behavior, untwisted bottom left, and with natural twist bottom right.](assets/images/computation-modeling_h3_comp-enzym-design.jpg)
Cell Wall Components
We use multi-scale modeling to investigate the properties of cellulose, hemicellulose, lignin, and pectin. Together they compose cell walls and are the source of biofuels and biomaterials. Our modeling investigates their properties for better conversion and material design.
![Illustration of core carbon metabolic model of Actinobacillus succinogenes that shows fluxes between central metabolic nodes. Production of succinate is limited by the redox balance of the cell, in which NADH consuming reactions (orange) must be balanced with NADH producing reactions (green).](assets/images/metabolic-modeling-capability.jpg)
Metabolic Modeling
To increase yields and redirect production to desired products, we model the existing pathways of microbial conversion of sugars to fuels and high-value products.
![Illustration showing how catalysts are designed to convert biomass to liquid fuels and chemicals. On the left is a stylized tree/forest labeled "Biomass"; partial hexagons and hexagonal lines are to the right with an arrow pointing right and labeled "Zeolites"; the arrow points to "Furans," a pentagon shape with four bars jutting out; a right arrow then points to stylized icons of "Plastic Resins," "Diesel," and "Gasoline" trucks and pumps.](assets/images/quantum-mechanical-models-capability.jpg)
Quantum Mechanical Models
NREL studies chemical and electronic properties and processes to reduce barriers in conversion and upgrading of renewable sources.
![Illustration of a computational tool for computing multi-phase adsorption on catalytic surfaces.](assets/images/computation-modeling_h3_quantum-mechanical.jpg)
Catalyst Science and Design
We use ab initio reaction energetics coupled with experimental kinetic models to predict catalyst performance under realistic reaction conditions.
![microscopic photos of poplar and pine particles and their respective 3D models](assets/images/computation-modeling_micro-photos-of-poplar.jpg)
Biomass Mesostructure and Continuum Multiphysics Modeling
We are developing realistic, species-specific particle models that predict how feedstocks and blendstocks will perform in various conversion processes.
![Illustration of a reactor simulation of a bubbling catalyst reactor: a three-dimensional cylinder displaying infrared temperatures in blue, green, yellow, orange, and red, with red and orange block-like structures forming a circle and coming out of the top of the cylinder.](assets/images/reactor-simulations-capability.jpg)
Reactor Simulations
We use multiphase (solid, liquid, and/or s + gas) flow simulations to describe, understand, and optimize reactors used in catalytic pyrolysis, enzymatic hydrolysis, and fermentation. Efforts also include coupling flow models with reaction kinetics and interphase transport. experiments.
Research Team
Principal Investigators
Peter Ciesielski
Research Scientist, Bioenergy and Biomaterials
Peter.Ciesielski@nrel.gov | 303-384-7691Lintao Bu
James Lischeske
Ambarish Nag
Hariswaran Sitaraman
Mike Sprague
Related and Integrated Programs
Biological and Catalytic Conversion of Sugars and Lignin
Biomass Deconstruction and Pretreatment
Enzyme and Microbial Development
Heterogeneous Catalysis for Thermochemical Conversion
Thermochemical Process Integration, Scale-Up, and Piloting
Collaborators
BioEnergy Science Center
C3Bio: Center for Direct Catalytic Conversion of Biomass to Biofuels
Computational Pyrolysis Consortium