Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Protein dynamics”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Asymmetric unwrapping of nucleosomal DNA propagates asymmetric opening and dissociation of the histone core

The nucleosome core particle (NCP) is the basic structural unit for genome packaging in eukaryotic cells and consists of DNA wound around a core of eight histone proteins. DNA access is modulated through dynamic processes of NCP disassembly. Partly disassembled structures, such as the hexasome (containing six histones) and the tetrasome (four histones), are important for transcription regulation in vivo. However, the pathways for their formation have been difficult to characterize. We combine time-resolved (TR) small-angle X-ray scattering and TR-FRET to correlate changes in the DNA conformations with composition of the histone core during salt-induced disassembly of canonical NCPs. We find that H2A–H2B histone dimers are released sequentially, with the first dimer being released after the DNA has formed an asymmetrically unwrapped, teardrop-shape DNA structure. This finding suggests that the octasome-to-hexasome transition is guided by the asymmetric unwrapping of the DNA. Furthermore, the link between DNA structure and histone composition suggests a potential mechanism for the action of proteins that alter nucleosome configurations such as histone chaperones and chromatin remodeling complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Sequence-specific backbone resonance assignments and microsecond timescale molecular dynamics simulation of human eosinophil-derived neurotoxin

Eight active canonical members of the pancreatic-like ribonuclease A (RNase A) superfamily have been identified in human. All structural homologs share similar RNA-degrading functions, while also cumulating other various biological activities in different tissues. Additionally, the functional homologs eosinophil-derived neurotoxin (EDN, or RNase 2) and eosinophil cationic protein (ECP, or RNase 3) are known to be expressed and secreted by eosinophils in response to infection, and have thus been postulated to play an important role in host defense and inflammatory response. We recently initiated the biophysical and dynamical investigation of several vertebrate RNase homologs and observed that clustering residue dynamics appear to be linked with the phylogeny and biological specificity of several members. Here we report the 1 H, 13 C and 15 N backbone resonance assignments of human EDN (RNase 2) and its molecular dynamics simulation on the microsecond timescale, providing means to pursue this comparative atomic-scale functional and dynamical analysis by NMR and computation over multiple time frames.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Kinematic Flexibility Analysis: Hydrogen Bonding Patterns Impart a Spatial Hierarchy of Protein Motion

Elastic network models (ENMs) and constraint-based, topological rigidity analysis are two distinct, coarse-grained approaches to study conformational flexibility of macromolecules. In the two decades since their introduction, both have contributed significantly to insights into protein molecular mechanisms and function. However, despite a shared purpose of these approaches, the topological nature of rigidity analysis, and thereby the absence of motion modes, has impeded a direct comparison. We present an alternative, kinematic approach to rigidity analysis, which circumvents these drawbacks. We introduce a novel protein hydrogen bond network spectral decomposition, which provides an orthonormal basis for collective motions modulated by noncovalent interactions, analogous to the eigenspectrum of normal modes. The zero modes decompose proteins into rigid clusters identical to those from topological rigidity, while nonzero modes rank protein motions by their hydrogen bond collective energy penalty. Our kinematic flexibility analysis bridges topological rigidity theory and ENM, enabling a detailed analysis of motion modes obtained from both approaches. Analysis of a large, structurally diverse data set revealed that collectivity of protein motions, reported by the Shannon entropy, is significantly reduced for rigidity theory compared to normal mode approaches. Strikingly, kinematic flexibility analysis suggests that the hydrogen bonding network encodes a protein-fold specific, spatial hierarchy of motions, which goes nearly undetected in ENM. This hierarchy reveals distinct motion regimes that rationalize experimental and simulated protein stiffness variations. Kinematic motion modes highly correlate with reported crystallographic B factors and molecular dynamics simulations of adenylate kinase. A formal expression for changes in free energy derived from the spectral decomposition indicates that motions across nearly 40% of modes obey enthalpy–entropy compensation. Taken together, our results suggest that hydrogen bond networks have evolved to modulate protein structure and dynamics, which can be efficiently probed by kinematic flexibility analysis.

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

