Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein folding”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Foldy: An open-source web application for interactive protein structure analysis

Foldy is a cloud-based application that allows non-computational biologists to easily utilize advanced AI-based structural biology tools, including AlphaFold and DiffDock. With many deployment options, it can be employed by individuals, labs, universities, and companies in the cloud without requiring hardware resources, but it can also be configured to utilize locally available computers. Foldy enables scientists to predict the structure of proteins and complexes up to 6000 amino acids with AlphaFold, visualize Pfam annotations, and dock ligands with AutoDock Vina and DiffDock. In our manuscript, we detail Foldy’s interface design, deployment strategies, and optimization for various user scenarios. We demonstrate its application through case studies including rational enzyme design and analyzing proteins with domains of unknown function. Furthermore, we compare Foldy’s interface and management capabilities with other open and closed source tools in the field, illustrating its practicality in managing complex data and computation tasks. Our manuscript underlines the benefits of Foldy as a day-to-day tool for life science researchers, and shows how Foldy can make modern tools more accessible and efficient.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

First comprehensive identification of cardiac proteins with putative increased O-GlcNAc levels during pressure overload hypertrophy

Protein posttranslational modifications (PTMs) by O-GlcNAc globally rise during pressure-overload hypertrophy (POH). However, a major knowledge gap exists on the specific proteins undergoing changes in O-GlcNAc levels during POH primarily because this PTM is low abundance and easily lost during standard mass spectrometry (MS) conditions used for protein identification. Methodologies have emerged to enrich samples for O-GlcNAcylated proteins prior to MS analysis. Accordingly, our goal was to identify the specific proteins undergoing changes in O-GlcNAc levels during POH. We used C57/Bl6 mice subjected to Sham or transverse aortic constriction (TAC) to create POH. From the hearts, we labelled the O-GlcNAc moiety with tetramethylrhodamine azide (TAMRA) before sample enrichment by TAMRA immunoprecipitation (IP). We used LC-MS/MS to identify and quantify the captured putative O-GlcNAcylated proteins. We identified a total of 700 putative O-GlcNAcylated proteins in Sham and POH. Two hundred thirty-three of these proteins had significantly increased enrichment in POH over Sham suggesting higher O-GlcNAc levels whereas no proteins were significantly decreased by POH. We examined two MS identified metabolic enzymes, CPT1B and the PDH complex, to validate by immunoprecipitation. We corroborated increased O-GlcNAc levels during POH for CPT1B and the PDH complex. Enzyme activity assays suggests higher O-GlcNAcylation increases CPT1 activity and decreases PDH activity during POH. In summary, we generated the first comprehensive list of proteins with putative changes in O-GlcNAc levels during POH. Our results demonstrate the large number of potential proteins and cellular processes affected by O-GlcNAc and serve as a guide for testing specific O-GlcNAc-regulated mechanisms during POH.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Genome Analysis of a Verrucomicrobial Endosymbiont With a Tiny Genome Discovered in an Antarctic Lake

