MOLS sampling and its applications in structural biophysics

This review describes the MOLS method and its applications. This computational method has been developed in our laboratory primarily to explore the conformational space of small peptides and identify features of interest, particularly the minima, i.e., the low energy conformations. A systematic “brute-force” search through the vast conformational space for such features faces the insurmountable problem of combinatorial explosion, whilst other techniques, e.g., Monte Carlo searches, are somewhat limited in their region of exploration and may be considered inexhaustive. The MOLS method, on the other hand, uses a sampling technique commonly employed in experimental design theory to identify a small sample of the conformational space that nevertheless retains information about the entire space. The information is extracted using a technique that is a variant of the self-consistent mean field technique, which has been used to identify, for example, the optimal set of side-chain conformations in a protein. Applications of the MOLS method to understand peptide structure, predict the structures of loops in proteins, predict three-dimensional structures of small proteins, and arrive at the best conformation, orientation, and positions of a small molecule ligand in a protein receptor site have all yielded satisfactory results.     

Structural analysis of coxsackievirus A7 reveals conformational changes associated with uncoating

Coxsackievirus A7 (CAV7) is a rarely detected and poorly characterized serotype of the Enterovirus species Human enterovirus A (HEV-A) within the Picornaviridae family. The CAV7-USSR strain has caused polio-like epidemics and was originally thought to represent the fourth poliovirus type, but later evidence linked this strain to the CAV7-Parker prototype. Another isolate, CAV7-275/58, was also serologically similar to Parker but was noninfectious in a mouse model. Sequencing of the genomic region encoding the capsid proteins of the USSR and 275/58 strains and subsequent comparison with the corresponding amino acid sequences of the Parker strain revealed that the Parker and USSR strains are nearly identical, while the 275/58 strain is more distant. Using electron cryomicroscopy and three-dimensional image reconstruction, the structures of the CAV7-USSR virion and empty capsid were resolved to 8.2-Å and 6.1-Å resolutions, respectively. This is one of the first detailed structural analyses of the HEV-A species. Using homology modeling, reconstruction segmentation, and flexible fitting, we constructed a pseudoatomic T = 1 (pseudo T = 3) model incorporating the three major capsid proteins (VP1 to VP3), addressed the conformational changes of the capsid and its constituent viral proteins occurring during RNA release, and mapped the capsid proteins'' variable regions to the structure. During uncoating, VP4 and RNA are released analogously to poliovirus 1, the interfaces of VP2 and VP3 are rearranged, and VP1 rotates. Variable regions in the capsid proteins were predicted to map mainly to the surface of VP1 and are thus likely to affect the tropism and pathogenicity of CAV7.     

Proteasomal dysfunction and ER stress triggers 2′-hydroxy-retrochalcone-induced paraptosis in cancer cells

Chalcones are biologically active class of compounds, known for their anticancer activities. Here we show for the first time that out of the six synthetic derivatives of chalcone tested, 2′-hydroxy-retrochalcone (HRC) was the most effective in inducing extensive cytoplasmic vacuolation mediated death called paraptosis in malignant breast and cervical cancer cells. The cell death by HRC is found to be nonapoptotic in nature due to the absence of DNA fragmentation, PARP cleavage, and phosphatidylserine externalization. It was also found to be nonautophagic as there was an increase in the levels of autophagic markers LC3I, LC3II and p62. Immunofluorescence with the endoplasmic reticulum (ER) marker protein calreticulin showed that the cytoplasmic vacuoles formed were derived from the ER. This ER dilation was due to ER stress as evidenced from the increase in polyubiquitinated proteins, Bip and CHOP. Docking studies revealed that HRC could bind to the Thr1 residue on the active site of the chymotrypsin-like subunit of the proteasome. The inhibition of proteasomal activity was further confirmed by the fluorescence based assay of the chymotrypsin-like subunit of the 26S proteasome. The cell death by HRC was also triggered by the collapse of mitochondrial membrane potential and depletion of ATP. Pretreatment with thiol antioxidants and cycloheximide were able to inhibit this programmed cell death. Thus our data suggest that HRC can effectively kill cancer cells via paraptosis, an alternative death pathway and can be a potential lead molecule for anticancer therapy.     

Crystal structure of NC1 domains. Structural basis for type IV collagen assembly in basement membranes

Type IV collagen, which is present in all metazoan, exists as a family of six homologous alpha(IV) chains, alpha1-alpha6, in mammals. The six chains assemble into three different triple helical protomers and self-associate as three distinct networks. The network underlies all epithelia as a component of basement membranes, which play important roles in cell adhesion, growth, differentiation, tissue repair and molecular ultrafiltration. The specificity of both protomer and network assembly is governed by amino acid sequences of the C-terminal noncollagenous (NC1) domain of each chain. In this study, the structural basis for protomer and network assembly was investigated by determining the crystal structure of the ubiquitous [(alpha1)(2).alpha2](2) NC1 hexamer of bovine lens capsule basement membrane at 2.0 A resolution. The NC1 monomer folds into a novel tertiary structure. The (alpha1)(2).alpha2 trimer is organized through the unique three-dimensional domain swapping interactions. The differences in the primary sequences of the hypervariable region manifest in different secondary structures, which determine the chain specificity at the monomer-monomer interfaces. The trimer-trimer interface is stabilized by the extensive hydrophobic and hydrophilic interactions without a need for disulfide cross-linking.     

Finding rigid bodies in protein structures: Application to flexible fitting into cryoEM maps

We present RIBFIND, a method for detecting flexibility in protein structures via the clustering of secondary structural elements (SSEs) into rigid bodies. To test the usefulness of the method in refining atomic structures within cryoEM density we incorporated it into our flexible fitting protocol (Flex-EM). Our benchmark includes 13 pairs of protein structures in two conformations each, one of which is represented by a corresponding cryoEM map. Refining the structures in simulated and experimental maps at the 5–15 Å resolution range using rigid bodies identified by RIBFIND shows a significant improvement over using individual SSEs as rigid bodies. For the 15 Å resolution simulated maps, using RIBFIND-based rigid bodies improves the initial fits by 40.64% on average, as compared to 26.52% when using individual SSEs. Furthermore, for some test cases we show that at the sub-nanometer resolution range the fits can be further improved by applying a two-stage refinement protocol (using RIBFIND-based refinement followed by an SSE-based refinement). The method is stand-alone and could serve as a general interactive tool for guiding flexible fitting into EM maps.