Structural and thermodynamic framework for PIEZO1 modulation by small molecules.

Mechanosensitive PIEZO channels constitute potential pharmacological targets for multiple clinical conditions, spurring the search for potent chemical PIEZO modulators. Among them is Yoda1, a widely used synthetic small molecule PIEZO1 activator discovered through cell-based high-throughput screening. Yoda1 is thought to bind to PIEZO1's mechanosensory arm domain, sandwiched between two transmembrane regions near the channel pore. However, how the binding of Yoda1 to this region promotes channel activation remains elusive. Here, we first demonstrate that cross-linking PIEZO1 repeats A and B with disulfide bridges reduces the effects of Yoda1 in a redox-dependent manner, suggesting that Yoda1 acts by perturbing the contact between these repeats. Using molecular dynamics-based absolute binding free energy simulations, we next show that Yoda1 preferentially occupies a deeper, amphipathic binding site with higher affinity in PIEZO1 open state. Using Yoda1's binding poses in open and closed states, relative binding free energy simulations were conducted in the membrane environment, recapitulating structure-activity relationships of known Yoda1 analogs. Through virtual screening of an 8 million-compound library using computed fragment maps of the Yoda1 binding site, we subsequently identified two chemical scaffolds with agonist activity toward PIEZO1. This study supports a pharmacological model in which Yoda1 activates PIEZO1 by wedging repeats A and B, providing a structural and thermodynamic framework for the rational design of PIEZO1 modulators. Beyond PIEZO channels, the three orthogonal computational approaches employed here represent a promising path toward drug discovery in highly heterogeneous membrane protein systems.

Effects of Xylanase A double mutation on substrate specificity and structural dynamics.

While protein activity is traditionally studied with a major focus on the active site, the activity of enzymes has been hypothesized to be linked to the flexibility of adjacent regions, warranting more exploration into how the dynamics in these regions affects catalytic turnover. One such enzyme is Xylanase A (XylA), which cleaves hemicellulose xylan polymers by hydrolysis at internal ß-1,4-xylosidic linkages. It contains a "thumb" region whose flexibility has been suggested to affect the activity. The double mutation D11F/R122D was previously found to affect activity and potentially bias the thumb region to a more open conformation. We find that the D11F/R122D double mutation shows substrate-dependent effects, increasing activity on the non-native substrate ONPX2 but decreasing activity on its native xylan substrate. To characterize how the double mutant causes these kinetics changes, nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations were used to probe structural and flexibility changes. NMR chemical shift perturbations revealed structural changes in the double mutant relative to the wild-type, specifically in the thumb and fingers regions. Increased slow-timescale dynamics in the fingers region was observed as intermediate-exchange line broadening. Lipari-Szabo order parameters show negligible changes in flexibility in the thumb region in the presence of the double mutation. To help understand if there is increased energetic accessibility to the open state upon mutation, alchemical free energy simulations were employed that indicated thumb opening is more favorable in the double mutant. These studies aid in further characterizing how flexibility in adjacent regions affects the function of XylA.

Synthesis and unexpected binding of monofluorinated N,N'-diacetylchitobiose and LacdiNAc to wheat germ agglutinin.

Fluorination of carbohydrate ligands of lectins is a useful approach to examine their binding profile, improve their metabolic stability and lipophilicity, and convert them into 19F NMR-active probes. However, monofluorination of monovalent carbohydrate ligands often leads to a decreased or completely lost affinity. By chemical glycosylation, we synthesized the full series of methyl ß-glycosides of N,N'-diacetylchitobiose (GlcNAcß(1-4)GlcNAcß1-OMe) and LacdiNAc (GalNAcß(1-4)GlcNAcß1-OMe) systematically monofluorinated at all hydroxyl positions. A competitive enzyme-linked lectin assay revealed that the fluorination at the 6'-position of chitobioside resulted in an unprecedented increase in affinity to wheat germ agglutinin (WGA) by one order of magnitude. For the first time, we have characterized the binding profile of a previously underexplored WGA ligand LacdiNAc. Surprisingly, 4'-fluoro-LacdiNAc bound WGA even stronger than unmodified LacdiNAc. These observations were interpreted using molecular dynamic calculations along with STD and transferred NOESY NMR techniques, which gave evidence for the strengthening of CH/π interactions after deoxyfluorination of the side chain of the non-reducing GlcNAc. These results highlight the potential of fluorinated glycomimetics as high-affinity ligands of lectins and 19F NMR-active probes.

