D614G substitution at the hinge region enhances the stability of trimeric SARS-CoV-2 spike protein

Mutations in the spike protein of SARS-CoV-2 are the major causes for the modulation of ongoing COVID-19 infection. Currently, the D614G substitution in the spike protein has become dominant worldwide. It is associated with higher infectivity than the ancestral (D614)variant. We demonstrate using Gaussian network model-based normal mode analysis that the D614G substitution occurs at the hinge region that facilitates domain-domain motions between receptor binding domain and S2 region of the spike protein. Computer-aided mutagenesis and inter-residue energy calculations reveal that contacts involving D614 are energetically frustrated. However, contacts involving G614 are energetically favourable, implying the substitution strengthens residue contacts that are formed within as well as between protomers. We also find that the free energy difference (ΔΔG) between two variants is -2.6 kcal/mol for closed and -2.0 kcal/mol for 1-RBD up conformation. Thus, the thermodynamic stability has increased upon D614G substitution. Whereas the reverse mutation in spike protein structures having G614 substitution has resulted in the free energy differences of 6.6 kcal/mol and 6.3 kcal/mol for closed and 1-RBD up conformations, respectively, indicating that the overall thermodynamic stability has decreased. These results suggest that the D614G substitution modulates the flexibility of spike protein and confers enhanced thermodynamic stability irrespective of conformational states. This data concurs with the known information demonstrating increased availability of the functional form of spikeprotein trimer upon D614G substitution.     

D614G mutation in the SARS-CoV-2 spike protein enhances viral fitness by desensitizing it to temperature-dependent denaturation

The D614G mutation in the spike protein of SARS-CoV-2 alters the fitness of the virus, leading to the dominant form observed in the COVID-19 pandemic. However, the molecular basis of the mechanism by which this mutation enhances fitness is not clear. Here we demonstrated by cryo-electron microscopy that the D614G mutation resulted in increased propensity of multiple receptor-binding domains (RBDs) in an upward conformation poised for host receptor binding. Multiple substates within the one RBD-up or two RBD-up conformational space were determined. According to negative staining electron microscopy, differential scanning calorimetry, and differential scanning fluorimetry, the most significant impact of the mutation lies in its ability to eliminate the unusual cold-induced unfolding characteristics and to significantly increase the thermal stability under physiological pH. The D614G spike variant also exhibited exceptional long-term stability when stored at 37 °C for up to 2 months. Our findings shed light on how the D614G mutation enhances the infectivity of SARS-CoV-2 through a stabilizing mutation and suggest an approach for better design of spike protein-based conjugates for vaccine development.     

Making sense of spike D614G in SARS-CoV-2 transmission

Severe acute respiratory syndrome coronavirus 2(SARS-CoV-2), the etiologic agent of the current coronavirus disease 2019(COVID-19) pandemic, has evolved to adapt to human host and transmission over the past 12 months. One prominent adaptive mutation is the asparagine-to-glycine substitution at amino acid position 614 in the viral spike protein(D614G), which has become dominant in the currently circulating virus strains. Since spike protein determines host ranges, tissue tropism, and pathogenesis through binding to the cellular receptor of angiotensin converting enzyme 2(ACE2), the D614G mutation is hypothesized to enhance viral fitness in human host, leading to increased transmission during the global pandemic. Here we summarize the recent progress on the role of the D614G mutation in viral replication, pathogenesis, transmission, and vaccine and therapeutic antibody development. These findings underscore the importance in closely monitoring viral evolution and defining their functions to ensure countermeasure efficacy against newly emerging variants.     

Engineered disulfide reveals structural dynamics of locked SARS-CoV-2 spike

The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct “S-R/x3” to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics.     

Elucidation of interactions regulating conformational stability and dynamics of SARS-CoV-2 S-protein

The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.     

Plasticity of 150-Loop in Influenza Neuraminidase Explored by Hamiltonian Replica Exchange Molecular Dynamics Simulations

