Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “backbone”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microstructure and heteroatom dictate the doping mechanism and thermoelectric properties of poly(alkyl-chalcogenophenes)

Heteroatom substitution can favorably alter electronic transport in conductive polymers to improve their thermoelectric performance. Here we report the spectroscopic, structural, and thermoelectric properties of poly(3–(3',7'-dimethyloctyl) chalcogenophenes) or P3RX doped with 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ), where the heteroatom [X = thiophene (T), selenophene (Se), tellurophene (Te)], the doping methodology, and extent of doping are systematically varied. Spectroscopic measurements reveal that while all P3RX polymers are appreciably doped, the doping mechanism is inherently different. Poly(3-hexylthiophene) (P3HT, used as a control) and poly(3–(3',7'-dimethyloctyl)tellurophene) (P3RTe) are doped primarily via integer charge transfer (ICT), whereas poly(3–(3',7'-dimethyloctyl)selenophene) (P3RSe) and poly(3–(3',7'-dimethyloctyl)thiophene) (P3RT) are doped via charge transfer complex (CTC) mechanisms. Despite these differences, all polymers saturate with roughly the same number of F4TCNQ counterions (1 dopant per 4 to 6 heterocycles), reinforcing the idea that the extent of charge transfer varies with the doping mechanism. Grazing incidence wide-angle x-ray scattering measurements provide insight into the structural driving forces behind different doping mechanisms—P3RT and P3RSe have similar microstructures in which F4TCNQ intercalates between the π-stacked backbones resulting in CTC doping (localized carriers), while P3HT and P3RTe have microstructures in which F4TCNQ intercalates in the alkyl side chain region, giving rise to ICT doping (delocalized carriers). These structural and spectroscopic observations shed light on why P3HT and P3RTe obtain electrical conductivities ca. 3 S/cm, while P3RT and P3RSe have conductivities <10–3 S/cm under the same thin film processing conditions. Ultimately, this work quantifies the effects of heteroatom, microstructural ordering, extent of doping, and doping mechanism, thereby providing rational guidance for designing future thermoelectric polymer-dopant systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-throughput screening of tribological properties of monolayer films using molecular dynamics and machine learning

Monolayer films have shown promise as a lubricating layer to reduce friction and wear of mechanical devices with separations on the nanoscale. These films have a vast design space with many tunable properties that can affect their tribological effectiveness. For example, terminal group chemistry, film composition, and backbone chemistry can all lead to films with significantly different tribological properties. This design space, however, is very difficult to explore without a combinatorial approach and an automatable, reproducible, and extensible workflow to screen for promising candidate films. Here, using the Molecular Simulation Design Framework (MoSDeF), a combinatorial screening study was performed to explore 9747 unique monolayer films (116 964 total simulations) and a machine learning (ML) model using a random forest regressor, an ensemble learning technique, to explore the role of terminal group chemistry and its effect on tribological effectiveness. The most promising films were found to contain small terminal groups such as cyano and ethylene. The ML model was subsequently applied to screen terminal group candidates identified from the ChEMBL small molecule library. Approximately 193 131 unique film candidates were screened with approximately a five order of magnitude speed-up in analysis compared to simulation alone. The ML model was thus able to be used as a predictive tool to greatly speed up the initial screening of promising candidate films for future simulation studies, suggesting that computational screening in combination with ML can greatly increase the throughput in combinatorial approaches to generate in silico data and then train ML models in a controlled, self-consistent fashion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Network stability of FKM rubber compounds under hydrogen pressure

Fluorocarbon-based rubber (FKM) compounds are typically utilized in the aeronautical industry due to its exceptional chemical and thermal resistance and recently has drawn significant interest from the hydrogen community for various hydrogen applications, but there exists little data evaluating its hydrogen compatibility. In this paper, a formulated FKM compound was exposed to ultrahigh hydrogen pressure (90 MPa) and analyzed by attenuated total reflectance Fourier-transform infrared (ATR-FTIR) spectroscopy and crosslinking densities by solvent swelling measurement. ATR-FTIR spectroscopic data demonstrated that high-pressure hydrogen not only caused cleavages in polymer backbone and cross linkages, but also induced additional degradation of FKM compounds by breaking intermolecular hydrogen bonding between carboxylic acids, which were produced during rubber molding and post curing processes. This hydrogen pressure-driven degradation was significantly higher towards the core of the material than at the surface due to its weaker bonding between polymers and fillers.