The mutualism effector MiSSP7 of Laccaria bicolor alters the interactions between the poplar JAZ6 protein and its associated proteins

Despite the pivotal role of jasmonic acid in the outcome of plant-microorganism interactions, JA-signaling components in roots of perennial trees like western balsam poplar (Populus trichocarpa) are poorly characterized. Here we decipher the poplar-root JA-perception complex centered on PtJAZ6, a co-repressor of JA-signaling targeted by the effector protein MiSSP7 from the ectomycorrhizal basidiomycete Laccaria bicolor during symbiotic development. Through protein–protein interaction studies in yeast we determined the poplar root proteins interacting with PtJAZ6. Moreover, we assessed via yeast triple-hybrid how the mutualistic effector MiSSP7 reshapes the association between PtJAZ6 and its partner proteins. In the absence of the symbiotic effector, PtJAZ6 interacts with the transcription factors PtMYC2s and PtJAM1.1. In addition, PtJAZ6 interacts with it-self and with other Populus JAZ proteins. Finally, MiSSP7 strengthens the binding of PtJAZ6 to PtMYC2.1 and antagonizes PtJAZ6 homo-/heterodimerization. We conclude that a symbiotic effector secreted by a mutualistic fungus may promote the symbiotic interaction through altered dynamics of a JA-signaling-associated protein–protein interaction network, maintaining the repression of PtMYC2.1-regulated genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Structure and dynamics of a constitutively active neurotensin receptor

Many G protein-coupled receptors show constitutive activity, resulting in the production of a second messenger in the absence of an agonist; and naturally occurring constitutively active mutations in receptors have been implicated in diseases. To gain insight into mechanistic aspects of constitutive activity, we report here the 3.3 Å crystal structure of a constitutively active, agonist-bound neurotensin receptor (NTSR1) and molecular dynamics simulations of agonist-occupied and ligand-free receptor. Comparison with the structure of a NTSR1 variant that has little constitutive activity reveals uncoupling of the ligand-binding domain from conserved connector residues, that effect conformational changes during GPCR activation. Furthermore, molecular dynamics simulations show strong contacts between connector residue side chains and increased flexibility at the intracellular receptor face as features that coincide with robust signalling in cells. In conclusion, the loss of correlation between the binding pocket and conserved connector residues, combined with altered receptor dynamics, possibly explains the reduced neurotensin efficacy in the constitutively active NTSR1 and a facilitated initial engagement with G protein in the absence of agonist.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: An accurate correction scheme for electrostatic finite-size effects

The calculation of a protein-ligand binding free energy based on molecular dynamics (MD) simulations generally relies on a thermodynamic cycle in which the ligand is alchemically inserted into the system, both in the solvated protein and free in solution. The corresponding ligand-insertion free energies are typically calculated in nanoscale computational boxes simulated under periodic boundary conditions and considering electrostatic interactions defined by a periodic lattice-sum. This is distinct from the ideal bulk situation of a system of macroscopic size simulated under non-periodic boundary conditions with Coulombic electrostatic interactions. This discrepancy results in finite-size effects, which affect primarily the charging component of the insertion free energy, are dependent on the box size, and can be large when the ligand bears a net charge, especially if the protein is charged as well. This article investigates finite-size effects on calculated charging free energies using as a test case the binding of the ligand 2-amino-5-methylthiazole (net charge +1 e) to a mutant form of yeast cytochrome c peroxidase in water. Considering different charge isoforms of the protein (net charges −5, 0, +3, or +9 e), either in the absence or the presence of neutralizing counter-ions, and sizes of the cubic computational box (edges ranging from 7.42 to 11.02 nm), the potentially large magnitude of finite-size effects on the raw charging free energies (up to 17.1 kJ mol{sup −1}) is demonstrated. Two correction schemes are then proposed to eliminate these effects, a numerical and an analytical one. Both schemes are based on a continuum-electrostatics analysis and require performing Poisson-Boltzmann (PB) calculations on the protein-ligand system. While the numerical scheme requires PB calculations under both non-periodic and periodic boundary conditions, the latter at the box size considered in the MD simulations, the analytical scheme only requires three non-periodic PB calculations for a given system, its dependence on the box size being analytical. The latter scheme also provides insight into the physical origin of the finite-size effects. These two schemes also encompass a correction for discrete solvent effects that persists even in the limit of infinite box sizes. Application of either scheme essentially eliminates the size dependence of the corrected charging free energies (maximal deviation of 1.5 kJ mol{sup −1}). Because it is simple to apply, the analytical correction scheme offers a general solution to the problem of finite-size effects in free-energy calculations involving charged solutes, as encountered in calculations concerning, e.g., protein-ligand binding, biomolecular association, residue mutation, pK{sub a} and redox potential estimation, substrate transformation, solvation, and solvent-solvent partitioning.

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

