Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein interaction”

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 1,531 records · Page 85

Identification of Small-Molecule Inhibitors of Fibroblast Growth Factor 23 Signaling via In Silico Hot Spot Prediction and Molecular Docking to α-Klotho

Fibroblast growth factor 23 (FGF23) is a therapeutic target for treating hereditary and acquired hypophosphatemic disorders, such as X-linked hypophosphatemic (XLH) rickets and tumor-induced osteomalacia (TIO), respectively. FGF23-induced hypophosphatemia is mediated by signaling through a ternary complex formed by FGF23, the FGF receptor (FGFR), and α-Klotho. Currently, disorders of excess FGF23 are treated with an FGF23-blocking antibody, burosumab. Small-molecule drugs that disrupt protein/protein interactions necessary for the ternary complex formation offer an alternative to disrupting FGF23 signaling. In this study, the FGF23:α-Klotho interface was targeted to identify small-molecule protein/protein interaction inhibitors since it was computationally predicted to have a large fraction of hot spots and two druggable residues on α-Klotho. We further identified Tyr433 on the KL1 domain of α-Klotho as a promising hot spot and α-Klotho as an appropriate drug-binding target at this interface. Subsequently, we performed in silico docking of ~5.5 million compounds from the ZINC database to the interface region of α-Klotho from the ternary crystal structure. Following docking, 24 and 20 compounds were in the final list based on the lowest binding free energies to α-Klotho and the largest number of contacts with Tyr433, respectively. Herein, five compounds were assessed experimentally by their FGF23-mediated extracellular signal-regulated kinase (ERK) activities in vitro, and two of these reduced activities significantly. Both these compounds were predicted to have favorable binding affinities to α-Klotho but not have a large number of contacts with the hot spot Tyr433. ZINC12409120 was found experimentally to disrupt FGF23:α-Klotho interaction to reduce FGF23-mediated ERK activities by 70% and have a half maximal inhibitory concentration (IC50) of 5.0 ± 0.23 μM. Molecular dynamics (MD) simulations of the ZINC12409120:α-Klotho complex starting from in silico docking poses reveal that the ligand exhibits contacts with residues on the KL1 domain, the KL1–KL2 linker, and the KL2 domain of α-Klotho simultaneously, thereby possibly disrupting the regular function of α-Klotho and impeding FGF23:α-Klotho interaction. ZINC12409120 is a candidate for lead optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ligand-induced transmembrane conformational coupling in monomeric $\mathrm{EGFR}$

Single pass cell surface receptors regulate cellular processes by transmitting ligand-encoded signals across the plasma membrane via changes to their extracellular and intracellular conformations. This transmembrane signaling is generally initiated by ligand binding to the receptors in their monomeric form. While subsequent receptor-receptor interactions are established as key aspects of transmembrane signaling, the contribution of monomeric receptors has been challenging to isolate due to the complexity and ligand-dependence of these interactions. By combining membrane nanodiscs produced with cell-free expression, single-molecule Förster Resonance Energy Transfer measurements, and molecular dynamics simulations, we report that ligand binding induces intracellular conformational changes within monomeric, full-length epidermal growth factor receptor (EGFR). Our observations establish the existence of extracellular/intracellular conformational coupling within a single receptor molecule. We implicate a series of electrostatic interactions in the conformational coupling and find the coupling is inhibited by targeted therapeutics and mutations that also inhibit phosphorylation in cells. Collectively, these results introduce a facile mechanism to link the extracellular and intracellular regions through the single transmembrane helix of monomeric EGFR, and raise the possibility that intramolecular transmembrane conformational changes upon ligand binding are common to single-pass membrane proteins.

59 BASIC BIOLOGICAL SCIENCES↗

Deactivation of SARS-CoV-2 by Boronic Acid-Functionalized Polymer

Boronic acid-modified polymers (BAMPs) can interact with glycoproteins and other glycosylated compounds through covalent binding of the boronic acid moieties to saccharide residues. As a first step toward evaluating the utility of BAMPs as SARS-CoV-2 antiviral agents, this COVID-19 rapid response LDRD was intended to examine the effect of BAMPs on SARS-CoV-2 spike glycoprotein and its subsequent binding with ACE2 receptor protein. Multiple different approaches were attempted in order to determine whether BAMPs based on poly(ethylene glycol) and poly(ethylenimine) bind the spike protein, but failed to produce a definitive answer. However, two different enzyme-linked immunosorbent assays clearly showed no discernable effect of boronic acid in inhibiting spike-ACE2 binding.