Polymers, hydrogen, stability, degradation, Charac↗
DOE OSTI2023

Sterically attenuated electronic communication in cobalt complexes of meridional isoquinoline-derived ligands for applications in electrocatalysis

The ability to synthetically tune the ligand frameworks of redox-active molecules is of critical importance to the economy of solar fuels because manipulating their redox properties can afford control over the operating potentials of sustained electrocatalytic or photoelectrocatalytic processes. The electronic and steric properties of 2,2′:6′,2″-terpyridine (Terpy) ligand frameworks can be tuned by functional group substitution on ligand backbones, and these correlate strongly to their Hammett parameters. The synthesis of a new series of tridentate meridional ligands of 2,4,6-trisubstituted pyridines that engineers the ability to finely tune the redox potentials of cobalt complexes to more positive potentials than that of their Terpy analogs is achieved by aryl-functionalizing at the four-position and by including isoquinoline at the two- and six-positions of pyridine (Aryl-DiQ). Their cobalt complex syntheses, their electronic properties, and their catalytic activity for carbon dioxide (CO2) reduction are reported and compared to their Terpy analogs. The cobalt derivatives generally experience a positive shift in their redox features relative to the Terpy-based analogs, covering a complementary potential range. Although those evaluated fail to produce any quantifiable products for the reduction of CO2 and suffer from long-term instability, these results suggest possible alternate strategies for stabilizing these compounds during catalysis. We speculate that lower equilibrium association constants to the cobalt center are intrinsic to these ligands, which originate from a steric interaction between protons on the pyridine and isoquinoline moieties. Nevertheless, the new Aryl-DiQ ligand framework has been engineered to selectively tune homoleptic cobalt complexes’ redox potentials.

Chemistry↗
DOE OSTI2024

Bond length alternation of π -conjugated polymers predicted by the Fermi–Löwdin orbital self-interaction correction method

π-conjugated polymers have been used in a wide range of practical applications, partly due to their unique properties that originate in the delocalization of electrons through the polymer backbone. The level of delocalization can be characterized by the induced bond length alternation (BLA), with shorter BLA connected with strong delocalization and vice versa. The accurate description of this structural parameter can be considered a benchmark for testing the capability of different electronic structure methods for self-interaction error (SIE) removal and electron correlation inclusion. Density functional theory (DFT), in its local or semi-local flavors, suffers from SIE and, thus, underestimates the BLA compared to self-interaction-free methods. In this work, we utilize the Fermi–Löwdin orbital self-interaction correction (FLOSIC) method for one-electron self-interaction removal to characterize the BLA of five oligomers with increasing length extrapolated to the polymeric limit. We compare the self-interaction-free BLA to several DFT approximations, Møller–Plesset second-order perturbation theory (MP2), and the BLA obtained with the domain based local pair natural orbital CCSD(T) [DLPNO-CCSD(T)] approximation. Our findings show that FLOSIC corrects for the small BLA given by (semi-)local DFT approximations, but it tends to overcorrect with respect to CAM-B3LYP, MP2, and DLPNO-CCSD(T).

Chemistry↗
DOE OSTI2024

Spectral broadening and vibronic dynamics of the S2 state of canthaxanthin in the orange carotenoid protein

