The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems.

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Large loop conformation sampling using the activation relaxation technique, ART-nouveau method.

We present an adaptation of the ART-nouveau energy surface sampling method to the problem of loop structure prediction. This method, previously used to study protein folding pathways and peptide aggregation, is well suited to the problem of sampling the conformation space of large loops by targeting probable folding pathways instead of sampling exhaustively that space. The number of sampled conformations needed by ART nouveau to find the global energy minimum for a loop was found to scale linearly with the sequence length of the loop for loops between 8 and about 20 amino acids. Considering the linear scaling dependence of the computation cost on the loop sequence length for sampling new conformations, we estimate the total computational cost of sampling larger loops to scale quadratically compared to the exponential scaling of exhaustive search methods.

Ultraviolet (280-400 nm)-induced DNA damage in the eggs and larvae of Calanus finmarchicus G. (Copepoda) and Atlantic cod (Gadus morhua).

In previous work, we evaluated the effects of ultraviolet (UV = 280-400 nm) radiation on the early life stages of a planktonic Calanoid copepod (Calanus finmarchicus Gunnerus) and of Atlantic cod (Gadus morhua). Both are key species in North Atlantic food webs. To further describe the potential impacts of UV exposure on the early life stages of these two species, we measured the wavelength-specific DNA damage (cyclobutane pyrimidine dimer [CPD] formation per megabase of DNA) induced under controlled experimental exposure to UV radiation. UV-induced DNA damage in C. finmarchicus and cod eggs was highest in the UV-B exposure treatments. Under the same spectral exposures, CPD loads in C. finmarchicus eggs were higher than those in cod eggs, and for both C. finmarchicus and cod embryos, CPD loads were generally lower in eggs than in larvae. Biological weighting functions (BWF) and exposure response curves that explain most of the variability in CPD production were derived from these data. Comparison of the BWF revealed significant differences in sensitivity to UV-B: C. finmarchicus is more sensitive than cod, and larvae are more sensitive than eggs. This is consistent with the raw CPD values. Shapes of the BWF were similar to each other and to a quantitative action spectrum for damage to T7 bacteriophage DNA that is unshielded by cellular material. The strong similarities in the shapes of the weighting functions are not consistent with photoprotection by UV-absorbing compounds, which would generate features in BWF corresponding to absorption bands. The BWF reported in this study were applied to assess the mortality that would result from accumulation of a given CPD load: for both C. finmarchicus and cod eggs, an increased load of 10 CPD Mb(-1) of DNA due to UV exposure would result in approximately 10% mortality.

Molecular dynamics simulations of the bacterial ABC transporter SAV1866 in the closed form.

The ATP binding cassette (ABC) transporter family of proteins contains members involved in ATP-mediated import or export of ligands at the cell membrane. For the case of exporters, the translocation mechanism involves a large-scale conformational change that involves a clothespin-like motion from an inward-facing open state, able to bind ligands and adenosine triphosphate (ATP), to an outward-facing closed state. Our work focuses on SAV1866, a bacterial member of the ABC transporter family for which the structure is known for the closed state. To evaluate the ability of this protein to undergo conformational changes at physiological temperature, we first performed conventional molecular dynamics (MD) on the cocrystallized adenosine diphosphate (ADP)-bound structure and on a nucleotide-free structure. With this assessment of SAV1866's stability, conformational changes were induced by steered molecular dynamics (SMD), in which the nucleotide binding domains (NBD) were pushed apart, simulating the ATP hydrolysis energy expenditure. We found that the transmembrane domain is not easily perturbed by large-scale motions of the NBDs.

Use of Umbrella Sampling to Calculate the Entrance/Exit Pathway for Z-Pro-Prolinal Inhibitor in Prolyl Oligopeptidase.

Prolyl oligopeptidase (POP), a member of the prolyl endopeptidase family, is known to play a role in several neurological disorders. Its primary function is to cleave a wide range of small oligopeptides, including neuroactive peptides. We have used force biased molecular dynamics simulation to study the binding mechanism of POP. We examined three possible binding pathways using Steered Molecular Dynamics (SMD) and Umbrella Sampling (US) on a crystal structure of porcine POP with bound Z-pro-prolinal (ZPP). Using SMD, an exit pathway between the first and seventh blade of the ß-propeller domain of POP was found to be a nonviable route. US on binding pathways through the ß-propeller tunnel and the TYR190-GLN208 flexible loop at the interface between both POP domains allowed us to isolate the flexible loop pathway as the most probable. Further analysis of that pathway suggests a long-range covariation of the interdomain H-bond network, which indicates the possibility of large-scale domain reorientation observed in bacterial homologues and hypothesized to also occur in human POP.

Properties of the membrane binding component of catechol-O-methyltransferase revealed by atomistic molecular dynamics simulations.

We used atomistic simulations to study the membrane-bound form of catechol-O-methyltransferase (MB-COMT). In particular we investigated the 26-residue transmembrane α-helical segment of MB-COMT together with the 24-residue fragment that links the transmembrane component to the main protein unit that was not included in our model. In numerous independent simulations we observed the formation of a salt bridge between ARG27 and GLU40. The salt bridge closed the flexible loop that formed in the linker and kept it in the vicinity of the membrane-water interface. All simulations supported this conclusion that the linker has a clear affinity for the interface and preferentially arranges its residues to reside next to the membrane, without a tendency to relocate into the water phase. Furthermore, an extensive analysis of databases for sequences of membrane proteins that have a single transmembrane helical segment brought about an interesting view that the flexible loop observed in our work can be a common structural element in these types of proteins. In the same spirit we close the article by discussing the role of salt bridges in the formation of three-dimensional structures of membrane proteins that exhibit a single transmembrane helix.

The complex folding pathways of protein A suggest a multiple-funnelled energy landscape.

Folding proteins into their native states requires the formation of both secondary and tertiary structures. Many questions remain, however, as to whether these form into a precise order, and various pictures have been proposed that place the emphasis on the first or the second level of structure in describing folding. One of the favorite test models for studying this question is the B domain of protein A, which has been characterized by numerous experiments and simulations. Using the activation-relaxation technique coupled with a generic energy model (optimized potential for efficient peptide structure prediction), we generate more than 50 folding trajectories for this 60-residue protein. While the folding pathways to the native state are fully consistent with the funnel-like description of the free energy landscape, we find a wide range of mechanisms in which secondary and tertiary structures form in various orders. Our nonbiased simulations also reveal the presence of a significant number of non-native beta and alpha conformations both on and off pathway, including the visit, for a non-negligible fraction of trajectories, of fully ordered structures resembling the native state of nonhomologous proteins.