Neuraminidase (NA) of influenza is a key target for antiviral inhibitors, and the 150-cavity in group-1 NA provides new insight in treating this disease. However, NA of 2009 pandemic influenza (09N1) was found lacking this cavity in a crystal structure. To address the issue of flexibility of the 150-loop, Hamiltonian replica exchange molecular dynamics simulations were performed on different groups of NAs. Free energy landscape calculated based on the volume of 150-cavity indicates that 09N1 prefers open forms of 150-loop. The turn A (residues 147–150) of the 150-loop is discovered as the most dynamical motif which induces the inter-conversion of this loop among different conformations. In the turn A, the backbone dynamic of residue 149 is highly related with the shape of 150-loop, thus can function as a marker for the conformation of 150-loop. As a contrast, the closed conformation of 150-loop is more energetically favorable in N2, one of group-2 NAs. The D147-H150 salt bridge is found having no correlation with the conformation of 150-loop. Instead the intimate salt bridge interaction between the 150 and 430 loops in N2 variant contributes the stabilizing factor for the closed form of 150-loop. The clustering analysis elaborates the structural plasticity of the loop. This enhanced sampling simulation provides more information in further structural-based drug discovery on influenza virus.     

SARS‐CoV‐2 Alpha,Beta, and Delta variants display enhanced Spike‐mediated syncytia formation

Severe COVID‐19 is characterized by lung abnormalities, including the presence of syncytial pneumocytes. Syncytia form when SARS‐CoV‐2 spike protein expressed on the surface of infected cells interacts with the ACE2 receptor on neighboring cells. The syncytia forming potential of spike variant proteins remain poorly characterized. Here, we first assessed Alpha (B.1.1.7) and Beta (B.1.351) spread and fusion in cell cultures, compared with the ancestral D614G strain. Alpha and Beta replicated similarly to D614G strain in Vero, Caco‐2, Calu‐3, and primary airway cells. However, Alpha and Beta formed larger and more numerous syncytia. Variant spike proteins displayed higher ACE2 affinity compared with D614G. Alpha, Beta, and D614G fusion was similarly inhibited by interferon‐induced transmembrane proteins (IFITMs). Individual mutations present in Alpha and Beta spikes modified fusogenicity, binding to ACE2 or recognition by monoclonal antibodies. We further show that Delta spike also triggers faster fusion relative to D614G. Thus, SARS‐CoV‐2 emerging variants display enhanced syncytia formation.     

Impairment of autophagosome-lysosome fusion in the buff mutant mice with the VPS33AD251E mutation

The HOPS (homotypic fusion and protein sorting) complex functions in endocytic and autophagic pathways in both lower eukaryotes and mammalian cells through its involvement in fusion events between endosomes and lysosomes or autophagosomes and lysosomes. However, the differential molecular mechanisms underlying these fusion processes are largely unknown. Buff (bf) is a mouse mutant that carries an Asp251-to-Glu point mutation (D251E) in the VPS33A protein, a tethering protein and a core subunit of the HOPS complex. Bf mice showed impaired spontaneous locomotor activity, motor learning, and autophagic activity. Although the gross anatomy of the brain was apparently normal, the number of Purkinje cells was significantly reduced. Furthermore, we found that fusion between autophagosomes and lysosomes was defective in bf cells without compromising the endocytic pathway. The direct association of mutant VPS33A^D251E with the autophagic SNARE complex, STX17 (syntaxin 17)-VAMP8-SNAP29, was enhanced. In addition, the VPS33A^D251E mutation enhanced interactions with other HOPS subunits, namely VPS41, VPS39, VPS18, and VPS11, except for VPS16. Reduction of the interactions between VPS33A^Y440D and several other HOPS subunits led to decreased association with STX17. These results suggest that the VPS33A^D251E mutation plays dual roles by increasing the HOPS complex assembly and its association with the autophagic SNARE complex, which selectively affects the autophagosome-lysosome fusion that impairs basal autophagic activity and induces Purkinje cell loss.     