We have performed a series of broadband multidimensional electronic spectroscopy experiments to probe the electronic and vibrational dynamics of the canthaxanthin chromophore of the Orange Carotenoid Protein (OCP) from Synechocystis sp. PCC 6803 in its photoactivated red state, OCPR. Cross-peaks observed below the diagonal of the two-dimensional electronic spectrum indicate that absorption transitions prepare the bright S2 state of the ketocarotenoid canthaxanthin near to a sequence of conical intersections, allowing passage to the dark S1 state via the Sx intermediate in <50 fs. Rapid damping of excited-state coherent wavepacket motions suggests that the branching coordinates of the conical intersections include out-of-plane deformation and C=C stretching coordinates of the π-conjugated isoprenoid backbone. The unusual proximity of the Franck–Condon S2 state structure to the conical intersections with Sx and S1 suggests that the protein surroundings of canthaxanthin prepare it to function as an excitation energy trap in the OCPR–phycobilisome complex. Numerical simulations using the multimode Brownian oscillator model demonstrate that the ground-state absorption spectrum of OCPR overlaps with the fluorescence emission spectrum of allophycocyanin due to spectral broadening derived especially from the intramolecular motions of the canthaxanthin chromophore in its binding site.

Yang, Kunwei↗
DOE OSTI2024

A modular and extensible CHARMM-compatible model for all-atom simulation of polypeptoids

Peptoids (N-substituted glycines) are a class of sequence-defined synthetic peptidomimetic polymers with applications including drug delivery, catalysis, and biomimicry. Classical molecular simulations have been used to predict and understand the conformational dynamics of single chains and their self-assembly into morphologies including sheets, tubes, spheres, and fibrils. The CGenFF-NTOID model based on the CHARMM General Force Field has demonstrated success in accurate all-atom molecular modeling of peptoid structure and thermodynamics. Extension of this force field to new peptoid side chains has historically required reparameterization of side chain bonded interactions against ab initio data. This fitting protocol improves the accuracy of the force field but is also burdensome and precludes modular extensibility of the model to arbitrary peptoid sequences. In this work, we develop and demonstrate a Modular Side Chain CGenFF-NTOID (MoSiC-CGenFF-NTOID) as an extension of CGenFF-NTOID employing a modular decomposition of the peptoid backbone and side chain parameterizations, wherein arbitrary side chains within the large family of substituted methyl groups (i.e., –CH 3 , –CH 2 R, –CHRR', and –CRR'R") are directly ported from CGenFF. We validate this approach against ab initio calculations and experimental data to develop a MoSiC-CGenFF-NTOID model for all 20 natural amino acid side chains along with 13 commonly used synthetic side chains and present an extensible paradigm to efficiently determine whether a novel side chain can be directly incorporated into the model or whether refitting of the CGenFF parameters is warranted. We make the model freely available to the community along with a tool to perform automated initial structure generation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Key Factors in Semi-Generic Coarse-Grained Modeling of Solid Polymer Electrolytes

Given the long length and time scales of interest and strong interactions present in solid polymer electrolytes, coarse-grained molecular models can be useful in understanding the molecular basis of their structure and transport properties. Rather than representing atoms individually, coarse-grained models use beads representing groups of atoms, making the system simpler and more efficient. Generic coarse-grained models, that are not built based on matching a particular atomistic system, can provide insight into the important physical considerations that apply across different chemical systems, though it can be unclear how to best match to experimental systems and capture relevant experimental behaviors with as simple of a model as possible. Here, we build on generic bead-spring models, but use stiff angle potentials, set different bead properties for different polymer types, and include additional ion parameters, creating a semi-generic coarse-grained model that can map more specifically to polymers and copolymers with different chain architectures, component glass transition temperatures, and ion solvation behaviors. We set most parameters with basic homopolymer data such as glass transition temperatures, Kuhn length, density, and dielectric constant, with further adjustment based on data from the polymer electrolyte system. We specifically model polystyrene-block-poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) with lithium triflate salt, considering the POEM to be made of a poly(methylmethacrylate) backbone with poly(ethylene oxide) side chains. Solvation of ions is accounted for by additional polymer-ion interactions of the form -S/r4, plus additional lithium-polymer Lennard-Jones interactions. We find these potentials have different effects and discuss strategies for setting these parameters.

Zhang, Yuanhao↗
DOE OSTI2020

Modular repeat protein sculpting using rigid helical junctions

