Graphene BioFET sensors for SARS-CoV-2 detection: a multiscale simulation approach.

Biological Field-Effect Transistors (BioFETs) have already demonstrated enormous potential for detecting minute amounts of ions and molecules. The use of two-dimensional (2D) materials has been shown to boost their performance and to enable the design of new applications. This combination deserves special interest in the current pandemic caused by the SARS-CoV-2 virus which demands fast, reliable and cheap detection methods. However, in spite of the experimental advances, there is a lack of a comprehensive and in-depth computational approach to capture the mechanisms underlying the sensor behaviour. Here, we present a multiscale platform that combines detailed atomic models of the molecules with mesoscopic device-level simulations. The fine-level description exploited in this approach accounts for the charge distribution of the receptor, its reconfiguration when the target binds to it, and the consequences in terms of sensitivity on the transduction mechanism. The results encourage the further exploration of improved sensor designs and 2D materials combined with diverse receptors selected to achieve the desired specificity.

Pulmonary Arterial Pruning and Longitudinal Change in Percent Emphysema and Lung Function: The Genetic Epidemiology of COPD Study.

BACKGROUND: Pulmonary endothelial damage has been shown to precede the development of emphysema in animals, and vascular changes in humans have been observed in COPD and emphysema. RESEARCH QUESTION: Is intraparenchymal vascular pruning associated with longitudinal progression of emphysema on CT imaging or decline in lung function over 5 years? STUDY DESIGN AND METHODS: The Genetic Epidemiology of COPD Study enrolled ever smokers with and without COPD from 2008 through 2011. The percentage of emphysema-like lung, or "percent emphysema," was assessed at baseline and after 5 years on noncontrast CT imaging as the percentage of lung voxels < -950 Hounsfield units. An automated CT imaging-based tool assessed and classified intrapulmonary arteries and veins. Spirometry measures are postbronchodilator. Pulmonary arterial pruning was defined as a lower ratio of small artery volume (< 5 mm2 cross-sectional area) to total lung artery volume. Mixed linear models included demographics, anthropomorphics, smoking, and COPD, with emphysema models also adjusting for CT imaging scanner and lung function models adjusting for clinical center and baseline percent emphysema. RESULTS: At baseline, the 4,227 participants were 60 ± 9 years of age, 50% were women, 28% were Black, 47% were current smokers, and 41% had COPD. Median percent emphysema was 2.1 (interquartile range, 0.6-6.3) and progressed 0.24 percentage points/y (95% CI, 0.22-0.26 percentage points/y) over 5.6 years. Mean FEV1 to FVC ratio was 68.5 ± 14.2% and declined 0.26%/y (95% CI, -0.30 to -0.23%/y). Greater pulmonary arterial pruning was associated with more rapid progression of percent emphysema (0.11 percentage points/y per 1-SD increase in arterial pruning; 95% CI, 0.09-0.16 percentage points/y), including after adjusting for baseline percent emphysema and FEV1. Arterial pruning also was associated with a faster decline in FEV1 to FVC ratio (-0.04%/y per 1-SD increase in arterial pruning; 95% CI, -0.008 to -0.001%/y). INTERPRETATION: Pulmonary arterial pruning was associated with faster progression of percent emphysema and more rapid decline in FEV1 to FVC ratio over 5 years in ever smokers, suggesting that pulmonary vascular differences may be relevant in disease progression. TRIAL REGISTRY: ClinicalTrials.gov; No.: NCT00608764; URL: www.clinicaltrials.gov.

Protein conformational flexibility modulates kinetics and thermodynamics of drug binding.

Structure-based drug design has often been restricted by the rather static picture of protein-ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery.

Three steps to gold: mechanism of protein adsorption revealed by Brownian and molecular dynamics simulations.

The addition of three N-terminal histidines to ß-lactamase inhibitor protein was shown experimentally to increase its binding potency to an Au(111) surface substantially but the binding mechanism was not resolved. Here, we propose a complete adsorption mechanism for this fusion protein by means of a multi-scale simulation approach and free energy calculations. We find that adsorption is a three-step process: (i) recognition of the surface predominantly by the histidine fusion peptide and formation of an encounter complex facilitated by a reduced dielectric screening of water in the interfacial region, (ii) adsorption of the protein on the surface and adoption of a specific binding orientation, and (iii) adaptation of the protein structure on the metal surface accompanied by induced fit. We anticipate that the mechanistic features of protein adsorption to an Au(111) surface revealed here can be extended to other inorganic surfaces and proteins and will therefore aid the design of specific protein-surface interactions.