Regulation of epitope exposure in the gp41 membrane-proximal external region through interactions at the apex of HIV-1 Env

Glycoprotein Env of human immunodeficiency virus type 1 (HIV-1) mediates viral entry through membrane fusion. Composed of gp120 and gp41 subunits arranged as a trimer-of-heterodimers, Env adopts a metastable, highly dynamic conformation on the virion surface. This structural plasticity limits the temporospatial exposure of many highly conserved, neutralizing epitopes, contributing to the difficulty in developing effective HIV-1 vaccines. Here, we employed antibody neutralization of HIV-1 infectivity to investigate how inter- and intra-gp120 interactions mediated by variable loops V1/V2 and V3 at the Env apex regulate accessibility of the gp41 membrane-proximal external region (MPER) at the Env base. Swapping the V3 loop from Env_SF162 into the Env_HXB2 background shifted MPER exposure from the prefusogenic state to a functional intermediate conformation that was distinct from the prehairpin-intermediate state sensitive to gp41-targeted fusion inhibitors. The V3-loop swap had a profound impact on global protein dynamics, biasing the equilibrium to a closed conformation resistant to most anti-gp120 antibodies, stabilizing the protein to both cold- and soluble CD4-induced Env inactivation, and increasing the CD4 requirements for viral entry. Further dissection of the Env_HXB2 V3 loop revealed that residue 306 uniquely modulated epitope exposure and trimer stability. The R306S substitution substantially decreased sensitivity to antibodies targeting the gp41 MPER and, surprisingly, the gp120 V3-loop crown (residues 312–315), but had only modest effects on exposure of intervening gp120 epitopes. Furthermore, the point mutation reduced soluble CD4-induced inactivation, but had no impact on cold inactivation. The residue appeared to exert its effects by electrostatically modifying the strength of intra-subunit interactions between the V1/V2 and V3 loops. The distinct patterns of neutralization and stability pointed to a novel prefusogenic Env conformation along the receptor activation pathway and suggested that apical Env-regulation of gp41 MPER exposure can be decoupled from much of the dynamics of gp120 subunits.     

Population Dynamics and Structural Effects at Short and Long Range Support the Hypothesis of the Selective Advantage of the G614 SARS-CoV-2 Spike Variant

SARS-CoV-2 epidemics quickly propagated worldwide, sorting virus genomic variants in newly established propagules of infections. Stochasticity in transmission within and between countries or an actual selective advantage could explain the global high frequency reached by some genomic variants. Using statistical analyses, demographic reconstructions, and molecular dynamics simulations, we show that the globally invasive G614 spike variant 1) underwent a significant demographic expansion in most countries explained neither by stochastic effects nor by overrepresentation in clinical samples, 2) increases the spike S1/S2 furin-like site conformational plasticity (short-range effect), and 3) modifies the internal motion of the receptor-binding domain affecting its cross-connection with other functional domains (long-range effect). Our results support the hypothesis of a selective advantage at the basis of the spread of the G614 variant, which we suggest may be due to structural modification of the spike protein at the S1/S2 proteolytic site, and provide structural information to guide the design of variant-specific drugs.     

Detection of the SARS-CoV-2 D614G mutation using engineered Cas12a guide RNA

Detection of pathogens with single-nucleotide variations is indispensable for the disease tracing, but remains technically challenging. The D614G mutation in the SARS-CoV-2 spike protein is known to markedly enhance viral infectivity but is difficult to detect. Here, we report an effective approach called “synthetic mismatch integrated crRNA guided Cas12a detection” (symRNA-Cas12a) to detect the D614 and G614 variants effectively. Using this method, we systemically screened a pool of crRNAs that contain all the possible nucleotide substitutions covering the -2 to +2 positions around the mutation and identify one crRNA that can efficiently increase the detection specificity by 13-fold over the ancestral crRNA. With this selected crRNA, the symRNA-Cas12a assay can detect as low as 10 copies of synthetic mutant RNA and the results are confirmed to be accurate by Sanger sequencing. Overall, we have developed the symRNA-Cas12a method to specifically, sensitively and rapidly detect the SARS-CoV-2 D614G mutation.     