The ability to precisely design large proteins with diverse shapes would enable applications ranging from the design of protein binders that wrap around their target to the positioning of multiple functional sites in specified orientations. We describe a protein backbone design method for generating a wide range of rigid fusions between helix-containing proteins and use it to design 75,000 structurally unique junctions between monomeric and homo-oligomeric de novo designed and ankyrin repeat proteins (RPs). Of the junction designs that were experimentally characterized, 82% have circular dichroism and solution small-angle X-ray scattering profiles consistent with the design models and are stable at 95 °C. Crystal structures of four designed junctions were in close agreement with the design models with rmsds ranging from 0.9 to 1.6 Å. Electron microscopic images of extended tetrameric structures and ~10-nm-diameter “L” and “V” shapes generated using the junctions are close to the design models, demonstrating the control the rigid junctions provide for protein shape sculpting over multiple nanometer length scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

DNA origami protection and molecular interfacing through engineered sequence-defined peptoids

DNA nanotechnology has established approaches for designing programmable and precisely controlled nanoscale architectures through specific Watson−Crick base-pairing, molecular plasticity, and intermolecular connectivity. In particular, superior control over DNA origami structures could be beneficial for biomedical applications, including biosensing, in vivo imaging, and drug and gene delivery. However, protecting DNA origami structures in complex biological fluids while preserving their structural characteristics remains a major challenge for enabling these applications. Here, we developed a class of structurally well-defined peptoids to protect DNA origamis in ionic and bioactive conditions and systematically explored the effects of peptoid architecture and sequence dependency on DNA origami stability. The applicability of this approach for drug delivery, bioimaging, and cell targeting was also demonstrated. A series of peptoids (PE1–9) with two types of architectures, termed as “brush” and “block,” were built from positively charged monomers and neutral oligo-ethyleneoxy monomers, where certain designs were found to greatly enhance the stability of DNA origami. Through experimental and molecular dynamics studies, we demonstrated the role of sequence-dependent electrostatic interactions of peptoids with the DNA backbone. We showed that octahedral DNA origamis coated with peptoid (PE2) can be used as carriers for anticancer drug and protein, where the peptoid modulated the rate of drug release and prolonged protein stability against proteolytic hydrolysis. Finally, we synthesized two alkyne-modified peptoids (PE8 and PE9), conjugated with fluorophore and antibody, to make stable DNA origamis with imaging and cell-targeting capabilities. Our results demonstrate an approach toward functional and physiologically stable DNA origami for biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The myosin interacting-heads motif present in live tarantula muscle explains tetanic and posttetanic phosphorylation mechanisms

Striated muscle contraction involves sliding of actin thin filaments along myosin thick filaments, controlled by calcium through thin filament activation. In relaxed muscle, the two heads of myosin interact with each other on the filament surface to form the interacting-heads motif (IHM). A key question is how both heads are released from the surface to approach actin and produce force. Here we used time-resolved synchrotron X-ray diffraction to study tarantula muscle before and after tetani. The patterns showed that the IHM is present in live relaxed muscle. Tetanic contraction produced only a very small backbone elongation, implying that mechanosensing—proposed in vertebrate muscle—is not of primary importance in tarantula. Rather, thick filament activation results from increases in myosin phosphorylation that release a fraction of heads to produce force, with the remainder staying in the ordered IHM configuration. After the tetanus, the released heads slowly recover toward the resting, helically ordered state. During this time the released heads remain close to actin and can quickly rebind, enhancing the force produced by posttetanic twitches, structurally explaining posttetanic potentiation. Taken together, these results suggest that, in addition to stretch activation in insects, two other mechanisms for thick filament activation have evolved to disrupt the interactions that establish the relaxed helices of IHMs: one in invertebrates, by either regulatory light-chain phosphorylation (as in arthropods) or Ca 2+ -binding (in mollusks, lacking phosphorylation), and another in vertebrates, by mechanosensing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Disruption of a key ligand-H-bond network drives dissociative properties in vamorolone for Duchenne muscular dystrophy treatment