Application of an alchemical free energy method for the prediction of thermostable DuraPETase variants.

Non-equilibrium (NEQ) alchemical free energy calculations are an emerging tool for accurately predicting changes in protein folding free energy resulting from amino acid mutations. In this study, this method in combination with the Rosetta ddg monomer tool was applied to predict more thermostable variants of the polyethylene terephthalate (PET) degrading enzyme DuraPETase. The Rosetta ddg monomer tool efficiently enriched promising mutations prior to more accurate prediction by NEQ alchemical free energy calculations. The relative change in folding free energy of 96 single amino acid mutations was calculated by NEQ alchemical free energy calculation. Experimental validation of ten of the highest scoring variants identified two mutations (DuraPETaseS61M and DuraPETaseS223Y) that increased the melting temperature (Tm) of the enzyme by up to 1 °C. The calculated relative change in folding free energy showed an excellent correlation with experimentally determined Tm resulting in a Pearson's correlation coefficient of r = - 0.84. Limitations in the prediction of strongly stabilizing mutations were, however, encountered and are discussed. Despite these challenges, this study demonstrates the practical applicability of NEQ alchemical free energy calculations in prospective enzyme engineering projects. KEY POINTS: • Rosetta ddg monomer enriches stabilizing mutations in a library of DuraPETase variants • NEQ free energy calculations accurately predict changes in Tm of DuraPETase • The DuraPETase variants S223Y, S42M, and S61M have increased Tm.

Molecular dynamics studies of the solubility behavior of carbon dioxide (CO2 ), difluoromethane (R-32), 1-chloro-3,3,3-trifluoropropene (R-1233zd(E)) and 2,3,3,3-tetrafluoro-1-propene (R-1234yf) in pentaerythritol tetra(2-ethylhexanoate) (PEB8), pentaerythritol tetrabutyrate (PEC4) and pentaerythritol tetraoctanoate (PEC8).

To reduce the climate impact of thermal engines such as heat pumps or refrigeration machines, refrigerants with a low global warming potential need to be paired with fitting lubricants. As the contamination of those liquid components influences the efficiency and lifetime of these machines, knowledge about their solubility behavior is of great interest. Molecular simulations offer mighty tools to investigate these solubilities while giving structural insight into the systems. Here the solubility behavior of CO2 , R-32, R-1233zd(E), and R-1234yf in PEB8, PEC4, and PEC8 is compared through the solvation free energy ∆GSolv obtained by molecular dynamics simulations. To derive ∆GSolv at low computational cost, an iterative method is used to find an optimal number and distribution of intermediate states. The resulting distributions are investigated with regard to different parameters of the employed softcore-potential. ∆GSolv values for the different refrigerant-lubricant pairings at different temperatures are provided, followed by a structural analysis.

NE-RDFE: A protocol and toolkit for computing relative dissociation free energies with GROMACS between dissimilar molecules using bidirectional nonequilibrium dual topology schemes.

We describe a step-by-step protocol and toolkit for the computation of the relative dissociation free energy (RDFE) with the GROMACS molecular dynamics package, based on a novel bidirectional nonequilibrium alchemical approach. The proposed methodology does not require any intervention on the code and allows computing with good accuracy the RDFE between small molecules with arbitrary differences in volume, charge, and chemical topology. The procedure is illustrated for the challenging SAMPL9 batch of host-guest pairs. The article is supplemented by a detailed online tutorial, available at https://procacci.github.io/vdssb_gromacs/NE-RDFE and by a public Zenodo repository available at https://zenodo.org/record/6982932.

Study of SQ109 analogs binding to mycobacterium MmpL3 transporter using MD simulations and alchemical relative binding free energy calculations.

N-geranyl-N΄-(2-adamantyl)ethane-1,2-diamine (SQ109) is a tuberculosis drug that has high potency against Mycobacterium tuberculosis (Mtb) and may function by blocking cell wall biosynthesis. After the crystal structure of MmpL3 from Mycobacterium smegmatis in complex with SQ109 became available, it was suggested that SQ109 inhibits Mmpl3 mycolic acid transporter. Here, we showed using molecular dynamics (MD) simulations that the binding profile of nine SQ109 analogs with inhibitory potency against Mtb and alkyl or aryl adducts at C-2 or C-1 adamantyl carbon to MmpL3 was consistent with the X-ray structure of MmpL3 - SQ109 complex. We showed that rotation of SQ109 around carbon-carbon bond in the monoprotonated ethylenediamine unit favors two gauche conformations as minima in water and lipophilic solvent using DFT calculations as well as inside the transporter's binding area using MD simulations. The binding assays in micelles suggested that the binding affinity of the SQ109 analogs was increased for the larger, more hydrophobic adducts, which was consistent with our results from MD simulations of the SQ109 analogues suggesting that sizeable C-2 adamantyl adducts of SQ109 can fill a lipophilic region between Y257, Y646, F260 and F649 in MmpL3. This was confirmed quantitatively by our calculations of the relative binding free energies using the thermodynamic integration coupled with MD simulations method with a mean assigned error of 0.74 kcal mol-1 compared to the experimental values.

