Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

Doxycycline reduces the migration of tuberous sclerosis complex‐2 null cells ‐ effects on RhoA‐GTPase and focal adhesion kinase

Lymphangioleiomyomatosis (LAM) is associated with dysfunction of the tuberous sclerosis complex (TSC) leading to enhanced cell proliferation and migration. This study aims to examine whether doxycycline, a tetracycline antibiotic, can inhibit the enhanced migration of TSC2‐deficient cells, identify signalling pathways through which doxycycline works and to assess the effectiveness of combining doxycycline with rapamycin (mammalian target of rapamycin complex 1 inhibitor) in controlling cell migration, proliferation and wound closure. TSC2‐positive and TSC2‐negative mouse embryonic fibroblasts (MEF), 323‐TSC2‐positive and 323‐TSC2‐null MEF and Eker rat uterine leiomyoma (ELT3) cells were treated with doxycycline or rapamycin alone, or in combination. Migration, wound closure and proliferation were assessed using a transwell migration assay, time‐lapse microscopy and manual cell counts respectively. RhoA‐GTPase activity, phosphorylation of p70S6 kinase (p70S6K) and focal adhesion kinase (FAK) in TSC2‐negative MEF treated with doxycycline were examined using ELISA and immunoblotting techniques. The enhanced migration of TSC2‐null cells was reduced by doxycycline at concentrations as low as 20 pM, while the rate of wound closure was reduced at 2–59 μM. Doxycycline decreased RhoA‐GTPase activity and phosphorylation of FAK in these cells but had no effect on the phosphorylation of p70S6K, ERK1/2 or AKT. Combining doxycycline with rapamycin significantly reduced the rate of wound closure at lower concentrations than achieved with either drug alone. This study shows that doxycycline inhibits TSC2‐null cell migration. Thus doxycycline has potential as an anti‐migratory agent in the treatment of diseases with TSC2 dysfunction.     

Exploring the conformational and binding properties of unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 through docking and molecular dynamics simulations

B‐cell lymphoma (Bcl‐2) is commonly associated with the progression and preservation of cancer and certain lymphomas; therefore, it is considered as a biological target against cancer. Nevertheless, evidence of all its structural binding sites has been hidden because of the lack of a complete Bcl‐2 model, given the presence of a flexible loop domain (FLD), which is responsible for its complex behavior. FLD region has been implicated in phosphorylation, homotrimerization, and heterodimerization associated with Bcl‐2 antiapoptotic function. In this contribution, homology modeling, molecular dynamics (MD) simulations in the microsecond (µs) time‐scale and docking calculations were combined to explore the conformational complexity of unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 systems. Conformational ensembles generated through MD simulations allowed for identifying the most populated unphosphorylated/phosphorylated monomeric conformations, which were used as starting models to obtain trimeric complexes through protein–protein docking calculations, also submitted to µs MD simulations. Principal component analysis showed that FLD represents the main contributor to total Bcl‐2 mobility, and is affected by phosphorylation and oligomerization. Subsequently, based on the most representative unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 conformations, docking studies were initiated to identify the ligand binding site of several known Bcl‐2 inhibitors to explain their influence in homo‐complex formation and phosphorylation. Docking studies showed that the different conformational states experienced by FLD, such as phosphorylation and oligomerization, play an essential role in the ability to make homo and hetero‐complexes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 393–413, 2016.     

The role of disulfide bonds and N‐terminus in the structural properties of hepcidins: Insights from molecular dynamics simulations

The main purpose of this work is to analyse, by means of molecular dynamics (MD) simulations both structural and mechanical‐dynamical differences between Hepcidin‐20 and Hepcidin‐25 in both oxidized and reduced states in aqueous solution. Results indicate that the presence of disulfide bonds is essential, in both peptides, for maintaining their β‐hairpin motif. As a matter of fact, the lack of this intra‐peptide covalent interactions produces an almost immediate deviation from the oxidized, plausibly active, structure in both the systems. Interestingly, reduced Hepcidin‐25 turns out to be characterized by a highly fluctuating structure which is found to rapidly span a large number of configurations at equilibrium. On the other hand, loss of disulfide bonds in the shorter peptide, results in a more compact and relatively rigid double‐turn structure. Comparison of mechanical–dynamical properties and sidechains–sidechains interactions in oxidized Hepcidin‐20 and Hepcidin‐25 strongly suggest also the key role of N‐terminus in the aggregation tendency of Hepcidin‐25. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 917–926, 2010.     

The β‐sheet breakers and π‐stacking

