Search NASA⌕ Search

SEARCH · Search NASA

Results for “protein design”

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 451 records · Page 25

Expanding the synthetic biology toolbox of Cupriavidus necator for establishing fatty acid production

Abstract The Gram-negative betaproteobacterium Cupriavidus necator is a chemolithotroph that can convert carbon dioxide into biomass. Cupriavidus necator has been engineered to produce a variety of high-value chemicals in the past. However, there is still a lack of a well-characterized toolbox for gene expression and genome engineering. Development and optimization of biosynthetic pathways in metabolically engineered microorganisms necessitates control of gene expression via functional genetic elements such as promoters, ribosome binding sites (RBSs), and codon optimization. In this work, a set of inducible and constitutive promoters were validated and characterized in C. necator, and a library of RBSs was designed and tested to show a 50-fold range of expression for green fluorescent protein (gfp). The effect of codon optimization on gene expression in C. necator was studied by expressing gfp and mCherry genes with varied codon-adaptation indices and was validated by expressing codon-optimized variants of a C12-specific fatty acid thioesterase to produce dodecanoic acid. We discuss further hurdles that will need to be overcome for C. necator to be widely used for biosynthetic processes.

59 BASIC BIOLOGICAL SCIENCES↗

Restoring Pre-Industrial CO2 Levels While Achieving Sustainable Development Goals

Unless humanity achieves United Nations Sustainable Development Goals (SDGs) by 2030 and restores the relatively stable climate of pre-industrial CO2 levels (as early as 2140), species extinctions, starvation, drought/floods, and violence will exacerbate mass migrations. This paper presents conceptual designs and techno-economic analyses to calculate sustainable limits for growing high-protein seafood and macroalgae-for-biofuel. We review the availability of wet solid waste and outline the mass balance of carbon and plant nutrients passing through a hydrothermal liquefaction process. The paper reviews the availability of dry solid waste and dry biomass for bioenergy with CO2 capture and storage (BECCS) while generating Allam Cycle electricity. Sufficient wet-waste biomass supports quickly building hydrothermal liquefaction facilities. Macroalgae-for-biofuel technology can be developed and straightforwardly implemented on SDG-achieving high protein seafood infrastructure. The analyses indicate a potential for (1) 0.5 billion tonnes/yr of seafood; (2) 20 million barrels/day of biofuel from solid waste; (3) more biocrude oil from macroalgae than current fossil oil; and (4) sequestration of 28 to 38 billion tonnes/yr of bio-CO2. Carbon dioxide removal (CDR) costs are between 25–33% of those for BECCS with pre-2019 technology or the projected cost of air-capture CDR.

54 ENVIRONMENTAL SCIENCES↗

A user-friendly plug-and-play cyclic olefin copolymer-based microfluidic chip for room-temperature, fixed-target serial crystallography

