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 397 records · Page 22

Machine Learning‐Augmented Molecular Dynamics Simulations (MD) Reveal Insights Into the Disconnect Between Affinity and Activation of ZTP Riboswitch Ligands

Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

Chemistry↗

Deciphering the Cofilin Oligomers via Intermolecular Disulfide Bond Formation: A Coarse-Grained Molecular Dynamics Approach to Understanding Cofilin’s Regulation on Actin Filaments

Cofilin, a key actin-binding protein, orchestrates the dynamics of the actomyosin network through its actin-severing activity and by promoting the recycling of actin monomers. Recent experimental work suggests that cofilin also forms functionally distinct oligomers through thiol post-translational modification (PTM) that encourages actin nucleation and assembly. Despite these advances, the structural conformations of cofilin oligomers that modulate actin activity remain elusive because there are combinatorial ways to oxidize thiols in cysteines to form disulfide bonds rapidly. This study employs molecular dynamics simulations to investigate human cofilin 1 as a case study for exploring cofilin dimers via disulfide bond formation. Using the free energy profiling, our simulations unveil a range of probable cofilin dimer structures not represented in current Protein Data Bank entries. These candidate dimers are characterized by their distinct population distributions and relative free energies. Of particular note is a dimer featuring an interface between cysteines 139 and 147 residues, which demonstrates stable free energy characteristics and intriguingly symmetrical geometry. In contrast, the experimentally proposed dimer structure exhibits a less stable free energy profile. Here, we also evaluate frustration quantification based on the energy landscape theory in the protein-protein interactions at the dimer interfaces. Notably, the 39-39 dimer configuration emerges as a promising candidate for forming cofilin tetramers, as substantiated by frustration analysis. Additionally, docking simulations with actin filaments further evaluate the stability of these cofilin dimer-actin complexes. Our findings thus offer a computational framework for understanding the role of thiol post-translational modification cofilin proteins in regulating oligomerization, and the subsequent cofilin-mediated actin dynamics in the actomyosin network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular basis for higher affinity of SARS-CoV-2 spike RBD for human ACE2 receptor

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused substantially more infections, deaths, and economic disruptions than the 2002-2003 SARSCoV. The key to understanding SARS-CoV-2's higher infectivity lies partly in its host receptor recognition mechanism. Experiments show that the human angiotensin converting enzyme 2 (ACE2) protein, which serves as the primary receptor for both CoVs, binds to the receptor binding domain (RBD) of CoV-2's spike protein stronger than SARS-CoV's spike RBD. The molecular basis for this difference in binding affinity, however, remains unexplained from X-ray structures. To go beyond insights gained from X-ray structures and investigate the role of thermal fluctuations in structure, we employ all-atom molecular dynamics simulations. Microseconds-long simulations reveal that while CoV and CoV-2 spike-ACE2 interfaces have similar conformational binding modes, CoV-2 spike interacts with ACE2 via a larger combinatorics of polar contacts, and on average, makes 45% more polar contacts. Correlation analysis and thermodynamic calculations indicate that these differences in the density and dynamics of polar contacts arise from differences in spatial arrangements of interfacial residues, and dynamical coupling between interfacial and non-interfacial residues. Furthermore, these results recommend that ongoing efforts to design spike-ACE2 peptide blockers will benefit from incorporating dynamical information as well as allosteric coupling effects.

59 BASIC BIOLOGICAL SCIENCES↗

Structural basis of transcription: RNA polymerase II substrate binding and metal coordination using a free-electron laser

Catalysis and translocation of multisubunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near-atomic resolution and precise arrangement of key active site components have been elusive. Here, we present the free-electron laser (FEL) structures of a matched ATP-bound Pol II and the hyperactive Rpb1 T834P bridge helix (BH) mutant at the highest resolution to date. The radiation-damage-free FEL structures reveal the full active site interaction network, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg 2+ , and, more importantly, a putative third (site C) Mg 2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structures indicate that the third Mg 2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/BH interactions induce conformational changes that could allow translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the T834P mutant reveal rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Science & Technology - Other Topics↗