Elevated O-GlcNAcylation stabilizes FOXM1 by its reduced degradation through GSK-3β inactivation in a human gastric carcinoma cell line, MKN45 cells

Highlights: • FOXM1 protein level and the cell proliferation were upregulated in a human gastric cancer cell line, MKN45 cells. • The ubiquitination and degradation of FOXM1 protein were both suppressed by high glucose and TMG treatment. • However, the O-GlcNAcylation was not detected on FOXM1 but not on GSK-3β. • The degradation of FOXM1 protein were suppressed by the treatment of GSK-3β inhibitors. • The higher risk of gastric cancer in diabetic patients could be due to O-GlcNAcylation-mediated FOXM1 stabilization. O-GlcNAcylation is a dynamic post-translational modification of cytonuclear proteins for intracellular signaling. Elevated O-GlcNAcylation is a general feature of cancer and contributes to cancer progression, and recent studies indicate the contribution to increasing incidence of various types of cancer in diabetic patients. However, the role of O-GlcNAcylation in tumor progression is not fully elucidated. Forkhead box M1 (FOXM1), a master mitotic transcription factor, has been implicated in all major hallmarks of cancer, and is wildly expressed in solid tumors. Given that FOXM1 expression was reported to be elevated in gastric cancer, we examined the effect of high glucose or an inhibitor of O-GlcNAc hydrolase, Thiamet G (TMG), on FOXM1 protein expression in a human gastric cancer cell line, MKN45 cells, and confirmed that FOXM1 protein level and the cell proliferation were upregulated. To investigate the molecular mechanisms by which FOXM1 protein expression is regulated by O-GlcNAcylation, the effect of high glucose and TMG on FOXM1 ubiquitination was examined in MKN45 cells. As a result, the ubiquitination and degradation of FOXM1 protein were both suppressed by high glucose and TMG treatment. However, the O-GlcNAcylation was not detected on FOXM1 but not on GSK-3β. High glucose and TMG treatment increased phospho-serine 9 GSK-3β, an inactive form, and the degradation of FOXM1 protein was suppressed by treatment of GSK-3β inhibitors in MKN45 cells. Taken together, we suggest that high glucose and elevated O-GlcNAcylation stabilize FOXM1 protein by its reduced degradation via GSK-3β inactivation in MKN45 cells, suggesting that the higher risk of gastric cancer in diabetic patients could be partially due to O-GlcNAcylation-mediated FOXM1 stabilization.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Conformational heterogeneity in closed and open states of the KcsA potassium channel in lipid bicelles

