Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular conformation”

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 163 records · Page 9

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↗

Continuous in vitro evolution of catalytic function

A population of RNA molecules that catalyze the template-directed ligation of RNA substrates was made to evolve in a continuous manner in the test tube. A simple serial transfer procedure was used to achieve approximately 300 successive rounds of catalysis and selective amplification in 52 hours. During this time, the population size was maintained against an overall dilution of 3 x 10(298). Both the catalytic rate and amplification rate of the RNAs improved substantially as a consequence of mutations that accumulated during the evolution process. Continuous in vitro evolution makes it possible to maintain laboratory "cultures" of catalytic molecules that can be perpetuated indefinitely.

Non-NASA Center↗

Dynamics and activation of membrane-bound B cell receptor assembly

B-cell receptor (BCR) complexes are expressed on the surface of a B-cell and are critical in antigen recognition and modulating the adaptive immune response. Even though the relevance of antibodies has been known for almost a hundred years, the antigen-dependent activation mechanism of B-cells has remained elusive. Several models have been proposed for BCR activation, including cross-linking, conformation-induced oligomerization, and dissociation activation models. Recently, the first cryo-EM structures of the human B-cell antigen receptor of the IgM and IgG isotypes have been published that validates the asymmetric organization of the BCR complex. Here, we carry out extensive molecular dynamics simulations to probe the conformational changes upon antigen binding and the influence of the membrane lipids. We identify two critical dynamical events that could be associated with antigen-dependent activation of BCR. First, antigen binding causes increased flexibility in regions distal to the antigen binding site. Second, antigen binding alters the rearrangement of IgM transmembrane helices, including the relative interaction of Igα/Igβ that mediates intracellular signaling. Furthermore, these transmembrane rearrangements lead to changes in localized lipid composition. Our work indirectly supports the conformational-change induced models of BCR activation and contributes to the understanding of the antigen-dependent activation mechanism of BCRs.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling study on the cleavage step of the self-splicing reaction in group I introns

A three-dimensional model of the Tetrahymena thermophila group I intron is used to further explore the catalytic mechanism of the transphosphorylation reaction of the cleavage step. Based on the coordinates of the catalytic core model proposed by Michel and Westhof (Michel, F., Westhof, E. J. Mol. Biol. 216, 585-610 (1990)), we first converted their ligation step model into a model of the cleavage step by the substitution of several bases and the removal of helix P9. Next, an attempt to place a trigonal bipyramidal transition state model in the active site revealed that this modified model for the cleavage step could not accommodate the transition state due to insufficient space. A lowering of P1 helix relative to surrounding helices provided the additional space required. Simultaneously, it provided a better starting geometry to model the molecular contacts proposed by Pyle et al. (Pyle, A. M., Murphy, F. L., Cech, T. R. Nature 358, 123-128. (1992)), based on mutational studies involving the J8/7 segment. Two hydrated Mg2+ complexes were placed in the active site of the ribozyme model, using the crystal structure of the functionally similar Klenow fragment (Beese, L.S., Steitz, T.A. EMBO J. 10, 25-33 (1991)) as a guide. The presence of two metal ions in the active site of the intron differs from previous models, which incorporate one metal ion in the catalytic site to fulfill the postulated roles of Mg2+ in catalysis. The reaction profile is simulated based on a trigonal bipyramidal transition state, and the role of the hydrated Mg2+ complexes in catalysis is further explored using molecular orbital calculations.

NASA Discipline Exobiology↗

Computer modeling of polymers

A Polymer Molecular Analysis Display System (p-MADS) was developed for computer modeling of polymers. This method of modeling allows for the theoretical calculation of molecular properties such as equilibrium geometries, conformational energies, heats of formations, crystal packing arrangements, and other properties. Furthermore, p-MADS has the following capabilities: constructing molecules from internal coordinates (bonds length, angles, and dihedral angles), Cartesian coordinates (such as X-ray structures), or from stick drawings; manipulating molecules using graphics and making hard copy representation of the molecules on a graphics printer; and performing geometry optimization calculations on molecules using the methods of molecular mechanics or molecular orbital theory.

Green, Terry J.↗

Molecular deformation and stress-strain behavior of poly(bisphenol-A-diphenyl sulfone)

