Search

Previously we established a family of macrocyclic peptide triazoles (cPTs) that inactivate the Env protein complex of HIV-1, and identified the pharmacophore that engages Env’s receptor binding pocket. Here, we examined the hypothesis that the side chains of both components of the triazole Pro - Trp segment of cPT pharmacophore work in tandem to make intimate contacts with two proximal subsites of the overall CD4 binding site of gp120 to stabilize binding and function. Variations of the triazole Pro R group, which previously had been significantly optimized, led to identification of a variant MG-II-20 that contains a pyrazole substitution. MG-II-20 has improved functional properties over previously examined variants, with Kd for gp120 in the nM range. In contrast, new variants of the Trp indole side chain, with either methyl- or bromo- components appended, had disruptive effects on gp120 binding, reflecting the sensitivity of function to changes in this component of the encounter complex. Plausible in silico models of cPT:gp120 complex structures were obtained that are consistent with the overall hypothesisof occupancy by the triazole Pro and Trp side chains, respectively, into the β20/21 and Phe43 sub-cavities. The overall results strengthen the definition of the cPT-Env inactivator binding site and provide a new lead composition (MG-II-20) as well as structure-function findings to guide future HIV-1 Env inactivator design.

Screening already approved drugs for activity against a novel pathogen can be an important part of global rapid-response strategies in pandemics. Such high-throughput repurposing screens have already identifed several existing drugs with potential to combat SARS-CoV-2. However, moving these hits forward for possible development into drugs specifcally against this pathogen requires unambiguous identifcation of their corresponding targets, something the high-throughput screens are not typically designed to reveal. We present here a new computational inverse-docking protocol that uses all-atom protein structures and a combination of docking methods to rank-order targets for each of several existing drugs for which a plurality of recent high-throughput screens detected anti-SARS-CoV-2 activity. We demonstrate validation of this method with known drug-target pairs, including both non-antiviral and antiviral compounds. We subjected 152 distinct drugs potentially suitable for repurposing to the inverse docking procedure. The most common preferential targets were the human enzymes TMPRSS2 and PIKfyve, followed by the viral enzymes Helicase and PLpro. All compounds that selected TMPRSS2 are known serine protease inhibitors, and those that selected PIKfyve are known tyrosine kinase inhibitors. Detailed structural analysis of the docking poses revealed important insights into why these selections arose, and could potentially lead to more rational design of new drugs against these targets. Fil: Ribone, Sergio Roman. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Unidad de Investigación y Desarrollo en Tecnología Farmacéutica. Universidad Nacional de Córdoba. Facultad de Ciencias Químicas. Unidad de Investigación y Desarrollo en Tecnología Farmacéutica; Argentina Fil: Paz, Sergio Alexis. Universidad Nacional de Córdoba. Facultad de Ciencias Químicas. Departamento de Química Teórica y Computacional; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigaciones en Físico-química de Córdoba. Universidad Nacional de Córdoba. Facultad de Ciencias Químicas. Instituto de Investigaciones en Físico-química de Córdoba; Argentina Fil: Abrams, Cameron F.. Drexel University; Estados Unidos Fil: Villarreal, Marcos Ariel. Universidad Nacional de Córdoba. Facultad de Ciencias Químicas. Departamento de Química Teórica y Computacional; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigaciones en Físico-química de Córdoba. Universidad Nacional de Córdoba. Facultad de Ciencias Químicas. Instituto de Investigaciones en Físico-química de Córdoba; Argentina

[Image: see text] The design and synthesis of butyl chain derivatives at the indane ring 3-position of our lead CD4-mimetic compound BNM-III-170 that inhibits human immunodeficiency virus (HIV-1) infection are reported. Optimization efforts were guided by crystallographic and computational analysis of the small-molecule ligands of the Phe43 cavity of the envelope glycoprotein gp120. Biological evaluation of 11–21 revealed that members of this series of CD4-mimetic compounds are able to inhibit HIV-1 viral entry into target cells more potently and with greater breadth compared to BNM-III-170. Crystallographic analysis of the binding pocket of 14, 16, and 17 revealed a novel hydrogen bonding interaction between His105 and a primary hydroxyl group on the butyl side chain. Further optimization of this interaction with the His105 residue holds the promise of more potent CD4-mimetic compounds.

Properties of epoxy thermosets can be varied broadly to suit design requirements by altering the chemistry of the component agents. Atomistically-detailed molecular dynamics simulations are well-suited for molecular insight into the structure-property relationship for a rational tailoring of the chemistry. Since the macroscopic properties of interest for applications emerge hierarchically from molecular-scale chemical interactions, seamless integration of experiment, computation, and theory is of great interest. Recently, a Specific Volume--Cooling Rate analysis protocol was successfully developed to quantitatively compare the volumetric properties of an epoxy network model with experimental results in the literature, in spite of the nine orders of magnitude mismatch in the accessible time-scales. Here, we extend the application of the method for two epoxy networks in the same class of chemistry but whose monomers have a higher number of repeating units compared to the previous one for validating the generality of our approach. We observed that atomistic simulations are able to predict the experimental temperature trend of the specific volume within 0.4% for both these networks. Using the William-Landel-Ferry equation to account for rate differences, we also see good agreement between the computational and experimental values of the glass transition temperature.