Acetylcholinesterase: Structure, dynamics, and interactions with organophosphorus compounds

Acetylcholinesterase (AChE) is an enzyme that hydrolyzes the neurotransmitter acetylcholine (ACh), removing it from the synaptic cleft after the transmission of an electrical signal, making it an essential component of chemical neurotransmission. AChE is a serine hydrolase, containing a catalytic triad of Ser/His/Glu. AChE is a prime target for pharmaceuticals treating a variety of neurological disorders. It is also the target of synthetic organophosphorus (OP) compounds that have been used as pesticides and chemical warfare agents. OP compounds contain a potent leaving group, such as fluorine, and act by forming a covalent adduct with the catalytic serine of the AChE active site. A wealth of structural information is available for AChE, including over 300 structures, including a subset of structures in complex with drugs as well as OP compounds. This review will highlight the interactions between OP compounds and AChE from a structural and computational perspective, with a discussion of access to the active site, as well as side reactions that lead to dealkylation of the OP-catalytic serine adduct, a process known as aging. We conclude that while the majority of the conformational changes needed to accommodate the OP compounds are localized to the acyl loop in the crystal structures, molecular dynamics simulations highlight the potential for a far more dynamic enzyme.

59 BASIC BIOLOGICAL SCIENCES↗

Polyampholyte physics: Liquid–liquid phase separation and biological condensates

Here we review the current understanding of molecular interactions that govern liquid-liquid phase separation (LLPS) of biological condensates. The connection between sequence, chain conformation, and phase separation of intrinsically disordered proteins (IDPs) and their model polyampholytes is discussed. In particular, we highlight how the charge pattern influences the conformation and phase behavior of natural IDPs. We then describe recent results from theoretical treatments of polyampholytes implementing random phase approximation, field-theoretic simulations, and transfer matrix theory that show an increase in charge segregation results in an increased tendency to phase separate.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of melt shearing on D -mannitol crystal twisting in the presence of small molecule and macromolecular additives

Here, many molecular crystals grow as twisted lamellae from the melt when the driving force is high, but the mechanisms governing the spontaneous formation of helicoidal crystallites are varied and those operative in particular cases are difficult to identify. Polyvinylpyrrolidone (PVP) induces the crystallization of D-mannitol as banded spherulites comprised of twisted fibrils emanating radially from spherulite nucleation centers. The molecular weight of the PVP phase affects the twisting pitch, P, (rotation of the crystallite by 180°), in D-mannitol thin films crystallized from the melt. At a relatively high crystallization temperature (T c =130°C), P was sensitive to the PVP molecular weight, ranging from 390±20 μm for D-mannitol films incorporating 15 wt% 10 kDa PVP to 20±3 μm for the same weight of 1300 kDa PVP. Magnitudes of complex viscosities, |η*|, of D-mannitol/PVP melts measured via small-amplitude oscillatory shearing were strongly dependent on the PVP molecular weight, but were not correlated to P(T). Instead, P was sensitive to the dynamic PVP chain conformation during D-mannitol crystallization. Under steady torsional shear, P decreased from~30 μm to~8μm with increasing shear rates from 0.01 to 100 s -1 for D-mannitol films crystallized at 130°C in the presence of 15 wt% 10 or 1300 kDa PVP. Shear forces decrease the entanglement density of polymer chains while orienting the chains along the stream lines of the viscometric flow, indicating that the conformations of macromolecular additives can affect the pitch of banded spherulites. By contrast, P in D-mannitol twisted by D-sorbitol was independent of shear rate.

36 MATERIALS SCIENCE↗

Machine Learning Reveals the Critical Interactions for SARS-CoV-2 Spike Protein Binding to ACE2