The strain-birefringence response of poly(bisphenol-A-diphenyl sulfone) is found to be independent of temperature at temperatures below -100 C; at higher temperatures, the response becomes slightly dependent on temperature, with lower birefringence seen at higher temperatures. The stress-strain behavior and the stress-birefringence response both depend on temperature over the entire testing temperature range (-179 C to 150 C) studied; this dependence, however, is not pronounced. The evidence is seen as suggesting that the polymer molecules respond to deformation by undergoing conformational rearrangements; the mode of the molecular deformation remains unchanged for temperatures of -100 C or lower. At higher temperatures, the average length of the chain segments involved in the rearrangement increases. The stress-strain response is attributed mainly to chain orientation. The entropic contribution deriving from chain orientation at temperatures below -100 C is still substantial. The modest temperature dependence of the stress-strain response suggests that the energy barriers for the chain segments involved in the rearrangement are relatively low.

Hong, S.-D.↗

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

3DRISM↗

Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor

Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.

59 BASIC BIOLOGICAL SCIENCES↗

Operando Investigation of the Molecular Origins of Dipole Switching in P(VDF‐TrFE‐CFE) Terpolymer for Large Adiabatic Temperature Change

Abstract Relaxor ferroelectric polymers exhibiting a giant electrocaloric effect (ECE) can potentially be used to create next‐generation solid‐state coolers. Under an electric field, poly(vinylidene fluoride‐trifluoroethylene‐chlorofluoroethylene) terpolymer goes through a large dipolar entropy change producing a high adiabatic temperature change (Δ T ECE ). This work resolves the molecular origins of the large entropy change behind the electric field‐induced dipole switching. A Fourier transform infrared spectroscopy equipped with a high voltage source is used to operandoly observe the characteristic molecular vibrational modes. A short‐range trans (T) conformation of the CF 2 ‐CH 2 dyads interrupted by a gauche (G) conformation, e.g., TTTG in the terpolymer chain, undergoes a dynamic transformation that leads to a corresponding Δ T ECE whenever an electric field is applied. The molecular dynamics simulation also proves that the energy barrier that the transformation from TTTGs into a long T sequence overcomes is smaller than that for all other conformations. A mixed solvent system is used to obtain T3G‐enriched terpolymer films exhibiting a 4.02 K Δ T ECE at 60 MV m −1 and these films are employed to manufacture a 2‐layer‐cascaded cooling device that achieves a 6.7 K temperature lift, the highest reported value for a 2‐layer cascaded device made of fluoropolymers.

Zhu, Yuan↗

Multi-scale computer-aided design and photo-controlled macromolecular synthesis boosting uranium harvesting from seawater

By integrating multi-scale computational simulation with photo-regulated macromolecular synthesis, this study presents a new paradigm for smart design while customizing polymeric adsorbents for uranium harvesting from seawater. A dissipative particle dynamics (DPD) approach, combined with a molecular dynamics (MD) study, is performed to simulate the conformational dynamics and adsorption process of a model uranium grabber, i.e., PAO m -b-PPEGMA n , suggesting that the maximum adsorption capacity with atomic economy can be achieved with a preferred block ratio of 0.18. The designed polymers are synthesized using the PET-RAFT polymerization in a microfluidic platform, exhibiting a record high adsorption capacity of uranium (11.4 ± 1.2 mg/g) in real seawater within 28 days. This study offers an integrated perspective to quantitatively assess adsorption phenomena of polymers, bridging metal-ligand interactions at the molecular level with their spatial conformations at the mesoscopic level. The established protocol is generally adaptable for target-oriented development of more advanced polymers for broadened applications.

97 MATHEMATICS AND COMPUTING↗

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural mechanisms for binding and activation of a contact-quenched fluorophore by RhoBAST

The fluorescent light-up aptamer RhoBAST, which binds and activates the fluorophore–quencher conjugate tetramethylrhodamine-dinitroaniline with high affinity, super high brightness, remarkable photostability, and fast exchange kinetics, exhibits excellent performance in super-resolution RNA imaging. Here we determine the co-crystal structure of RhoBAST in complex with tetramethylrhodamine-dinitroaniline to elucidate the molecular basis for ligand binding and fluorescence activation. The structure exhibits an asymmetric “A”-like architecture for RhoBAST with a semi-open binding pocket harboring the xanthene of tetramethylrhodamine at the tip, while the dinitroaniline quencher stacks over the phenyl of tetramethylrhodamine instead of being fully released. Molecular dynamics simulations show highly heterogeneous conformational ensembles with the contact-but-unstacked fluorophore–quencher conformation for both free and bound tetramethylrhodamine-dinitroaniline being predominant. The simulations also show that, upon RNA binding, the fraction of xanthene-dinitroaniline stacked conformation significantly decreases in free tetramethylrhodamine-dinitroaniline. This highlights the importance of releasing dinitroaniline from xanthene tetramethylrhodamine to unquench the RhoBAST–tetramethylrhodamine-dinitroaniline complex. Using SAXS and ITC, we characterized the magnesium dependency of the folding and binding mode of RhoBAST in solution and indicated its strong structural robustness. The structures and binding modes of relevant fluorescent light-up aptamers are compared, providing mechanistic insights for rational design and optimization of this important fluorescent light-up aptamer-ligand system.