We previously discovered a class of recombinant lectin conjugates, denoted lectin DLIs (‘dual-acting lytic inhibitors’) that bind to the HIV-1 envelope (Env) protein trimer and cause both lytic inactivation of HIV-1 virions and cytotoxicity of Env-expressing cells. To facilitate mechanistic investigation of DLI function, we derived the simplified prototype microvirin (MVN)-DLI, containing an MVN domain that binds high-mannose glycans in Env, connected to a DKWASLWNW sequence (denoted ‘Trp3’) derived from the membrane-associated region of gp41. The relatively much stronger affinity of the lectin component than Trp3 argues that the lectin functions to capture Env to enable Trp3 engagement and consequent Env membrane disruption and virolysis. The relatively simplified engagement pattern of MVN with Env opened up the opportunity, pursued here, to use recombinant glycan knockout gp120 variants to identify the precise Env binding site for MVN that drives DLI engagement and lysis. Using mutagenesis combined with a series of biophysical and virological experiments, we identified a restricted set of residues, N262, N332 and N448, all localized in a cluster on the outer domain of gp120, as the essential epitope for MVN binding. By generating these mutations in the corresponding HIV-1 virus, we established that the engagement of this glycan cluster with the lectin domain of MVN*-DLI is the trigger for DLI-derived virus and cell inactivation. Beyond defining the initial encounter step for lytic inactivation, this study provides a guide to further elucidate DLI mechanism, including the stoichiometry of Env trimer required for function, and downstream DLI optimization.

Sizing emulsions used as glass fiber surface treatments in composites manufacturing are aqueous suspensions of hydrophobic film formers, surface coupling agents, and surfactants. We employ all-atom molecular dynamics simulations to characterize droplet structures in several aqueous blends of the film-former diglycidyl ether of bisphenol A, coupling agent glycidoxypropyl trimethoxysilane, and a triblock copolymer surfactant (Pluronic L35 PEO/PPO copolymer). We show that the quasi-equilibrium states of emulsion droplets are invariant to different initial configurations. We examine the role of the surfactant in determining coupling agent partitioning between the droplet shell and corona and coupling agent cluster size distributions. This work takes a step toward systematic understanding of the sizing chemistry to optimize the interface between the glass and the resin in commercially relevant composites.

A modified single sweep approach is proposed for generating free energy landscapes. The approach replaces the use of temperature-accelerated molecular dynamics (TAMD) to generate centres in collect...

[Image: see text] With approximately 37 million people living with HIV worldwide and an estimated 2 million new infections reported each year, the need to derive novel strategies aimed at eradicating HIV-1 infection remains a critical worldwide challenge. One potential strategy would involve eliminating infected cells via antibody-dependent cellular cytotoxicity (ADCC). HIV-1 has evolved sophisticated mechanisms to conceal epitopes located in its envelope glycoprotein (Env) that are recognized by ADCC-mediating antibodies present in sera from HIV-1 infected individuals. Our aim is to circumvent this evasion via the development of small molecules that expose relevant anti-Env epitopes and sensitize HIV-1 infected cells to ADCC. Rapid elaboration of an initial screening hit using parallel synthesis and structure-based optimization has led to the development of potent small molecules that elicit this humoral response. Efforts to increase the ADCC activity of this class of small molecules with the aim of increasing their therapeutic potential was based on our recent cocrystal structures with gp120 core.

The HIV-1 envelope glycoprotein (Env) trimer mediates cell entry and is conformationally dynamic1-8. Imaging by single-molecule fluorescence resonance energy transfer (smFRET) has revealed that, on the surface of intact virions, mature pre-fusion Env transitions from a pre-triggered conformation (state 1) through a default intermediate conformation (state 2) to a conformation in which it is bound to three CD4 receptor molecules (state 3)8-10. It is currently unclear how these states relate to known structures. Breakthroughs in the structural characterization of the HIV-1 Env trimer have previously been achieved by generating soluble and proteolytically cleaved trimers of gp140 Env that are stabilized by a disulfide bond, an isoleucine-to-proline substitution at residue 559 and a truncation at residue 664 (SOSIP.664 trimers)5,11-18. Cryo-electron microscopy studies have been performed with C-terminally truncated Env of the HIV-1JR-FL strain in complex with the antibody PGT15119. Both approaches have revealed similar structures for Env. Although these structures have been presumed to represent the pre-triggered state 1 of HIV-1 Env, this hypothesis has never directly been tested. Here we use smFRET to compare the conformational states of Env trimers used for structural studies with native Env on intact virus. We find that the constructs upon which extant high-resolution structures are based predominantly occupy downstream conformations that represent states 2 and 3. Therefore, the structure of the pre-triggered state-1 conformation of viral Env that has been identified by smFRET and that is preferentially stabilized by many broadly neutralizing antibodies-and thus of interest for the design of immunogens-remains unknown.

The effect of alkyl chain in the amidoamine crosslinker on the yield and post-yield behavior using non-equilibrium molecular dynamics simulation is studied in this work. Specifically the deformation behavior of two thermosets with and without the alkyl chain are simulated. The yield stress and strain hardening modulus were predicted. It was found that the thermoset with alkyl chain had lower yield stress and strain hardening modulus compared to the one without. Stress partitioning revealed that intermolecular Lennard-Jones (LJ) and covalent interactions contribute significantly to the yield stress and strain hardening modulus in these thermosets while the electrostatic interactions have no significant effect on both. The effect of methylene interactions on the yield and post-yield behavior were also investigated. It was found that the lower yield stress in the thermoset with alkyl chain is due to the dilution effect and the decrease in strain hardening modulus is due to both dilution effect and intramolecular LJ interactions of methylenes. Simulations reveal a trade off between dilution effect and covalent interactions of the system which can be utilized to tune the mechanical properties of thermosets. The major contribution of this work is the computational approach that can be used to discover the effect of additives or other molecular constituents on the deformation characteristics of a thermoset regardless of the deformation mode tensile or compressive.

The seasonal nature in the outbreaks of respiratory viral infections with increased transmission during low temperatures has been well established. The current COVID-19 pandemic makes no exception, and temperature has been suggested to play a role on the viability and transmissibility of SARS-CoV-2. The receptor binding domain (RBD) of the Spike glycoprotein binds to the angiotensin-converting enzyme 2 (ACE2) to initiate viral fusion. Studying the effect of temperature on the receptor-Spike interaction, we observed a significant and stepwise increase in RBD-ACE2 affinity at low temperatures, resulting in slower dissociation kinetics. This translated into enhanced interaction of the full Spike to ACE2 receptor and higher viral attachment at low temperatures. Interestingly, the RBD N501Y mutation, present in emerging variants of concern (VOCs) that are fueling the pandemic worldwide, bypassed this requirement. This data suggests that the acquisition of N501Y reflects an adaptation to warmer climates, a hypothesis that remains to be tested.