59 BASIC BIOLOGICAL SCIENCES↗

An active role of the amyloplasts and nuclei of root statocytes in graviperception

Three main phases are discerned in the gravitropic reaction: perception of a gravitational stimulus, its transduction, and fixation of the reaction resulting in bending of an organ. According to the starch-statolith hypothesis of Nemec and Haberlandt, amyloplasts in the structurally and functionally specialized graviperceptive cells (statocytes) sediment in the direction of a gravitational vector in the distal part of a cell while a nucleus is in the proximal one. If amyloplasts appear to act as gravity sensors, the receptors, which interact with sedimented amyloplasts, and next signaling are still unclear. An analysis of the structural-functional organization of cells in different root cap layers of such higher plants as pea, Arabidopsis thaliana, and Brassica rapa grown under 1 g, on the clinostats, and in microgravity, allows us to support the hypothesis that amyloplasts function as statoliths in statocytes, but they may not be only the passive statolithic mass. We propose that amyloplasts fulfill a more complex function by interacting with a receptor, which is a nucleus, in transduction of some signal to it. Gravity-induced statolith movement in certain order leads to a new functional connection between gravity susceptors--amyloplasts and a receptor--a nucleus receiving some signal presumedly of a mechanical or biochemical nature from the amyloplasts. During gravitropism, sugar signaling could induce expression of genes encoding auxin transport proteins in a nucleus giving the nucleus an intermediate role in signal trunsduction following perception. c 2001 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

short duration↗

GeneLab for High Schools: Data Mining for the Next Generation

Modern biological sciences have become increasingly based on molecular biology and high-throughput molecular techniques, such as genomics, transcriptomics, and proteomics. NASA Scientists and the NASA Space Biology Program have aimed to examine the fundamental building blocks of life (RNA, DNA and protein) in order to understand the response of living organisms to space and aid in fundamental research discoveries on Earth. In an effort to enable NASA funded science to be available to everyone, NASA has collected the data from omics studies and curated them in a data system called GeneLab. Whilst most college-level interns, academics and other scientists have had some interaction with omics data sets and analysis tools, high school students often have not. Therefore, the Space Biology Program is implementing a new Summer Program for high-school students that aims to inspire the next generation of scientists to learn about and get involved in space research using GeneLabs Data System. The program consists of three main components core learning modules, focused on developing students knowledge on the Space Biology Program and Space Biology research, Genelab and the data system, and previous research conducted on model organisms in space; networking and team work, enabling students to interact with guest lecturers from local universities and their fellow peers, and also enabling them to visit local universities and genomics centers around the Bay area; and finally an independent learning project, whereby students will be required to form small groups, analyze a dataset on the Genelab platform, generate a hypothesis and develop a research plan to test their hypothesis. This program will not only help inspire high-school students to become involved in space-based research but will also help them develop key critical thinking and bioinformatics skills required for most college degrees and furthermore, will enable them to establish networks with their peers and connections with university Professors that may help them achieve their educational goals.

genelab↗

De Novo Design of High‐Affinity Miniprotein Binders Targeting Francisella Tularensis Virulence Factor

Abstract Francisella tularensis poses considerable public health risk due to its high infectivity and potential for bioterrorism. Francisella‐like lipoprotein (Flpp3), a key virulence factor unique to Francisella, plays critical roles in infection and immune evasion, making it a promising target for therapeutic development. However, the lack of well‐defined binding pockets and structural information on native interactions has hindered structure‐guided ligand discovery against Flpp3. Here, we used a combination of physics‐based and deep‐learning methods to design high‐affinity miniprotein binders targeting two distinct sites on Flpp3. We identified four binders for site I with binding affinities ranging between 24–110 nM. For the second site, an initial binder showed a dissociation constant ( K D ) of 81 nM, and subsequent site saturation mutagenesis yielded variants with sub‐nanomolar affinities. Circular dichroism confirmed the topology of designed miniproteins. The X‐ray crystal structure of Flpp3 in complex with a site I binder is nearly identical to the design model (Cα root‐mean‐square deviation (RMSD): 0.9 Å). These designed miniproteins provide research tools to explore the roles of Flpp3 in tularemia and should enable the development of new therapeutic candidates.