SARS-CoV and SARS-CoV-2 bind to the human ACE2 receptor in practically identical conformations, although several residues of the receptor-binding domain (RBD) differ between them. Herein, we have used molecular dynamics (MD) simulations, machine learning (ML), and free-energy perturbation (FEP) calculations to elucidate the differences in binding by the two viruses. Although only subtle differences were observed from the initial MD simulations of the two RBD–ACE2 complexes, ML identified the individual residues with the most distinctive ACE2 interactions, many of which have been highlighted in previous experimental studies. FEP calculations quantified the corresponding differences in binding free energies to ACE2, and examination of MD trajectories provided structural explanations for these differences. Lastly, the energetics of emerging SARS-CoV-2 mutations were studied, showing that the affinity of the RBD for ACE2 is increased by N501Y and E484K mutations but is slightly decreased by K417N.

60 APPLIED LIFE SCIENCES↗

Molecular simulation of glycerol-derived triether podands for lithium ion solvation

Solvate ionic liquids (ILs) are promising candidates for several applications due to their stability, high coulombic efficiency, and low volatility. In this work, we investigate the solvation of lithium-bistriflimide by different glycerol-derived triether solvents, using molecular dynamics simulations. Very strong interactions between Li + and the solvent oxygen sites are found, leading to significant conformational changes in the solvent. Here, by comparing the conformation of the neat solvents with their IL mixtures at different concentrations and temperatures, we find that the presence of Li + induces a distinct crown-like structure in the solvent molecules. The Li + cations and the surrounding solvent form a podand complex, which is stable even at elevated temperatures. These glycerol-derived solvents exhibit distinct interactions with Li + cations which may be exploited in electrolytic applications or lithium recovery processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gating choreography and mechanism of the human proton-activated chloride channel $\mathrm{ASOR}$

The proton-activated chloride channel ASOR (TMEM206/PAC) permeates anions across cellular membranes in response to acidification, thereby enhancing acid-induced cell death and regulating endocytosis. The molecular mechanisms of pH-dependent control are not understood, in part because structural information for an activated conformation of ASOR is lacking. Here, we reconstitute function from purified protein and present a 3.1-Å-resolution cryo–electron microscopy structure of human ASOR at acidic pH in an activated conformation. The work contextualizes a previous acidic pH structure as a desensitized conformation. Combined with electrophysiological studies and high-resolution structures of resting and desensitized states, the work reveals mechanisms of proton sensing and ion pore gating. Clusters of extracellular acidic residues function as pH sensors and coalesce when protonated. Ensuing conformational changes induce metamorphosis of transmembrane helices to fashion an ion conduction pathway unique to the activated conformation. The studies identify a new paradigm of channel gating in this ubiquitous ion channel.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular structure of r/GCG/d/TATACGC/ - A DNA-RNA hybrid helix joined to double helical DNA

The molecule r(GCG)d(TATACGC) is self-complementary and forms two DNA-RNA hybrid segments surrounding a central region of double helical DNA; its molecular structure has been solved by X-ray analysis. All three parts of the molecule adopt a conformation which is close to that seen in the 11-fold RNA double helix. The conformation of the ribonucleotides is partly determined by water molecules bridging between the ribose O2' hydroxyl group and cytosine O2. The hybrid-DNA duplex junction contains no structural discontinuities. However, the central DNA TATA sequence has some structural irregularities.

Wang, A. H.-J.↗

An algorithm for converting a virtual-bond chain into a complete polypeptide backbone chain

A systematic analysis is presented of the algorithm for converting a virtual-bond chain, defined by the coordinates of the alpha-carbons of a given protein, into a complete polypeptide backbone. An alternative algorithm, based upon the same set of geometric parameters used in the Purisima-Scheraga algorithm but with a different "linkage map" of the algorithmic procedures, is proposed. The global virtual-bond chain geometric constraints are more easily separable from the loal peptide geometric and energetic constraints derived from, for example, the Ramachandran criterion, within the framework of this approach.

NASA Discipline Exobiology↗

Lessons from an evolving rRNA: 16S and 23S rRNA structures from a comparative perspective

