Single-strand conformation polymorphism (SSCP) for the analysis of genetic variation

The accurate analysis of genetic variation has major implications in many areas of biomedical research, including the identification of infectious agents (such as parasites), the diagnosis of infections, and the detection of unknown or known disease-causing mutations. Mutation scanning methods, including PCR-coupled single-strand conformation polymorphism (SSCP), have significant advantages over many other nucleic acid techniques for the accurate analysis of allelic and mutational sequence variation. The present protocol describes the SSCP method of analysis, including all steps from the small-scale isolation of genomic DNA and PCR amplification of target sequences, through to the gel-based separation of amplicons and scanning for mutations by SSCP (either by the analysis of radiolabeled amplicons in mutation detection enhancement (MDE) gels or by non-isotopic SSCP using precast GMA gels). The subsequent sequence analysis of polymorphic bands isolated from gels is also detailed. The SSCP protocol can readily detect point mutations for amplicon sizes of up to 450-500 bp, and usually takes 1-2 days to carry out. This user-friendly, low-cost, potentially high-throughput platform has demonstrated the utility to study a wide range of pathogens and diseases, and has the potential to be applied to any gene of any organism.     

Denaturing gradient electrophoresis (DGE) and single-strand conformation polymorphism (SSCP) molecular fingerprintings revisited by simulation and used as a tool to measure microbial diversity

The exact extent of microbial diversity remains unknowable. Nevertheless, fingerprinting patterns [denaturing gradient electrophoresis (DGE), single-strand conformation polymorphism (SSCP)] provide an image of a microbial ecosystem and contain diversity data. We generated numerical simulation fingerprinting patterns based on three types of distribution (uniform, geometric and lognormal) with a range of units from 10 to 500,000. First, simulated patterns containing a diversity of around 1000 units or more gave patterns similar to those obtained in experiments. Second, the number of bands or peaks saturated quickly to about 35 and were unrelated to the degree of diversity. Finally, assuming lognormal distribution, we used an estimator of diversity on in silico and experimental fingerprinting patterns. Results on in silico patterns corresponded to the simulation inputs. Diversity results in experimental patterns were in the same range as those obtained from the same DNA sample in molecular inventories. Thus, fingerprinting patterns contain extractable data about diversity although not on the basis of a number of bands or peaks, as is generally assumed to be the case.     

Orientation and conformation of a lipase at an interface studied by molecular dynamics simulations

Electron density profiles calculated from molecular dynamics trajectories are used to deduce the orientation and conformation of Thermomyces lanuginosa lipase and a mutant adsorbed at an air-water interface. It is demonstrated that the profiles display distinct fine structures, which uniquely characterize enzyme orientation and conformation. The density profiles are, on the nanosecond timescale, determined by the average enzyme conformation. We outline a computational scheme that from a single molecular dynamics trajectory allows for extraction of electron density profiles referring to different orientations of the lipase relative to an implicit interface. Profiles calculated for the inactive and active conformations of the lipase are compared with experimental electron density profiles measured by x-ray reflectivity for the lipase adsorbed at an air-water interface. The experimental profiles contain less fine structural information than the calculated profiles because the resolution of the experiment is limited by the intrinsic surface roughness of water. Least squares fits of the calculated profiles to the experimental profiles provide areas per adsorbed enzyme and suggest that Thermomyces lanuginosa lipase adsorbs to the air-water interface in a semiopen conformation with the lid oriented away from the interface.     

Direct comparison of single-strand conformation polymorphism (SSCP) and denaturing gradient gel electrophoresis (DGGE) to characterize a microbial community on the basis of 16S rRNA gene fragments

Characterization of microbial communities using single-strand conformation polymorphism (SSCP) was compared with that using denaturing gradient gel electrophoresis (DGGE). This comparison was based on the V3-4 region (Escherichia coli positions: 341-806) of 16S rRNA gene of bacterial or archaeal communities obtained from a methanogenic bioreactor. Significant differences in the bacterial banding profiles were observed while attempting to detect the diversity of the community and its succession during the reactor operation. The SSCP produced a number of sharp bands and differentiated the bacterial community structures to which the DGGE gave an identical pattern. On the other hand, the SSCP and DGGE provided similar succession patterns for archaeal community.     

Aeromonas detection and characterization using genus-specific PCR and single-strand conformation polymorphism (SSCP)

