Atomic insight into the CD4 binding-induced conformational changes in HIV-1 gp120

The entry of HIV-1 into a target cell requires gp120 and receptor CD4 as well as coreceptor CCR5/CXCR4 recognition events associated with conformational changes of the involved proteins. The binding of CD4 to gp120 is the initiation step of the whole process involving structural rearrangements that are crucial for subsequent pathways. Despite the wealth of knowledge about the gp120/CD4 interactions, details of the conformational changes occurring at this stage remain elusive. We have performed molecular dynamics simulations in explicit solvent based on the gp120/CD4/CD4i crystal structure in conjunction with modeled V3 and V4 loops to gain insight into the dynamics of the binding process. Three differentiated interaction modes between CD4 and gp120 were found, which involve electrostatics, hydrogen bond and van der Waals networks. A "binding funnel" model is proposed based on the dynamical nature of the binding interface together with a CD4-attraction gradient centered in gp120 at the CD4-Phe43-binding cavity. Distinct dynamical behaviors of free and CD4-bound gp120 were monitored, which likely represent the ground and pre-fusogenic states, respectively. The transition between these states revealed concerted motions in gp120 leading to: i) loop contractions around the CD4-Phe43-insertion cavity; ii) stabilization of the four-stranded "bridging sheet" structure; and iii) translocation and clustering of the V3 loop and the bridging sheet leading to the formation of the coreceptor binding site. Our results provide new insight into the dynamic of the underlying molecular recognition mechanism that complements the biochemical and structural studies.     

Domain swapping of CD4 upon dimerization

It has recently been shown that disulfide bond Cys130-Cys159 in domain 2 of monomeric CD4 is involved in the formation of CD4 disulfide-bonded dimers on cell surfaces and that it can influence the permissiveness of cells to HIV infection. Because this disulfide bond is buried in the monomer, a large conformational change must take place in order to allow for such disulfide exchange. Using standard optimization techniques, whose efficiency was first checked in the well-documented CD2 case, we have shown that 3D domain swapping is a likely candidate for the conformational change, the hinge loop, or linker, being loop E-F. Indeed, as a consequence of domain swapping, because Cys130 and Cys159 belong to beta-strands C and F, respectively, two disulfide bonds become established between Cys130 in one monomer and Cys159 in the other one. Such a disulfide exchange has already been observed when the nuclear magnetic resonance (NMR) structure of the prion protein was compared to the crystallographic, dimeric one. In both cases, domain swapping implies disulfide exchange because the linker is located in the sequence between two disulfide-bonded cysteines. As in the CD2 case, the proposed configuration of the CD4 dimer is found as a pair of neighboring monomers in the crystallographic unit cell. Moreover, because in this configuration the epitope of monoclonal antibody MT151, which does not compete with Gp120 for CD4 binding, is in the cleft between the pair of CD4 monomers, it is suggested that MT151 achieves its HIV-blocking activity by interfering with the formation of CD4 domain-swapped dimers on cell surface.     

Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.     

Morphological and molecular characterisation and differentiation of <Emphasis Type="Italic">Sargassum baccularia</Emphasis> and <Emphasis Type="Italic">S. polycystum</Emphasis> (Phaeophyta)

Two Sargassum species (S. baccularia and S. polycystum) collected from Teluk Kemang and Cape Rachado, Port Dickson, Negeri Sembilan, Malaysia, which are alike in morphology except for the rhizoidal system and vesicles, were characterised using random amplified polymorphism DNA (RAPD). The genomic DNA of both species was isolated from the leaves using a modified CTAB method. Four random primers, that is, OPA2, OPA3, OPA4 and OPA13, successfully amplified the DNA. The polymorphisms generated by these four primers were analysed using the Dice Coefficient of Similarity and cluster analysis was carried out using GelCompar II Version 2.0 (Applied Maths, Kortrijk, Belgium) based on UPGMA. DNA analysis showed that three primers were able to differentiate the two species. Morphological analysis using Principal Components Analysis (PCA) and discriminant function analysis supported the molecular data. Both species are characterised by heavily muricate main branches, oblong-lanceolate leaves with dentate margins and discoid holdfasts and spherical vesicles; both are dioecious. The only difference is that S. polycystum has secondary holdfasts transformed into stolons. This last characteristic is therefore a very important criterion and may contribute to the difference shown by DNA analysis.     

Novel Techniques for Weak Alignment of Proteins in Solution Using Chemical Tags Coordinating Lanthanide Ions