Fibrillation of β‐amyloid is recognized as a key process leading to the development of Alzheimer's disease. Small peptides called β‐sheet breakers were found to inhibit the process of β‐amyloid fibrillation and to dissolve amyloid fibrils in vitro, in vivo, and in cell culture studies [1,2]. The mechanism by which peptide inhibition takes place remains elusive and a detailed model needs to be established. Here, we present new insights into the possible role of consecutive Phe residues, present in the structure of β‐sheet breakers, supported by the results obtained by means of MD simulations. We performed a 30‐ns MD of two β‐sheet breakers: iAβ5 (LPFFD) and iAβ6 (LPFFFD) which have two and three consecutive Phe residues, respectively. We have found that Phe rings in these peptides tend to form stacked conformations. For one of the peptides – iAβ6 – the calculated electrostatic contribution to free energy of one of the conformers with three rings stacked (c2) is significantly lower than that corresponding to the unstacked one (c1), two rings stacked (c0) and second conformer with three rings stacked (c3). This may favor the interaction of the c2 conformer with the target on amyloid fibril. We hypothesize that the mechanism of inhibition of amyloidogenesis by β‐sheet breaker involves competition among π‐stacked Phe residues of the inhibitor and π‐stacking within the β‐amyloid fibril. iAβ6 may be a promising candidate for a lead compound of amyloidogenesis inhibitors. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Mechanism of allosteric propagation across a β‐sheet structure investigated by molecular dynamics simulations

The bacterial adhesin FimH consists of an allosterically regulated mannose‐binding lectin domain and a covalently linked inhibitory pilin domain. Under normal conditions, the two domains are bound to each other, and FimH interacts weakly with mannose. However, under tensile force, the domains separate and the lectin domain undergoes conformational changes that strengthen its bond with mannose. Comparison of the crystallographic structures of the low and the high affinity state of the lectin domain reveals conformational changes mainly in the regulatory inter‐domain region, the mannose binding site and a large β sheet that connects the two distally located regions. Here, molecular dynamics simulations investigated how conformational changes are propagated within and between different regions of the lectin domain. It was found that the inter‐domain region moves towards the high affinity conformation as it becomes more compact and buries exposed hydrophobic surface after separation of the pilin domain. The mannose binding site was more rigid in the high affinity state, which prevented water penetration into the pocket. The large central β sheet demonstrated a soft spring‐like twisting. Its twisting motion was moderately correlated to fluctuations in both the regulatory and the binding region, whereas a weak correlation was seen in a direct comparison of these two distal sites. The results suggest a so called “population shift” model whereby binding of the lectin domain to either the pilin domain or mannose locks the β sheet in a rather twisted or flat conformation, stabilizing the low or the high affinity state, respectively. Proteins 2016; 84:990–1008. © 2016 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Functional roles of magnesium binding to extracellular signal-regulated kinase 2 explored by molecular dynamics simulations and principal component analysis

Molecular dynamics (MD) simulations coupled with principal component (PC) analysis were carried out to study functional roles of Mg²⁺ binding to extracellular signal-regulated kinase 2 (ERK2). The results suggest that Mg²⁺ binding heavily decreases eigenvalue of the first principal component and totally inhibits motion strength of ERK2, which favors stabilization of ERK2 structure. Binding free energy predictions indicate that Mg²⁺ binding produces an important effect on binding ability of adenosine triphosphate (ATP) to ERK2 and strengthens the ATP binding. The calculations of residue-based free energy decomposition show that lack of Mg²⁺ weakens interactions between the hydrophobic rings of ATP and five residues I29, V37, A50, L105, and L154. Hydrogen bond analyses also prove that Mg²⁺ binding increases occupancies of hydrogen bonds formed between ATP and residues K52, Q103, D104, and M106. We expect that this study can provide a significant theoretical hint for designs of anticancer drugs targeting ERK2.     

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2‐3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM‐PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics. Copyright © 2014 John Wiley & Sons, Ltd.     

All-atom molecular dynamics simulations of an artificial sodium channel in a lipid bilayer: the effect of water solvation/desolvation of the sodium ion

All-atom molecular dynamics is used to investigate the transport of Na⁺ across a 1,2-dioleoyl-sn-glycero-3-phosphocholine lipid bilayer facilitated by a diazacrown hydraphile. Specifically, the free energy of Na⁺ passing through the bilayer is calculated using the adaptive biasing force method to study the free energy associated with the increase in Na⁺ transport in the presence of the hydraphile molecule. The results show that water interaction greatly influences Na⁺ transport through the lipid bilayer as water is pulled through the bilayer with Na⁺ forming a water channel. The hydraphile causes a reduction in the free energy barrier for the transport of Na⁺ through the head group part of the lipid bilayer since it complexes the Na⁺ reducing the necessity for water to be complexed and, therefore, dragged through with Na⁺, an energetically unfavorable process. The free energy associated with Na⁺ being desolvated within the bilayer is significantly decreased in the presence of the hydraphile molecule; the hydraphile increases the number of solvation states of Na⁺ that can be adopted, and this increase in the number of available configurations provides an entropic explanation for the success of the hydraphile.     