NMR structure of the second intracellular loop of the alpha 2A adrenergic receptor: evidence for a novel cytoplasmic helix

A major, unresolved question in signal transduction by G protein coupled receptors (GPCRs) is to understand how, at atomic resolution, a GPCR activates a G protein. A step toward answering this question was made with the determination of the high-resolution structure of rhodopsin; we now know the intramolecular interactions that characterize the resting conformation of a GPCR. To what degree does this structure represent a structural paradigm for other GPCRs, especially at the cytoplasmic surface where GPCR-G protein interaction occurs and where the sequence homology is low among GPCRs? To address this question, we performed NMR studies on approximately 35-residue-long peptides including the critical second intracellular loop (i2) of the alpha 2A adrenergic receptor (AR) and of rhodopsin. To stabilize the secondary structure of the peptide termini, 4-12 residues from the adjacent transmembrane helices were included and structures determined in dodecylphosphocholine micelles. We also characterized the effects on an alpha 2A AR peptide of a D130I mutation in the conserved DRY motif. Our results show that in contrast to the L-shaped loop in the i2 of rhodopsin, the i2 of the alpha 2A AR is predominantly helical, supporting the hypothesis that there is structural diversity within GPCR intracellular loops. The D130I mutation subtly modulates the helical structure. The spacing of nonpolar residues in i2 with helical periodicity is a predictor of helical versus loop structure. These data should lead to more accurate models of the intracellular surface of GPCRs and of receptor-mediated G protein activation.     

Distinct dynamics and interaction patterns in H‐ and K‐Ras oncogenic P‐loop mutants

Despite years of study, the structural or dynamical basis for the differential reactivity and oncogenicity of Ras isoforms and mutants remains unclear. In this study, we investigated the effects of amino acid variations on the structure and dynamics of wild type and oncogenic mutants G12D, G12V, and G13D of H‐ and K‐Ras proteins. Based on data from µs‐scale molecular dynamics simulations, we show that the overall structure of the proteins remains similar but there are important differences in dynamics and interaction networks. We identified differences in residue interaction patterns around the canonical switch and distal loop regions, and persistent sodium ion binding near the GTP particularly in the G13D mutants. Our results also suggest that different Ras variants have distinct local structural features and interactions with the GTP, variations that have the potential to affect GTP release and hydrolysis. Furthermore, we found that H‐Ras proteins and particularly the G12V and G13D variants are significantly more flexible than their K‐Ras counterparts. Finally, while most of the simulated proteins sampled the effector‐interacting state 2 conformational state, G12V and G13D H‐Ras adopted an open switch state 1 conformation that is defective in effector interaction. These differences have implications for Ras GTPase activity, effector or exchange factor binding, dimerization and membrane interaction. Proteins 2017; 85:1618–1632. © 2017 Wiley Periodicals, Inc.     

The Val-210-Ile pathogenic Creutzfeldt–Jakob disease mutation increases both the helical and aggregation propensities of a sequence corresponding to helix-3 of PrPC