Predicting antibiotic resistance in complex protein targets using alchemical free energy methods.

Drug resistant Mycobacterium tuberculosis, which mostly results from single nucleotide polymorphisms in antibiotic target genes, poses a major threat to tuberculosis treatment outcomes. Relative binding free energy (RBFE) calculations can rapidly predict the effects of mutations, but this approach has not been tested on large, complex proteins. We use RBFE calculations to predict the effects of M. tuberculosis RNA polymerase and DNA gyrase mutations on rifampicin and moxifloxacin susceptibility respectively. These mutations encompass a range of amino acid substitutions with known effects and include large steric perturbations and charged moieties. We find that moderate numbers (n = 3-15) of short RBFE calculations can predict resistance in cases where the mutation results in a large change in the binding free energy. We show that the method lacks discrimination in cases with either a small change in energy or that involve charged amino acids, and we investigate how these calculation errors may be decreased.

Alchemical Design of Pharmacological Chaperones with Higher Affinity for Phenylalanine Hydroxylase.

Phenylketonuria (PKU) is a rare metabolic disease caused by variations in a human gene, PAH, encoding phenylalanine hydroxylase (PAH), and the enzyme converting the essential amino acid phenylalanine into tyrosine. Many PKU-causing variations compromise the conformational stability of the encoded enzyme, decreasing or abolishing its catalytic activity, and leading to an elevated concentration of phenylalanine in the blood, which is neurotoxic. Several therapeutic approaches have been developed to treat the more severe manifestations of the disorder, but they are either not entirely effective or difficult to adhere to throughout life. In a search for novel pharmacological chaperones to treat PKU, a lead compound was discovered (compound IV) that exhibited promising in vitro and in vivo chaperoning activity on PAH. The structure of the PAH-IV complex has been reported. Here, using alchemical free energy calculations (AFEC) on the structure of the PAH-IV complex, we design a new generation of compound IV-analogues with a higher affinity for the enzyme. Seventeen novel analogues were synthesized, and thermal shift and isothermal titration calorimetry (ITC) assays were performed to experimentally evaluate their stabilizing effect and their affinity for the enzyme. Most of the new derivatives bind to PAH tighter than lead compound IV and induce a greater thermostabilization of the enzyme upon binding. Importantly, the correspondence between the calculated alchemical binding free energies and the experimentally determined ΔΔGb values is excellent, which supports the use of AFEC to design pharmacological chaperones to treat PKU using the X-ray structure of their complexes with the target PAH enzyme.

SAMPL7 TrimerTrip host-guest binding affinities from extensive alchemical and end-point free energy calculations.

The prediction of host-guest binding affinities with computational modelling is still a challenging task. In the 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge, a new host named TrimerTrip was synthesized and the thermodynamic parameters of 16 structurally diverse guests binding to the host were characterized. In the TrimerTrip-guest challenge, only structures of the host and the guests are provided, which indicates that the predictions of both the binding poses and the binding affinities are under assessment. In this work, starting from the binding poses obtained from our previous enhanced sampling simulations in the configurational space, we perform extensive alchemical and end-point free energy calculations to calculate the host-guest binding affinities retrospectively. The alchemical predictions with two widely accepted charge schemes (i.e. AM1-BCC and RESP) are in good agreement with the experimental reference, while the end-point estimates perform poorly in reproducing the experimental binding affinities. Aside from the absolute value of the binding affinity, the rank of binding free energies is also crucial in drug design. Surprisingly, the end-point MM/PBSA method seems very powerful in reproducing the experimental rank of binding affinities. Although the length of our simulations is long and the intermediate spacing is dense, the convergence behavior is not very good, which may arise from the flexibility of the host molecule. Enhanced sampling techniques in the configurational space may be required to obtain fully converged sampling. Further, as the length of sampling in alchemical free energy calculations already achieves several hundred ns, performing direct simulations of the binding/unbinding event in the physical space could be more useful and insightful. More details about the binding pathway and mechanism could be obtained in this way. The nonequilibrium method could also be a nice choice if one insists to use the alchemical method, as the intermediate sampling is avoided to some extent.