Duchenne muscular dystrophy is a genetic disorder that shows chronic and progressive damage to skeletal and cardiac muscle leading to premature death. Antiinflammatory corticosteroids targeting the glucocorticoid receptor (GR) are the current standard of care but drive adverse side effects such as deleterious bone loss. Through subtle modification to a steroidal backbone, a recently developed drug, vamorolone, appears to preserve beneficial efficacy but with significantly reduced side effects. We use combined structural, biophysical, and biochemical approaches to show that loss of a receptor-ligand hydrogen bond drives these remarkable therapeutic effects. Moreover, vamorolone uniformly weakens coactivator associations but not corepressor associations, implicating partial agonism as the main driver of its dissociative properties. Additionally, we identify a critical and evolutionarily conserved intramolecular network connecting the ligand to the coregulator binding surface. Interruption of this allosteric network by vamorolone selectively reduces GR-driven transactivation while leaving transrepression intact. Our results establish a mechanistic understanding of how vamorolone reduces side effects, guiding the future design of partial agonists as selective GR modulators with an improved therapeutic index.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structural basis for differential recognition of phosphohistidine-containing peptides by 1-pHis and 3-pHis monoclonal antibodies

In 2015, monoclonal antibodies (mAbs) that selectively recognize the 1-pHis or 3-pHis isoforms of phosphohistidine were developed by immunizing rabbits with degenerate Ala/Gly peptides containing the nonhydrolyzable phosphohistidine (pHis) analog- phosphotriazolylalanine (pTza). Here, we report structures of five rabbit mAbs bound to cognate pTza peptides: SC1-1 and SC50-3 that recognize 1-pHis, and their 3-pHis–specific counterparts, SC39-4, SC44-8, and SC56-2. These cocrystal structures provide insights into the binding modes of the pTza phosphate group that are distinct for the 1- and 3-pHis mAbs with the selectivity arising from specific contacts with the phosphate group and triazolyl ring. The mode of phosphate recognition in the 3-pHis mAbs recapitulates the Walker A motif, as present in kinases. The complementarity-determining regions (CDRs) of four of the Fabs interact with the peptide backbone rather than peptide side chains, thus conferring sequence independence, whereas SC44-8 shows a proclivity for binding a GpHAGA motif mediated by a sterically complementary CDRL3 loop. Specific hydrogen bonding with the triazolyl ring precludes recognition of pTyr and other phosphoamino acids by these mAbs. Kinetic binding experiments reveal that the affinity of pHis mAbs for pHis and pTza peptides is submicromolar. Bound pHis mAbs also shield the pHis peptides from rapid dephosphorylation. The epitope–paratope interactions illustrate how these anti-pHis antibodies are useful for a wide range of research techniques and this structural information can be utilized to improve the specificity and affinity of these antibodies toward a variety of pHis substrates to understand the role of histidine phosphorylation in healthy and diseased states.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Computationally designed peptide macrocycle inhibitors of New Delhi metallo-β-lactamase 1

The rise of antibiotic resistance calls for new therapeutics targeting resistance factors such as the New Delhi metallo-β-lactamase 1 (NDM-1), a bacterial enzyme that degrades β-lactam antibiotics. We present structure-guided computational methods for designing peptide macrocycles built from mixtures of l - and d -amino acids that are able to bind to and inhibit targets of therapeutic interest. Our methods explicitly consider the propensity of a peptide to favor a binding-competent conformation, which we found to predict rank order of experimentally observed IC 50 values across seven designed NDM-1- inhibiting peptides. We were able to determine X-ray crystal structures of three of the designed inhibitors in complex with NDM-1, and in all three the conformation of the peptide is very close to the computationally designed model. In two of the three structures, the binding mode with NDM-1 is also very similar to the design model, while in the third, we observed an alternative binding mode likely arising from internal symmetry in the shape of the design combined with flexibility of the target. Although challenges remain in robustly predicting target backbone changes, binding mode, and the effects of mutations on binding affinity, our methods for designing ordered, binding-competent macrocycles should have broad applicability to a wide range of therapeutic targets.

42 ENGINEERING↗
DOE OSTI2021

Transferrin receptor targeting by de novo sheet extension