Gokce‐Alpkilic, Gizem [Molecular Engineering and S↗

Preparation and Characterization of Fluorescent Derivatives of Chicken Egg White Lysozyme

Fluorescence is one of the most versatile and powerful tools for the study of macromolecules. While most proteins are intrinsically fluorescent, working at crystallization concentrations require the use of covalently prepared derivatives added as tracers. This approach requires derivatives that do not markedly affect the crystal packing. We have prepared a number of fluorescent derivatives of chicken egg white lysozyme with probes bound to one of two different sites on the protein molecule. Lucifer yellow, cascade blue, and 5-(2-aminoethyl)aminonapthalene-l-sulfonic acid (EDANS) have been attached to the side chain carboxyl of asp101 using a carbodiimide coupling procedure. asp101 lies within the active site cleft, and it is believed that the probes are at least partially "buried" within that cleft. Lucifer yellow and MANS probes with iodoacetamide reactive groups have been bound to hisl5, located on the "back side" of the molecule relative to the active site. The fluorescently labeled protein is readily purified from the starting material by cation exchange chromatography. All the derivatives fluoresce in both the solution and the crystalline states. Fluorescence characterization has focused on determining the bound probe quantum yields, lifetimes, absorption and emission spectra, and quenching by added solutes in comparison to the free probe. No appreciable changes are found in the lifetimes of any of the probes except for cascade blue, where Tau(sub free) = 3.52 ns vrs Tau(sub bound) = 2.8 ns. Spectral shifts are found in most cases. Particularly strong quenching upon binding is found in the case of the cascade blue derivative, likely due to probe interactions with the active site cleft. While none of the asp101 bound probes are well quenched by commonly employed solutes, such as potassium and sodium iodide, acrylamide, primuline, the chloride salts of manganese, cesium, and cobalt, trifluoroacetamide, trichloroethanol, and thallium iodide, in those cases where quenching is observed the bound probe is less efficiently quenched relative to the free probe. This indicates that the bound probes are less accessible to the bulk solution, an expected finding for attachment within the active site cleft. Attempts have been made to bind other molecules to these sites, with varying success. Interestingly, all three probes contain one or more sulfonate ((Ar-S03)-) groups. We have not been successful in binding analogous probes without sulfate groups such as pyrene, or with derivatized sulfonate groups such as dansyl type probes, analogous to MANS but where the sulfonate group is derivatized, Ar-S02-N2C2H7. None of the probes is rigidly bound to the protein, i.e., they all have a probe motion superimposed on that of the protein.

Sumida, John↗

Mutations of Omicron Variant at the Interface of the Receptor Domain Motif and Human Angiotensin-Converting Enzyme-2

The most recent Omicron variant of SARS-CoV-2 has caused global concern and anxiety. The only thing certain about this strain, with a large number of mutations in the spike protein, is that it spreads quickly, seems to evade immune defense, and mitigates the benefits of existing vaccines. Based on the ultra-large-scale ab initio computational modeling of the receptor binding motif (RBM) and the human angiotensin-converting enzyme-2 (ACE2) interface, we provide the details of the effect of Omicron mutations at the fundamental atomic scale level. In-depth analysis anchored in the novel concept of amino acid-amino acid bond pair units (AABPU) indicates that mutations in the Omicron variant are connected with (i) significant changes in the shape and structure of AABPU components, together with (ii) significant increase in the positive partial charge, which facilitates the interaction with ACE2. We have identified changes in bonding due to mutations in the RBM. The calculated bond order, based on AABPU, reveals that the Omicron mutations increase the binding strength of RBM to ACE2. Our findings correlate with and are instrumental to explain the current observations and can contribute to the prediction of next potential new variant of concern.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

EPR studies of ferredoxin in spinach and cyanobacterial thylakoids related to photosystem I-driven NADP + reduction

Photosynthetic light-dependent reactions occur in thylakoid membranes where embedded proteins capture light energy and convert it to chemical energy in the form of ATP and NADPH for use in carbon fixation. One of these integral membrane proteins is Photosystem I (PSI). PSI catalyzes light-driven transmembrane electron transfer from plastocyanin (Pc) to oxidized ferredoxin (Fd). Electrons from reduced Fd are used by the enzyme ferredoxin-NADP + reductase (FNR) for the reduction of NADP + to NADPH. Fd and Pc are both small soluble proteins whereas the larger FNR enzyme is associated with the membrane. To investigate electron shuttling between these diffusible and embedded proteins, thylakoid photoreduction of NADP + was studied. As isolated, both spinach and cyanobacterial thylakoids generate NADPH upon illumination without extraneous addition of Fd. These findings indicate that isolated thylakoids either (i) retain a “pool” of Fd which diffuses between PSI and membrane bound FNR or (ii) that a fraction of PSI is associated with Fd, with the membrane environment facilitating PSI-Fd-FNR interactions which enable multiple turnovers of the complex with a single Fd. To explore the functional association of Fd with PSI in thylakoids, electron paramagnetic resonance (EPR) spectroscopic methodologies were developed to distinguish the signals for the reduced Fe-S clusters of PSI and Fd. Temperature-dependent EPR studies show that the EPR signals of the terminal [4Fe-4S] cluster of PSI can be distinguished from the [2Fe-2S] cluster of Fd at > 30 K. At 50 K, the cw X-band EPR spectra of cyanobacterial and spinach thylakoids reduced with dithionite exhibit EPR signals of a [2Fe-2S] cluster with g-values g x = 2.05, g y = 1.96, and g z = 1.89, confirming that Fd is present in thylakoid preparations capable of NADP + photoreduction. Quantitation of the EPR signals of P 700 + and dithionite reduced Fd reveal that Fd is present at a ratio of ~ 1 Fd per PSI monomer in both spinach and cyanobacterial thylakoids. Light-driven electron transfer from PSI to Fd in thylakoids confirms Fd is functionally associated (< 0.4 Fd/PSI) with the acceptor end of PSI in isolated cyanobacterial thylakoids. Furthermore, these EPR experiments provide a benchmark for future spectroscopic characterization of Fd interactions involved in multistep relay of electrons following PSI charge separation in the context of photosynthetic thylakoid microenvironments.

Electron paramagnetic resonance↗

An in vitro BRAF activation assay elucidates molecular mechanisms driving disassembly of the autoinhibited BRAF state

The RAF kinases (ARAF, BRAF, and CRAF) are essential components of the RAS-ERK signaling pathway, which controls vital cellular processes and is frequently dysregulated in human disease. Notably, mutations that alter BRAF function are prominent drivers of human cancer and certain RASopathy disorders, making BRAF an important target for therapeutic intervention. Despite extensive research, several aspects of BRAF regulation remain unclear. In this study, we developed an in vitro BRAF activation assay using purified autoinhibited BRAF:14-3-3 2 :MEK complexes. Our results show that fully processed, active-state KRAS alone can promote dimer-dependent BRAF activation. Moreover, we found that phosphatidylserine (PS)-containing liposomes synergized with KRAS to promote BRAF activation, achieving activity levels comparable to those observed with BRAF proteins that constitutively dimerize. In contrast, the SMP phosphatase complex had only a minimal effect on BRAF catalytic activity in this system but mediated the dephosphorylation of the negative regulatory pS365 14-3-3 binding site in a manner that was accelerated by the presence of KRAS alone or KRAS and 30% PS liposomes. Finally, we show that inhibitors blocking the BRAF RBD:KRAS interaction were able to suppress the in vitro activation of BRAF, underscoring the critical role of RAS binding in initiating the disassembly of the BRAF autoinhibited state. Thus, this assay provides valuable insights into the steps required for BRAF activation and can serve as an effective screening tool for identifying compounds that may inhibit this process and have therapeutic potential.

BRAF↗

Structure of the rabies virus glycoprotein trimer bound to a prefusion-specific neutralizing antibody

Rabies infection is nearly 100% lethal if untreated and kills more than 50,000 people annually, many of them children. Existing rabies vaccines target the rabies virus glycoprotein (RABV-G) but generate short-lived immune responses, likely because the protein is heterogeneous under physiological conditions. Here, we report the 3.39 Å cryo–electron microscopy structure of trimeric, prefusion RABV-G complexed with RVA122, a potently neutralizing human antibody. RVA122 binds to a quaternary epitope at the top of RABV-G, bridging domains and stabilizing RABV-G protomers in a prefusion state. RABV-G trimerization involves side-to-side interactions between the central α helix and adjacent loops, rather than contacts between central helices, and interactions among the fusion loops at the glycoprotein base. These results provide a basis from which to develop improved rabies vaccines based on RABV-G stabilized in the prefusion conformation.

59 BASIC BIOLOGICAL SCIENCES↗

Broad betacoronavirus neutralization by a stem helix–specific human antibody

Targeting a range of betacoranaviruses In the past 20 years, three highly pathogenic β-coronaviruses have crossed from animals to humans, including the most recent: severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). A spike protein that decorates these viruses has an S1 domain that binds host cell receptors and an S2 domain that fuses the viral and cell membranes to allow cell entry. The S1 domain is the target of many neutralizing antibodies but is more genetically variable than S2, and antibodies can exert selective pressure, leading to resistant variants. Pinto et al . identified five monoclonal antibodies that interact with a helix in the S2 domain. The most broadly neutralizing antibody inhibited all β-coronavirus subgenera and reduced viral burden in hamsters infected with SARS-CoV-2. —VV

