Disulfide folding pathways of cystine knot proteins: Tying the knot within the circular backbone of the cyclotides

The plant cyclotides are a fascinating family of circular proteins that contain a cyclic cystine knot motif. The knotted topology and cyclic nature of the cyclotides pose interesting questions about folding mechanisms and how the knotted arrangement of disulfide bonds is formed. In the current study we have examined the oxidative refolding and reductive unfolding of the prototypic cyclotide, kalata B1. A stable two-disulfide intermediate accumulated during oxidative refolding but not in reductive unfolding. Mass spectrometry and NMR spectroscopy were used to show that the intermediate contained a native-like structure with two native disulfide bonds topologically similar to the intermediate isolated for the related cystine knot protein EETI-II (LeNguyen, D., Heitz, A., Chiche, L., El Hajji, M., and Castro B. (1993) Protein Sci. 2, 165-174). However, the folding intermediate observed for kalata B1 is not the immediate precursor of the three-disulfide native peptide and does not accumulate in the reductive unfolding process, in contrast to the intermediate observed for EETI-II. These alternative pathways of linear and cyclic cystine knot proteins appear to be related to the constraints imposed by the cyclic backbone of kalata B1 and the different ring size of the cystine knot. The three-dimensional structure of a synthetic version of the two-disulfide intermediate of kalata B1 in which Ala residues replace the reduced Cys residues provides a structural insight into why the two-disulfide intermediate is a kinetic trap on the folding pathway.

Kalata B8, a novel antiviral circular protein, exhibits conformational flexibility in the cystine knot motif

The cyclotides are a family of circular proteins with a range of biological activities and potential pharmaceutical and agricultural applications. The biosynthetic mechanism of cyclization is unknown and the discovery of novel sequences may assist in achieving this goal. In the present study, we have isolated a new cyclotide from Oldenlandia affinis, kalata B8, which appears to be a hybrid of the two major subfamilies (Mobius and bracelet) of currently known cyclotides. We have determined the three-dimensional structure of kalata B8 and observed broadening of resonances directly involved in the cystine knot motif, suggesting flexibility in this region despite it being the core structural element of the cyclotides. The cystine knot motif is widespread throughout Nature and inherently stable, making this apparent flexibility a surprising result. Further-more, there appears to be isomerization of the peptide backbone at an Asp-Gly sequence in the region involved in the cyclization process. Interestingly, such isomerization has been previously characterized in related cyclic knottins from Momordica cochinchinensis that have no sequence similarity to kalata B8 apart from the six conserved cysteine residues and may result from a common mechanism of cyclization. Kalata B8 also provides insight into the structure-activity relationships of cyclotides as it displays anti-HIV activity but lacks haemolytic activity. The 'uncoupling' of these two activities has not previously been observed for the cyclotides and may be related to the unusual hydrophilic nature of the peptide.

Molecular dynamics studies on nucleoside 2',3'-cyclic phosphates

2',3'-cyclic nucleotides are intermediates and substrates of Ribonuclease (RNase)-catalysed reactions. The characterization of the equilibrium conformation as well as the flexibility inherent in these molecules helps in understanding the enzymatic action of RNases. The present study explores parameters like phase angle, glycosydic torsion angle and hydrogen bond to find possible interrelationship between them through Molecular Dynamics (MD) simulations on 3'-GMP, 3'-UMP, A>p, G>p, U>p, C>p, GpA>p and UpA>p. Interesting results of the effect of cyclisation and other constraints such as hydrogen bond between certain groups on the equilibrium ribose conformation have emerged from this study.

Torsional properties of bamboo-like structured Cu nanowires

Recently, researchers reported that nanowires (NWs) are often polycrystalline, which contain grain or twin boundaries that transect the whole NW normal to its axial direction into a bamboo like structure. In this work, large-scale molecular dynamics simulation is employed to investigate the torsional behaviours of bamboo-like structured Cu NWs. The existence of grain boundaries is found to induce a considerably large reduction to the critical angle, and the more of grain boundaries the less reduction appears, whereas, the presence of twin boundaries only results in a relatively smaller reduction to the critical angle. The introduction of grain boundaries reduces the torsional rigidity of the NW, whereas, the twin boundaries exert insignificant influence to the torsional rigidity. NWs with grain boundaries are inclined to produce a local HCP structure during loading, and the plastic deformation is usually evenly distributed along the axial axis of the NW. The plastic deformation of both perfect NW and NWs with twin boundaries is dominated by the nucleation and propagation of parallel intrinsic stacking faults. This study will enrich the current understanding of the mechanical properties of NWs, which will eventually shed lights on their applications.

Psi-angle dependence of geminal proton–carbonyl-13C coupling constants in peptides. A theoretical investigation†

Theoretical calculations of the geminal carbonyl-13C- proton coupling constant, 2J(C′H), in α-amino acids have been carried out using Dirac Vector model and Penney-Dirac bond order formulations. The results indicate that the couplings are dependent on the backbone torsion angle psi (ψ) of the amino acid residues in peptides. The meagre available experimental data seem to support the theoretical findings.