Molecular simulations reveal that the long range fluctuations of human DPP III change upon ligand binding.

The experimentally determined structures of human dipeptidyl peptidase III (DPP III) for the wild-type protein and for the complex of its E451A mutant with the peptide substrate, tynorphin, differ significantly in their overall shape. The two domains of the enzyme are separated by a wide cleft in the structure of the ligand-free enzyme, while in the ligand-bound mutant they are very close to each other, and the protein structure is extremely compact. Here, we applied a range of molecular dynamics simulation techniques to investigate the DPP III conformational landscape and the influence of ligand binding on the protein structure and dynamics. We used conventional, accelerated and steered methods to simulate DPP III and its complexes with tynorphin and with the preferred, synthetic, substrate Arg-Arg-2-naphthylamide. We found that DPP III can adopt a number of different forms in solution. The compact forms are more stable, but the open and partially closed states, spanning a wide range of conformations, can more effectively recognize the substrate which preferentially binds to the five-stranded ß-core of the lower DPP III domain. The simulations indicated the existence of a dynamic equilibrium between open and semi-closed states and revealed two ways that the protein can close, leading to two distinct compact structures. The way in which the protein closes depends on the presence of the ligand.

Structural basis for the varying propensities of different amino acids to adopt the collagen conformation.

Although previous experimental studies have shown the positional preference of different amino acids (AAs) to form a stable triple helical collagen motif, the structural basis for the variations in the sequence and the positional propensity has not been systematically investigated. Thus, we have here probed the origin of the structural stability offered by the 20 naturally occurring AAs to collagen by means of classical molecular dynamics (MD) simulation. Simulations were carried out on 39 collagen-like peptides employing a host-guest approach. The results show that the propensity of the different AAs to adopt collagen-like conformations depends primarily on their ϕ and ψ angle preferences. Changes in these angles upon substitution of different AAs in the X(AA) and Y(AA) positions in the canonical ((Gly-X(AA)-Y(AA))(7))(3) motif dictate the formation of interchain hydrogen bonds, solvent interactions, and puckering of neighboring imino acids and, thus, the structural stability of the collagen. The role of solvent-mediated hydrogen bonds in the stabilization of collagen has also been elucidated from the MD simulations. In addition to the conventional hydrogen bonds known to be present in collagen, a hitherto unidentified direct interchain hydrogen bond, between the X(AA) N-H group and the Hyp O-H group of the neighboring chain, was observed during the simulations. Its occupancy was ∼36% when Leu was present at the X(AA) position.

Homodimeric enzymes as drug targets.