59 BASIC BIOLOGICAL SCIENCES↗

Measurement report: The Fifth International Workshop on Ice Nucleation phase 1 (FIN-01): intercomparison of single-particle mass spectrometers

Knowledge of the chemical composition and mixing state of aerosols at a single-particle level is critical for gaining insights into atmospheric processes. One common tool to make these measurements is single-particle mass spectrometry. There remains a need to compare the performance of different single-particle mass spectrometers (SPMSs). An intercomparison of SPMSs was conducted at the Aerosol Interaction and Dynamics in the Atmosphere (AIDA) chamber at the Karlsruhe Institute of Technology (KIT) in November 2014, as part of the first phase of the Fifth International Workshop on Ice Nucleation (FIN-01). In this paper we compare size distributions and mass spectra of atmospherically relevant particle types measured by five SPMSs. These include different minerals, desert and soil dusts, soot, bioaerosol (Snomax; protein granule), secondary organic aerosol (SOA), and SOA-coated mineral particles. Most SPMSs reported similar vacuum aerodynamic diameter (d va ) within typical instrumental ranges from ~100–200 nm (lower limit) to ~2–3 µm (upper limit). In general, all SPMSs exhibited a wide dynamic range (up to ~10 3 ) and high signal-to-noise ratio (up to ~10 4 ) in mass spectra. Common spectral features with small diversities in mass spectra were found with high average Pearson's correlation coefficients, i.e., for average positive spectra r avg-pos =0.74 ± 0.12 and average negative spectra r avg-neg =0.67 ± 0.22. We found that instrument-specific detection efficiency (DE) was more dependent on particle size than particle type, and particle identification favored the use of bipolar, rather than monopolar, instruments. Particle classification from “blind experiments” showed that all instruments differentiated SOA, soot, and soil dust and detected subtle changes in the particle internal mixing but had difficulties differentiating among specific mineral types and dusts. This study helps to further understand the capabilities and limitations of the single-particle mass spectrometry technique in general and the specific performance of the instrument in characterizing atmospheric aerosol particles.