Obtaining accurate and reproducible free energies from molecular simulations is somewhat tricky due to incomplete knowledge of crucial slow degrees of freedom leading to hidden barriers that can stymie sampling. Employing a sufficiently large number of collective variables (CV) and ensuring ergodic sampling in orthogonal CV space, perhaps via tempering methods, can reduce these issues to some extent. For complex systems with high-dimensional free energy landscapes, both these approaches become computationally expensive. For high-dimensional landscapes, efficient exploration can be enabled by using temperature-accelerated MD (TAMD) and identification and characterization of minimum free energy pathways connecting minima can be found by using the string method (SM). Both TAMD and SM use mean-force estimates from finite MD simulations and are thus susceptible to sampling restrictions from hidden variables. A recent development in parallel tempering methods, "generalized replica exchange solute tempering" (gREST), can enhance sampling at a reasonable computational cost with its flexibility to target very specific "solutes" which can include arbitrary independent variables. Considering the advantages of both methods, we implement gREST-enabled TAMD and SM. By considering two different collective variable representations of the pentapeptide neurotransmitter met-enkephalin, we show that both gREST-enabled TAMD and SM yield more accurate and reproducible free energy predictions than TAMD and SM alone. Given the moderate computational cost of gREST compared with other replica-exchange methods, gREST-enabled SM represents a more attractive method for characterizing free energy minima and pathways among them for a large variety of systems.

A strategy has been established for the synthesis of a family of bifunctional HIV-1 inhibitor covalent conjugates with the potential to bind simultaneously to both the gp120 and gp41 subunits of the HIV-1 envelope glycoprotein trimeric complex (Env). One component of the conjugates is derived from BNM-III-170, a small-molecule CD4 mimic that binds to gp120. The second component, comprised of the peptide DKWASLWNW ("Trp3"), was derived from the N-terminus of the HIV-1 gp41 Membrane Proximal External Region (MPER) and found previously to bind to the gp41 subunit of Env. The resulting bifunctional conjugates were shown to inhibit virus cell infection with low micromolar potency and to induce lysis of the HIV-1 virion. Crucially, virolysis was found to be dependent on the covalent linkage of the BNM-III-170 and Trp3 domains, as coadministration of a mixture of the un-cross-linked components proved to be nonlytic. However, a significant magnitude of lytic activity was observed in Env-negative and other control pseudoviruses, suggesting parallel mechanisms of action of the conjugates involving Env interaction and direct membrane disruption. Computational modeling suggested strong membrane-binding activity of BNM-III-170, which may underly the nonspecific virolytic effects of the conjugates. To investigate the scope of the membrane effect, cell-based cytotoxicity and membrane permeability assays were performed employing flow cytometry. Here, we observed a dose-dependent and specific cytotoxic effect on HIV-1 Env-expressing cells by the small-molecule bifunctional inhibitor. Most importantly, Env-negative cells were not susceptible to the cytotoxic effect upon exposure to this construct at concentrations where cell-killing effects were observed for Env-positive cells. Computational structural modeling supports a mechanism in which the bifunctional inhibitors bind to the gp120 and gp41 subunits in tandem in open-state Env trimers and induce relative motion of the gp120 subunits consistent with models of Env inactivation. This observation supports the idea that the cell-killing effect of the small-molecule bifunctional inhibitor is due to specific Env conformational triggering. This work lays important groundwork to advance a small-molecule bifunctional inhibitor approach for eliminating Env-expressing infected cells and the eradication of HIV-1.

A novel severe acute respiratory (SARS)-like coronavirus (SARS-CoV-2) is responsible for the current global coronavirus disease 2019 (COVID-19) pandemic, infecting millions of people and causing hundreds of thousands of deaths. The viral entry of SARS-CoV-2 depends on an interaction between the receptor binding domain of its trimeric Spike glycoprotein and the human angiotensin converting enzyme 2 (ACE2) receptor. A better understanding of the Spike/ACE2 interaction is still required to design anti-SARS-CoV-2 therapeutics. Here, we investigated the degree of cooperativity of ACE2 within both the SARS-CoV-2 and the closely related SARS-CoV-1 membrane-bound S glycoproteins. We show that there exist differential inter-protomer conformational transitions between both Spike trimers. Interestingly, the SARS-CoV-2 spike exhibits a positive cooperativity for monomeric soluble ACE2 binding when compared to the SARS-CoV-1 spike, which might have more structural restrains. Our findings can be of importance in the development of therapeutics that block the Spike/ACE2 interaction.

The Dual‐Acting Virolytic Entry Inhibitors, or DAVEI's, are a class of recombinant chimera fusion proteins consisting of a lectin, a flexible polypeptide linker, and a fragment of the membrane‐proximal external region (MPER) of HIV‐1 gp41. DAVEIs trigger virolysis of HIV‐1 virions through interactions with the trimeric envelope glycoprotein complex (Env), though the details of these interactions are not fully determined as yet. The purpose of this work was to use structural modeling to rationalize a dependence of DAVEI potency on the molecular length of the linker connecting the two components. We used temperature accelerated molecular dynamics and on‐the‐fly parameterization to compute free energy versus end‐to‐end distance for two different linker lengths, DAVEI L0 (His(6)) and DAVEI L2 ([Gly(4)Ser](2)His(6)). Additionally, an envelope model was created based on a cryo‐electron microscopy‐derived structure of a cleaved, soluble Env construct, with high‐mannose glycans added which served as putative docking locations for the lectin, along with MPER added that served as a putative docking location for the MPER region of DAVEI (MPER(DAVEI)). Using MD simulation, distances between the lectin C‐terminus and Env gp41 MPER were measured. We determined that none of the glycans were close enough to gp41 MPER to allow DAVEI L0 to function, while one, N448, will allow DAVEI L2 to function. These findings are consistent with the previously determined dependence of lytic function on DAVEI linker lengths. This supports the hypothesis that DAVEI's engage Env at both glycans and the Env MPER in causing membrane poration and lysis.