Binding thermodynamics and interaction patterns of human purine nucleoside phosphorylase-inhibitor complexes from extensive free energy calculations.

Human purine nucleoside phosphorylase (hPNP) plays a significant role in the catabolism of deoxyguanosine. The trimeric protein is an important target in the treatment of T-cell cancers and autoimmune disorders. Experimental studies on the inhibition of the hPNP observe that the first ligand bound to one of three subunits effectively inhibits the protein, while the binding of more ligands to the subsequent sites shows negative cooperativities. In this work, we performed extensive end-point and alchemical free energy calculations to determine the binding thermodynamics of the trimeric protein-ligand system. 13 Immucillin inhibitors with experimental results are under calculation. Two widely accepted charge schemes for small molecules including AM1-BCC and RESP are adopted for ligands. The results of RESP are in better agreement with the experimental reference. Further investigations of the interaction networks in the protein-ligand complexes reveal that several residues play significant roles in stabilizing the complex structure. The most commonly observed ones include PHE200, GLU201, MET219, and ASN243. The conformations of the protein in different protein-ligand complexes are observed to be similar. We expect these insights to aid the development of potent drugs targeting hPNP.

Interaction of Thymine DNA Glycosylase with Oxidised 5-Methyl-cytosines in Their Amino- and Imino-Forms.

Thymine DNA Glycosylase (TDG) is an enzyme of the base excision repair mechanism and removes damaged or mispaired bases from DNA via hydrolysis of the glycosidic bond. Specificity is of high importance for such a glycosylase, so as to avoid the damage of intact DNA. Among the substrates reported for TDG are mispaired uracil and thymine but also formyl-cytosine and carboxyl-cytosine. Methyl-cytosine and hydroxylmethyl-cytosine are, in contrast, not processed by the TDG enzyme. We have in this work employed molecular dynamics simulations to explore the conformational dynamics of DNA carrying a formyl-cytosine or carboxyl-cytosine and compared those to DNA with the non-cognate bases methyl-cytosine and hydroxylmethyl-cytosine, as amino and imino tautomers. Whereas for the mispairs a wobble conformation is likely decisive for recognition, all amino tautomers of formyl-cytosine and carboxyl-cytosine exhibit the same Watson-Crick conformation, but all imino tautomers indeed form wobble pairs. The conformational dynamics of the amino tautomers in free DNA do not exhibit differences that could be exploited for recognition, and also complexation to the TDG enzyme does not induce any alteration that would indicate preferable binding to one or the other oxidised methyl-cytosine. The imino tautomers, in contrast, undergo a shift in the equilibrium between a closed and a more open, partially flipped state, towards the more open form upon complexation to the TDG enzyme. This stabilisation of the more open conformation is most pronounced for the non-cognate bases methyl-cytosine and hydroxyl-cytosine and is thus not a likely mode for recognition. Moreover, calculated binding affinities for the different forms indicate the imino forms to be less likely in the complexed DNA. These findings, together with the low probability of imino tautomers in free DNA and the indifference of the complexed amino tautomers, suggest that discrimination of the oxidised methyl-cytosines does not take place in the initial complex formation.

Non-active site mutants of HIV-1 protease influence resistance and sensitisation towards protease inhibitors.

BACKGROUND: HIV-1 can develop resistance to antiretroviral drugs, mainly through mutations within the target regions of the drugs. In HIV-1 protease, a majority of resistance-associated mutations that develop in response to therapy with protease inhibitors are found in the protease's active site that serves also as a binding pocket for the protease inhibitors, thus directly impacting the protease-inhibitor interactions. Some resistance-associated mutations, however, are found in more distant regions, and the exact mechanisms how these mutations affect protease-inhibitor interactions are unclear. Furthermore, some of these mutations, e.g. N88S and L76V, do not only induce resistance to the currently administered drugs, but contrarily induce sensitivity towards other drugs. In this study, mutations N88S and L76V, along with three other resistance-associated mutations, M46I, I50L, and I84V, are analysed by means of molecular dynamics simulations to investigate their role in complexes of the protease with different inhibitors and in different background sequence contexts. RESULTS: Using these simulations for alchemical calculations to estimate the effects of mutations M46I, I50L, I84V, N88S, and L76V on binding free energies shows they are in general in line with the mutations' effect on [Formula: see text] values. For the primary mutation L76V, however, the presence of a background mutation M46I in our analysis influences whether the unfavourable effect of L76V on inhibitor binding is sufficient to outweigh the accompanying reduction in catalytic activity of the protease. Finally, we show that L76V and N88S changes the hydrogen bond stability of these residues with residues D30/K45 and D30/T31/T74, respectively. CONCLUSIONS: We demonstrate that estimating the effect of both binding pocket and distant mutations on inhibitor binding free energy using alchemical calculations can reproduce their effect on the experimentally measured [Formula: see text] values. We show that distant site mutations L76V and N88S affect the hydrogen bond network in the protease's active site, which offers an explanation for the indirect effect of these mutations on inhibitor binding. This work thus provides valuable insights on interplay between primary and background mutations and mechanisms how they affect inhibitor binding.