The process of ion channel gating—opening and closing—involves local and global structural changes in the channel in response to external stimuli. Conformational changes depend on the energetic landscape that underlies the transition between closed and open states, which plays a key role in ion channel gating. For the prokaryotic, pH-gated potassium channel KcsA, closed and open states have been extensively studied using structural and functional methods, but the dynamics within each of these functional states as well as the transition between them is not as well understood. In this study, we used solution nuclear magnetic resonance (NMR) spectroscopy to investigate the conformational transitions within specific functional states of KcsA. We incorporated KcsA channels into lipid bicelles and stabilized them into a closed state by using either phosphatidylcholine lipids, known to favor the closed channel, or mutations designed to trap the channel shut by disulfide cross-linking. A distinct state, consistent with an open channel, was uncovered by the addition of cardiolipin lipids. Using selective amino acid labeling at locations within the channel that are known to move during gating, we observed at least two different slowly interconverting conformational states for both closed and open channels. The pH dependence of these conformations and the predictable disruptions to this dependence observed in mutant channels with altered pH sensing highlight the importance of conformational heterogeneity for KcsA gating.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Mechanism of the association between Na + binding and conformations at the intracellular gate in neurotransmitter:sodium symporters

Neurotransmitter:sodium symporters (NSSs) terminate neurotransmission by Na + -dependent reuptake of released neurotransmitters. Previous studies suggested that Na + -binding reconfigures dynamically coupled structural elements in an allosteric interaction network (AIN) responsible for function-related conformational changes, but the intramolecular pathway of this mechanism has remained uncharted. Here we describe a new approach for the modeling and analysis of intramolecular dynamics in the bacterial NSS homolog LeuT. From microsecond-scale molecular dynamics simulations and cognate experimental verifications in both LeuT and human dopamine transporter (hDAT), we apply the novel method to identify the composition and the dynamic properties of their conserved AIN. In LeuT, two different perturbations disrupting Na+ binding and transport ( i.e. replacing Na + with Li + or the Y268A mutation at the intracellular gate) affect the AIN in strikingly similar ways. In contrast, other mutations that affect the intracellular gate (i.e. R5A and D369A) do not significantly impair Na + cooperativity and transport. Our analysis shows these perturbations to have much lesser effects on the AIN, underscoring the sensitivity of this novel method to the mechanistic nature of the perturbation. Notably, this set of observations holds as well for hDAT, where the aligned Y335A, R60A, and D436A mutations also produce different impacts on Na + dependence. Furthermore, the detailed AIN generated from our method is shown to connect Na + binding with global conformational changes that are critical for the transport mechanism. Lastly, that the AIN between the Na + binding sites and the intracellular gate in bacterial LeuT resembles that in eukaryotic hDAT highlights the conservation of allosteric pathways underlying NSS function.

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

Phosphoryl transfer reaction snapshots in crystals: Insights into the mechanism of protein kinase a catalytic subunit

To study the catalytic mechanism of phosphorylation catalyzed by cAMP-dependent protein kinase (PKA) a structure of the enzyme-substrate complex representing the Michaelis complex is of specific interest as it can shed light on the structure of the transition state. However, all previous crystal structures of the Michaelis complex mimics of the PKA catalytic subunit (PKAc) were obtained with either peptide inhibitors or ATP analogs. Here we utilized Ca 2+ ions and sulfur in place of the nucleophilic oxygen in a 20-residue pseudo-substrate peptide (CP20) and ATP to produce a close mimic of the Michaelis complex. In the ternary reactant complex, the thiol group of Cys-21 of the peptide is facing Asp-166 and the sulfur atom is positioned for an in-line phosphoryl transfer. Replacement of Ca 2+ cations with Mg 2+ ions resulted in a complex with trapped products of ATP hydrolysis: phosphate ion and ADP. As a result, the present structural results in combination with the previously reported structures of the transition state mimic and phosphorylated product complexes complete the snapshots of the phosphoryl transfer reaction by PKAc, providing us with the most thorough picture of the catalytic mechanism to date.

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

Linearly scaling and almost Hamiltonian dielectric continuum molecular dynamics simulations through fast multipole expansions

Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved—up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages.

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

Effects of Pressure and Temperature on the Atomic Fluctuations of Dihydrofolate Reductase from a Psychropiezophile and a Mesophile