Over the past two decades, serial X-ray crystallography has enabled the structure determination of a wide range of proteins. With the advent of X-ray free-electron lasers (XFELs), ever-smaller crystals have yielded high-resolution diffraction and structure determination. A crucial need to continue advancement is the efficient delivery of fragile and micrometre-sized crystals to the X-ray beam intersection. This paper presents an improved design of an all-polymer microfluidic `chip' for room-temperature fixed-target serial crystallography that can be tailored to broadly meet the needs of users at either synchrotron or XFEL light sources. The chips are designed to be customized around different types of crystals and offer users a friendly, quick, convenient, ultra-low-cost and robust sample-delivery platform. Compared with the previous iteration of the chip [Gilbile et al. (2021), Lab Chip , 21 , 4831–4845], the new design eliminates cleanroom fabrication. It has a larger imaging area to volume, while maintaining crystal hydration stability for both in situ crystallization or direct crystal slurry loading. Crystals of two model proteins, lysozyme and thaumatin, were used to validate the effectiveness of the design at both synchrotron (lysozyme and thaumatin) and XFEL (lysozyme only) facilities, yielding complete data sets with resolutions of 1.42, 1.48 and 1.70 Å, respectively. Overall, the improved chip design, ease of fabrication and high modifiability create a powerful, all-around sample-delivery tool that structural biologists can quickly adopt, especially in cases of limited sample volume and small, fragile crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design Principles for Using Amphiphilic Polymers To Create Microporous Water

Aqueous dispersions of microporous nanocrystals with dry, gas-accessible pores—referred to as “microporous water”—enable high densities of gas molecules to be transported through water. For many applications of microporous water, generalizable strategies are required to functionalize the external surface of microporous particles to control their dispersibility, stability, and interactions with other solution-phase components—including catalysts, proteins, and cells—while retaining as much of their internal pore volume as possible. Here, in this work, we establish design principles for the noncovalent surface functionalization of hydrophobic metal–organic frameworks with amphiphilic polymers that render the particles dispersible in water and enhance their hydrolytic stability. Specifically, we show that block co-polymers with persistence lengths that exceed the micropore aperture size of zeolitic imidazolate frameworks (ZIFs) can dramatically enhance ZIF particle dispersibility and stability while preserving porosity and >80% of the theoretical O 2 carrying capacity. Moreover, enhancements in hydrolytic stability are greatest when the polymer can form strong bonds to exposed metal sites on the external particle surface. More broadly, our insights provide guidelines for controlling the interface between polymers and metal–organic framework particles in aqueous environments to augment the properties of microporous water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recombinant expression of computationally designed peptide-bundlemers in Escherichia coli

Computational design of fully artificial peptides is extensively researched by material scientists and engineers for the construction of novel nanostructures and biomaterials. Such design has yielded a peptide-based building block or bundlemer, a coiled coil peptide assembly that undergoes further physical-covalent interactions to form 1D, 2D and, potentially, 3D hierarchical assemblies and displays targeted and biomimetic material properties. Recombinant expression is a convenient, flexible tool to synthesize such artificial and modified peptides in large quantities while also enabling economical synthesis of isotopically labeled peptides and longer protein-like artificial peptides. This report describes the protocol for recombinant expression of a 31-amino acid, computationally designed bundlemer-forming peptide in Escherichia coli. Peptide yields of 10 mgs per liter of media were achieved which highlights complementary advantages of recombinant expression technique relative to conventional laboratory-scale synthesis, such as solid-phase peptide synthesis.

36 MATERIALS SCIENCE↗

Decoding the protein–ligand interactions using parallel graph neural networks

Abstract Protein–ligand interactions (PLIs) are essential for biochemical functionality and their identification is crucial for estimating biophysical properties for rational therapeutic design. Currently, experimental characterization of these properties is the most accurate method, however, this is very time-consuming and labor-intensive. A number of computational methods have been developed in this context but most of the existing PLI prediction heavily depends on 2D protein sequence data. Here, we present a novel parallel graph neural network (GNN) to integrate knowledge representation and reasoning for PLI prediction to perform deep learning guided by expert knowledge and informed by 3D structural data. We develop two distinct GNN architectures: $$\hbox {GNN}_{\mathrm{F}}$$ GNN F is the base implementation that employs distinct featurization to enhance domain-awareness, while $$\hbox {GNN}_{\mathrm{P}}$$ GNN P is a novel implementation that can predict with no prior knowledge of the intermolecular interactions. The comprehensive evaluation demonstrated that GNN can successfully capture the binary interactions between ligand and protein’s 3D structure with 0.979 test accuracy for $$\hbox {GNN}_{\mathrm{F}}$$ GNN F and 0.958 for $$\hbox {GNN}_{\mathrm{P}}$$ GNN P for predicting activity of a protein–ligand complex. These models are further adapted for regression tasks to predict experimental binding affinities and $$\hbox {pIC}_{\mathrm{50}}$$ pIC 50 crucial for compound’s potency and efficacy. We achieve a Pearson correlation coefficient of 0.66 and 0.65 on experimental affinity and 0.50 and 0.51 on $$\hbox {pIC}_{\mathrm{50}}$$ pIC 50 with $$\hbox {GNN}_{\mathrm{F}}$$ GNN F and $$\hbox {GNN}_{\mathrm{P}}$$ GNN P , respectively, outperforming similar 2D sequence based models. Our method can serve as an interpretable and explainable artificial intelligence (AI) tool for predicted activity, potency, and biophysical properties of lead candidates. To this end, we show the utility of $$\hbox {GNN}_{\mathrm{P}}$$ GNN P on SARS-Cov-2 protein targets by screening a large compound library and comparing the prediction with the experimentally measured data.

59 BASIC BIOLOGICAL SCIENCES↗

Enhanced Binding of SARS-CoV-2 Spike Protein to Receptor by Distal Polybasic Cleavage Sites

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein plays a crucial role in binding the human cell receptor ACE2 that is required for viral entry. Many studies have been conducted to target the structures of RBD–ACE2 binding and to design RBD-targeting vaccines and drugs. Nevertheless, mutations distal from the SARS-CoV-2 RBD also impact its transmissibility and antibody can target non-RBD regions, suggesting the incomplete role of the RBD region in the spike protein–ACE2 binding. Here, in order to elucidate distant binding mechanisms, we analyze complexes of ACE2 with the wild-type spike protein and with key mutants via large-scale all-atom explicit solvent molecular dynamics simulations. We find that though distributed approximately 10 nm away from the RBD, the SARS-CoV-2 polybasic cleavage sites enhance, via electrostatic interactions and hydration, the RBD–ACE2 binding affinity. A negatively charged tetrapeptide (GluGluLeuGlu) is then designed to neutralize the positively charged arginine on the polybasic cleavage sites. We find that the tetrapeptide GluGluLeuGlu binds to one of the three polybasic cleavage sites of the SARS-CoV-2 spike protein lessening by 34% the RBD–ACE2 binding strength. This significant binding energy reduction demonstrates the feasibility to neutralize RBD–ACE2 binding by targeting this specific polybasic cleavage site. Our work enhances understanding of the binding mechanism of SARS-CoV-2 to ACE2, which may aid the design of therapeutics for COVID-19 infection.

60 APPLIED LIFE SCIENCES↗

Using mechanistic models and machine learning to design single-color multiplexed nascent chain tracking experiments

mRNA translation is the ubiquitous cellular process of reading messenger-RNA strands into functional proteins. Over the past decade, large strides in microscopy techniques have allowed observation of mRNA translation at a single-molecule resolution for self-consistent time-series measurements in live cells. Dubbed Nascent chain tracking (NCT), these methods have explored many temporal dynamics in mRNA translation uncaptured by other experimental methods such as ribosomal profiling, smFISH, pSILAC, BONCAT, or FUNCAT-PLA. However, NCT is currently restricted to the observation of one or two mRNA species at a time due to limits in the number of resolvable fluorescent tags. In this work, we propose a hybrid computational pipeline, where detailed mechanistic simulations produce realistic NCT videos, and machine learning is used to assess potential experimental designs for their ability to resolve multiple mRNA species using a single fluorescent color for all species. Our simulation results show that with careful application this hybrid design strategy could in principle be used to extend the number of mRNA species that could be watched simultaneously within the same cell. We present a simulated example NCT experiment with seven different mRNA species within the same simulated cell and use our ML labeling to identify these spots with 90% accuracy using only two distinct fluorescent tags. We conclude that the proposed extension to the NCT color palette should allow experimentalists to access a plethora of new experimental design possibilities, especially for cell Signaling applications requiring simultaneous study of multiple mRNAs.

59 BASIC BIOLOGICAL SCIENCES↗

Synthesis and in silico and in vitro evaluation of trimethoxy-benzamides designed as anti-prion derivatives

Transmissible spongiform encephalopathies (TSEs), also known as prion diseases, are neurodegenerative disorders which affect mammals, including the human species, and arise after the conversion of the monomeric cellular prion protein (PrP{sup C}) into the aggregated scrapie form (PrP{sup Sc}). There is no therapy to treat TSEs and the identification of compounds that bind PrP{sup C}, preventing its conversion into PrP{sup Sc}, is a viable therapeutic strategy. We designed and synthesized six novel trimethoxy-benzamide compounds as anti-prion drug candidates. Molecular docking analyses predicted that all the derivatives bind to a hotspot region located in the PrP globular domain with very similar spatial orientation and interaction mode. Although none of the analogs inhibited in vitro-aggregation of recombinant PrP (rPrP) in a cell-free conversion assay, the RT-QuIC, compound 8a accelerated rPrP conversion into PrP{sup Sc}-like species. STD-NMR and ITC analyses indicated that both 8a and 8b bind to rPrP{sup 90–231}. These analogs were toxic to PrP{sup Sc}-infected cell lines, hence we could not assess their anti-prion activity by using this cellular approach, although this toxicity was cell line-dependent. These results point out that the 4-amino-quinoline trimethoxy-benzamide scaffold described herein represents a novel chemical pattern useful as a starting point for future structural optimization in the design of PrP ligands with improved affinity and safety profiles.

60 APPLIED LIFE SCIENCES↗

A protective and broadly binding antibody class engages the influenza virus hemagglutinin head at its stem interface

ABSTRACT Influenza infection and vaccination impart strain-specific immunity that protects against neither seasonal antigenic variants nor the next pandemic. However, antibodies directed to conserved sites can confer broad protection. Here, we identify and characterize a class of human antibodies that engage a previously undescribed, conserved epitope on the influenza hemagglutinin (HA) protein. Prototype antibody S8V1-157 binds at the normally occluded interface between the HA head and stem. Antibodies to this HA head–stem interface epitope are non-neutralizingin vitrobut protect against lethal influenza infection in mice. These antibodies bind to most influenza A subtypes and seasonal human variants, and are present at low frequencies in the memory B cell populations of multiple human donors. Vaccines designed to elicit these antibodies might contribute to “universal” influenza immunity. IMPORTANCE Antibodies to the influenza virus hemagglutinin (HA) protein confer the strongest protection against infection. Human antibodies elicited by infection and/or vaccination fail to protect against antigenically novel animal, pandemic, or human seasonal viruses. Improved vaccines are needed. We identify a novel class of antibodies that bind most divergent HA subtypes and all seasonal human HA antigenic variants tested. These antibodies confer protection from lethal influenza challenge in animal models. The corresponding epitope on the HA head is occluded by its interaction with the stem and is inaccessible in the well-resolved prefusion state. The immunogenicity of this head–stem interface indicates that poorly understood conformations of HA presenting widely conserved surfaces are explored in biochemical, cell-based, andin vivoassays. Head–stem interface antibodies warrant further investigation as an avenue to improve influenza vaccines and therapeutics.

Microbiology↗

A neutralizing antibody target in early HIV-1 infection was recapitulated in rhesus macaques immunized with the transmitted/founder envelope sequence

Transmitted/founder (T/F) HIV-1 envelope proteins (Envs) from infected individuals that developed neutralization breadth are likely to possess inherent features desirable for vaccine immunogen design. To explore this premise, we conducted an immunization study in rhesus macaques (RM) using T/F Env sequences from two human subjects, one of whom developed potent and broad neutralizing antibodies (Z1800M) while the other developed little to no neutralizing antibody responses (R66M) during HIV-1 infection. Using a DNA/MVA/protein immunization protocol, 10 RM were immunized with each T/F Env. Within each T/F Env group, the protein boosts were administered as either monomeric gp120 or stabilized trimeric gp140 protein. All vaccination regimens elicited high titers of antigen-specific IgG, and two animals that received monomeric Z1800M Env gp120 developed autologous neutralizing activity. Using early Env escape variants isolated from subject Z1800M as guides, the serum neutralizing activity of the two immunized RM was found to be dependent on the gp120 V5 region. Interestingly, the exact same residues of V5 were also targeted by a neutralizing monoclonal antibody (nmAb) isolated from the subject Z1800M early in infection. Glycan profiling and computational modeling of the Z1800M Env gp120 immunogen provided further evidence that the V5 loop is exposed in this T/F Env and was a dominant feature that drove neutralizing antibody targeting during infection and immunization. An expanded B cell clonotype was isolated from one of the neutralization-positive RM and nmAbs corresponding to this group demonstrated V5-dependent neutralization similar to both the RM serum and the human Z1800M nmAb. The results demonstrate that neutralizing antibody responses elicited by the Z1800M T/F Env in RM converged with those in the HIV-1 infected human subject, illustrating the potential of using immunogens based on this or other T/F Envs with well-defined immunogenicity as a starting point to drive breadth.

59 BASIC BIOLOGICAL SCIENCES↗

Drug discovery efforts at George Mason University

With over 39,000 students, and research expenditures in excess of $200 million, George Mason University (GMU) is the largest R1 (Carnegie Classification of very high research activity) university in Virginia. Mason scientists have been involved in the discovery and development of novel diagnostics and therapeutics in areas as diverse as infectious diseases and cancer. Below are highlights of the efforts being led by Mason researchers in the drug discovery arena. To enable targeted cellular delivery, and non-biomedical applications, Veneziano and colleagues have developed a synthesis strategy that enables the design of self-assembling DNA nanoparticles (DNA origami) with prescribed shape and size in the 10 to 100 nm range. The nanoparticles can be loaded with molecules of interest such as drugs, proteins and peptides, and are a promising new addition to the drug delivery platforms currently in use. The investigators also recently used the DNA origami nanoparticles to fine tune the spatial presentation of immunogens to study the impact on B cell activation. These studies are an important step towards the rational design of vaccines for a variety of infectious agents. To elucidate the parameters for optimizing the delivery efficiency of lipid nanoparticles (LNPs), Buschmann, Paige and colleagues have devised methods for predicting and experimentally validating the pKa of LNPs based on the structure of the ionizable lipids used to formulate the LNPs. These studies may pave the way for the development of new LNP delivery vehicles that have reduced systemic distribution and improved endosomal release of their cargo post administration. To better understand protein-protein interactions and identify potential drug targets that disrupt such interactions, Luchini and colleagues have developed a methodology that identifies contact points between proteins using small molecule dyes. The dye molecules noncovalently bind to the accessible surfaces of a protein complex with very high affinity, but are excluded from contact regions. When the complex is denatured and digested with trypsin, the exposed regions covered by the dye do not get cleaved by the enzyme, whereas the contact points are digested. The resulting fragments can then be identified using mass spectrometry. The data generated can serve as the basis for designing small molecules and peptides that can disrupt the formation of protein complexes involved in disease processes. For example, using peptides based on the interleukin 1 receptor accessory protein (IL-1RAcP), Luchini, Liotta, Paige and colleagues disrupted the formation of IL-1/IL-R/IL-1RAcP complex and demonstrated that the inhibition of complex formation reduced the inflammatory response to IL-1B. Working on the discovery of novel antimicrobial agents, Bishop, van Hoek and colleagues have discovered a number of antimicrobial peptides from reptiles and other species. DRGN-1, is a synthetic peptide based on a histone H1-derived peptide that they had identified from Komodo Dragon plasma. DRGN-1 was shown to disrupt bacterial biofilms and promote wound healing in an animal model. The peptide, along with others, is being developed and tested in preclinical studies. Other research by van Hoek and colleagues focuses on in silico antimicrobial peptide discovery, screening of small molecules for antibacterial properties, as well as assessment of diffusible signal factors (DFS) as future therapeutics. The above examples provide insight into the cutting-edge studies undertaken by GMU scientists to develop novel methodologies and platform technologies important to drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗

Hallucination of closed repeat proteins containing central pockets

In pseudocyclic proteins, such as TIM barrels, β barrels, and some helical transmembrane channels, a single subunit is repeated in a cyclic pattern, giving rise to a central cavity that can serve as a pocket for ligand binding or enzymatic activity. Inspired by these proteins, we devised a deep-learning-based approach to broadly exploring the space of closed repeat proteins starting from only a specification of the repeat number and length. Biophysical data for 38 structurally diverse pseudocyclic designs produced in Escherichia coli are consistent with the design models, and the three crystal structures we were able to obtain are very close to the designed structures. Docking studies suggest the diversity of folds and central pockets provide effective starting points for designing small-molecule binders and enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Three-dimensional structure-guided evolution of a ribosome with tethered subunits

We report RNA-based macromolecular machines, such as the ribosome, have functional parts reliant on structural interactions spanning sequence-distant regions. These features limit evolutionary exploration of mutant libraries and confound three-dimensional structure-guided design. To address these challenges, we describe Evolink (evolution and linkage), a method that enables high-throughput evolution of sequence-distant regions in large macromolecular machines, and library design guided by computational RNA modeling to enable exploration of structurally stable designs. Using Evolink, we evolved a tethered ribosome with a 58% increased activity in orthogonal protein translation and a 97% improvement in doubling times in SQ171 cells compared to a previously developed tethered ribosome, and reveal new permissible sequences in a pair of ribosomal helices with previously explored biological function. The Evolink approach may enable enhanced engineering of macromolecular machines for new and improved functions for synthetic biology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Engineering sensitivity and specificity of AraC-based biosensors responsive to triacetic acid lactone and orsellinic acid

Abstract We previously described the design of triacetic acid lactone (TAL) biosensor ‘AraC-TAL1’, based on the AraC regulatory protein. Although useful as a tool to screen for enhanced TAL biosynthesis, this variant shows elevated background (leaky) expression, poor sensitivity and relaxed inducer specificity, including responsiveness to orsellinic acid (OA). More sensitive biosensors specific to either TAL or OA can aid in the study and engineering of polyketide synthases that produce these and similar compounds. In this work, we employed a TetA-based dual-selection to isolate new TAL-responsive AraC variants showing reduced background expression and improved TAL sensitivity. To improve TAL specificity, OA was included as a ‘decoy’ ligand during negative selection, resulting in the isolation of a TAL biosensor that is inhibited by OA. Finally, to engineer OA-specific AraC variants, the iterative protein redesign and optimization computational framework was employed, followed by 2 rounds of directed evolution, resulting in a biosensor with 24-fold improved OA/TAL specificity, relative to AraC-TAL1.

Wang, Zhiqing↗

Enhancing chemical bioproduction with rational control of bacterial post-translational modifications

Efficient conversion of inexpensive feedstocks to valuable chemicals by microbes is critical for a robust bioeconomy, but the ability to rationally design bacteria is hampered by insufficient knowledge of how post translational modifications (PTMs) control bacterial protein function and thus bioproduction phenotypes. Our study will focus on the lysine acetylation, a ubiquitous bacterial PTM that can affect the function of enzymes in central metabolism that are often critical for bioproduction processes, disrupt transcriptional regulation, and reduce translation. However, most lysine acetylation data is observational, which means that we do not know when, how, and what specific acetylated residues affect protein function and bacterial physiology. For our model host, we will use a Pseudomonas putida strain that we previously engineered to convert lignocellulosic feedstocks into chemicals such as itaconic acid (ITA). With this strain, we use a dynamic two-stage bioproduction process in which ITA is produced during a non-growth associated production phase. Production is highest during growth stages when lysine acetylation is low in other organisms (early stationary phase) and stalls in conditions where acetylation is highest (late stationary phase). The switch from high to stalled ITA production is also correlated with an unexpected increase in acetate levels – the precursor to non-enzymatic lysine acetylation. As such, we predict that lysine acetylation plays a substantial role in regulating the metabolic pathways required for ITA production. We will develop a generalizable approach that combines high-throughput genetic screens and cutting-edge genome engineering with state-of-the-art proteomics, metabolomics, and genetic code expansion methods to identify and modulate lysine acetylation patterns in bacteria. Ultimately, these strategies aim to manipulate protein expression and acetylation patterns to enhance bioproduction phenotypes (e.g., sustained ITA production in late stationary phase).

60 APPLIED LIFE SCIENCES↗

spammR: an R package designed for analysis and integration of spatial multi-omic measurements

Spatial omics is a young and evolving field and as such shows rapid development of novel technologies and analysis methods to measure transcripts, proteins, metabolites, and post-translational modifications at high spatial resolution. These advances in technology have enabled the simultaneous generation of abundance profiles for multiple different omics types and associated microscopy imaging data, as well as their analysis in a spatial context. However, most analytical tools are designed for spatial transcriptomics platforms and are challenging to use in other contexts such as mass spectrometry-based measurements or metagenomics. To this end we present spammR (spatial analysis of multi-omics measurements in R), an R package that enables end-to-end analysis with a specific focus on mass-spectrometry derived spatial omics datasets with (1) smaller sample sizes and spatial sparsity of samples, (2) considerable missingness, and (3) no a-priori knowledge about proteins or genes of interest, relying on a fully data-driven approach.

spammR↗

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