Absolute Protein Binding Free Energy Simulations for Ligands with Multiple Poses, a Thermodynamic Path That Avoids Exhaustive Enumeration of the Poses.

We propose a free energy calculation method for receptor-ligand binding, which have multiple binding poses that avoids exhaustive enumeration of the poses. For systems with multiple binding poses, the standard procedure is to enumerate orientations of the binding poses, restrain the ligand to each orientation, and then, calculate the binding free energies for each binding pose. In this study, we modify a part of the thermodynamic cycle in order to sample a broader conformational space of the ligand in the binding site. This modification leads to more accurate free energy calculation without performing separate free energy simulations for each binding pose. We applied our modification to simple model host-guest systems as a test, which have only two binding poses, by using a single decoupling method (SDM) in implicit solvent. The results showed that the binding free energies obtained from our method without knowing the two binding poses were in good agreement with the benchmark results obtained by explicit enumeration of the binding poses. Our method is applicable to other alchemical binding free energy calculation methods such as the double decoupling method (DDM) in explicit solvent. We performed a calculation for a protein-ligand system with explicit solvent using our modified thermodynamic path. The results of the free energy simulation along our modified path were in good agreement with the results of conventional DDM, which requires a separate binding free energy calculation for each of the binding poses of the example of phenol binding to T4 lysozyme in explicit solvent. © 2019 Wiley Periodicals, Inc.

Alchemical free energy calculations via metadynamics: Application to the theophylline-RNA aptamer complex.

We propose a computational workflow for robust and accurate prediction of both binding poses and their affinities at early stage in designing drug candidates. Small, rigid ligands with few intramolecular degrees of freedom, for example, fragment-like molecules, have multiple binding poses, even at a single binding site, and their affinities are often close to each other. We explore various structures of ligand binding to a target through metadynamics using a small number of collective variables, followed by reweighting to obtain the atomic coordinates. After identifying each binding pose by cluster analysis, we perform alchemical free energy calculations on each structure to obtain the overall value. We applied this protocol in computing free energy of binding for the theophylline-RNA aptamer complex. Of the six (meta)stable structures found, the most favorable binding structure is consistent with the structure obtained by NMR. The overall free energy of binding reproduces the experimental values very well.

Multi-phase Boltzmann weighting: accounting for local inhomogeneity in molecular simulations of water-octanol partition coefficients in the SAMPL6 challenge.

Accurately computing partition coefficients is a pivotal part of drug discovery. Specifically, octanol-water partition coefficients can provide information into hydrophobicity of drug-like molecules, as well as a de facto representation of membrane permeability. However, one challenge facing the computation of partition coefficients is the need to encapsulate various microscopic environments. These include areas of largely bulk solvent (i.e., either water or octanol) or regions where octanol is saturated with water or areas of higher salt concentration. Also, tautomeric effects require consideration. Thus, we present a Boltzmann weighting approach that incorporates transfer free energies across varying microscopic media, as well as varying tautomeric state, to compute partition coefficients in the SAMPL6 challenge.

Computational Methods Used to Explore Transport Events in Biological Systems.

Transport of various molecules facilitated with membrane proteins is necessary for maintaining homeostasis in living cells. In humans, dysfunction of these proteins leads to many diseases. Thus, understanding how the membrane proteins function may help using them as therapeutic targets. To successfully investigate the mechanistic aspects of transport, the choice of appropriate methods is crucial. We review the computational methods that have proven most effective in investigating transport events, specifically, deterministic time-dependent classical molecular dynamics and its enhanced sampling variants, as well as methods based on Brownian dynamics. We describe technical aspects of these methods and examples of their novel variants or combinations that have been recently and successfully applied in the transport studies. We also discuss the difficulties related to these methods and provide possible solutions to avoid them.