During human immunodeficiency virus type 1 (HIV-1) entry into cells, the viral envelope glycoprotein (Env) trimer [(gp120/gp41)(3)] binds the receptors CD4 and CCR5 and fuses the viral and cell membranes. CD4 binding changes Env from a pretriggered (state-1) conformation to more open downstream conformations. BMS-378806 (here called BMS-806) blocks CD4-induced conformational changes in Env important for entry and is hypothesized to stabilize a state-1-like Env conformation, a key vaccine target. Here, we evaluated the effects of BMS-806 on the conformation of Env on the surface of cells and virus-like particles. BMS-806 strengthened the labile, noncovalent interaction of gp120 with the Env trimer, enhanced or maintained the binding of most broadly neutralizing antibodies, and decreased the binding of poorly neutralizing antibodies. Thus, in the presence of BMS-806, the cleaved Env on the surface of cells and virus-like particles exhibits an antigenic profile consistent with a state-1 conformation. We designed novel BMS-806 analogues that stabilized the Env conformation for several weeks after a single application. These long-acting BMS-806 analogues may facilitate enrichment of the metastable state-1 Env conformation for structural characterization and presentation to the immune system. IMPORTANCE The envelope glycoprotein (Env) spike on the surface of human immunodeficiency virus type 1 (HIV-1) mediates the entry of the virus into host cells and is also the target for antibodies. During virus entry, Env needs to change shape. Env flexibility also contributes to the ability of HIV-1 to evade the host immune response; many shapes of Env raise antibodies that cannot recognize the functional Env and therefore do not block virus infection. We found that an HIV-1 entry inhibitor, BMS-806, stabilizes the functional shape of Env. We developed new variants of BMS-806 that stabilize Env in its natural state for long periods of time. The availability of such long-acting stabilizers of Env shape will allow the natural Env conformation to be characterized and tested for efficacy as a vaccine.

The Phe43 cavity of HIV-1 envelope glycoproteins (Env) is an attractive druggable target. New promising compounds, including small CD4 mimetics (CD4mc), were shown to insert deeply into this cavity. Here, we identify a new network of residues that helps to shape this highly conserved CD4 binding pocket and characterize the structural determinants responsible for Env sensitivity to small CD4 mimetics., The HIV-1 envelope glycoproteins (Env) undergo conformational changes upon interaction of the gp120 exterior glycoprotein with the CD4 receptor. The gp120 inner domain topological layers facilitate the transition of Env to the CD4-bound conformation. CD4 engages gp120 by introducing its phenylalanine 43 (Phe43) in a cavity (“the Phe43 cavity”) located at the interface between the inner and outer gp120 domains. Small CD4-mimetic compounds (CD4mc) can bind within the Phe43 cavity and trigger conformational changes similar to those induced by CD4. Crystal structures of CD4mc in complex with a modified CRF01_AE gp120 core revealed the importance of these gp120 inner domain layers in stabilizing the Phe43 cavity and shaping the CD4 binding site. Our studies reveal a complex interplay between the gp120 inner domain and the Phe43 cavity and generate useful information for the development of more-potent CD4mc.

Molecular dynamics simulations of a solvent-free coarse-grained lipid model are used to characterize the mechanisms by which lipid-bilayer hemifusion diaphragm (HD) intermediates relax, across a range of global compositions of negative intrinsic curvature (NIC) lipids and neutral-curvature lipids. At low concentrations of NIC lipids, rapid fission produces a double bilayer end state through a lateral diffusion-based mechanism enabled by spontaneous rim-pore defects. At moderately higher NIC lipid concentrations, rim pores are absent and stable leaflet three-junctions persist, revealing an HD relaxation mechanism entirely reliant on lipid flip-flop, and end states that are either stable fusion pores or stable HD's. These fusogenic systems exhibit dynamics highly dependent on NIC lipid concentration via an underlying sensitivity of flip-flop rates for neutral lipids on NIC lipid concentration. This work illustrates that HD dynamics may be altered through regulation of lipid composition in the immediate three-junction region. This work further highlights the potential role of flippases in biological fusion and the importance of lipid composition on fusion dynamics.

The effects of pendant alkyl chain length n ranging from n = 0 to n = 10 in amidoamine crosslinkers on mechanical properties of epoxy thermosets are explored in this work. The glassy-state Young’s modulus was estimated using non-equilibrium molecular dynamics (MD) simulations and compared with experiments. Both simulations and experiments showed that Young’s modulus decreases with increase in n. Stress partitioning based on molecular interaction types showed that both Lennard-Jones and covalent bond interactions were responsible for this sensitivity, with Coulombic interactions playing no significant role. The dependence of Young’s modulus on n was strain-rate dependent in the simulations, with moderate to high strain rates showing no sensitivity. A strong correlation was observed between Young’s modulus from non-equilibrium MD and volume fraction of methylenes estimated from equilibrium MD. Poisson’s ratios of all systems predicted from the simulations were insensitive to n, indicating a lack of anisotropy. The information revealed here on the roles of various intermolecular interactions on the mechanical properties of these thermosets could be useful for design of crosslinkers.