A peptide corresponding to the third helical region within the PrP^C protein, from residues 198 to 218 (helix-3), was synthesised with and without the familial 210-Val to Ile Creutzfeldt–Jakob disease mutation. The NMR structure of PrP^C predicts no global variation in stability for this mutation, indicating that local sequence rather than global structural factors are involved in the pathological effects of this mutation. ¹H NMR analysis of peptides with and without this mutation indicated that it had no significant effect on local helical structure. Temperature denaturation studies monitored by CD showed that the mutation increased the helical content within this region (helical propensity), but did not stabilise the helix toward denaturation (helical stability). Aggregation data indicated that, in addition to increasing helical propensity, this mutation increased the aggregation propensity of this sequence. CD and NMR data indicate that helical interactions, stabilised by the Val-210-Ile mutation, may precede the formation of β-sheet aggregates in this peptide sequence. Therefore, this pathological mutation probably does not facilitate PrP^C to PrP^Sc conversion by directly destabilising the helical structure of PrP^C, but may preferentially stabilise PrP^Sc by facilitating β-sheet formation within this sequence region of PrP. In addition, helical interactions between helix-3 in two or more PrP^C molecules may promote conversion to PrP^Sc.     

Mutational sensitivity of D614G in spike protein of SARS-CoV-2 in Jordan

BackgroundSpike protein is the surface glycoprotein of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) necessary for the entry of the virus via the transmembrane receptors of the human respiratory cells causing COVID-19 disease.AimHere, we aimed to predict the three-dimensional monomer structure of spike protein of SARS-CoV-2 from 20 Jordanian nasopharyngeal samples and to determine the percentage of single amino acid variants (SAV) in the spike protein of SARS-CoV-2.MethodsThe output of the Protein Homology/analogY Recognition Engine V 2.0 (Phyre2) found four single amino acid variants in the spike gene.ResultsThe first variant represented by 5% of samples that showed tyrosine deletion at Y144 located in the N terminal domain. The second and the dominant variant, represented by 62%, showed aspartate a coil amino acid substitution to glycine an extracellular amino acid at D614G located in the spike recognition binding site. The third variant, represented by 5%, showed aspartate substitution to tyrosine at D1139Y, and the fourth variant, represented by 5% glycine substitution to serine at G1167S.ConclusionOur results have shown low mutational sensitivity in all variants except to D614G the one with the most likely neutral mutational sensitivity that all variants might not explicitly affect the function of spike glycoprotein. However, D614G might change the viral conformational plasticity and hence a potential viral fitness gain but one must be cautious about drawing any concrete conclusions about the severity of symptoms and viral transmission from genomic data only.General significanceStudying mutations such as D614G in deep is essential to control the pandemic in terms of immune systems, antibodies, or even vaccines.     

Modelling conformational state dynamics and its role on infection for SARS-CoV-2 Spike protein variants

The SARS-CoV-2 Spike protein needs to be in an open-state conformation to interact with ACE2 to initiate viral entry. We utilise coarse-grained normal mode analysis to model the dynamics of Spike and calculate transition probabilities between states for 17081 variants including experimentally observed variants. Our results correctly model an increase in open-state occupancy for the more infectious D614G via an increase in flexibility of the closed-state and decrease of flexibility of the open-state. We predict the same effect for several mutations on glycine residues (404, 416, 504, 252) as well as residues K417, D467 and N501, including the N501Y mutation recently observed within the B.1.1.7, 501.V2 and P1 strains. This is, to our knowledge, the first use of normal mode analysis to model conformational state transitions and the effect of mutations on such transitions. The specific mutations of Spike identified here may guide future studies to increase our understanding of SARS-CoV-2 infection mechanisms and guide public health in their surveillance efforts.     

Molecular mechanism of HIV-1 gp120 mutations that reduce CD4 binding affinity