Refinement of X-ray data on dual cosubstrate specificity of CK2 kinase by free energy calculations based on molecular dynamics simulation

Free energy differences of binding of adenosine triphosphate (ATP) and guanine triphosphate (GTP) to the protein kinase CK2 (casein kinase 2) were calculated, using molecular dynamics (MD) simulations and the thermodynamic cycle approach. Good agreement with experimental data was obtained. Simulations confirm observations based on crystallographic data that specifically interacting water molecules in the binding site region of CK2 kinase play a key role in its ability to use ATP or GTP as equally efficient phosphate donors. We point out that to obtain quantitatively reasonable results, it was necessary to modify original X-ray data by assuming the presence of an additional water molecule in the CK2 binding site structure with GTP.     

Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

Essential function of the N‐termini tails of the proteasome for the gating mechanism revealed by molecular dynamics simulations

Proteasome is involved in the degradation of proteins. Proteasome activators bind to the proteasome core particle (CP) and facilitate opening a gate of the CP, where Tyr8 and Asp9 in the N‐termini tails of the CP form the ordered open gate. In a double mutant (Tyr8Gly/Asp9Gly), the N‐termini tails are disordered and the stabilized open‐gate conformation cannot be formed. To understand the gating mechanism of the CP for the translocation of the substrate, four different molecular dynamics simulations were carried out: ordered‐ and Tyr8Gly/Asp9Gly disordered‐gate models of the CP complexed with an ATP‐independent PA26 and ordered‐ and disordered‐gate models of the CP complexed with an ATP‐dependent PAN‐like activator. The free‐energies of the translocation of a polypeptide substrate moving through the gate were estimated. In the ordered‐gate models, the substrate in the activator was more stable than that in the CP. The conformational entropy of the N‐termini tails of the CP was larger when the substrate was in the activator than in the CP. In the disordered‐gate models, the substrate in the activator was more destabilized than in the ordered‐gate models. The mutated N‐termini tails became randomized and their increased conformational entropy could no longer increase further even when the substrate was in the activator, meaning the randomized N‐termini tails had lost the ability to stabilize the substrate in the activator. Thus, it was concluded that the dynamics of the N‐termini tails entropically play a key role in the translocation of the substrate. Proteins 2014; 82:1985–1999. © 2014 Wiley Periodicals, Inc.     

Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: Insights from molecular dynamics simulations

Protein aggregation into insoluble fibrillar structures known as amyloid characterizes several neurodegenerative diseases, including Alzheimer's, Huntington's and Creutzfeldt‐Jakob. Transthyretin (TTR), a homotetrameric plasma protein, is known to be the causative agent of amyloid pathologies such as FAP (familial amyloid polyneuropathy), FAC (familial amyloid cardiomiopathy) and SSA (senile systemic amyloidosis). It is generally accepted that TTR tetramer dissociation and monomer partial unfolding precedes amyloid fibril formation. To explore the TTR unfolding landscape and to identify potential intermediate conformations with high tendency for amyloid formation, we have performed molecular dynamics unfolding simulations of WT‐TTR and L55P‐TTR, a highly amyloidogenic TTR variant. Our simulations in explicit water allow the identification of events that clearly discriminate the unfolding behavior of WT and L55P‐TTR. Analysis of the simulation trajectories show that (i) the L55P monomers unfold earlier and to a larger extent than the WT; (ii) the single α‐helix in the TTR monomer completely unfolds in most of the L55P simulations while remain folded in WT simulations; (iii) L55P forms, early in the simulations, aggregation‐prone conformations characterized by full displacement of strands C and D from the main β‐sandwich core of the monomer; (iv) L55P shows, late in the simulations, severe loss of the H‐bond network and consequent destabilization of the CBEF β‐sheet of the β‐sandwich; (v) WT forms aggregation‐compatible conformations only late in the simulations and upon extensive unfolding of the monomer. These results clearly show that, in comparison with WT, L55P‐TTR does present a much higher probability of forming transient conformations compatible with aggregation and amyloid formation.     