The de novo design of polar protein–protein interactions is challenging because of the thermodynamic cost of stripping water away from the polar groups. Here, we describe a general approach for designing proteins which complement exposed polar backbone groups at the edge of beta sheets with geometrically matched beta strands. We used this approach to computationally design small proteins that bind to an exposed beta sheet on the human transferrin receptor (hTfR), which shuttles interacting proteins across the blood–brain barrier (BBB), opening up avenues for drug delivery into the brain. We describe a design which binds hTfR with a 20 nM K d , is hyperstable, and crosses an in vitro microfluidic organ-on-a-chip model of the human BBB. Our design approach provides a general strategy for creating binders to protein targets with exposed surface beta edge strands.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural OFF/ON transitions of myosin in relaxed porcine myocardium predict calcium-activated force

Contraction in striated muscle is initiated by calcium binding to troponin complexes, but it is now understood that dynamic transition of myosin between resting, ordered OFF states on thick filaments and active, disordered ON states that can bind to thin filaments is critical in regulating muscle contractility. These structural OFF to ON transitions of myosin are widely assumed to correspond to transitions from the biochemically defined, energy-sparing, super-relaxed (SRX) state to the higher ATPase disordered-relaxed (DRX) state. Here, in this study, we examined the effect of 2’-deoxy-ATP (dATP), a naturally occurring energy substrate for myosin, on the structural OFF to ON transitions of myosin motors in porcine cardiac muscle thick filaments. Small-angle X-ray diffraction revealed that titrating dATP in relaxation solutions progressively moves the myosin heads from ordered OFF states on the thick filament backbone to disordered ON states closer to thin filaments. Importantly, we found that the structural OFF to ON transitions are not equivalent to the biochemically defined SRX to DRX transitions and that the dATP-induced structural OFF to ON transitions of myosin motors in relaxed muscle are strongly correlated with submaximal force augmentation by dATP. These results indicate that structural OFF to ON transitions of myosin in relaxed muscle can predict the level of force attained in calcium-activated cardiac muscle. Computational modeling and stiffness measurements suggest a final step in the OFF to ON transition may involve a subset of DRX myosins that form weakly bound cross-bridges prior to becoming active force-producing cross-bridges.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

De novo design of small beta barrel proteins

Small beta barrel proteins are attractive targets for computational design because of their considerable functional diversity despite their very small size (<70 amino acids). However, there are considerable challenges to designing such structures, and there has been little success thus far. Because of the small size, the hydrophobic core stabilizing the fold is necessarily very small, and the conformational strain of barrel closure can oppose folding; also intermolecular aggregation through free beta strand edges can compete with proper monomer folding. Here, we explore the de novo design of small beta barrel topologies using both Rosetta energy–based methods and deep learning approaches to design four small beta barrel folds: Src homology 3 (SH3) and oligonucleotide/oligosaccharide-binding (OB) topologies found in nature and five and six up-and-down-stranded barrels rarely if ever seen in nature. Both approaches yielded successful designs with high thermal stability and experimentally determined structures with less than 2.4 Å rmsd from the designed models. Using deep learning for backbone generation and Rosetta for sequence design yielded higher design success rates and increased structural diversity than Rosetta alone. The ability to design a large and structurally diverse set of small beta barrel proteins greatly increases the protein shape space available for designing binders to protein targets of interest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Exact reaction coordinates for flap opening in HIV-1 protease

The primary goal of protein science is to understand how proteins function, which requires understanding the functional dynamics responsible for transitions between different functional structures of a protein. A central concept is the exact reaction coordinates that can determine the value of committor for any protein configuration, which provide the optimal description of functional dynamics. Despite intensive efforts, identifying the exact reaction coordinates (RCs) in complex molecules remains a formidable challenge. Using the recently developed generalized work functional, we report the discovery of the exact RCs for an important functional process—the flap opening of HIV-1 protease. Our results show that this process has six RCs, each one is a linear combination of ~240 backbone dihedrals, providing the precise definition of collectivity and cooperativity in the functional dynamics of a protein. Applying bias potentials along each RC can accelerate flap opening by 10 3 to 10 4 folds. The success in identifying the RCs of a protein with 198 residues represents a significant progress beyond that of the alanine dipeptide, currently the only other complex molecule for which the exact RCs for its conformational changes are known. Our results suggest that the generalized work functional (GWF) might be the fundamental operator of mechanics that controls protein dynamics.

Science & Technology - Other Topics↗