Structural characterization of the POTRA domains from A. baumannii reveals new conformations in BamA
Not Available
Cottom, Claire Overly (ORCID:0000000349192114)↗
Engineering topics
Gumbart, James C.
Publications and source records attributed to Gumbart, James C..
Generative $β$-hairpin design using a residue-based physicochemical property landscape
De novo peptide design is a new frontier that has broad application potential in the biological and biomedical fields. Most existing models for de novo peptide design are largely based on sequence homology that can be restricted based on evolutionarily derived protein sequences and lack the physicochemical context essential in protein folding. Generative machine learning for de novo peptide design is a promising way to synthesize theoretical data that are based on, but unique from, the observable universe. In this study, we created and tested a custom peptide generative adversarial network intended to design peptide sequences that can fold into the -hairpin secondary structure. This deep neural network model is designed to establish a preliminary foundation of the generative approach based on physicochemical and conformational properties of 20 canonical amino acids, for example, hydrophobicity and residue volume, using extant structure-specific sequence data from the PDB. The beta generative adversarial network model robustly distinguishes secondary structures of hairpin from α helix and intrinsically disordered peptides with an accuracy of up to 96% and generates artificial -hairpin peptide sequences with minimum sequence identities around 31% and 50% when compared against the current NCBI PDB and nonredundant databases, respectively. These results highlight the potential of generative models specifically anchored by physicochemical and conformational property features of amino acids to expand the sequence-to-structure landscape of proteins beyond evolutionary limits.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Competition between inside-out unfolding and pathogenic aggregation in an amyloid-forming β-propeller
Studies of folded-to-misfolded transitions using model protein systems reveal a range of unfolding needed for exposure of amyloid-prone regions for subsequent fibrillization. Here, we probe the relationship between unfolding and aggregation for glaucoma-associated myocilin. Mutations within the olfactomedin domain of myocilin (OLF) cause a gain-of-function, namely cytotoxic intracellular aggregation, which hastens disease progression. Aggregation by wild-type OLF (OLF WT ) competes with its chemical unfolding, but only below the threshold where OLF loses tertiary structure. Representative moderate (OLF D380A ) and severe (OLF I499F ) disease variants aggregate differently, with rates comparable to OLF WT in initial stages of unfolding, and variants adopt distinct partially folded structures seen along the OLF WT urea-unfolding pathway. Whether initiated with mutation or chemical perturbation, unfolding propagates outward to the propeller surface. In sum, for this large protein prone to amyloid formation, the requirement for a conformational change to promote amyloid fibrillization leads to direct competition between unfolding and aggregation.
59 BASIC BIOLOGICAL SCIENCES↗
Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor
Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ~5-50 cm -1 . This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effect of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal/mol -1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Lpp positions peptidoglycan at the AcrA-TolC interface in the AcrAB-TolC multidrug efflux pump
The multidrug efflux pumps of Gram-negative bacteria are a class of complexes that span the periplasm, coupling both the inner and outer membranes to expel toxic molecules. The best-characterized example of these tripartite pumps is the AcrAB-TolC complex of Escherichia coli. However, how the complex interacts with the peptidoglycan (PG) cell wall, which is anchored to the outer membrane (OM) by Braun’s lipoprotein (Lpp), is still largely unknown. In this work, we present molecular dynamics simulations of a complete, atomistic model of the AcrAB-TolC complex with the inner membrane, OM, and PG layers all present. We find that the PG localizes to the junction of AcrA and TolC, in agreement with recent cryo-tomography data. Free-energy calculations reveal that the positioning of PG is determined by the length and conformation of multiple Lpp copies anchoring it to the OM. The distance between the PG and OM measured in cryo-electron microscopy images of wild-type E. coli also agrees with the simulation-derived spacing. Finally, sequence analysis of AcrA suggests a conserved role for interactions with PG in the assembly and stabilization of efflux pumps, one that may extend to other trans-envelope complexes as well.
59 BASIC BIOLOGICAL SCIENCES↗
Machine Learning Reveals the Critical Interactions for SARS-CoV-2 Spike Protein Binding to ACE2
SARS-CoV and SARS-CoV-2 bind to the human ACE2 receptor in practically identical conformations, although several residues of the receptor-binding domain (RBD) differ between them. Herein, we have used molecular dynamics (MD) simulations, machine learning (ML), and free-energy perturbation (FEP) calculations to elucidate the differences in binding by the two viruses. Although only subtle differences were observed from the initial MD simulations of the two RBD–ACE2 complexes, ML identified the individual residues with the most distinctive ACE2 interactions, many of which have been highlighted in previous experimental studies. FEP calculations quantified the corresponding differences in binding free energies to ACE2, and examination of MD trajectories provided structural explanations for these differences. Lastly, the energetics of emerging SARS-CoV-2 mutations were studied, showing that the affinity of the RBD for ACE2 is increased by N501Y and E484K mutations but is slightly decreased by K417N.
60 APPLIED LIFE SCIENCES↗
Tuning Proton Transfer Thermodynamics in SARS-CoV-2 Main Protease: Implications for Catalysis and Inhibitor Design
The catalytic reaction in SARS-CoV-2 main protease is activated by a proton transfer (PT) from Cys145 to His41. The same PT is likely also required for the covalent binding of some inhibitors. Here we use a multiscale computational approach to investigate the PT thermodynamics in the apo enzyme and in complex with two potent inhibitors, N3 and the α-ketoamide 13b. We show that with the inhibitors the free energy cost to reach the charge-separated state of the active-site dyad is lower, with N3 inducing the most significant reduction. We also show that a few key sites (including specific water molecules) significantly enhance or reduce the thermodynamic feasibility of the PT reaction, with selective desolvation of the active site playing a crucial role. The approach presented is a cost-effective procedure to identify the enzyme regions that control the activation of the catalytic reaction and is thus also useful to guide the design of inhibitors.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Bifunctional Janus Particles as Multivalent Synthetic Nanoparticle Antibodies (SNAbs) for Selective Depletion of Target Cells
Monoclonal antibodies (mAb) have had a transformative impact on treating cancers and immune disorders. SNAbs However, their use is limited by high development time and monetary cost, manufacturing complexities, suboptimal pharmacokinetics, and availability of disease-specific targets. To address some of these challenges, we developed an entirely synthetic, multivalent, Janus nanotherapeutic platform, called Synthetic Nanoparticle Antibodies (SNAbs). SNAbs, with phage-display-identified cell-targeting ligands on one "face" and Fc-mimicking ligands on the opposite "face", were synthesized using a custom, multistep, solid-phase chemistry method. SNAbs efficiently targeted and depleted myeloid-derived immune-suppressor cells (MDSCs) from mouse-tumor and rat-trauma models, ex vivo. Systemic injection of MDSC-targeting SNAbs efficiently depleted circulating MDSCs in a mouse triple-negative breast cancer model, enabling enhanced T cell and Natural Killer cell infiltration into tumors. Finally, our results demonstrate that SNAbs are a versatile and effective functional alternative to mAbs, with advantages of a plug- and-play, cell-free manufacturing process, and high-throughput screening (HTS)-enabled library of potential targeting ligands.
60 APPLIED LIFE SCIENCES↗
Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease
The main protease (M pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics.
59 BASIC BIOLOGICAL SCIENCES↗
Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19
In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.
60 APPLIED LIFE SCIENCES↗
β-Barrel proteins tether the outer membrane in many Gram-negative bacteria
Gram-negative bacteria have a cell envelope that comprises an outer membrane (OM), a peptidoglycan (PG) layer and an inner membrane (IM). The OM and PG are load-bearing, selectively permeable structures that are stabilized by cooperative interactions between IM and OM proteins. In Escherichia coli , Braun’s lipoprotein (Lpp) forms the only covalent tether between the OM and PG and is crucial for cell envelope stability; however, most other Gram-negative bacteria lack Lpp so it has been assumed that alternative mechanisms of OM stabilization are present. In this work we used a glycoproteomic analysis of PG to show that β-barrel OM proteins are covalently attached to PG in several Gram-negative species, including Coxiella burnetii, Agrobacterium tumefaciens and Legionella pneumophila . In C. burnetii , we found that four different types of covalent attachments occur between OM proteins and PG, with tethering of the β-barrel OM protein BbpA becoming most abundant in the stationary phase and tethering of the lipoprotein LimB similar throughout the cell cycle. Using a genetic approach, we demonstrate that the cell cycle-dependent tethering of BbpA is partly dependent on a developmentally regulated L,D-transpeptidase (Ldt). We use our findings to propose a model of Gram-negative cell envelope stabilization that includes cell cycle control and an expanded role for Ldts in covalently attaching surface proteins to PG.
59 BASIC BIOLOGICAL SCIENCES↗