Effect of DNA on the conformational dynamics of the endonucleases I‐DmoI as provided by molecular dynamics simulations

The conformational behavior of the wild‐type endonucleases I‐DmoI and two of its mutants has been studied in the presence and in the absence of DNA target sequences by means of extended molecular dynamics simulations. Our results show that in the absence of DNA, the three protein forms explore a similar essential conformational space, whereas when bound to the same DNA target sequence of 25 base pairs, they diversify and restrain the subspace explored. In addition, the differences in the essential subspaces explored by the residues near the catalytic site for both the bound and unbound forms are discussed in background of the experimental protein activity.     

Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation

High level of hematopoietic cell kinase (Hck) is associated with drug resistance in chronic myeloid leukemia. Additionally, Hck activity has also been connected with the pathogenesis of HIV-1 and chronic obstructive pulmonary disease. In this study, three-dimensional (3D) QSAR pharmacophore models were generated for Hck based on experimentally known inhibitors. A best pharmacophore model, Hypo1, was developed with high correlation coefficient (0.975), Low RMS deviation (0.60) and large cost difference (49.31), containing three ring aromatic and one hydrophobic aliphatic feature. It was further validated by the test set (r?=?0.96) and Fisher’s randomization method (95%). Hypo 1 was used as a 3D query for screening the chemical databases, and the hits were further screened by applying Lipinski’s rule of five and ADMET properties. Selected hit compounds were subjected to molecular docking to identify binding conformations in the active site. Finally, the appropriate binding modes of final hit compounds were revealed by molecular dynamics (MD) simulations and free energy calculation studies. Hence, we propose the final three hit compounds as virtual candidates for Hck inhibitors.     

Global fold of human cannabinoid type 2 receptor probed by solid‐state 13C‐, 15N‐MAS NMR and molecular dynamics simulations

The global fold of human cannabinoid type 2 (CB₂) receptor in the agonist‐bound active state in lipid bilayers was investigated by solid‐state ¹³C‐ and ¹⁵N magic‐angle spinning (MAS) NMR, in combination with chemical‐shift prediction from a structural model of the receptor obtained by microsecond‐long molecular dynamics (MD) simulations. Uniformly ¹³C‐ and ¹⁵N‐labeled CB₂ receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. ¹³C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of C_α, C_β, and C?O bands of superimposed resonances with predictions from protein structures generated by MD. The experimental NMR spectra matched the calculated spectra reasonably well indicating agreement of the global fold of the protein between experiment and simulations. In particular, the ¹³C chemical shift distribution of C_α resonances was shown to be very sensitive to both the primary amino acid sequence and the secondary structure of CB₂. Thus the shape of the C_α band can be used as an indicator of CB₂ global fold. The prediction from MD simulations indicated that upon receptor activation a rather limited number of amino acid residues, mainly located in the extracellular Loop 2 and the second half of intracellular Loop 3, change their chemical shifts significantly (≥1.5 ppm for carbons and ≥5.0 ppm for nitrogens). Simulated two‐dimensional ¹³C_α(i)? ¹³C?O(i) and ¹³C?O(i)? ¹⁵NH(i + 1) dipolar‐interaction correlation spectra provide guidance for selective amino acid labeling and signal assignment schemes to study the molecular mechanism of activation of CB₂ by solid‐state MAS NMR. Proteins 2014; 82:452–465. © 2013 Wiley Periodicals, Inc.     

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 MPER₃. We found that a W666A mutation had a large positive for all three protomers, while W672A had a large positive for only two of the three protomers, with the other tryptophans remaining unimportant contributors to MPER₃ 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 MPER₃. 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.     

Conformational changes of the p53-binding cleft of MDM2 revealed by molecular dynamics simulations

Two 35-ns molecular dynamics simulations of both ligated [mouse double minute protein 2 (MDM2(p53))] and unligated (MDM2(apo)) structures of human MDM2 bound to the N-terminal domain of the tumor suppressor p53 have been performed. Analysis of the dynamics revealed that the most flexible region of MDM2 was the p53-binding cleft. When MDM2 was bound to p53, a wider and more stable topology of the cleft was obtained, while unligated MDM2 showed a narrower and highly flexible cleft. It was also found that the dynamics involved in the opening/closing motions were due to the movement of different domains of the protein, which is in agreement with recent experimental data. Considering our results, a mechanism in which p53 might be recognized and attached to MDM2 is proposed, and some implications on future directions for in silico anticancer drug design efforts are discussed. In summary, the observations made here would be very useful not only for better understanding of the biological implications of the MDM2 dynamics, but also for future efforts in anticancer drug design and discovery.