Organic Lake in Antarctica is a marine-derived, cold (–13°C), stratified (oxic-anoxic), hypersaline (>200 gl –1 ) system with unusual chemistry (very high levels of dimethylsulfide) that supports the growth of phylogenetically and metabolically diverse microorganisms. Symbionts are not well characterized in Antarctica. However, unicellular eukaryotes are often present in Antarctic lakes and theoretically could harbor endosymbionts. Here, we describe Candidatus Organicella extenuata, a member of the Verrucomicrobia with a highly reduced genome, recovered as a metagenome-assembled genome with genetic code 4 (UGA-to-Trp recoding) from Organic Lake. It is closely related to Candidatus Pinguicocccus supinus (163,218 bp, 205 genes), a newly described cytoplasmic endosymbiont of the freshwater ciliate Euplotes vanleeuwenhoeki (Serra et al., 2020). At 158,228 bp (encoding 194 genes), the genome of Ca. Organicella extenuata is among the smallest known bacterial genomes and similar to the genome of Ca. Pinguicoccus supinus (163,218 bp, 205 genes). Ca. Organicella extenuata retains a capacity for replication, transcription, translation, and protein-folding while lacking any capacity for the biosynthesis of amino acids or vitamins. Notably, the endosymbiont retains a capacity for fatty acid synthesis (type II) and iron–sulfur (Fe-S) cluster assembly. Metagenomic analysis of 150 new metagenomes from Organic Lake and more than 70 other Antarctic aquatic locations revealed a strong correlation in abundance between Ca. Organicella extenuata and a novel ciliate of the genus Euplotes. Like Ca. Pinguicoccus supinus, we infer that Ca. Organicella extenuata is an endosymbiont of Euplotes and hypothesize that both Ca. Organicella extenuata and Ca. Pinguicocccus supinus provide fatty acids and Fe-S clusters to their Euplotes host as the foundation of a mutualistic symbiosis. The discovery of Ca. Organicella extenuata as possessing genetic code 4 illustrates that in addition to identifying endosymbionts by sequencing known symbiotic communities and searching metagenome data using reference endosymbiont genomes, the potential exists to identify novel endosymbionts by searching for unusual coding parameters.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Plant root associated chitinases: structures and functions

Chitinases degrade chitin, a linear homopolymer of β-1,4-linked N-acetyl-D-glucosamine (GlcNAc) residues found in the cell walls of fungi and the exoskeletons of arthropods. They are secreted by the roots into the rhizosphere, a complex and dynamic environment where intense nutrient exchange occurs between plants and microbes. Here we modeled, expressed, purified, and characterized Zea mays and Oryza sativa root chitinases, and the chitinase of a symbiotic bacterium, Chitinophaga oryzae 1303 for their activities with chitin, di-, tri-, and tetra-saccharides and Aspergillus niger, with the goal of determining their role(s) in the rhizosphere and better understanding the molecular mechanisms underlying plant-microbe interactions. We show that Zea mays basic endochitinase (ZmChi19A) and Oryza sativa chitinase (OsChi19A) are from the GH19 chitinase family. The Chitinophaga oryzae 1303 chitinase (CspCh18A) belongs to the GH18 family. The three enzymes have similar apparent K M values of (20-40 µM) for the substrate 4-MU-GlcNAc 3 . They vary in their pH and temperature optima with OsChi19A activity optimal between pH 5–7 and 30–40°C while ZmChi19A and CspCh18A activities were optimal at pH 7-9 and 50–60°C. Modeling and site-directed mutation of ZmChi19A identified the catalytic cleft and the active residues E147 and E169 strategically positioned at ~8.6Å from each other in the folded protein. Cleavage of 4-MU-GlcNAc 3 was unaffected by the absence of the CBD but diminished in the absence of the flexible C-terminal domain. However, unlike for the soluble substrate, the CBD and the newly identified flexible C-terminal domain were vital for inhibiting Aspergillus niger growth. The results are consistent with the involvement of the plant chitinases in defense against pathogens like fungi that have chitin exoskeletons. In summary, we have characterized the functional features and structural domains necessary for the activity of two plant root chitinases that are believed to be involved in plant defense and a bacterial chitinase that, along with the plant chitinases, may participate in nutrient recycling in the rhizosphere.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Mechanical Properties of a Solvated Biomolecule: RGD (1FUV) Peptide

The mechanical properties of proteins/peptides play an essential role in their functionalities and implications, as well as their structure and dynamic properties. Understanding mechanical properties is pivotal to our knowledge of protein folding and the molecular basis of diverse cellular processes. Herein, we present a computational approach using ab initio quantum mechanical calculations to determine the mechanical properties—such as bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio—of a solvated Arg-Gly-Asp (RGD) peptide model. Since this peptide serves as the RGD-directed integrin recognition site and may participate in cellular adhesion, it is considered a promising small peptide for medicinal applications. This successful approach paves the way for investigating larger and more complex biomolecules.