The 16S and 23S rRNA higher-order structures inferred from comparative analysis are now quite refined. The models presented here differ from their immediate predecessors only in minor detail. Thus, it is safe to assert that all of the standard secondary-structure elements in (prokaryotic) rRNAs have been identified, with approximately 90% of the individual base pairs in each molecule having independent comparative support, and that at least some of the tertiary interactions have been revealed. It is interesting to compare the rRNAs in this respect with tRNA, whose higher-order structure is known in detail from its crystal structure (36) (Table 2). It can be seen that rRNAs have as great a fraction of their sequence in established secondary-structure elements as does tRNA. However, the fact that the former show a much lower fraction of identified tertiary interactions and a greater fraction of unpaired nucleotides than the latter implies that many of the rRNA tertiary interactions remain to be located. (Alternatively, the ribosome might involve protein-rRNA rather than intramolecular rRNA interactions to stabilize three-dimensional structure.) Experimental studies on rRNA are consistent to a first approximation with the structures proposed here, confirming the basic assumption of comparative analysis, i.e., that bases whose compositions strictly covary are physically interacting. In the exhaustive study of Moazed et al. (45) on protection of the bases in the small-subunit rRNA against chemical modification, the vast majority of bases inferred to pair by covariation are found to be protected from chemical modification, both in isolated small-subunit rRNA and in the 30S subunit. The majority of the tertiary interactions are reflected in the chemical protection data as well (45). On the other hand, many of the bases not shown as paired in Fig. 1 are accessible to chemical attack (45). However, in this case a sizeable fraction of them are also protected against chemical modification (in the isolated rRNA), which suggests that considerable higher-order structure remains to be found (although all of it may not involve base-base interactions and so may not be detectable by comparative analysis). The agreement between the higher-order structure of the small-subunit rRNA and protection against chemical modification is not perfect, however; some bases shown to covary canonically are accessible to chemical modification (45).(ABSTRACT TRUNCATED AT 400 WORDS).

Review↗

Zuotin, a putative Z-DNA binding protein in Saccharomyces cerevisiae

A putative Z-DNA binding protein, named zuotin, was purified from a yeast nuclear extract by means of a Z-DNA binding assay using [32P]poly(dG-m5dC) and [32P]oligo(dG-Br5dC)22 in the presence of B-DNA competitor. Poly(dG-Br5dC) in the Z-form competed well for the binding of a zuotin containing fraction, but salmon sperm DNA, poly(dG-dC) and poly(dA-dT) were not effective. Negatively supercoiled plasmid pUC19 did not compete, whereas an otherwise identical plasmid pUC19(CG), which contained a (dG-dC)7 segment in the Z-form was an excellent competitor. A Southwestern blot using [32P]poly(dG-m5dC) as a probe in the presence of MgCl2 identified a protein having a molecular weight of 51 kDa. The 51 kDa zuotin was partially sequenced at the N-terminal and the gene, ZUO1, was cloned, sequenced and expressed in Escherichia coli; the expressed zuotin showed similar Z-DNA binding activity, but with lower affinity than zuotin that had been partially purified from yeast. Zuotin was deduced to have a number of potential phosphorylation sites including two CDC28 (homologous to the human and Schizosaccharomyces pombe cdc2) phosphorylation sites. The hexapeptide motif KYHPDK was found in zuotin as well as in several yeast proteins, DnaJ of E.coli, csp29 and csp32 proteins of Drosophila and the small t and large T antigens of the polyoma virus. A 60 amino acid segment of zuotin has similarity to several histone H1 sequences. Disruption of ZUO1 in yeast resulted in a slow growth phenotype.

NASA Discipline Exobiology↗

Structural comparison of allophycocyanin variants reveals the molecular basis for their spectral differences