Determining the effects of extreme conditions on proteins from “extremophilic” and mesophilic microbes is important for understanding how life adapts to living at extremes as well as how extreme conditions can be used for sterilization and food preservation. Previous molecular dynamics simulations of dihydrofolate reductase (DHFR) from a psychropiezophile (cold- and pressure-loving), Moritella profunda (Mp), and a mesophile, Escherichia coli (Ec), at various pressures and temperatures indicate that atomic fluctuations, which are important for enzyme function, increase with both temperature and pressure. Here, the factors that cause increases in atomic fluctuations in the simulations are examined. The fluctuations increase with temperature not only because of greater thermal energy and thermal expansion of the protein but also because hydrogen bonds between protein atoms are weakened. However, the increase in fluctuations with pressure cannot be due to thermal energy, which remains constant, nor the compressive effects of pressure, but instead, the hydrogen bonds are also weakened. In addition, increased temperature causes larger increases in fluctuations of the loop regions of MpDHFR than EcDHFR, and increased pressure causes both increases and decreases in fluctuations of the loops, which differ between the two.

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

General trends of dihedral conformational transitions in a globular protein

In this paper, dihedral conformational transitions are analyzed systematically in a model globular protein, cytochrome P450cam, to examine their structural and chemical dependences through combined conventional molecular dynamics (cMD), accelerated molecular dynamics (aMD) and adaptive biasing force (ABF) simulations. The aMD simulations are performed at two acceleration levels, using dihedral and dual boost, respectively. In comparison with cMD, aMD samples protein dihedral transitions approximately two times faster on average using dihedral boost, and ~3.5 times faster using dual boost. In the protein backbone, significantly higher dihedral transition rates are observed in the bend, coil, and turn flexible regions, followed by the β bridge and β sheet, and then the helices. Moreover, protein side chains of greater length exhibit higher transition rates on average in the aMD-enhanced sampling. Side chains of the same length (particularly Nχ = 2) exhibit decreasing transition rates with residues when going from hydrophobic to polar, then charged and aromatic chemical types. The reduction of dihedral transition rates is found to be correlated with increasing energy barriers as identified through ABF free energy calculations. In conclusion, these general trends of dihedral conformational transitions provide important insights into the hierarchical dynamics and complex free energy landscapes of functional proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Polymorphisms in Fibronectin Binding Protein A of Staphylococcus Aureus are Associated with Infection of Cardiovascular Devices

Medical implants, like cardiovascular devices, improve the quality of life for countless individuals but may become infected with bacteria like Staphylococcus aureus. Such infections take the form of a bio-film, a structured community of bacterial cells adherent to the surface of a solid substrate. Every bio-film begins with an attractive force or bond between bacterium and substratum. We used atomic force microscopy to probe experimentally forces between a fibronectin-coated surface (i.e., proxy for an implanted cardiac device) and fibronectin-binding receptors on the surface of individual living bacteria from each of 80 clinical isolates of S. aureus. These isolates originated from humans with infected cardiac devices (CDI; n = 26), uninfected cardiac devices (n = 20), and the anterior nares of asymptomatic subjects (n = 34). CDI isolates exhibited a distinct bindingforce signature and had speci!c single amino acid polymorphisms in fibronectin-binding protein A corresponding to E652D, H782Q, and K786N. In silico molecular dynamics simulations demonstrate that residues D652, Q782, and N786 in fibronectin-binding protein A form extra hydrogen bonds with fibronectin, complementing the higher binding force and energy measured by atomic force microscopy for the CDI isolates. This study is significant, because it links pathogenic bacteria biofilms from the length scale of bonds acting across a nanometer-scale space to the clinical presentation of disease at the human dimension.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2016

Effect of ethanol concentrations on temperature driven structural changes of chymotrypsin inhibitor 2