A molecule with an anisotropic magnetic susceptibility is spontaneously aligned in a static magnetic field. Alignment of such a molecule yields residual dipolar couplings and pseudocontact shifts. Lanthanide ions have recently been successfully used to provide an anisotropic magnetic susceptibility in target molecules either by replacing a calcium ion with a lanthanide ion in calcium-binding proteins or by attaching an EDTA derivative to a cysteine residue via a disulfide bond. Here we describe a novel enantiomerically pure EDTA derived tag that aligns stronger due to its shorter linker and does not suffer from stereochemical diversity upon lanthanide complexation. We observed residual (15)N,(1)H-dipolar couplings of up to 8 Hz at 800 MHz induced by a single alignment tensor from this tag.     

Rotational-echo double-resonance NMR-restrained model of the ternary complex of 5-enolpyruvylshikimate-3-phosphate synthase

The 46-kD enzyme 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase catalyzes the condensation of shikimate-3-phosphate (S3P) and phosphoenolpyruvate to form EPSP. The reaction is inhibited by N-(phosphonomethyl)-glycine (Glp), which, in the presence of S3P, binds to EPSP synthase to form a stable ternary complex. We have used solid-state NMR and molecular modeling to characterize the EPSP synthase-S3P-Glp ternary complex. Modeling began with the crystal coordinates of the unliganded protein, published distance restraints, and information from the chemical modification and mutagenesis literature on EPSP synthase. New inter-ligand and ligand-protein distances were obtained. These measurements utilized the native (31)P in S3P and Glp, biosynthetically (13)C-labeled S3P, specifically (13)C and (15)N labeled Glp, and a variety of protein-(15)N labels. Several models were investigated and tested for accuracy using the results of both new and previously published rotational-echo double resonance (REDOR) NMR experiments. The REDOR model is compared with the recently published X-ray crystal structure of the ternary complex, PDB code 1G6S. There is general agreement between the REDOR model and the crystal structure with respect to the global folding of the two domains of EPSP synthase and the relative positioning of S3P and Glp in the binding pocket. However, some of the REDOR data are in disagreement with predictions based on the coordinates of 1G6S, particularly those of the five arginines lining the binding site. We attribute these discrepancies to substantive differences in sample preparation for REDOR and X-ray crystallography. We applied the REDOR restraints to the 1G6S coordinates and created a REDOR-refined xray structure that agrees with the NMR results.     

Uniform and Residue-specific 15N-labeling of Proteins on a Highly Deuterated Background

A general method for stable-isotope labeling of large proteins is introduced and applied for studies of the E. coli GroE chaperone proteins by solution NMR. In addition to enabling the residue-specific (15)N-labeling of proteins on a highly deuterated background, it is also an efficient approach for uniform labeling. The method meets the requirements of high-level deuteration, minimal cross-labeling and high protein yield, which are crucial for NMR studies of structures with sizes above 150 kDa. The results obtained with the new protocol are compared to other strategies for protein labeling, and evaluated with regard to the influence of external factors on the resulting isotope labeling patterns. Applications with the GroE system show that these strategies are efficient tools for studies of structure, dynamics and intermolecular interactions in large supramolecular complexes, when combined with TROSY- and CRINEPT-based experimental NMR schemes.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Diffusion ordered spectroscopy as a complement to size exclusion chromatography in oligosaccharide analysis

A series of N-acetyl-chitooligosaccharides (GlcNAc)(1-6) have been studied by a nuclear magnetic resonance (NMR) method, diffusion ordered spectroscopy (DOSY). DOSY has also been applied to two additional synthetic related oligosaccharides [GlcNH(2)-(GlcNAc)(4) and GlcNH(2)-(GlcNAc)(2)-GlcNAcSO(3)Na]. A plot of the log of the determined diffusion coefficients (logD) of (GlcNAc)(n) versus the log of molecular weight was linear (6 points, R(2) = 0.995). The molecular weights of the two synthetic chitin derivatives could be estimated to within 10% error. The processed NMR data of all the chitooligosaccharides was also plotted in a polyacrylamide gel-like format to aid visual interpretation. Moreover, the logD value of the NMR signal resonances of a chitin-binding protein (hevein) changed as a function of a given titrated ligand, (GlcNAc)(6). Evidence for a 2:1 hevein:(GlcNAc)(6) complex is detected by DOSY at high hevein:(GlcNAc)(6) ratios. This data is consistent with published analytical ultracentrifugation and isothermal titration calorimetry data. A 1:1 complex is preferred at higher ligand concentrations. DOSY can complement size exclusion chromatography in carbohydrate research with the advantage that oligosaccharides are more readily detected by NMR.