The seasonal nature of outbreaks of respiratory viral infections with increased transmission during low temperatures has been well established. Accordingly, temperature has been suggested to play a role on the viability and transmissibility of SARS-CoV-2, the virus responsible for the COVID-19 pandemic. The receptor-binding domain (RBD) of the Spike glycoprotein is known to bind to its host receptor angiotensin-converting enzyme 2 (ACE2) to initiate viral fusion. Using biochemical, biophysical, and functional assays to dissect the effect of temperature on the receptor-Spike interaction, we observed a significant and stepwise increase in RBD-ACE2 affinity at low temperatures, resulting in slower dissociation kinetics. This translated into enhanced interaction of the full Spike glycoprotein with the ACE2 receptor and higher viral attachment at low temperatures. Interestingly, the RBD N501Y mutation, present in emerging variants of concern (VOCs) that are fueling the pandemic worldwide (including the B.1.1.7 (α) lineage), bypassed this requirement. This data suggests that the acquisition of N501Y reflects an adaptation to warmer climates, a hypothesis that remains to be tested.

12-Hydroxystearic acid (12HSA) and its derivatives are well-known organogelators, and they play critical roles in a variety of applications. The overall aggregate structure of 12HSA is sensitive to the chirality at the 12th carbon, but a fundamental understanding of this dependence is lacking. In this study, molecular dynamics simulations were conducted on microsecond long time scales for (1) (R)-12HSA, (2) (S)-12HSA, and (3) a 50/50 racemic mixture, each solvated at 12.5 wt % in explicit hexane. Self-assembly was accelerated by turning off alkyl chain dihedral gauche states and forcing the molecules to adopt an all-trans conformation. The stability of the resulting aggregates was tested by quenching them with access to gauche states restored. Ordered aggregates produced from optically pure (R)- and (S)-12HSA remained stable for at least 1 μs. The characteristic ordered structure observed is termed a "ring-of-rings" motif, and it contains two twisted six-membered ringlike bundles connected through acetic acid dimerization and surrounded by six satellite bundles. The chirality at the 12th carbon dictates the overall twist of the rings and thereby the handedness of the aggregates. Racemic mixtures did not produce stable ordered aggregates likely due to insufficient enantiomerically pure ring formation. The most prevalent finite-size effect observed was the stochastic formation of percolating aggregates, which were later avoided by using solvent-permeable, solute-impermeable, confining walls. The resulting ordered aggregates were in all important ways identical to those produced in unconfined systems. The combination of cycling off and on gauche states and the semipermeable walls may be an important new way to study the self-assembly underlying aggregation at industrially relevant concentrations of surfactants in organic solvents.

Polyamine crosslinkers functionalized with pendant fatty acids are commercially used in epoxy coatings in part to enhance water barrier properties for corrosion prevention. However, a systematic understanding of the links between monomer molecular structure and material properties of such systems remains elusive, which limits our ability to design newer and more environmentally friendly coating systems. In this work, the effect of n-alkyl chain length in fatty-acid-functionalized polyamines on crosslinked epoxy-amidoamine systems is studied using a combined experimental and simulation-based approach. Both experiments and simulations show that density decreases and volume-expansion coefficients increase with increasing pendant chain length. Furthermore, it is shown that the trends in density and coefficient of volume thermal expansion with chain length obtained from simulations are consistent with an ideal mixing approximation. Interestingly, however, the glass-transition temperature T g is found to be insensitive to chain length. Molecular simulations reveal that increasing the alkyl chain length from four to ten carbons does not introduce new flexibility mechanisms to the dense thermosets, which explains the T g insensitivity. This work demonstrates a new way to significantly decrease the density of a thermoset polymer without compromising desirable properties such as high T g and low coefficient of thermal expansion, and therefore may provide sounder rationale for molecular-based design of epoxy-based coatings.

Peptide triazole (PT) antagonists interact with gp120 subunits of HIV-1 Env trimers to block host cell receptor interactions, trigger gp120 shedding, irreversibly inactivate virus and inhibit infection. Despite these enticing functions, understanding the structural mechanism of PT-Env trimer encounter has been limited. In this work, we combined competition interaction analysis and computational simulation to demonstrate PT binding to the recombinant soluble trimer, BG505 SOSIP.664, a stable variant that resembles native virus spikes in binding to CD4 receptor as well as known conformationally-dependent Env antibodies. Binding specificity and computational modeling fit with encounter through complementary PT pharmacophore Ile-triazolePro-Trp interaction with a 2-subsite cavity in the Env gp120 subunit of SOSIP trimer similar to that in monomeric gp120. These findings argue that PTs are able to recognize and bind a closed prefusion state of Env trimer upon HIV-1 encounter. The results provide a structural model of how PTs exert their function on virion trimeric spike protein and a platform to inform future antagonist design. Proteins 2017; 85:843-851. © 2016 Wiley Periodicals, Inc.

Dual-acting virucidal entry inhibitors (DAVEIs) have previously been shown to cause irreversible inactivation of HIV-1 Env-presenting pseudovirus by lytic membrane transformation. This study examined whether this transformation could be generalized to include membranes of Env-presenting cells. Flow cytometry was used to analyze HEK293T cells transiently transfected with increasing amounts of DNA encoding JRFL Env, loaded with calcein dye, and treated with serial dilutions of microvirin (Q831K/M83R)-DAVEI. Comparing calcein retention against intact Env expression (via Ab 35O22) on individual cells revealed effects proportional to Env expression. "Low-Env" cells experienced transient poration and calcein leakage, while "high-Env" cells were killed. The cell-killing effect was confirmed with an independent mitochondrial activity-based cell viability assay, showing dose-dependent cytotoxicity in response to DAVEI treatment. Transfection with increasing quantities of Env DNA showed further shifts toward "High-Env" expression and cytotoxicity, further reinforcing the Env dependence of the observed effect. Controls with unlinked DAVEI components showed no effect on calcein leakage or cell viability, confirming a requirement for covalently linked DAVEI compounds to achieve Env transformation. These data demonstrate that the metastability of Env is an intrinsic property of the transmembrane protein complex and can be perturbed to cause membrane disruption in both virus and cell contexts.