A series of atomistic molecular dynamics (MD) simulations of a small enzymatic protein Chymotrypsin Inhibitor 2 (CI2) in water-ethanol mixed solutions were carried out to explore the underlying mechanism of ethanol driven conformational changes of the protein. Efforts have been made to probe the influence of ethanol concentrations ranging from 0% to 75% (v/v) at ambient condition (300 K (T1)) and at elevated temperatures (375 K (T2) and 450 K (T3)) to investigate the temperature induced conformational changes of the protein further. Our study showed that the effect of varying ethanol concentrations on protein’s structure is almost insignificant at T1 and T2 temperatures whereas at T3 temperature, partial unfolding of CI2 in 10% ethanol solution followed by full unfolding of the protein at ethanol concentrations above 25% occurs. However, interestingly, at T3 temperature CI2’s native structure was found to be retained in pure water (0% ethanol solution) indicating that the cosolvent ethanol do play an important role in thermal denaturation of CI2. Such observations were quantified in the light of root-mean-square deviations (RMSDs) and radius of gyration. Although higher RMSD values of β-sheet over α-helix indicate complete destruction of the β-structure of CI2 at high ethanol concentrations, the associated time scale showed that the faster melting of α-helix happens over β-sheet. Around 60%-80% of initial native contacts of the protein were found broken with the separation of hydrophobic core consisting eleven residues at ethanol concentrations greater than 25%. This leads protein to expand with the increase in solvent accessible surface area. The interactions between protein and solvent molecules showed that protein’s solvation shell preferred to accommodate ethanol molecules as compared to water thereby excluded water molecules from CI2’s surface. Further, concentration dependent differential self-aggregation behavior of ethanol is likely to regulate the replacement of relatively fast diffused water by low diffused ethanol molecules from protein’s surface during the unfolding process.

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

Generative $β$-hairpin design using a residue-based physicochemical property landscape

De novo peptide design is a new frontier that has broad application potential in the biological and biomedical fields. Most existing models for de novo peptide design are largely based on sequence homology that can be restricted based on evolutionarily derived protein sequences and lack the physicochemical context essential in protein folding. Generative machine learning for de novo peptide design is a promising way to synthesize theoretical data that are based on, but unique from, the observable universe. In this study, we created and tested a custom peptide generative adversarial network intended to design peptide sequences that can fold into the -hairpin secondary structure. This deep neural network model is designed to establish a preliminary foundation of the generative approach based on physicochemical and conformational properties of 20 canonical amino acids, for example, hydrophobicity and residue volume, using extant structure-specific sequence data from the PDB. The beta generative adversarial network model robustly distinguishes secondary structures of hairpin from α helix and intrinsically disordered peptides with an accuracy of up to 96% and generates artificial -hairpin peptide sequences with minimum sequence identities around 31% and 50% when compared against the current NCBI PDB and nonredundant databases, respectively. These results highlight the potential of generative models specifically anchored by physicochemical and conformational property features of amino acids to expand the sequence-to-structure landscape of proteins beyond evolutionary limits.

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

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

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

Synthetic and systems biology for microbial production of commodity chemicals

The combination of synthetic and systems biology is a powerful framework to study fundamental questions in biology and produce chemicals of immediate practical application such as biofuels, polymers, or therapeutics. However, we cannot yet engineer biological systems as easily and precisely as we engineer physical systems. In this review, we describe the path from the choice of target molecule to scaling production up to commercial volumes. We present and explain some of the current challenges and gaps in our knowledge that must be overcome in order to bring our bioengineering capabilities to the level of other engineering disciplines. Challenges start at molecule selection, where a difficult balance between economic potential and biological feasibility must be struck. Pathway design and construction have recently been revolutionized by next-generation sequencing and exponentially improving DNA synthesis capabilities. Although pathway optimization can be significantly aided by enzyme expression characterization through proteomics, choosing optimal relative protein expression levels for maximum production is still the subject of heuristic, non-systematic approaches. Toxic metabolic intermediates and proteins can significantly affect production, and dynamic pathway regulation emerges as a powerful but yet immature tool to prevent it. Host engineering arises as a much needed complement to pathway engineering for high bioproduct yields; and systems biology approaches such as stoichiometric modeling or growth coupling strategies are required. A final, and often underestimated, challenge is the successful scale up of processes to commercial volumes. Sustained efforts in improving reproducibility and predictability are needed for further development of bioengineering.

59 BASIC BIOLOGICAL SCIENCES↗