Based on sequence alignment, oligonucleotide primers targeting the Aeromonas extracellular lipase gene were developed for PCR detection of member of the genus. A pair of primers designed for conserved regions of the gene amplified a 276?bp sequence in all Aeromonas species and tested strains, but did not have a positive result with other Gram-positive and Gram-negative bacteria, showing high specificity and sensitivity. Selective enrichment in alkaline peptone water, followed by centrifugation, and direct usage of cells suspension as template, detected initial populations of 10?c.f.u.?ml(-1). Single-strand conformation polymorphism analysis of the PCR products allowed the characterization of Aeromonas strains with a high discriminatory power (Simpson's index?=?0.988). The method presented here provides a useful tool for the rapid detection of Aeromonas and the characterization of Aeromonas isolates.     

Single-stranded conformation polymorphism (SSCP)-driven indirect sequencing in detection of short deletion

To seek for novel rare and/or sporadic mutations within genomic neighborhood of common −344 C>T (rs2427827) FCER1A proximal promoter polymorphism, which by its prevalence could have masked the presence of less frequent genetic variants in our previous single-stranded conformational polymorphism (SSCP) mutational search study, SSCP screening was repeated using primers fixing −344 C>T variant. The genomic region surrounding −344 C>T polymorphism was confirmed to be fairly conservative and only a single sporadic mutation (present in 1 out of 196 chromosomes), a 6-bp deletion −262 to 257 del CTAGAC, was found. From the methodological point of view, we demonstrated a successful detection of a short-length polymorphism using SSCP screening in a population, in which the same genomic region contained frequent single-nucleotide polymorphic variant. In conjunction with subsequent cloning of aberrantly migrating PCR products and SSCP-driven indirect sequencing of the clones, this method offers a superior (to direct sequencing of PCR products) detection of rare mutations. The cost-effective method applied by us enables a reliable characterization of infrequent (thus present only in heterozygotic form) short-length variance of the sequence which are difficult to interpret by direct sequencing.     

Homology modeling and molecular dynamics simulations of Dengue virus NS2B/NS3 protease: insight into molecular interaction

The pathogenic West Nile virus (WNV) and Dengue virus (DV) are growing global threats for which there are no specific treatments. Both viruses possess a two component NS2B/NS3 protease which cleaves viral precursor proteins. Whereas for the WNV protease two crystal structures in complex with an inhibitor have been solved recently, no such information is available for the DV protease. Here, we report the generation of a homology model of DV NS2B/NS3 protease. Since it is known from the related WNV protease that it adopts a distinct conformation in free and in inhibitor‐complexed form, a special emphasis was given to the analysis of the protease flexibility. Therefore, several models of DV NS2B/NS3 protease complexed with the peptidic inhibitor (Bz‐Nle(P4)‐Lys(P3)‐Arg(P2)‐Arg(P1)‐H) were generated. The first DV protease model (DV‐1) was constructed using the available crystal structure of the apo DV NS2B/NS3 protease. The second model (DV‐2) was built taking the WNV NS3/NS2B protease in the inhibitor‐complexed form as the template structure. Molecular dynamics simulations which were carried out for the WNV crystal structures as well as for the DV models provided an understanding of the role of NS2B for maintaining the protease in the active conformation. It was also demonstrated that NS2B is not only important for maintaining NS3 in the active form, but is also essential for establishing the interaction between residues from the S2 pocket and the peptidic inhibitor. The DV NS2B/NS3 model in the productive conformation can now be used for structure‐based design purposes. Copyright © 2009 John Wiley & Sons, Ltd.     

Tyrosine aminotransferase: biochemical and structural properties and molecular dynamics simulations

Tyrosine aminotransferase (TAT) catalyzes the transamination of tyrosine and other aromatic amino acids. The enzyme is thought to play a role in tyrosinemia type II, hepatitis and hepatic carcinoma recovery. The objective of this study is to investigate its biochemical and structural characteristics and substrate specificity in order to provide insight regarding its involvement in these diseases. Mouse TAT (mTAT) was cloned from a mouse cDNA library, and its recombinant protein was produced using Escherichia coli cells and purified using various chromatographic techniques. The recombinant mTAT is able to catalyze the transamination of tyrosine using α-ketoglutaric acid as an amino group acceptor at neutral pH. The enzyme also can use glutamate and phenylalanine as amino group donors and p-hydroxyphenylpyruvate, phenylpyruvate and alpha-ketocaproic acid as amino group acceptors. Through macromolecular crystallography we have determined the mTAT crystal structure at 2.9 ? resolution. The crystal structure revealed the interaction between the pyridoxal-5′-phosphate cofactor and the enzyme, as well as the formation of a disulphide bond. The detection of disulphide bond provides some rational explanation regarding previously observed TAT inactivation under oxidative conditions and reactivation of the inactive TAT in the presence of a reducing agent. Molecular dynamics simulations using the crystal structures of Trypanosoma cruzi TAT and human TAT provided further insight regarding the substrate-enzyme interactions and substrate specificity. The biochemical and structural properties of TAT and the binding of its cofactor and the substrate may help in elucidation of the mechanism of TAT inhibition and activation.     

Beta-hairpin conformation of fibrillogenic peptides: structure and alpha-beta transition mechanism revealed by molecular dynamics simulations

Understanding the conformational transitions that trigger the aggregation and amyloidogenesis of otherwise soluble peptides at atomic resolution is of fundamental relevance for the design of effective therapeutic agents against amyloid-related disorders. In the present study the transition from ideal alpha-helical to beta-hairpin conformations is revealed by long timescale molecular dynamics simulations in explicit water solvent, for two well-known amyloidogenic peptides: the H1 peptide from prion protein and the Abeta(12-28) fragment from the Abeta(1-42) peptide responsible for Alzheimer's disease. The simulations highlight the unfolding of alpha-helices, followed by the formation of bent conformations and a final convergence to ordered in register beta-hairpin conformations. The beta-hairpins observed, despite different sequences, exhibit a common dynamic behavior and the presence of a peculiar pattern of the hydrophobic side-chains, in particular in the region of the turns. These observations hint at a possible common aggregation mechanism for the onset of different amyloid diseases and a common mechanism in the transition to the beta-hairpin structures. Furthermore the simulations presented herein evidence the stabilization of the alpha-helical conformations induced by the presence of an organic fluorinated cosolvent. The results of MD simulation in 2,2,2-trifluoroethanol (TFE)/water mixture provide further evidence that the peptide coating effect of TFE molecules is responsible for the stabilization of the soluble helical conformation.     

Demonstration of structural polymorphism among MB3 light chains by two-dimensional gel electrophoresis

The heavy and light chain subunits of MB3 molecules were isolated from KT2 (DKT2, DR4, MB3 homozygous), ER (Dw4, DR4, MB3 homozygous), JMe (Dw5, DR5, MB3 homozygous), EBV-Sh (DSh, DRw6.2, MB3 homozygous), and EBV-Ky (DKy, DRw9, MB3 homozygous) cells and were compared with one another by two-dimensional gel electrophoresis. The MB3 light chains from KT2, ER, and EBV-Ky cells were clearly different in terms of their isoelectric points, whereas those from ER, JMe, and EBV-Sh cells were indistinguishable. No differences in charge or m.w. were noted for the MB3 heavy chains from the five cell lines. Thus, three out of the five MB3-positive, D/DR-disparate cell lines were found to express structurally distinct MB3 molecules, demonstrating that MB3 is a public serologic specificity shared by at least three structurally distinct MB (human I-A-like) molecules. Because the DR light chain subunits isolated from EBV-Wa, KT2, ER, JMe, EBV-Sh, and EBV-Ky cells differed from one another in their isoelectric points, the DR light chains were apparently more polymorphic than the MB3 light chains.     

Learning the deformation mechanism of poly(vinylidine fluoride-co-chlorotrifluoroethylene): an insight into strain-induced microstructure evolution via molecular dynamics

Learning the micro-mechanisms of fluorinated polymers during mechanical response is more difficult than that of common polymers due to the unique intrinsic characteristics of the fluorine element. In this paper, we applied molecular dynamics simulations to study deformation mechanisms of poly(vinylidine fluoride-co-chlorotrifluoroethylene) during uniaxial tension. We analyzed the variations of individual energy components and structural distribution curves versus strain in addition to the commonly used stress-strain curves and microstructure evolutions during stretching. The elastic limit is ??=?0.02, ??=?0.06 is the yield point, ??=?0.24 is the termination of the softening, necking occurs at 0.24?<???<?0.5, strain hardening occurs at 0.5?<???<?2.6, and ??=?2.6 is the damage or brake point. The elastic behavior of the material does not rely on strain rate, the obvious effect of strain rate can be seen at the yield region and strain softening region, and the stress values are not influenced by strain rates at the softening and hardening stages. Overall, total potential energy is mainly correlated with non-bonded energy, and the proportion of ΔE_coul overwhelms all the others. The energy components are ordered: ΔE_coul?>?ΔE_vdwl >?>?ΔE_angle?>?ΔE_dihed?>?ΔE_bond. The chain conformation at yield point is almost unchanged compared with the pre-stretching conformation. The chain conformations at the end of strain softening changes more obviously than that at yield point. The molecular chains maintain random coil structure before strain hardening, and switch into a stretch chain conformation gradually during strain hardening. The maximum change in bond angle during the stretching process is F-C-H, the largest change in bond length is the C-Cl bond, and the largest change in dihedral angle is H-C-C-H. The change of non-bonded interaction in the poly(VDF-co-CTFE) system is much larger than the bonding interaction, and the main factor affecting bonding interaction is the change of angles.

Open image in new window

Graphical abstract Poly(vinylidine fluoride-co-chlorotrifluoroethylene) during tension?

     

Characterisation of Trypanosoma cruzi populations by dna polymorphism of the cruzipain gene detected by single-stranded dna conformation polymorphism (SSCP) and direct sequencing

Fifty fresh isolates of Trypanosoma cruzi from Triatoma dimidiata vectors and 31 from patients with Chagas disease were analysed for DNA polymorphisms within the 432-bp core region of the cruzipain gene which encodes the active site of cathepsin L-like cystein proteinase. The cruzipain gene showed signs of polymorphism consisting of four different DNA sequences in Central and South American isolates of T. cruzi. The PCR fragments of Guatemalan isolates could be divided into three groups, Groups 1, 2 and 3, based on different patterns of single-stranded DNA conformation polymorphism. All of the strains isolated from Brazil, Chile, and Paraguay, except for the CL strain, showed a Group 4 pattern. Two to four isolates from each group were analysed by cloning and sequencing. A silent mutation occurred between Groups 1 and 2, and five nucleotides and two aa substitutions were detected between Groups 1 and 3. The DNA sequence of Group 4 contained five nucleotides and one aa substitution from Group 1. All of the DNA sequences corresponded well with the single-stranded DNA conformation polymorphism. The Group 1 isolates, the majority in the Guatemalan population (70/81, 86.4%), were isolated from both triatomines and humans, but Group 3 were isolated only from humans. Moreover, the Group 2 isolates were detected only in triatomine vectors (9/50; 18%), but never in humans (0/32, P<0.05) suggesting that this group has an independent life-cycle in sylvatic animals and is maintained by reservoir hosts other than humans.     

Unfolding of hen egg lysozyme by molecular dynamics simulations at 300K: insight into the role of the interdomain interface

We present the results of two 1.2 ns molecular dynamics (MD) unfolding simulations on hen egg lysozyme in water at 300K, performed using a new procedure called PEDC (Path Exploration With Distance Constraints). This procedure allows exploration of low energy structures as a function of increasing RMSD from the native structure, and offers especially the possibility of extensive exploration of the conformational space during the initial unfolding stages. The two independent MD simulations gave similar chronology of unfolding events: disruption of the active site, kinking of helix C, partial unfolding of the three-stranded beta-sheet to a two-stranded sheet (during which the helices A, B, and D remain to a great extent native), and finally unfolding of the beta-domain and partial unfolding of the alpha-domain in which hydrophobic clusters persist. We show particularly that the loss of hydrophobic contacts between the beta-sheet turn residues Leu55 and Ile56 and the hydrobic patch of the alpha-domain destabilizes the beta-domain and leads to its unfolding, suggesting that the correct embedding of these residues in the alpha-beta interface may constitute the rate limiting step in folding. These results are in accord with experimental observations on the folding/unfolding behavior of hen egg lysozyme at room temperature. They would also explain the loss of stability and the tendency to aggregation observed for the mutant Leu55Thr, and the slow refolding kinetics observed in the analogous amyloidogenic variant of human lysozyme.     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.