Human immunodeficiency virus (HIV-1) entry into cells is mediated by the viral envelope glycoprotein (Env) trimer, which consists of three gp120 exterior glycoproteins and three gp41 transmembrane glycoproteins. When gp120 binds sequentially to the receptors CD4 and CCR5 on the target cell, the metastable Env trimer is triggered to undergo entry-related conformational changes. PF-68742 is a small molecule that inhibits the infection of a subset of HIV-1 strains by interfering with an Env function other than receptor binding. Determinants of HIV-1 resistance to PF-68742 map to the disulfide loop and fusion peptide of gp41. Of the four possible PF-68742 stereoisomers, only one, MF275, inhibited the infection of CD4-positive CCR5-positive cells by some HIV-1 strains. MF275 inhibition of these HIV-1 strains occurred after CD4 binding but before the formation of the gp41 six-helix bundle. Unexpectedly, MF275 activated the infection of CD4-negative CCR5-positive cells by several HIV-1 strains resistant to the inhibitory effects of the compound in CD4-positive target cells. In contrast to CD4 complementation by CD4-mimetic compounds, activation of CD4-independent infection by MF275 did not depend upon the availability of the gp120 Phe 43 cavity. Sensitivity to inhibitors indicates that MF275-activated virus entry requires formation/exposure of the gp41 heptad repeat (HR1) as well as CCR5 binding. MF275 apparently activates a virus entry pathway parallel to that triggered by CD4 and CD4-mimetic compounds. Strain-dependent divergence in Env conformational transitions allows different outcomes, inhibition or activation, in response to MF275. Understanding the mechanisms of MF275 activity should assist efforts to optimize its utility. IMPORTANCE Envelope glycoprotein (Env) spikes on the surface of human immunodeficiency virus (HIV-1) bind target cell receptors, triggering changes in the shape of Env. We studied a small molecule, MF275, that also induced shape changes in Env. The consequences of MF275 interaction with Env depended on the HIV-1 strain, with infection by some viruses inhibited and infection by other viruses enhanced. These studies reveal the strain-dependent diversity of HIV-1 Envs as they undergo shape changes in proceeding down the entry pathway. Appreciation of this diversity will assist attempts to develop broadly active inhibitors of HIV-1 entry.

HIV-1 protease (PR) is the viral protein responsible for virion maturation, and its mechanisms of action remain incompletely understood. PR is dimeric and contains two flexible, symmetry-related flaps, which act as a gate to inhibit access to the binding pocket and hold the polypeptide substrate in the binding pocket once bound. Wide flap opening, a conformational change assumed to be necessary for substrate binding, is a rare event in the closed and bound form. In this study, we use molecular dynamics (MD) simulations and advanced MD techniques including temperature acceleration and string method in collective variables to study the conformational changes associated with substrate unbinding of both wild-type and F99Y mutant PR. The F99Y mutation is shown via MD to decouple the closing of previously unrecognized distal pockets from substrate unbinding. To determine whether or not the F99Y mutation affects the energetic cost of wide flap opening, we use string method in collective variables to determine the minimum free-energy mechanism for wide flap opening in concert with distal pocket closing. The results indicate that the major energetic cost in flap opening is disengagement of the two flap-tip Ile50 residues from each other and is not affected by the F99Y mutation.

Locating saddle points on free energy surfaces is key in characterizing multistate transition events in complicated molecular-scale systems. Because these saddle points represent transition states, determining minimum free energy pathways to these saddles and measuring their free energies relative to their connected minima are further necessary, for instance, to estimate transition rates. In this work, we propose a new multistring version of the climbing string method in collective variables to locate all saddles and corresponding pathways on free energy surfaces. The method uses dynamic strings to locate saddles and static strings to keep a history of prior strings converged to saddles. Interaction of the dynamic strings with the static strings is used to avoid the convergence to already-identified saddles. Additionally, because the strings approximate curves in collective-variable space, and we can measure free energy along each curve, identification of any saddle’s two connected minima is guaranteed. We demonstrate this method to map the network of stationary points in the 2D and 4D free energy surfaces of alanine dipeptide and alanine tripeptide, respectively.

The HIV-1 envelope glycoprotein (Env) trimer mediates virus entry into cells. The “closed” conformation of Env is resistant to nonneutralizing antibodies (nnAbs). These antibodies mostly recognize occluded epitopes that can be exposed upon binding of CD4 or small-molecule CD4 mimetics (CD4mc). Here, we describe a new family of small molecules that expose Env to nnAbs and sensitize infected cells to antibody-dependent cellular cytotoxicity (ADCC). These compounds have a limited capacity to inhibit virus infection directly but are able to sensitize viral particles to neutralization by otherwise nonneutralizing antibodies. Structural analysis shows that some analogs of this family of CD4mc engage the gp120 Phe43 cavity by contacting the highly conserved D368 residue, making them attractive scaffolds for drug development. IMPORTANCE HIV-1 has evolved multiple strategies to avoid humoral responses. One efficient mechanism is to keep its envelope glycoprotein (Env) in its “closed” conformation. Here, we report on a new family of small molecules that are able to “open up” Env, thus exposing vulnerable epitopes. This new family of molecules binds in the Phe43 cavity and contacts the highly conserved D368 residue. The structural and biological attributes of molecules of this family make them good candidates for drug development.

The design, synthesis and validation of a family of small molecule “Dual-Action Virucidal EntryInhibitors” (DAVEIs) has been achieved that result in irreversible lytic inactivation of HIV-1 virions. These constructs contained two functional components that endow the capacity to bindsimultaneously to both the gp120 and gp41 subunits of the HIV-1 Envelope glycoprotein (Env). One component is derived from BNM-III-170, a small molecule CD4 mimic warhead that binds togp120. The second component, a Trp3 peptide, is a 9-amino acid segment based on the gp41 Membrane Proximal External Region (MPER) that has been proposed to bind to the gp41 MPERdomain of the Env. The resulting smDAVEIs both inhibit infection with low micromolar potency and induce lysis of the HIV-1 virion. The lytic activity was selective for functional HIV-1 virions. Crucially, virolysis was found to be dependent on covalent tethering of the BNM-III-170 and Trp3 domains with various spacers, as coadministration of the un-crosslinked components proved not to be lytic. Computational modeling supports a mechanism in which DAVEIs bind to open-state Env trimers and induce relative motion of gp120 subunits that further opens the trimers. Overall, this work represents a promising new step toward the use of small-molecule DAVEIs for eradication of HIV.