The interaction of the HIV-1 fusion protein gp120 with its cellular receptor CD4 represents a crucial step of the viral infection process, thus rendering gp120 a promising target for the intervention with anti-HIV drugs. Naturally occurring mutations of gp120, however, can decrease its affinity for anti-infective ligands like therapeutic antibodies or soluble CD4. To understand this phenomenon on a structural level, we performed molecular dynamics simulations of two gp120 variants (termed gp120_3-2 and gp120_2-1), which exhibit a significantly decreased binding of soluble CD4. In both variants, the exchange of a nonpolar residue byglutamate was identified as an important determinant for reduced binding. However, those glutamates are located at different sequence positions and affect different steps of the recognition process: E471 in gp120_3-2 predominantly affects the CD4-bound conformation, whereas E372 in gp120_2-1 mainly modulates the conformational sampling of free gp120. Despite these differences, there exists an interesting similarity between the two variants: both glutamates exert their function by modulating the conformation and interactions of glycine-rich motifs (G366–G367, G471–G473) resulting in an accumulation of binding incompetent gp120 conformations or a loss of intermolecular gp120–CD4 hydrogen bonds. Thus, the present data suggests that interference with the structure and dynamics of glycine-rich stretches might represent a more widespread mechanism, by which gp120 mutations reduce binding affinity. This knowledge should be helpful to predict the resistance of novel gp120 mutations or to design gp120–ligands with improved binding properties.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:41     

Elbow Flexibility of the kt38 RNA Kink-Turn Motif Investigated by Free-Energy Molecular Dynamics Simulations

Kink-turns (K-turns) are common structural motifs that can introduce sharp kinks into double-stranded RNA, and have been proposed to mediate large-scale motions in the ribosome. K-turns consist of a bulge loop region flanked by trans sugar-Hoogsteen G:A pairs, and the sharp kink conformation is stabilized by A-minor interactions (adenine contacting a G:C basepair in the minor groove). Umbrella-sampling molecular dynamics simulations were used to disrupt an A-minor interaction in the ribosomal kt38 turn and to calculate the associated free-energy change. Coupling of umbrella sampling with replica exchanges between neighboring umbrella-sampling intervals could further improve the convergence of the free-energy calculations. The simulations revealed a coupled A-minor disruption and global opening of the K-turn motif, and allowed us to characterize several intermediate A-minor conformations. The calculated free-energy profile indicated a meta-stable, semi-open structure of slightly higher free energy (∼1 kcal mol⁻¹), separated by a small free-energy barrier (∼1.5 kcal mol⁻¹) from the closed (highly kinked) form. Both K-turn states are stabilized by distinct variants of the A-minor interaction. Further opening of the K-turn structure required significantly larger free-energy changes. The semi-open form had a reduced kink angle compatible with experimental data on K-turn solution structures, and opening was coupled to a continuous global unwinding of the K-turn motif. The range of free-energy changes associated with kt38 opening and unwinding are compatible with the idea that K-turns may facilitate biologically relevant motions during large-scale ribosome dynamics.     

Stabilization of a strained protein loop conformation through protein engineering.

Staphylococcal nuclease is found in two folded conformations that differ in the isomerization of the Lys 116-Pro 117 peptide bond, resulting in two different conformations of the residue 112-117 loop. The cis form is favored over the trans with an occupancy of 90%. Previous mutagenesis studies have shown that when Lys 116 is replaced by glycine, a trans conformation is stabilized relative to the cis conformation by the release of steric strain in the trans form. However, when Lys 116 is replaced with alanine, the resulting variant protein is identical to the wild-type protein in its structure and in the dominance of the cis configuration. The results of these studies suggested that any nuclease variant with a non-glycine residue at position 116 should also favor the cis form because of steric requirements of the beta-carbon at this position. In this report, we present a structural analysis of four nuclease variants with substitutions at position 116. Two variants, K116E and K116M, follow the "beta-carbon" hypothesis by favoring the cis form. Furthermore, the crystal structure of K116E is nearly identical to that of the wild-type protein. Two additional variants, K116D and K116N, provide exceptions to this simple "beta-carbon" rule in that the trans conformation is stabilized relative to the cis configuration by these substitutions. Crystallographic data indicate that this stabilization is effected through the addition of tertiary interactions between the side chain of position 116 with the surrounding protein and water structure. The detailed trans conformation of the K116D variant appears to be similar to the trans conformation observed in the K116G variant, suggesting that these two mutations stabilize the same conformation but through different mechanisms.