Biochemistry & Molecular Biology↗
NASA NTRS1993

The molecular matching problem

Molecular chemistry contains many difficult optimization problems that have begun to attract the attention of optimizers in the Operations Research community. Problems including protein folding, molecular conformation, molecular similarity, and molecular matching have been addressed. Minimum energy conformations for simple molecular structures such as water clusters, Lennard-Jones microclusters, and short polypeptides have dominated the literature to date. However, a variety of interesting problems exist and we focus here on a molecular structure matching (MSM) problem.

Kincaid, Rex K.↗
NASA NTRS2001

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗
NASA NTRS1998

Acid-growth response and alpha-expansins in suspension cultures of bright yellow 2 tobacco

The possibility that Bright Yellow 2 (BY2) tobacco (Nicotiana tabacum L.) suspension-cultured cells possess an expansin-mediated acid-growth mechanism was examined by multiple approaches. BY2 cells grew three times faster upon treatment with fusicoccin, which induces an acidification of the cell wall. Exogenous expansins likewise stimulated BY2 cell growth 3-fold. Protein extracted from BY2 cell walls possessed the expansin-like ability to induce extension of isolated walls. In western-blot analysis of BY2 wall protein, one band of 29 kD was recognized by anti-expansin antibody. Six different classes of alpha-expansin mRNA were identified in a BY2 cDNA library. Northern-blot analysis indicated moderate to low abundance of multiple alpha-expansin mRNAs in BY2 cells. From these results we conclude that BY2 suspension-cultured cells have the necessary components for expansin-mediated cell wall enlargement.

Non-NASA Center↗
NASA NTRS1995

Finding the global minimum: a fuzzy end elimination implementation

The 'fuzzy end elimination theorem' (FEE) is a mathematically proven theorem that identifies rotameric states in proteins which are incompatible with the global minimum energy conformation. While implementing the FEE we noticed two different aspects that directly affected the final results at convergence. First, the identification of a single dead-ending rotameric state can trigger a 'domino effect' that initiates the identification of additional rotameric states which become dead-ending. A recursive check for dead-ending rotameric states is therefore necessary every time a dead-ending rotameric state is identified. It is shown that, if the recursive check is omitted, it is possible to miss the identification of some dead-ending rotameric states causing a premature termination of the elimination process. Second, we examined the effects of removing dead-ending rotameric states from further considerations at different moments of time. Two different methods of rotameric state removal were examined for an order dependence. In one case, each rotamer found to be incompatible with the global minimum energy conformation was removed immediately following its identification. In the other, dead-ending rotamers were marked for deletion but retained during the search, so that they influenced the evaluation of other rotameric states. When the search was completed, all marked rotamers were removed simultaneously. In addition, to expand further the usefulness of the FEE, a novel method is presented that allows for further reduction in the remaining set of conformations at the FEE convergence. In this method, called a tree-based search, each dead-ending pair of rotamers which does not lead to the direct removal of either rotameric state is used to reduce significantly the number of remaining conformations. In the future this method can also be expanded to triplet and quadruplet sets of rotameric states. We tested our implementation of the FEE by exhaustively searching ten protein segments and found that the FEE identified the global minimum every time. For each segment, the global minimum was exhaustively searched in two different environments: (i) the segments were extracted from the protein and exhaustively searched in the absence of the surrounding residues; (ii) the segments were exhaustively searched in the presence of the remaining residues fixed at crystal structure conformations. We also evaluated the performance of the method for accurately predicting side chain conformations. We examined the influence of factors such as type and accuracy of backbone template used, and the restrictions imposed by the choice of potential function, parameterization and rotamer database. Conclusions are drawn on these results and future prospects are given.

NASA Program Exobiology↗
NASA NTRS1993

The 'tubulin-like' S1 protein of Spirochaeta is a member of the hsp65 stress protein family