Abstract Allophycocyanins are phycobiliproteins that absorb red light and transfer the energy to the reaction centers of oxygenic photosynthesis in cyanobacteria and red algae. Recently, it was shown that some allophycocyanins absorb far-red light and that one subset of these allophycocyanins, comprising subunits from the ApcD4 and ApcB3 subfamilies (FRL-AP), form helical nanotubes. The lowest energy absorbance maximum of the oligomeric ApcD4-ApcB3 complexes occurs at 709 nm, which is unlike allophycocyanin (AP; ApcA-ApcB) and allophycocyanin B (AP-B; ApcD-ApcB) trimers that absorb maximally at ~ 650 nm and ~ 670 nm, respectively. The molecular bases of the different spectra of AP variants are presently unclear. To address this, we structurally compared FRL-AP with AP and AP-B, performed spectroscopic analyses on FRL-AP, and leveraged computational approaches. We show that among AP variants, the α-subunit constrains pyrrole ring A of its phycocyanobilin chromophore to different extents, and the coplanarity of ring A with rings B and C sets a baseline for the absorbance maximum of the chromophore. Upon oligomerization, the α-chromophores of all AP variants exhibit a red shift of the absorbance maximum of ~ 25 to 30 nm and band narrowing. We exclude excitonic coupling in FRL-AP as the basis for this red shift and extend the results to discuss AP and AP-B. Instead, we attribute these spectral changes to a conformational alteration of pyrrole ring D, which becomes more coplanar with rings B and C upon oligomerization. This study expands the molecular understanding of light-harvesting attributes of phycobiliproteins and will aid in designing phycobiliproteins for biotechnological applications.

Gisriel, Christopher J.↗

Synthesis, aphicidal activity and conformation of novel insect kinin analogues as potential eco‐friendly insecticides

Abstract BACKGROUND The discovery of ecofriendly insecticides through a new strategy for aphid control is important because of the substantial resistance and unexpected eco‐toxicity to honeybees caused by traditional insecticides. The insect kinins, a class of multifunctional insect neuropeptides, are considered for potential application in pest control. In our previous work we developed several series of insect kinin analogues and found a promising lead II‐1 with good aphicidal activity. To seek further eco‐friendly aphicides, the optimization of II‐1 is carried out in this study. RESULTS Fifteen novel Yaa 3 modified analogues based on the lead II‐1 were synthesized. The aphicidal tests indicated that IV‐3 , IV‐5 and IV‐10 exhibited significant activity against the soybean aphid Aphis glycines with LC 50 values of 0.0029, 0.0072 and 0.0086 mmol L –1 , respectively, higher than that of lead II‐1 and the commercial Pymetrozine. The molecular modeling results showed that analogues II‐1 , IV‐3 , IV‐5 , IV‐7 and IV‐10 formed a β‐turn‐like conformation, while the conformation of analogues IV‐1 , IV‐2 and IV‐9 seemed to be linear. Some structural elements favorable for the activity were proposed based on the conformation‐activity relationship of the analogues. CONCLUSION Insect kinin analogues derived from lead II‐1 by modifying the hydrolysis site Yaa 3 with natural, sterically hindered α‐ and β‐amino acids showed great potential as eco‐friendly insecticides. Inspiringly, the most active analogue IV‐3 can be a candidate for further development. The β‐turn‐like conformation and the orientation of the aromatic rings of the side chain of Phe 2 and Trp 4 may be critical factors beneficial to activity. © 2019 Society of Chemical Industry

Zhang, Chuanliang↗

LongBondEliminator: A Molecular Simulation Tool to Remove Ring Penetrations in Biomolecular Simulation Systems

We develop a workflow, implemented as a plugin to the molecular visualization program VMD, that can fix ring penetrations with minimal user input. LongBondEliminator, detects ring piercing artifacts by the long, strained bonds that are the local minimum energy conformation during minimization for some assembled simulation system. The LongBondEliminator tool then automatically treats regions near these long bonds using multiple biases applied through NAMD. By combining biases implemented through the collective variables module, density-based forces, and alchemical techniques in NAMD, LongBondEliminator will iteratively alleviate long bonds found within molecular simulation systems. Through three concrete examples with increasing complexity, a lignin polymer, an viral capsid assembly, and a large, highly glycosylated protein aggrecan, we demonstrate the utility for this method in eliminating ring penetrations from classical MD simulation systems. The tool is available via gitlab as a VMD plugin, and has been developed to be generically useful across a variety of biomolecular simulations.

59 BASIC BIOLOGICAL SCIENCES↗