KR13, a peptide triazole thiol previously established to inhibit HIV-1 infection and cause virus lysis, was evaluated by flow cytometry against JRFL Env-presenting cells to characterize induced Env and membrane transformations leading to irreversible inactivation. Transiently transfected HEK293T cells were preloaded with calcein dye, treated with KR13 or its thiol-blocked analogue KR13b, fixed, and stained for gp120 (35O22), MPER (10E8), 6-helix-bundle (NC-1), immunodominant loop (50-69), and fusion peptide (VRC34.01). KR13 induced dose-dependent transformations of Env and membrane characterized by transient poration, MPER exposure, and 6-helix-bundle formation (analogous to native fusion events), but also reduced immunodominant loop and fusion peptide exposure. Using a fusion peptide mutant (V504E), we found that KR13 transformation does not require functional fusion peptide for poration. In contrast, simultaneous treatment with fusion inhibitor T20 alongside KR13 prevented membrane poration and MPER exposure, showing that these events require 6-helix-bundle formation. Based on these results, we formulated a model for PTT-induced Env transformation portraying how, in the absence of CD4/co-receptor signaling, PTT may provide alternate means of perturbing the metastable Env-membrane complex, and inducing fusion-like transformation. In turn, the results show that such transformations are intrinsic to Env and can be diverted for irreversible inactivation of the protein complex.

We study the thermodynamic stability of the native state of the human prion protein using a new free-energy method, replica-exchange on-the-fly parameterization. This method is designed to overcome hidden-variable sampling limitations to yield nearly error-free free-energy profiles along a conformational coordinate. We confirm that all four (M129V, D178N) polymorphs have a ground-state conformation with three intact β-sheet hydrogen bonds. Additionally, they are observed to have distinct metastabilities determined by the side-chain at position 129. We rationalize these findings with reference to the prion “strain” hypothesis, which links the variety of transmissible spongiform encephalopathy phenotypes to conformationally distinct infectious prion forms and classifies distinct phenotypes of sporadic Creutzfeldt-Jakob disease based solely on the 129 polymorphism. Because such metastable structures are not easily observed in structural experiments, our approach could potentially provide new insights into the conformational origins of prion diseases and other pathologies arising from protein misfolding and aggregation. Fil: Paz, Sergio Alexis. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Drexel University; Estados Unidos Fil: Vanden-Eijnden, Eric. University of New York; Estados Unidos Fil: Abrams, Cameron F.. Drexel University; Estados Unidos

Molecular dynamics simulations are used to study highly cross-linked epoxy networks comprised of furanyl epoxy monomer, 2,5-bis[(2-oxiranylmethoxy)methyl]-furan (BOF), that is cross-linked by two furanyl amine hardeners, 5,5'-methylenedifurfurylamine (DFDA) and 5,5'-ethylidenedifurfirylamine (CH3-DFDA). Important properties of these fully furan-based systems, including room temperature density, glass transition temperature, and Young's modulus are found to agree with previous experimental results. We also compare the simulated and experimental values of four fully furan-based thermosetting materials to those using the conventional resin diglycidyl ether of bisphenol A (DGEBA) cured with the two furanyl hardeners. Our simulation results predict a slight decrease in density and Young's modulus, but no impact on the glass transition temperature, upon adding the methyl group in DFDA. Detailed analyses of the MD trajectories reveal the underlying mechanisms responsible for the observed structure/property relations, which center on the lack of collinear covalent bonds in the BOF molecular structure. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016.

Small-molecule CD4 mimics (SMCM’s) bind to the gp120 subunit of the HIV-1 envelope glycoprotein (Env) and have been optimized to block cell infection in vitro. The lack of the V1/2 and V3 loops and the presence of the β2/3 and β20/21 strands (bridging sheet) in the available structures of the monomeric gp120 core may limit its applicability as a target for further synthetic optimization of SMCM potency and/or breadth. Here, we employ a combination of binding-site search, docking, estimation of protein–ligand interaction energy, all-atom molecular dynamics, and ELISA-based CD4-binding competition assays to create, characterize, and rationalize models of first- and second-generation of SMCM’s bound to the distinct, trimeric BG505 SOSIP.664 structures 4NCO and 4TVP containing V1/2 and V3 loops with no bridging sheet. We demonstrate that the in silico neutralization of the highly conserved D368 is necessary to obtain the correct orientation of SMCM in their binding site when docking against the monomeric gp120 core. The computational results correlate with IC50’s measured in CD4 binding competition ELISA and with KD’s measured on gp120 core monomer. This supports the hypothesis that the 4NCO trimeric structure represents a viable target for further SMCM’s optimization with advantages over both the 4TVP trimer and gp120 core monomer. Finally, the docking protocol has been optimized to screen compounds that can clearly interact with the highly conserved residue D368, increasing the likelihood of future optimizations to arrive at SMCM’s with a broader spectrum of activity.

12-Hydroxystearic acid (12HSA) is a well-known organogelator, and its metal salts and derivatives find roles in many important applications. The structures of aggregates of 12-hydroxysteric acid and its salts depend sensitively on cation type, but a fundamental understanding of this phenomenon is lacking. In this study, molecular dynamics simulations were conducted on the microsecond long time scales for (1) 12HSA and (2) its lithium salt, each at 12.5 wt % in explicit hexane solvent. Self-assembly was accelerated by using a modified potential to prohibit alkane chain dihedral gauche states (all-trans-12HSA) and then verified by continuation using standard force-field parameters. In three independent simulation, acceleration using "gauche-less" potentials resulted in self-assembled pseudocrystalline aggregates through formation of polarized five- and six-membered rings between inter-12-hydroxyl groups and head-to-head carboxylic acid dimerization. When subjected to the unmodified dihedral potential, two of the three structures remained stable after 1 μs of MD. Stable structures exhibited a "ring-of-rings" motif, composed of two six-membered acetic acid-dimerized ring bundles with six satellite rings, while the unstable structure did not. In strong contrast, the lithium salt produced a network of fibrils that spanned the volume of the sample. When lithium ions were substituted for carboxylic acid protons in the stable acid structures, they remained intact but lost their chiral nature. Both the acid and lithium structures displayed scattering peaks that agreed with experiment. Taken together, our results suggest that this ring-of-rings structure could be a primary feature of the self-assembly of 12HSA in organic solvents.