A 65-kDa protein (called S1) from Spirochaeta bajacaliforniensis was identified as 'tubulin-like' because it cross-reacted with at least four different antisera raised against tubulin and was isolated, with a co-polymerizing 45-kDa protein, by warm-cold cycling procedures used to purify tubulin from mammalian brain. Furthermore, at least three genera of non-cultivable symbiotic spirochetes (Pillotina, Diplocalyx, and Hollandina) that contain conspicuous 24-nm cytoplasmic tubules displayed a strong fluorescence in situ when treated with polyclonal antisera raised against tubulin. Here we summarize results that lead to the conclusion that this 65-kDa protein has no homology to tubulin. S1 is an hsp65 stress protein homologue. Hsp65 is a highly immunogenic family of hsp60 proteins which includes the 65-kDa antigens of Mycobacterium tuberculosis (an active component of Freund's complete adjuvant), Borrelia, Treponema, Chlamydia, Legionella, and Salmonella. The hsp60s, also known as chaperonins, include E. coli GroEL, mitochondrial and chloroplast chaperonins, the pea aphid 'symbionin' and many other proteins involved in protein folding and the stress response.

NASA Discipline Exobiology↗
NASA NTRS2006

Water as a matrix for life

"Follow the water" is the canonical strategy in searching for life in the universe. Conventionally, discussion of this topic is focused on how solvent supports organic chemistry sufficiently rich to seed life. Perhaps more importantly, solvent must promote self-organization of organic matter into functional structures capable of responding to environmental changes. This process is based on non-covalent interactions. They are constantly formed and broken in response to internal and external stimuli. This requires that their strength must be properly tuned. If they were too weak, the system would exhibit undesired, uncontrolled response to natural fluctuations of physical and chemical parameters. If they were too strong kinetics of biological processes would be slow and energetics costly. Non-covalent interactions are strongly mediated by the solvent. Specifically, high dielectric solvents for life are needed for solubility of polar species and flexibility of biological structures stabilized by electrostatic interactions. Water exhibits a remarkable trait that it promotes solvophobic interactions between non-polar species, which are responsible for self-organization phenomena such as the formation of cellular boundary structures, and protein folding and aggregation. Unusual temperature dependence of hydrophobic interactions - they often become stronger as temperature increases - is a consequence of the temperature insensitivity of properties of the liquid water. This contributes to the existence of robust life over a wide temperature range. Water is not the only liquid with favorable properties for supporting life. Other pure liquids or their mixtures that have high dielectric constants and simultaneously support some level of self-organization will be discussed.

Pohorille, Andrew↗
NASA NTRS2012

Is Water a Universal Solvent for Life?

There are strong reasons to believe that the laws, principles and constraints of physics and chemistry are universal. It is much less clear how this universality translates into our understanding of the origins of life. Conventionally, discussions of this topic focus on chemistry that must be sufficiently rich to seed life. Although this is clearly a prerequisite for the emergence of living systems, I propose to focus instead on self-organization of matter into functional structures capable of reproduction, evolution and responding to environmental changes. In biology, most essential functions are largely mediated by noncovalent interactions (interactions that do not involve making or breaking chemical bonds). Forming chemical bonds is only a small part of what living systems do. There are specific implications of this point of view for universality. I will concentrate on one of these implications. Strength of non-covalent interactions must be properly tuned. If they were too weak, the system would exhibit undesired, uncontrolled response to natural fluctuations of physical and chemical parameters. If they were too strong kinetics of biological processes would be slow and energetics costly. This balance, however, is not a natural property of complex chemical systems. Instead, it has to be achieved with the aid of an appropriate solvent for life. In particular, potential solvents for life must be characterized by a high dielectric constant to ensure solubility of polar species and sufficient flexibility of biological structures stabilized by electrostatic interactions. Among these solvents, water exhibits a remarkable trait that it also promotes solvophobic (hydrophobic) interactions between non-polar species, typically manifested by a tendency of these species to aggregate and minimize their contacts with the aqueous solvent. Hydrophobic interactions are responsible, at least in part, for many self-organization phenomena in biological systems, such as the formation of cellular boundary structures or protein folding. Strengths of electrostatic and hydrophobic interactions are similar and can be balanced over a wide range of temperatures, which considerably increases the repertoire of interactions that can be used to modulate biological functions. Some properties of water, e.g. its chemical activity against polymerization reactions, are considered as unfavorable to life. In actuality, this might be a favorable trait because life requires a balance between constructive and destructive processes. For example, molecules synthesized in response to specific conditions must be degraded once these conditions change. Otherwise regulation of biological processes would be virtually impossible. Water might not be the only liquid with favorable properties for supporting life. It has been proposed that formamide, which might be present elsewhere in the universe in sufficient quantities to warrant interest, could be a potential alternative to water for the origin of life. However, this will remain highly hypothetical until it is demonstrated in further studies on its physical, chemical and biological properties it is capable of mediating self-organization of matter and providing proper balance between different types of non-covalent interactions.

Pohorill, Andrew↗
NASA NTRS2018

Is Water Necessary for Life?

The universality of water as the solvent for life is usually justified by its role in supporting the rich organic chemistry. It has been pointed out, however, that even richer synthetic chemistry is possible in other organic solvents. Does it mean that water is not necessary for life? Here, other, essential criteria for solvent for life that have not been sufficiently considered are discussed. In biological systems, complex molecules are not only constantly synthesized but also degraded. Solvent-mediated degradation is essential for regulating cell content, preventing overcrowding and allowing for recycling organic material. Achieving a balance between synthetic and degradative processes is facile in water, but not in many other organic liquids. Thus, the so-called “water paradox” according to which water is both necessary to life and toxic to biopolymer synthesis might not be paradoxical at all. The machinery of life is based on non-covalent interactions that do not involve making or breaking chemical bonds. Their strength needs to be properly tuned. If they are too weak, there might be undesired response to natural fluctuations of physical or chemical parameters. If they are too strong, the kinetics and energetics of cellular processes could adversely influenced. The solvent must allow for balancing these interactions, which provides strong, universal constraints on the medium for life. Water influences non-covalent interactions mainly by two mechanisms. First, it reduces strong, electrostatic interactions between molecules, chemical groups or atoms carrying electric charge or dipole. Second, it induces the hydrophobic effect, the tendency to remove non-polar (hydrophobic) molecules and groups from direct contact with aqueous solution and, instead, interact with each other. In living systems, the hydrophobic effect is largely responsible for self-organization of molecules to more complex structures, such as aggregation of lipid molecules to form biological membranes and protein folding. Water exists as stable liquid in a large temperature range, and the hydrophobic effects are a consequence of the temperature insensitivity of essential properties of its liquid state. In summary, water accomplishes an amazing feat – it reduces strong interactions between dissolved species and simultaneously increases the strength of weak interactions, bringing all of them to the right range. Once we consider not only synthetic capabilities but also other required traits of the solvent for life, no viable alternative to water is currently known.

Life↗
DOE OSTI2023

Population-based heteropolymer design to mimic protein mixtures

Biological fluids, the most complex blends, have compositions that constantly vary and cannot be molecularly defined. Despite these uncertainties, proteins fluctuate, fold, function and evolve as programmed. We propose that in addition to the known monomeric sequence requirements, protein sequences encode multi-pair interactions at the segmental level to navigate random encounters; synthetic heteropolymers capable of emulating such interactions can replicate how proteins behave in biological fluids individually and collectively. Here, we extracted the chemical characteristics and sequential arrangement along a protein chain at the segmental level from natural protein libraries and used the information to design heteropolymer ensembles as mixtures of disordered, partially folded and folded proteins. For each heteropolymer ensemble, the level of segmental similarity to that of natural proteins determines its ability to replicate many functions of biological fluids including assisting protein folding during translation, preserving the viability of fetal bovine serum without refrigeration, enhancing the thermal stability of proteins and behaving like synthetic cytosol under biologically relevant conditions. Molecular studies further translated protein sequence information at the segmental level into intermolecular interactions with a defined range, degree of diversity and temporal and spatial availability. This framework provides valuable guiding principles to synthetically realize protein properties, engineer bio/abiotic hybrid materials and, ultimately, realize matter-to-life transformations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

De novo design of knotted tandem repeat proteins

De novo protein design methods can create proteins with folds not yet seen in nature. These methods largely focus on optimizing the compatibility between the designed sequence and the intended conformation, without explicit consideration of protein folding pathways. Deeply knotted proteins, whose topologies may introduce substantial barriers to folding, thus represent an interesting test case for protein design. Here we report our attempts to design proteins with trefoil (3 1 ) and pentafoil (5 1 ) knotted topologies. We extended previously described algorithms for tandem repeat protein design in order to construct deeply knotted backbones and matching designed repeat sequences (N = 3 repeats for the trefoil and N = 5 for the pentafoil). We confirmed the intended conformation for the trefoil design by X ray crystallography, and we report here on this protein’s structure, stability, and folding behaviour. The pentafoil design misfolded into an asymmetric structure (despite a 5-fold symmetric sequence); two of the four repeat-repeat units matched the designed backbone while the other two diverged to form local contacts, leading to a trefoil rather than pentafoil knotted topology. Our results also provide insights into the folding of knotted proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Controlling Amphiphilic Polymer Folding beyond the Primary Structure with Protein-Mimetic Di(Phenylalanine)

While methods for polymer synthesis have proliferated, their functionality pales in comparison to natural biopolymers–strategies are limited for building the intricate network of noncovalent interactions necessary to elicit complex, protein-like functions. Using a bioinspired di(phenylalanine) acrylamide (FF) monomer, we explored the impact of various noncovalent interactions in generating ordered assembled structures. Amphiphilic copolymers were synthesized that exhibit β-sheet-like local structure upon collapsing into single-chain assemblies in aqueous environments. Here, the systematic analysis of a series of amphiphilic copolymers illustrated that the global collapse is primarily driven by hydrophobic forces. Hydrogen-bonding and aromatic interactions stabilize local structure, as β-sheet-like interactions were identified via circular dichroism and thioflavin T fluorescence. Similar analysis of phenylalanine (F) and alanine-phenylalanine acrylamide (AF) copolymers found that distancing the aromatic residue from the polymer backbone is sufficient to induce β-sheet-like local structure akin to the FF copolymers; however, the interactions between AF subunits are less stable than those formed by FF. Further, hydrogen-bond donating hydrophilic monomers disrupt internal structure formed by FF within collapsed assemblies. Collectively, these results illuminate design principles for the facile incorporation of multiple facets of protein-mimetic, higher-order structure within folded synthetic polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗
DOE OSTI2021

The Promise of Mutation Resistant Drugs for SARS-CoV-2 That Interdict in the Folding of the Spike Protein Receptor Binding Domain

In recent work, we proposed that effective therapeutic drugs aimed at treating the SARS-CoV-2 infection could be developed based on interdicting in the early steps of the folding pathway of key viral proteins, including the receptor binding domain (RBD) of the spike protein. In order to provide for a drug target on the protein, the earliest contact-formation event along the dominant folding pathway of the RBD spike protein was predicted employing the Sequential Collapse Model (SCM). The segments involved in the predicted earliest contact were suggested to provide optimal folding interdiction target regions (FITRs) for potential therapeutic drugs, with a focus on folding interdicting peptides (FIPs). In this paper, we extend our analysis to include 13 known single mutations of the RBD spike protein as well as the triple mutation B1.351 and the recent double mutation B1.617.2. The results show that the location of the FITR does not change in any of the 15 studied mutations, providing for a mutation-resistant drug design strategy for the RBD-spike protein.

Bergasa-Caceres, Fernando↗