59 BASIC BIOLOGICAL SCIENCES↗

A mini-TGA protein modulates gene expression through heterogeneous association with transcription factors

Abstract TGA (TGACG-binding) transcription factors, which bind their target DNA through a conserved basic region leucine zipper (bZIP) domain, are vital regulators of gene expression in salicylic acid (SA)-mediated plant immunity. Here, we investigated the role of StTGA2.1, a potato (Solanum tuberosum) TGA lacking the full bZIP, which we named a mini-TGA. Such truncated proteins have been widely assigned as loss-of-function mutants. We, however, confirmed that StTGA2.1 overexpression compensates for SA-deficiency, indicating a distinct mechanism of action compared with model plant species. To understand the underlying mechanisms, we showed that StTGA2.1 can physically interact with StTGA2.2 and StTGA2.3, while its interaction with DNA was not detected. We investigated the changes in transcriptional regulation due to StTGA2.1 overexpression, identifying direct and indirect target genes. Using in planta transactivation assays, we confirmed that StTGA2.1 interacts with StTGA2.3 to activate StPRX07, a member of class III peroxidases (StPRX), which are known to play role in immune response. Finally, via structural modeling and molecular dynamics simulations, we hypothesized that the compact molecular architecture of StTGA2.1 distorts DNA conformation upon heterodimer binding to enable transcriptional activation. This study demonstrates how protein truncation can lead to distinct functions and that such events should be studied carefully in other protein families.

59 BASIC BIOLOGICAL SCIENCES↗

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA’s diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

99 GENERAL AND MISCELLANEOUS↗

Engineered living materials for capture, conversion, and recycling technologies

Continued dependence on crude oil and natural gas resources for fossil fuels has caused global atmospheric carbon dioxide (CO 2 ) emissions to increase to record-setting proportions. There is an urgent need for efficient and inexpensive carbon sequestration systems to mitigate large-scale CO 2 emissions from industrial flue gas. Carbonic anhydrase (CA) has shown high potential for enhanced CO 2 capture applications compared to conventional absorption-based methods currently utilized in various industrial settings. This study aims to understand structural aspects that contribute to the stability of CA enzymes critical for their applications in industrial processes, which require the ability to withstand conditions different from their native environments. Here, we evaluated the thermostability and enzyme activity of mesophilic and thermophilic CA variants at different temperature conditions and in the presence of atmospheric gas pollutants like nitrogen oxides (NO x ) and sulphur oxides (SO x ). Based on our enzyme activity assays and molecular dynamics simulations, we see increased conformational stability and CA activity levels in thermostable CA variants incubated week-long at different temperature conditions. The thermostable CA variants also retained high levels of CA activity despite changes in solution pH due to increasing NO x and SO x concentrations. Furthermore, a loss of CA activity was observed only at high concentrations of NO x /SO x that possibly can be minimized with appropriate buffered solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sugars as the optimal biosynthetic carbon substrate of aqueous life throughout the universe

Our previous analysis of the energetics of metabolism showed that both the biosynthesis of amino acids and lipids from sugars, and the fermentation of organic substrates, were energetically driven by electron transfer reactions resulting in carbon redox disproportionation (Weber, 1997). Redox disproportionation--the spontaneous (energetically favorable) direction of carbon group transformation in biosynthesis--is brought about and driven by the energetically downhill transfer of electron pairs from more oxidized carbon groups (with lower half-cell reduction potentials) to more reduced carbon groups (with higher half-cell reduction potentials). In this report, we compare the redox and kinetic properties of carbon groups in order to evaluate the relative biosynthetic capability of organic substrates, and to identify the optimal biosubstrate. This analysis revealed that sugars (monocarbonyl alditols) are the optimal biosynthetic substrate because they contain the maximum number of biosynthetically useful high energy electrons/carbon atom while still containing a single carbonyl group needed to kinetically facilitate their conversion to useful biosynthetic intermediates. This conclusion applies to aqueous life throughout the Universe because it is based on invariant aqueous carbon chemistry--primarily, the universal reduction potentials of carbon groups.

NASA Center ARC↗