Blinded prediction of protein-ligand binding affinity using Amber thermodynamic integration for the 2018 D3R grand challenge 4.

In the framework of the 2018 Drug Design Data Resource grand challenge 4, blinded predictions on relative binding free energy were performed for a set of 39 ligands of the Cathepsin S protein. We leveraged the GPU-accelerated thermodynamic integration of Amber 18 to advance our computational prediction. When our entry was compared to experimental results, a good correlation was observed (Kendall's τ: 0.62, Spearman's ρ: 0.80 and Pearson's R: 0.82). We designed a parallelized transformation map that placed ligands into several groups based on common alchemical substructures; TI transformations were carried out for each ligand to the relevant substructure, and between substructures. Our calculations were all conducted using the linear potential scaling scheme in Amber TI because we believe the softcore potential/dual-topology approach as implemented in current Amber TI is highly fault-prone for some transformations. The issue is illustrated by using two examples in which typical preparation for the dual-topology approach of Amber TI fails. Overall, the high accuracy of our prediction is a result of recent advances in force fields (ff14SB and GAFF), as well as rapid calculation of ensemble averages enabled by the GPU implementation of Amber. The success shown here in a blinded prediction strongly suggests that alchemical free energy calculation in Amber is a promising tool for future commercial drug design.

Calculate protein-ligand binding affinities with the extended linear interaction energy method: application on the Cathepsin S set in the D3R Grand Challenge 3.

We participated in the Cathepsin S (CatS) sub-challenge of the Drug Design Data Resource (D3R) Grand Challenge 3 (GC3) in 2017 to blindly predict the binding poses of 24 CatS-bound ligands, the binding affinity ranking of 136 ligands, and the binding free energies of a subset of 33 ligands in Stage 1A and Stage 2. Our submitted predictions ranked relatively well compared to the submissions from other participants. Here we present our methodologies used in the challenge. For the binding pose prediction, we employed the Glide module in the Schrodinger Suite 2017 and AutoDock Vina. For the binding affinity/free energy prediction, we carried out molecular dynamics simulations of the complexes in explicit water solvent with counter ions, and then estimated the binding free energies with our newly developed model of extended linear interaction energy (ELIE), which is inspired by two other popular end-point approaches: the linear interaction energy (LIE) method, and the molecular mechanics with Poisson-Boltzmann surface area solvation method (MM/PBSA). Our studies suggest that ELIE is a good trade-off between efficiency and accuracy, and it is appropriate for filling the gap between the high-throughput docking and scoring methods and the rigorous but much more computationally demanding methods like free energy perturbation (FEP) or thermodynamics integration (TI) in computer-aided drug design (CADD) projects.

Roles of conserved tryptophans in trimerization of HIV-1 membrane-proximal external regions: Implications for virucidal design via alchemical free-energy molecular simulations.

The Dual-Action Virolytic Entry Inhibitors, or "DAVEI's," are a class of recombinant fusions of a lectin, a linker polypeptide, and a 15-residue fragment from the membrane-proximal external region (MPER) of HIV-1 gp41. DAVEI's trigger rupture of HIV-1 virions, and the interaction site between DAVEI MPER and HIV-1 lies in the gp41 component of the envelope glycoprotein Env. Here, we explore the hypothesis that DAVEI MPER engages Env gp41 in a mode structurally similar to a crystallographic MPER trimer. We used alchemical free-energy perturbation to assess the thermodynamic roles of each of the four conserved tryptophan residues on each protomer of MPER3 . We found that a W666A mutation had a large positive ΔΔG for all three protomers, while W672A had a large positive ΔΔG for only two of the three protomers, with the other tryptophans remaining unimportant contributors to MPER3 stability. The protomer for which W672 is not important is unique in the placement of its W666 sidechain between the other two protomers. We show that the unique orientation of this W666 sidechain azimuthally rotates its protomer away from the orientation it would have if the trimer were symmetric, resulting in the diminished interaction of this W672 with the rest of MPER3 . Our findings are consistent with our previous experimental study of W-to-A mutants of DAVEI. This suggests that DAVEI MPER may engage HIV-1 Env to form a mixed trimer state in which one DAVEI MPER forms a trimer by displacing a more weakly interacting protomer of the endogenous Env MPER trimer.