54 ENVIRONMENTAL SCIENCES↗

Photoinduced Electron Transfer Informs on Pathway Coupling in Flavin-Based Electron Bifurcation

Flavin-based electron bifurcation (FBEB) is an enzymatic mechanism that generates extremely high-energy electrons to drive unfavorable chemical reactions. It is utilized by the NADHdependent ferredoxin:NADP + -oxidoreductase (Nfn) enzyme in hyperthermophile Pyrococcus furiosus to bifurcate electrons from NADPH into the coupled low-potential (endergonic) and highpotential (exergonic) pathways. This process enables P. furiosus to live in harsh and uninviting environments. Despite its biological importance, the mechanisms used by Nfn to facilitate exceptional directional control over short-lived, high-energy electrons and to prevent undesired transfer, particularly along the low-potential pathway, are still not well understood. To elucidate how the protein environment contributes to electronic control in the lowpotential pathway, new techniques must be utilized to probe these unstable intermediates. In this study, we have adapted lowtemperature photoexcitation combined with electron paramagnetic resonance (EPR) to accumulate the short-lived intermediate and place it in the context of the other cofactors involved in the low-potential pathway of Nfn. We observed coincident growth of both the radical intermediate and its nearby [4Fe-4S] cluster over 4.5 h of illumination with NADPH at cryogenic temperatures. The photogenerated paramagnetic species were stable in LN 2 storage indefinitely and recombined when warmed to higher temperatures. The results provide insights into the electron transfer steps and cofactor interactions along the low potential pathway, facilitating a more robust mechanistic understanding of the high-energy events of electron bifurcation. Furthermore, through comparison of cryogenic and room temperature experiments, a potential gating step involving the movement of key residues important for the reversibility of electron flow along this pathway is suggested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure-Function and Nearest Neighbor Interactions of the Cytochrome b 6 f Complex