We developed coarse-grained (CG) molecular representations of mixtures of diglycidyl ether of bisphenol-A (DGEBA) and poly(oxypropylene) diamine (POP-DA) for use in CG molecular dynamics (MD) simulations. In the CG representation, DGEBA is comprised of three beads of two types and POP-DA also by three beads of two types. Atomistic MD of liquid systems was performed to derive intra- and inter-bead potentials via Boltzmann inversion. While the bonded potentials, composed of bond stretching and angle bending, were parameterized directly from the distribution functions of all atomistic molecular dynamics trajectories, the non-bonded potentials were derived from the iterative Boltzmann Inversion with a given set of coarse-grained interactions. CG systems correctly reproduced liquid and crosslinked densities. Under uniaxial tension, the Young's modulus of the CG systems was much lower than the experimental value, and we show this arises from the assumed form of the extrapolated regions of the CG potentials. By stiffening these regions, we increased the CG Young's modulus of the crosslinked systems without sacrificing the correct prediction of density. This suggests that transferrable CG potentials can be optimized for use in non-equilibrium MD for property estimation.

Curing epoxies with a mixture of low- and high-T g diamines has been proposed as a way to increase thermoset toughness. We seek here to understand the origins of toughness enhancement in systems comprised of the diamines poly(oxypropylene)diamine (POPDA) and diethyltoluenediamine (DETDA) together with the epoxy resin diglycidyl ether of bisphenol A (DGEBA) via control of network isomerization. Two curing protocols at constant overall DGEBA/DETDA/POPDA 2:1 amine:epoxy stoichiometric composition are compared: (i) curing a liquid mixture of DGEBA, DETDA and POPDA, and (ii) partially curing DGEBA with POPDA (60% of amines reacted), then adding DETDA and more DGEBA to continue to a fully cured stoichiometric sample; the latter is referred to as the “partially reacted substructure” (PRS) method. PRS samples are 50% tougher than the compositionally-identical mixed samples yet have higher T g 's than the mixed samples. We show here that MD simulations of model systems provide a molecular-level rationale for this observation. First, MD yields reasonably accurate densities and T g 's. Lower T g 's in the mixed systems are correlated to larger thermal fluctuations in positions of monomer centers enabled by more uniform dispersion of the POPDA molecules. Furthermore, the onset of crosslink bond stretching under steady uniaxial tensile strain occurs at lower strains in the mixed samples, which correlates to their lower experimental ductility. This behavior is shown to arise from POPDA molecules in the PRS system more easily deforming from their unstrained conformations than they can in the mixed systems. These findings provide further guidance in the use of control over network isomerization at constant composition to enhance toughness of thermoset systems.

Salts of 12-hydroxystearate are important organogelators and grease thickeners, but a structural rationale for their rheological properties remains elusive. We use quantum mechanical calculations and molecular dynamics (MD) simulations to analyze aggregate structures for (1) ( R)-12-hydroxystearic acid (( R)-12HSA), (2) lithium ( R)-12-hydroxystearate (( R)-Li12HS), and (3) sodium ( R)-12-hydroxystearate (( R)-Na12HS). First, quantum mechanical calculations were used to establish the structure and complexation energies of dimers of acetic acid, lithium acetate, and sodium acetate. The expected acetic acid dimer is predicted, and both the lithium acetate and sodium acetate dimer formed a C2 h-symmetric structure. All dimers were sufficiently stable to allow modeling them as pseudocovalent complexes in all-atom, explicit solvent MD. After microsecond-long MD, all systems produced strong ringlike ordered nuclei. The C2 h lithium salt molecules produced aggregates that had the most efficient packing at the head group and a higher frequency of hydroxyl hydrogen bonding compared to the sodium salt. This ordering propensity explains the high melting temperature of ( R)-Li12HS. Also, the higher frequency of hydrogen bonding leads to fewer solvent-exposed hydrogen bond partners. This explains why lithium is a common counterion in high-temperature and water-resistant greases.

We introduce a two-dimensional version of the method called on-the-fly free energy parametrization (OTFP) to reconstruct free-energy surfaces using Molecular Dynamics simulations, which we name OTFP-2D. We first test the new method by reconstructing the well-known dihedral angles free energy surface of solvated alanine dipeptide. Then, we use it to investigate the process of K + ions translocation inside the Kv1.2 channel. By comparing a series of two-dimensional free energy surfaces for ion movement calculated with different conditions on the intercalated water molecules, we first recapitulate the widely accepted knock-on mechanism for ion translocation and then confirm that permeation occurs with water molecules alternated among the ions, in accordance with the latest experimental findings. From a methodological standpoint, our new OTFP-2D algorithm demonstrates the excellent sampling acceleration of temperature-accelerated molecular dynamics and the ability to efficiently compute 2D free-energy surfaces. It will therefore be useful in large variety complex biomacromolecular simulations. Fil: Paz, Sergio Alexis. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigaciones en Físico-química de Córdoba. Universidad Nacional de Córdoba. Facultad de Ciencias Químicas. Instituto de Investigaciones en Físico-química de Córdoba; Argentina Fil: Maragliano, Luca. Fondazione Istituto Italiano Di Tecnologia; Italia Fil: Abrams, Cameron F.. Drexel University; Estados Unidos