HELANAL: A program to characterise helix geometry in proteins,

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C-alpha atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.

Structural Polymorphism in d(T).d(A)*d(T) Triple Helices

DNA triple helices containing two thymine strands and one adenine strand have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that the structure with a parallel orientation of the third strand, in which the third strand base forms Hoogsteen hydrogen bonds with the adenine base in the Watson-Crick duplex, is energetically the most favourable. An antiparallel orientation of the third strand is also possible, in which the third strand base hydrogen bonds to both the bases in the Watson-Crick duplex. This structure is energetically comparable to the parallel structure. For the parallel triplex a 200ps molecular dynamics simulation starting from two different starting structures indicates that at 300K significant structural heterogeneity exists in this tripler structure. The results are compared with existing structural data on this class of triplexes derived from theoretical and NMR techniques.

Computer modelling studies of ribonuclease A-pyrimidine nucleotide complexes

Different modes of binding of pyrimidine monophosphates 2'-UMP, 3'-UMP, 2'-CMP and 3'-CMP to ribonuclease (RNase) A are studied by energy minimization in torsion angle and subsequently in Cartesian coordinate space. The results are analysed in the light of primary binding sites. The hydrogen bonding pattern brings out roles for amino acids such as Asn44 and Ser123 apart from the well known active site residues viz., His12,Lys41,Thr45 and His119. Amino acid segments 43-45 and 119-121 seem to be guiding the ligand binding by forming a pocket. Many of the active site charged residues display considerable movement upon nucleotide binding.

Extracting Ramachandran torsional angle distributions from 2D NMR data using Tikhonov regularization /

Solid state nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for studying structural and dynamical properties of disordered and partially ordered materials, such as glasses, polymers, liquid crystals, and biological materials. In particular, twodimensional( 2D) NMR methods such as ^^C-^^C correlation spectroscopy under the magicangle- spinning (MAS) conditions have been used to measure structural constraints on the secondary structure of proteins and polypeptides. Amyloid fibrils implicated in a broad class of diseases such as Alzheimer's are known to contain a particular repeating structural motif, called a /5-sheet. However, the details of such structures are poorly understood, primarily because the structural constraints extracted from the 2D NMR data in the form of the so-called Ramachandran (backbone torsion) angle distributions, g{^,'4)), are strongly model-dependent. Inverse theory methods are used to extract Ramachandran angle distributions from a set of 2D MAS and constant-time double-quantum-filtered dipolar recoupling (CTDQFD) data. This is a vastly underdetermined problem, and the stability of the inverse mapping is problematic. Tikhonov regularization is a well-known method of improving the stability of the inverse; in this work it is extended to use a new regularization functional based on the Laplacian rather than on the norm of the function itself. In this way, one makes use of the inherently two-dimensional nature of the underlying Ramachandran maps. In addition, a modification of the existing numerical procedure is performed, as appropriate for an underdetermined inverse problem. Stability of the algorithm with respect to the signal-to-noise (S/N) ratio is examined using a simulated data set. The results show excellent convergence to the true angle distribution function g{(j),ii) for the S/N ratio above 100.

Interpolating scattering amplitudes between the instant form and the front form of relativistic dynamics

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Cross-validated maximum likelihood enhances crystallographic simulated annealing refinement

Recently, the target function for crystallographic refinement has been improved through a maximum likelihood analysis, which makes proper allowance for the effects of data quality, model errors, and incompleteness. The maximum likelihood target reduces the significance of false local minima during the refinement process, but it does not completely eliminate them, necessitating the use of stochastic optimization methods such as simulated annealing for poor initial models. It is shown that the combination of maximum likelihood with cross-validation, which reduces overfitting, and simulated annealing by torsion angle molecular dynamics, which simplifies the conformational search problem, results in a major improvement of the radius of convergence of refinement and the accuracy of the refined structure. Torsion angle molecular dynamics and the maximum likelihood target function interact synergistically, the combination of both methods being significantly more powerful than each method individually. This is demonstrated in realistic test cases at two typical minimum Bragg spacings (dmin = 2.0 and 2.8 Å, respectively), illustrating the broad applicability of the combined method. In an application to the refinement of a new crystal structure, the combined method automatically corrected a mistraced loop in a poor initial model, moving the backbone by 4 Å.

Crystal structure of ammonium (3,5-dichlorophenoxy)acetate hemihydrate

In the structure of the title hydrated salt, NH4+·C8H5Cl2O3-·0.5H2O, where the anion derives from (3,5-di­chloro­phen­oxy)acetic acid, the ammonium cation is involved in extensive N-H...O hydrogen bonding with both carboxyl­ate and ether O-atom acceptors giving sheet structures lying parallel to (100). The water mol­ecule of solvation lies on a crystallographic twofold rotation axis and is involved in intra-sheet O-H...Ocarboxyl­ate hydrogen-bonding inter­actions. In the anion, the oxoacetate side chain assumes an antiperiplanar conformation with the defining C-O-C-C torsion angle = -171.33 (15)°.