Test prediction from our (i) high resolution (2.5 Å) structure of cyanobacterial cytochrome b 6 f complex (Structure, 22, 1008-1015, 2014), and (ii) analysis of quinone/quinol traffic within the b 6 f complex (Biophys. J., 107, 1620-1628, 2014) concerning the rate-limiting step in the linear (non-cyclic) electron transport chain. It was concluded that it arises from access of plastoquinol, the electron/proton donor, to the high potential electron carriers, i. e., the Fe 2 S 2 and cytochrome f proteins. The directed motion of the quinol proceeds via a narrow channel on the lumen or electrochemically positive, p- (lumen) side of the b 6 f complex which was defined by crystal structures of the b 6 f complex with quinol analogue inhibitors (J. Mol. Biol., 370, 39-52, 2007).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A conserved epitope III on hepatitis C virus E2 protein has alternate conformations facilitating cell binding or virus neutralization

Epitope III, a highly conserved amino acid motif of 524 APTYSW 529 on the hepatitis C virus (HCV) E2 glycoprotein, resides in the critical loop that binds to the host receptor CD81, thus making it one of the most important antibody targets for blocking HCV infections. Here, we have determined the X-ray crystal structure of epitope III at a 2.0-Å resolution when it was captured by a site-specific neutralizing antibody, monoclonal antibody 1H8 (mAb1H8). The snapshot of this complex revealed that epitope III has a relatively rigid structure when confined in the binding grooves of mAb1H8, which confers the residue specificity at both ends of the epitope. Such a high shape complementarity is reminiscent of the “lock and key” mode of action, which is reinforced by the incompatibility of an antibody binding with an epitope bearing specific mutations. By subtly positioning the side chains on the three residues of Tyr 527 , Ser 528 , and Trp 529 while preserving the spatial rigidity of the rest, epitope III in this cocrystal complex adopts a unique conformation that is different from previously described E2 structures. With further analyses of molecular docking and phage display–based peptide interactions, we recognized that it is the arrangements of two separate sets of residues within epitope III that create these discrete conformations for the epitope to interact selectively with either mAb1H8 or CD81. These observations thus raise the possibility that local epitope III conformational dynamics, in conjunction with sequence variations, may act as a regulatory mechanism to coordinate “mAb1H8-like” antibody-mediated immune defenses with CD81-initiated HCV infections.

59 BASIC BIOLOGICAL SCIENCES↗

Trichoderma harzianum transcriptome in response to the nematode Pratylenchus brachyurus

The root-lesion nematode Pratylenchus brachyurus causes extensive damage in several crops of economic importance. Fungi of Trichoderma genus have been highlighted as biopesticide agents in the control of several plant diseases. Although it is already widely used in agriculture, there are few studies, especially at the molecular level, that evaluate T. harzianum in the control of P. brachyurus . The aim of the present study was to investigate how interaction with the nematode P. brachyurus influences gene expression of T. harzianum , by using RNA-Seq analysis. Of the 13,932 predicted genes in the T. harzianum genome, 2,922 (21%) were differentially expressed in the presence of P. brachyurus , in relation to the absence of the nematode. Among the differentially expressed genes, we found genes encoding Carbohydrate Active EnZymes (CAZy), MEROPS peptidases, and proteins related to secondary metabolite synthesis. 118 pathways were identified as related to the biosynthesis of secondary metabolites. Additionally, the analysis identified 136 metabolic pathways related to these differentially expressed genes, among which we highlight: aminobenzoate degradation, xenobiotic metabolism by cytochrome P450, and sesquiterpenoid and triterpenoid biosynthesis. Our results contribute to a better understanding of the response of T. harzianum to the nematode P. bracyurus , the potential for biocontrol by this fungus.

59 BASIC BIOLOGICAL SCIENCES↗

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single–amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single– and multiple–amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.

60 APPLIED LIFE SCIENCES↗