Many enzymes and proteins are regulated by their quaternary structure and/or by their association in homo- and/or hetero-oligomer complexes. Thus, these protein-protein interactions can be good targets for blocking or modulating protein function therapeutically. The large number of oligomeric structures in the Protein Data Bank (http://www.rcsb.org/) reflects growing interest in proteins that function as multimeric complexes. In this review, we consider the particular case of homodimeric enzymes as drug targets. There is intense interest in drugs that inhibit dimerization of a functionally obligate homodimeric enzyme. Because amino acid conservation within enzyme interfaces is often low compared to conservation in active sites, it may be easier to achieve drugs that target protein interfaces selectively and specifically. Two main types of dimerization inhibitors have been developed: peptides or peptidomimetics based on sequences involved in protein-protein interactions, and small molecules that act at hot spots in protein-protein interfaces. Examples include inhibitors of HIV protease and HIV integrase. Studying the mechanisms of action and locating the binding sites of such inhibitors requires different techniques for different proteins. For some enzymes, ligand binding is only detectable in vivo or after unfolding of the complexes. Here, we review the structural features of dimeric enzymes and give examples of inhibition through interference in dimer stability. Several techniques for studying these complex phenomena will be presented.

Design and characterization of a mutation outside the active site of human thymidylate synthase that affects ligand binding.

Owing to its central role in DNA synthesis, human thymidylate synthase (hTS) is a well-established target for chemotherapeutic agents, such as fluoropyrimidines. The use of hTS inhibitors in cancer therapy is limited by their toxicity and the development of cellular drug resistance. Here, with the aim of shedding light on the structural role of the A-helix in fluoropyrimidine resistance, we have created a fluoropyrimidine-resistant mutant by making a single point mutation, Glu30Trp. We postulated that residue 30, which is located in the A-helix, close to but outside the enzyme active site, could have a long-range effect on inhibitor binding. The mutant shows 100 times lower specific activity with respect to the wild-type hTS and is resistant to the classical inhibitor, FdUMP, as shown by a 6-fold higher inhibition constant. Circular dichroism experiments show that the mutant is folded. The results of molecular modeling and simulation suggest that the Glu30Trp mutation gives rise to resistance by altering the hydrogen-bond network between residue 30 and the active site.

ProSAT2--Protein Structure Annotation Server.

ProSAT2 is a server to facilitate interactive visualization of sequence-based, residue-specific annotations mapped onto 3D protein structures. As the successor of ProSAT (Protein Structure Annotation Tool), it includes its features for visualizing SwissProt and PROSITE functional annotations. Currently, the ProSAT2 server can perform automated mapping of information on variants and mutations from the UniProt KnowledgeBase and the BRENDA enzyme information system onto protein structures. It also accepts and maps user-prepared annotations. By means of an annotation selector, the user can interactively select and group residue-based information according to criteria such as whether a mutation affects enzyme activity. The visualization of the protein structures is based on the WebMol Java molecular viewer and permits simultaneous highlighting of annotated residues and viewing of the corresponding descriptive texts. ProSAT2 is available at http://projects.villa-bosch.de/mcm/database/prosat2/.

Issues in high-throughput comparative modelling: a case study using the ubiquitin E2 conjugating enzymes.

Sequences of the ubiquitin-conjugating enzyme (UBC or E2) family were used as a test set to investigate issues associated with the high-throughput comparative modelling of protein structures. A semi-automatic method was initially developed with particular emphasis on producing models of a quality suitable for structural comparison. Structural and sequence features of the E2 family were used to improve the sequence alignment and the quality of the structural templates. Initially, failure to correct for subtle structural inconsistencies between templates lead to problems in the comparative analysis of the UBC electrostatic potentials. Modelling of known UBC structures using Modeller 4.0 showed that multiple templates produced, on average, no better models than the use of just one template, as judged by the root-mean-squared deviation between the comparative model and crystal structure backbones. Using four different quality-checking methods, for a given target sequence, it was not possible to distinguish the model most similar to the experimental structure. The UBC models were thus finally modelled using only the crystal structure template with the highest sequence identity to the target to be modelled, and producing only one model solution. Quality checking was used to reject models with obvious structural anomalies (e.g., bad side-chain packing). The resulting models have been used for a comparison of UBC structural features and of their electrostatic potentials. The work was extended through the development of a fully automated pipeline that identifies E2 sequences in the sequence databases, aligns and models them, and calculates the associated electrostatic potential.

Computational approaches to structural and functional analysis of plastocyanin and other blue copper proteins.

Computational techniques are becoming increasingly important in structural and functional biology, in particular as tools to aid the interpretation of experimental results and the design of new systems. This review reports on recent studies employing a variety of computational approaches to unravel the microscopic details of the structure-function relationships in plastocyanin and other proteins belonging to the blue copper superfamily. Aspects covered include protein recognition, electron transfer and protein-solvent interaction properties of the blue copper protein family. The relevance of integrating diverse computational approaches to address the analysis of a complex protein system, such as a cupredoxin metalloprotein, is emphasized.

Electrostatic analysis and Brownian dynamics simulation of the association of plastocyanin and cytochrome f.

The oxidation of cytochrome f by the soluble cupredoxin plastocyanin is a central reaction in the photosynthetic electron transfer chain of all oxygenic organisms. Here, two different computational approaches are used to gain new insights into the role of molecular recognition and protein-protein association processes in this redox reaction. First, a comparative analysis of the computed molecular electrostatic potentials of seven single and multiple point mutants of spinach plastocyanin (D42N, E43K, E43N, E43Q/D44N, E59K/E60Q, E59K/E60Q/E43N, Q88E) and the wt protein was carried out. The experimentally determined relative rates (k(2)) for the set of plastocyanin mutants are found to correlate well (r(2) = 0.90 - 0.97) with the computed measure of the similarity of the plastocyanin electrostatic potentials. Second, the effects on the plastocyanin/cytochrome f association rate of these mutations in the plastocyanin "eastern site" were evaluated by simulating the association of the wild type and mutant plastocyanins with cytochrome f by Brownian dynamics. Good agreement between the computed and experimental relative rates (k(2)) (r(2) = 0.89 - 0.92) was achieved for the plastocyanin mutants. The results obtained by applying both computational techniques provide support for the fundamental role of the acidic residues at the plastocyanin eastern site in the association with cytochrome f and in the overall electron-transfer process.

Comparative binding energy analysis of the substrate specificity of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10.

Comparative binding energy (COMBINE) analysis was conducted for 18 substrates of the haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 (DhlA): 1-chlorobutane, 1-chlorohexane, dichloromethane, 1,2-dichloroethane, 1,2-dichloropropane, 2-chloroethanol, epichlorohydrine, 2-chloroacetonitrile, 2-chloroacetamide, and their brominated analogues. The purpose of the COMBINE analysis was to identify the amino acid residues determining the substrate specificity of the haloalkane dehalogenase. This knowledge is essential for the tailoring of this enzyme for biotechnological applications. Complexes of the enzyme with these substrates were modeled and then refined by molecular mechanics energy minimization. The intermolecular enzyme-substrate energy was decomposed into residue-wise van der Waals and electrostatic contributions and complemented by surface area dependent and electrostatic desolvation terms. Partial least-squares projection to latent structures analysis was then used to establish relationships between the energy contributions and the experimental apparent dissociation constants. A model containing van der Waals and electrostatic intermolecular interaction energy contributions calculated using the AMBER force field explained 91% (73% cross-validated) of the quantitative variance in the apparent dissociation constants. A model based on van der Waals intermolecular contributions from AMBER and electrostatic interactions derived from the Poisson-Boltzmann equation explained 93% (74% cross-validated) of the quantitative variance. COMBINE models predicted correctly the change in apparent dissociation constants upon single-point mutation of DhlA for six enzyme-substrate complexes. The amino acid residues contributing most significantly to the substrate specificity of DhlA were identified; they include Asp124, Trp125, Phe164, Phe172, Trp175, Phe222, Pro223, and Leu263. These residues are suitable targets for modification by site-directed mutagenesis.

Comparative binding energy (COMBINE) analysis of influenza neuraminidase-inhibitor complexes.

Neuraminidase is a surface glycoprotein of influenza viruses that cleaves terminal sialic acids from carbohydrates. It is critical for viral release from infected cells and facilitates viral spread in the respiratory tract. The catalytic active site of neuraminidase is highly conserved in all type A and B influenza viruses, making it an excellent target for antiinfluenza drug design. Indeed, neuraminidase inhibitors have recently become available in the clinic for the treatment of influenza. Here, we describe the use of 3D structures of neuraminidase-inhibitor complexes to derive quantitative structure-activity relationships (QSARs) to aid understanding of the mechanism of inhibition and the discovery of new inhibitors. Crystal structures of neuraminidase-inhibitor complexes were used alongside modeled complexes to derive QSAR models by COMparative BINding Energy (COMBINE) analysis (Ortiz, A. R.; Pisabarro, M. T.; Gago, F.; Wade, R. C. J. Med. Chem. 1995, 38, 2681-2691). The neuraminidase proteins studied include type A subtypes N2 and N9 (which have ca. 50% sequence identity) and an active site mutant of the N9 subtype. The inhibitors include sialic acid and benzoic acid analogues with diverse frameworks and substitution groups. By considering the contributions of the protein residues and a key water molecule to the electrostatic and van der Waals intermolecular interaction energies, a predictive and robust QSAR model for binding to type A neuraminidase was obtained. In this QSAR model, 12 protein residues and 1 bound water molecule are highlighted as particularly important for inhibitory activity. This QSAR model provides guidelines for structural modification of current inhibitors and the design of novel inhibitors in order to optimize inhibitory activity.

Protein-protein association: investigation of factors influencing association rates by brownian dynamics simulations.

The rate of protein-protein association limits the response time due to protein-protein interactions. The bimolecular association rate may be diffusion-controlled or influenced, and in such cases, Brownian dynamics simulations of protein-protein diffusional association may be used to compute association rates. Here, we report Brownian dynamics simulations of the diffusional association of five different protein-protein pairs: barnase and barstar, acetylcholinesterase and fasciculin-2, cytochrome c peroxidase and cytochrome c, the HyHEL-5 antibody and hen egg lysozyme (HEL), and the HyHEL-10 antibody and HEL. The same protocol was used to compute the diffusional association rates for all the protein pairs in order to assess, by comparison to experimentally measured rates, whether the association of these proteins can be explained solely on the basis of diffusional encounter. The simulation protocol is similar to those previously derived for simulation of the association of barnase and barstar, and of acetylcholinesterase and fasciculin-2; these produced results in excellent agreement with experimental data for these protein pairs, with changes in association rate due to mutations reproduced within the limits of expected computational and modeling errors. Here, we find that for all protein pairs, the effects of mutations can be well reproduced by the simulations, even though the degree of the electrostatic translational and orientational steering varies widely between the cases. However, the absolute values of association rates for the acetylcholinesterase: fasciculin-2 and HyHEL-10 antibody: HEL pairs are overestimated. Comparison of bound and unbound protein structures shows that this may be due to gating resulting from protein flexibility in some of the proteins. This may lower the association rates compared to their bimolecular diffusional encounter rates.

How do substrates enter and products exit the buried active site of cytochrome P450cam? 1. Random expulsion molecular dynamics investigation of ligand access channels and mechanisms.

Cytochrome P450s form a ubiquitous protein family with functions including the synthesis and degradation of many physiologically important compounds and the degradation of xenobiotics. Cytochrome P450cam from Pseudomonas putida has provided a paradigm for the structural understanding of cytochrome P450s. However, the mechanism by which camphor, the natural substrate of cytochrome P450cam, accesses the buried active site is a long-standing puzzle. While there is recent crystallographic and simulation evidence for opening of a substrate-access channel in cytochrome P450BM-3, for cytochrome P450cam, no such conformational changes have been observed either in different crystal structures or by standard molecular dynamics simulations. Here, a novel simulation method, random expulsion molecular dynamics, is presented, in which substrate-exit channels from the buried active site are found by imposing an artificial randomly oriented force on the substrate, in addition to the standard molecular dynamics force field. The random expulsion molecular dynamics method was tested in simulations of the substrate-bound structure of cytochrome P450BM-3, and then applied to complexes of cytochrome P450cam with different substrates and with product. Three pathways were identified, one of which corresponds to a channel proposed earlier on the basis of crystallographic and site-directed mutagenesis data. Exit via the water-filled channel, which was previously suggested to be a product exit channel, was not observed. The pathways obtained by the random expulsion molecular dynamics method match well with thermal motion pathways obtained by an analysis of crystallographic B-factors. In contrast to large backbone motions (up to 4 A) observed in cytochrome P450BM-3 for the exit of palmitoleic acid, passage of camphor through cytochrome P450cam only requires small backbone motions (less than 2.4 A) in conjunction with side-chain rotations. Concomitantly, in almost all the exit trajectories, salt-links that have been proposed to act as ionic tethers between secondary structure elements of the protein, are perturbed.

How do substrates enter and products exit the buried active site of cytochrome P450cam? 2. Steered molecular dynamics and adiabatic mapping of substrate pathways.

Three possible channels by which substrates and products can exit from the buried active site of cytochrome P450cam have been identified by means of random expulsion molecular dynamics simulations. In the investigation described here, we computed estimates of the relative probabilities of ligand passage through the three channels using steered molecular dynamics and adiabatic mapping. For comparison, the same techniques are also applied to investigate substrate egress from cytochrome P450-BM3. The channel in cytochrome P450cam, for which there is the most supporting evidence from experiments (which we name pathway 2a), is computed to be the most probable ligand exit channel. It has the smallest computed unbinding work and force. For this channel, the ligand exits between the F/G loop and the B' helix. Two mechanistically distinct, but energetically similar routes through this channel were observed, showing that multiple pathways along one channel are possible. The probability of ligand exit via the next most probable channel (pathway 3), which is located between the I helix and the F and G helices, is estimated to be less than 1/10 of the probability of exit along pathway 2a. Low-frequency modes of the protein extracted from an essential dynamics analysis of a 1 ns duration molecular dynamics simulation of cytochrome P450cam with camphor bound, support the opening of pathway 2a on a longer timescale. On longer timescales, it is therefore expected that this pathway becomes more dominant than estimated from the present computations.

Multicopy molecular dynamics simulations suggest how to reconcile crystallographic and product formation data for camphor enantiomers bound to cytochrome P-450cam.

Multiple ligand binding modes are possible in many enzyme active sites; their presence in cytochrome P450cam (P450cam) is evident from crystallographic studies of the binding of thiocamphor and phenylimidazoles. Here, we use multicopy molecular dynamics simulations to compare the binding modes of (1R)- and (1S)-camphor in the active site of P450cam. Simulations with (1R)-camphor, the natural substrate, serve to calibrate our protocol: 19 out of 20 copies of (1R)-camphor converged to coordinates very close to those observed for (1R)-camphor in its crystallographic complex with P450cam during the simulations. Simulations with the (1S)-camphor enantiomer showed greater mobility of the substrate, consistent with spectroscopic data, and resulted in 3 major binding modes. One of these is similar to the major conformation (of the two conformations assigned) in a recently determined crystal structure, but this conformation is not correctly oriented for regiospecific hydroxylation at C-5. The simulations, however, provide evidence for reorientation of (1S)-camphor upon formation of the reactive Fe-O intermediate to an orientation suitable for hydroxylation. The simulations thus permit rationalisation of the apparent inconsistency between the crystal structure and the reaction products.

Blue copper proteins: a comparative analysis of their molecular interaction properties.

Blue copper proteins are type-I copper-containing redox proteins whose role is to shuttle electrons from an electron donor to an electron acceptor in bacteria and plants. A large amount of experimental data is available on blue copper proteins; however, their functional characterization is hindered by the complexity of redox processes in biological systems. We describe here the application of a semiquantitative method based on a comparative analysis of molecular interaction fields to gain insights into the recognition properties of blue copper proteins. Molecular electrostatic and hydrophobic potentials were computed and compared for a set of 33 experimentally-determined structures of proteins from seven blue copper subfamilies, and the results were quantified by means of similarity indices. The analysis provides a classification of the blue copper proteins and shows that (I) comparison of the molecular electrostatic potentials provides useful information complementary to that highlighted by sequence analysis; (2) similarities in recognition properties can be detected for proteins belonging to different subfamilies, such as amicyanins and pseudoazurins, that may be isofunctional proteins; (3) dissimilarities in interaction properties, consistent with experimentally different binding specificities, may be observed between proteins belonging to the same subfamily, such as cyanobacterial and eukaryotic plastocyanins; (4) proteins with low sequence identity, such as azurins and pseudoazurins, can have sufficient similarity to bind to similar electron donors and acceptors while having different binding specificity profiles.

The cisproline(i - 1)-aromatic(i) interaction: folding of the Ala-cisPro-Tyr peptide characterized by NMR and theoretical approaches.

Cisproline(i - 1)-aromatic(i) interactions have been detected in several short peptides in aqueous solution by analysis of anomalous chemical shifts measured by 1H-NMR spectroscopy. This formation of local structure is of importance for protein folding and binding properties. To obtain an atomic-detail characterisation of the cisproline(i - 1)-aromatic(i) interaction in terms of structure, energetics and dynamics, we studied the minimal peptide unit, blocked Ala-cisPro-Tyr, using computational and experimental techniques. Structural database analyses and a systematic search revealed two groups of conformations displaying a cisproline(i - 1)-aromatic(i) interaction. These conformations were taken as seeds for molecular dynamics simulations in explicit solvent at 278 K. During a total of 33.6 ns of simulation, all the 'folded' conformations and some 'unfolded' states were sampled. 1H- and 13C-chemical shifts and 3J-coupling constants were measured for the Ala-Pro-Tyr peptide. Excellent agreement was found between all the measured and computed NMR properties, showing the good quality of the force field. We find that under the experimental and simulation conditions, the Ala-cisPro-Tyr peptide is folded 90% of the time and displays two types of folded conformation which we denote 'a' and 'b'. The type a conformations are twice as populated as the type b conformations. The former have the tyrosine ring interacting with the alanine alpha proton and are enthalpically stabilised. The latter have the aromatic ring interacting with the proline side chain and are entropically stabilised. The combined and complementary use of computational and experimental techniques permitted derivation of a detailed scenario of the 'folding' of this peptide.