Insights into recently fragmented vole populations from combined genetic and demographic data

We combined demographic and genetic data to evaluate the effects of habitat fragmentation on the population structure of the California red-backed vole (Clethrionomys californicus). We analysed variation in the mitochondrial DNA (mtDNA) control region and five nuclear microsatellite loci in small samples collected from two forest fragments and an unfragmented control site in 1990-91. We intensively sampled the same forest fragments and two different control sites in 1998 and 1999. Vole abundances fluctuated greatly at sizes below 50 individuals per fragment. Fragment populations had significantly lower mtDNA allelic diversity than controls, but not nuclear heterozygosity or numbers of alleles. The use of only trapping and/or mtDNA marker data would imply that fragment populations are at least partially isolated and vulnerable to inbreeding depression. In contrast, the abundance estimates combined with microsatellite data show that small fragment populations must be linked to nearby forests by high rates of migration. These results provide evidence for the usefulness of combining genetic and demographic data to understand nonequilibrium population structure in recently fragmented habitats.     

Response of SCP-2L domain of human MFE-2 to ligand removal: binding site closure and burial of peroxisomal targeting signal

In the study of the structure and function relationship of human MFE-2, we have investigated the dynamics of human MFE-2SCP-2L (hSCP-2L) and its response to ligand removal. A comparison was made with homologous rabbit SCP-2. Breathing and a closing motion are found, identifiable with an adjustment in size and a closing off of the binding pocket. Crucial residues for structural integrity have been identified. Particularly mobile areas of the protein are loop 1 that is connecting helices A and C in space, and helix D, next to the entrance of the pocket. In hSCP-2L, the binding pocket gets occupied by Phe93, which is making a tight hydrophobic contact with Trp36. In addition, it is found that the C-terminal peroxisomal targeting signal (PTS1) that is solvent exposed in the complexed structure becomes buried when no ligand is present. Moreover, an anti-correlation exists between burial of PTS1 and the size of the binding pocket. The results are in accordance with plant nsLTPs, where a similar accommodation of binding pocket size was found after ligand binding/removal. Furthermore, the calculations support the suggestion of a ligand-assisted targeting mechanism.     

Structure, dynamics and binding characteristics of the second PDZ domain of PTP-BL

The PDZ domains of the protein tyrosine phosphatase PTP-BL mediate interactions by binding to specific amino acid sequences in target proteins. The solution structure of the second PDZ domain of PTP-BL, PDZ2, displays a compact fold with six β strands and two α-helices. A unique feature of this domain compared to the canonical PDZ fold is an extended flexible loop at the base of the binding pocket, termed L1, that folds back onto the protein backbone, a feature that is shared by both the murine and human orthologues. The structure of PDZ2 differs significantly from the orthologous human structure. A comparison of structural quality indicators clearly demonstrates that the PDZ2 ensemble is statistically more reasonable than that of the human orthologue. The analysis of ¹⁵N relaxation data for PDZ2 shows a normal pattern, with more rigid secondary structures and more flexible loop structures. Close to the binding pocket, Leu85 and Thr88 display greater mobility when compared to surrounding residues. Peptide binding studies demonstrated a lack of interaction between murine PDZ2 and the C terminus of the murine Fas/CD95 receptor, suggesting that the Fas/CD95 receptor is not an in vivo target for PDZ2. In addition, PDZ2 specifically binds the C termini of both human Fas/CD95 receptor and the RIL protein, despite RIL containing a non-canonical PDZ-interacting sequence of E-x-V. A model of PDZ2 with the RIL peptide reveals that the PDZ2 binding pocket is able to accommodate the bulkier side-chain of glutamic acid while maintaining crucial protein to peptide hydrogen bond interactions.     

PvuII-endonuclease induces structural alterations at the scissile phosphate group of its cognate DNA

We investigated the PvuII endonuclease with its cognate DNA by means of molecular dynamics simulations. Comparing the complexed DNA with a reference simulation of free DNA, we saw structural changes at the scissile phosphodiester bond. At this GpC step, the enzyme induces the highest twist and axial rise, inclination is increased and the minor groove widened. The distance between the scissile phosphate group and the phosphate group of the following thymine base is shortened significantly, indicating a substrate-assisted catalysis. A feasible reason for this vicinity is the catalytically important amino acid residue lysine 70, which bridges the free oxygen atoms of the successive phosphate groups. Due to this geometry, a compact reaction pocket is formed where a water molecule can be held, thus bringing the reaction partners for hydrolysis into contact. The O1-P-O2 angle of the scissile nucleotide is decreased, probably due to a complexation of the negative oxygen atoms through protein and solvent contacts.     

NMR solution structure and dynamics of the peptidyl-prolyl cis-trans isomerase domain of the trigger factor from Mycoplasma genitalium compared to FK506-binding protein

We have solved the solution structure of the peptidyl-prolyl cis-trans isomerase (PPIase) domain of the trigger factor from Mycoplasma genitalium by homo- and heteronuclear NMR spectroscopy. Our results lead to a well-defined structure with a backbone rmsd of 0.23 A. As predicted, the PPIase domain of the trigger factor adopts the FK506 binding protein (FKBP) fold. Furthermore, our NMR relaxation data indicate that the dynamic behavior of the trigger factor PPIase domain and of FKBP are similar. Structural variations when compared to FKBP exist in the flap region and within the bulges of strand 5 of the beta sheet. Although the active-site crevice is similar to that of FKBP, subtle steric variations in this region can explain why FK506 does not bind to the trigger factor. Sequence variability (27% identity) between trigger factor and FKBP results in significant differences in surface charge distribution and the absence of the first strand of the central beta sheet. Our data indicate, however, that this strand may be partially structured as "nascent" beta strand. This makes the trigger factor PPIase domain the most minimal representative of the FKBP like protein family of PPIases.     

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.     

Calcium-dependent homoassociation of E-cadherin by NMR spectroscopy: changes in mobility,conformation and mapping of contact regions

Cadherins are calcium-dependent cell surface proteins that mediate homophilic cellular adhesion. The calcium-induced oligomerization of the N-terminal two domains of epithelial cadherin (ECAD12) was followed by NMR spectroscopy in solution over a large range of protein (10 microM-5 mM) and calcium (0-5 mM) concentrations. Several spectrally distinct states could be distinguished that correspond to a calcium-free monomeric form, a calcium-bound monomeric form, and to calcium-bound higher oligomeric forms. Chemical shift changes between these different states define calcium-binding residues as well as oligomerization contacts. Information about the relative orientation and mobility of the ECAD12 domains in the various states was obtained from weak alignment and 15N relaxation experiments. The data indicate that the calcium-free ECAD12 monomer adopts a flexible, kinked conformation that occludes the dimer interface observed in the ECAD12 crystal structure. In contrast, the calcium-bound monomer is already in a straight, non-flexible conformation where this interface is accessible. This mechanism provides a rational for the calcium-induced adhesiveness. Oligomerization induces chemical shift changes in an area of domain CAD1 that is centered at residue Trp-2. These shift changes extend to almost the entire surface of domain CAD1 at high (5 mM) protein concentrations. Smaller additional clusters of shift perturbations are observed around residue A80 in CAD1 and K160 in CAD2. According to weak alignment and relaxation data, the symmetry of a predominantly dimeric solution aggregate at 0.6 mM ECAD12 differs from the approximate C2-symmetry of the crystalline dimer.     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.